Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,915,089 G→T H38N (CAT→AAT)  glmU ← fused N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,915,0890GT100.0% 233.1 / NA 60H38N (CAT→AAT) glmUfused N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (32/28);  total (32/28)

TGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATGCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTA  >  NC_000913/3914941‑3915233
                                                                                                                                                    |                                                                                                                                                
tgcattgcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCt                                                                                                                                                >  2:33296/1‑151 (MQ=255)
  cattgcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTg                                                                                                                                               >  2:822931/1‑150 (MQ=255)
      gcatGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAAcc                                                                                                                                           >  1:531427/1‑150 (MQ=255)
          gACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCg                                                                                                                                       <  2:274802/150‑1 (MQ=255)
          gACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGc                                                                                                                                      >  1:19800/1‑151 (MQ=255)
            cccGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCtt                                                                                                                                    <  1:446135/151‑1 (MQ=255)
              cGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTc                                                                                                                                  <  1:452934/151‑1 (MQ=255)
                 aCCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTccc                                                                                                                                <  2:118965/150‑1 (MQ=255)
                 aCCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCg                                                                                                                               >  1:168818/1‑151 (MQ=255)
                   ccAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGc                                                                                                                             <  1:258799/151‑1 (MQ=255)
                    cAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGGACCATCGCTTTCCCGGCa                                                                                                                            <  1:919753/151‑1 (MQ=255)
                    cAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCa                                                                                                                            <  2:118126/151‑1 (MQ=255)
                          ctctGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGt                                                                                                                       >  1:970314/1‑150 (MQ=255)
                            ctGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTa                                                                                                                      <  1:540743/149‑1 (MQ=255)
                              gCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTAt                                                                                                                     <  2:556865/148‑1 (MQ=255)
                              gCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATgca                                                                                                                  <  1:677030/151‑1 (MQ=255)
                                cTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATgcagc                                                                                                                >  1:306480/1‑151 (MQ=255)
                                       aCCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCg                                                                                                         >  2:311649/1‑151 (MQ=255)
                                            gTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGa                                                                                                    <  2:717696/151‑1 (MQ=255)
                                            gTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGa                                                                                                    <  2:964730/151‑1 (MQ=255)
                                              tAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGa                                                                                                       >  1:130927/1‑146 (MQ=255)
                                              tAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGAt                                                                                                   <  1:867371/150‑1 (MQ=255)
                                               aaGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGa                                                                                                       >  2:386035/1‑145 (MQ=255)
                                               aaGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATc                                                                                                  <  1:130006/150‑1 (MQ=255)
                                                         tcTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATgc                                                                                       >  1:566705/1‑151 (MQ=255)
                                                          cTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATgcg                                                                                      >  1:60183/1‑151 (MQ=255)
                                                            ttCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATgcgc                                                                                     <  1:172485/150‑1 (MQ=255)
                                                               aGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGc                                                                                 >  2:954246/1‑151 (MQ=255)
                                                               aGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGc                                                                                 >  2:696234/1‑151 (MQ=255)
                                                                     tGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGc                                                                           >  1:729211/1‑151 (MQ=255)
                                                                        tttAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCt                                                                         <  2:878264/150‑1 (MQ=255)
                                                                           aGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCgg                                                                     <  1:226947/151‑1 (MQ=255)
                                                                              aGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCa                                                                   >  2:229652/1‑150 (MQ=255)
                                                                               gATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCa                                                                   >  1:120908/1‑149 (MQ=255)
                                                                                aTCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAgg                                                                >  1:135935/1‑151 (MQ=255)
                                                                                aTCGCCGCCGTGACCGTACACCAGGTCAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAg                                                                 <  1:756784/150‑1 (MQ=255)
                                                                                 tCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAgg                                                                >  1:745254/1‑150 (MQ=255)
                                                                                  cgccgcCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGAt                                                              >  2:913290/1‑151 (MQ=255)
                                                                                     cgccgTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCAc                                                           >  2:193012/1‑151 (MQ=255)
                                                                                      gccgTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACt                                                          >  2:668384/1‑151 (MQ=255)
                                                                                        cgTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTAc                                                        <  1:694282/151‑1 (MQ=255)
                                                                                        cgTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTAc                                                        >  1:840701/1‑151 (MQ=255)
                                                                                             ccGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTc                                                    >  2:15983/1‑150 (MQ=255)
                                                                                                tACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCAt                                                  >  2:160689/1‑149 (MQ=255)
                                                                                                   acCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCat                                             >  1:295404/1‑151 (MQ=255)
                                                                                                              gTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACAt                                   >  1:677289/1‑150 (MQ=255)
                                                                                                                  gctgcgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATAcgcg                              <  1:241567/151‑1 (MQ=255)
                                                                                                                  gctgcgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATAcgcg                              <  1:721291/151‑1 (MQ=255)
                                                                                                                   ctgcgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGt                             >  2:162783/1‑151 (MQ=255)
                                                                                                                    tgcgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTc                            >  1:840325/1‑151 (MQ=255)
                                                                                                                     gcgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTc                            <  1:63570/150‑1 (MQ=255)
                                                                                                                     gcgcCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCACGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTc                            <  2:351948/150‑1 (MQ=255)
                                                                                                                        ccTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTg                         >  2:867885/1‑150 (MQ=255)
                                                                                                                              tcattcGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTa                   >  1:893150/1‑150 (MQ=255)
                                                                                                                               cattcGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTaa                  <  1:620891/150‑1 (MQ=255)
                                                                                                                                  tcGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTaa                  <  2:155390/147‑1 (MQ=255)
                                                                                                                                   cGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTg              <  2:168832/150‑1 (MQ=255)
                                                                                                                                          tCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGa       >  1:219785/1‑150 (MQ=255)
                                                                                                                                            aaTGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAAtt    <  1:229633/151‑1 (MQ=255)
                                                                                                                                               gACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTa  <  1:635640/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATGCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTTA  >  NC_000913/3914941‑3915233

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: