Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,317,992 C→T 5.4% intergenic (‑97/+14) phnL ← / ← phnK ribophosphonate triphosphate synthase subunit; putative ABC transporter‑related ATPase/carbon‑phosphorus lyase complex subunit, putative ATP transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,317,9920CT 5.4% 417.0 / 7.9 130intergenic (‑97/+14)phnL/phnKribophosphonate triphosphate synthase subunit; putative ABC transporter‑related ATPase/carbon‑phosphorus lyase complex subunit, putative ATP transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (65/57);  new base T (2/5);  total (67/63)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 135 total aligned reads displayed.

CCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGG  >  NC_000913/4317843‑4318139
                                                                                                                                                     |                                                                                                                                                   
ccGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTAcgccgc                                                                                                                                                     <  1:479825/150‑1 (MQ=255)
    tttGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                                 <  1:409421/150‑1 (MQ=38)
    tttGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                                 <  1:1234872/150‑1 (MQ=38)
     ttGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCACATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                                >  1:624678/1‑150 (MQ=255)
        cTGGTGCATGATTAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTGCGGGTGTCGGATGCTACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                             <  2:1090477/150‑2 (MQ=25)
        cTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcc                                                                                                                                             >  1:590170/1‑149 (MQ=255)
        cTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                             >  2:671168/1‑150 (MQ=255)
         ttgTGCAGGATGAATTTTTTACTGACGATTTGTACGTTAATCATTTTTCGCCTCGTGTTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAc                                                                                                                                            <  2:860730/148‑1 (MQ=255)
          ggTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAcc                                                                                                                                           <  1:196849/150‑1 (MQ=255)
          ggTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAcc                                                                                                                                           <  2:929332/150‑1 (MQ=255)
            tGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCgg                                                                                                                                         >  2:1635490/1‑150 (MQ=255)
            tGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTCATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCgg                                                                                                                                         >  1:867219/1‑150 (MQ=255)
              cAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCt                                                                                                                                       >  1:699259/1‑150 (MQ=255)
                ggATGAAGGTTTTACTGTCGGTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGGCGGATGCGACGCTGACGCGTCGTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCa                                                                                                                                     <  2:959108/150‑1 (MQ=255)
                  ctGAAGGTTTTACTGCCGATTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGTTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAAt                                                                                                                                   <  1:188274/149‑1 (MQ=255)
                  aTGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGGCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAAt                                                                                                                                   <  2:558707/150‑1 (MQ=255)
                      aGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTg                                                                                                                               <  1:42686/150‑1 (MQ=255)
                      aGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCACATCCGGCACCGGCTCAATTCTg                                                                                                                               <  2:1325019/150‑1 (MQ=255)
                       ggTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGTTGTCGGATTCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTATGCCGGATAATTCGTTTACGCCGCATCCGGCACCGGCTCAATTcatc                                                                                                                              >  1:1667332/1‑147 (MQ=37)
                       ggTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGc                                                                                                                              <  2:437567/150‑1 (MQ=255)
                       ggTTTTACTGACGTTTTGTACGTAAATCATTTTTCGCCTCGTGGTTCGGGGGTCGGCTGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTGGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGAACCGGCTCAATTCTGc                                                                                                                              <  1:1588896/150‑1 (MQ=37)
                              cTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCg                                                                                                                       >  1:118081/1‑150 (MQ=255)
                              cTGACGATTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGTAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCg                                                                                                                       >  1:881655/1‑150 (MQ=255)
                               tGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGa                                                                                                                      <  2:1154660/150‑1 (MQ=255)
                                gACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGAt                                                                                                                     <  1:60082/150‑1 (MQ=255)
                                 aCGTTTTGTACGTTATTCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATg                                                                                                                    <  1:329905/150‑1 (MQ=255)
                                      ttGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTTGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACAcc                                                                                                               <  2:832552/150‑1 (MQ=255)
                                        gTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACcag                                                                                                             >  1:1422393/1‑150 (MQ=255)
                                         tACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACcagc                                                                                                            >  1:238803/1‑150 (MQ=255)
                                         tACGTTAATCATTTTTCGCCTCGTCGTTCGGGTGTCGGCTGCGAAGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGCCACaagc                                                                                                            >  1:1230836/1‑150 (MQ=25)
                                          aCGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACcagca                                                                                                           <  1:899531/150‑1 (MQ=255)
                                          aCGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACcagca                                                                                                           >  1:1134145/1‑150 (MQ=255)
                                          aCGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACcagca                                                                                                           >  2:410946/1‑150 (MQ=255)
                                          aCGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACcagca                                                                                                           >  2:1798870/1‑150 (MQ=255)
                                           cGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACcagcag                                                                                                          >  2:132087/1‑150 (MQ=255)
                                           cGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACcagcag                                                                                                          >  1:1741185/1‑150 (MQ=255)
                                              tAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCCGCtg                                                                                                       >  2:112106/1‑150 (MQ=255)
                                              tAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCtg                                                                                                       <  1:1712225/150‑1 (MQ=255)
                                               aaTCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCtgt                                                                                                      >  2:1467126/1‑150 (MQ=255)
                                                aTCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCtgtg                                                                                                     <  1:1276364/150‑1 (MQ=255)
                                                aTCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCtgtg                                                                                                     <  1:1754159/150‑1 (MQ=255)
                                                 tCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTgtgt                                                                                                    <  2:395555/150‑1 (MQ=255)
                                                  cATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCCCCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGta                                                                                                   >  1:1360666/1‑150 (MQ=255)
                                                  cATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGta                                                                                                   <  1:1586191/150‑1 (MQ=255)
                                                  cATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATACGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGta                                                                                                   >  1:1368302/1‑150 (MQ=255)
                                                   aTTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGtat                                                                                                  >  1:357621/1‑150 (MQ=255)
                                                    tttttCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGGGtata                                                                                                 >  2:1326579/1‑150 (MQ=255)
                                                       ttCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACgg                                                                                              <  1:1121013/150‑1 (MQ=255)
                                                       ttCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACgg                                                                                              <  1:1126940/150‑1 (MQ=255)
                                                       ttCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCAGTTGTAGGCCTGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACgg                                                                                              <  2:1090629/150‑1 (MQ=255)
                                                         cGCCTCGTGGTTCGGGTGGCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCAATGACACCAGCAGCTGTGTATACGgat                                                                                            <  2:1794556/150‑1 (MQ=255)
                                                           ccTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGgatga                                                                                          >  1:1789449/1‑150 (MQ=255)
                                                                 ggTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCggg                                                                                    >  1:1045977/1‑150 (MQ=255)
                                                                 ggTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCggg                                                                                    <  2:210608/150‑1 (MQ=255)
                                                                            cGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCAc                                                                         <  1:783760/150‑1 (MQ=255)
                                                                              gATGCGACGCTGACGCGTCTTATCCGACCGACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTGTACGCCGCATCCGGCACCGGCTCAATTCTGCAAACCCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCAcgc                                                                       <  1:272006/150‑1 (MQ=255)
                                                                                  cGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGcggtc                                                                   >  2:1056982/1‑150 (MQ=255)
                                                                                  cGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGGCGAGCACGcggtc                                                                   >  1:1306330/1‑150 (MQ=255)
                                                                                        ggACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTaa                                                             <  1:150542/149‑1 (MQ=255)
                                                                                         gACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAATCCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAAc                                                            >  1:810051/1‑150 (MQ=255)
                                                                                         gACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAAc                                                            >  1:1160236/1‑150 (MQ=255)
                                                                                          aCGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAAcc                                                           <  1:376148/150‑1 (MQ=255)
                                                                                             cgTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCAc                                                        <  2:1714987/150‑1 (MQ=255)
                                                                                              gTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAAATCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCAct                                                       >  1:1036778/1‑150 (MQ=255)
                                                                                               tCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGGTAACCCActg                                                      >  2:349186/1‑149 (MQ=255)
                                                                                                cTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCActct                                                     <  2:789717/150‑1 (MQ=255)
                                                                                                 ttATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACtctc                                                    <  1:776887/150‑1 (MQ=255)
                                                                                                  tATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTcc                                                   >  1:1637331/1‑150 (MQ=255)
                                                                                                  tATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCAACTCTcc                                                   >  1:787938/1‑150 (MQ=255)
                                                                                                   aTCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTcca                                                  >  1:29334/1‑150 (MQ=255)
                                                                                                    tCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGCTGCGGGTCGTCGAGCACGCGGTCGGTTAACCCCCTCTccac                                                 >  2:1077292/1‑150 (MQ=255)
                                                                                                       gACCTACGGGGACCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAACCCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGGTAAGCCACTCTccaccgc                                              >  2:158218/1‑148 (MQ=255)
                                                                                                         ccTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACtt                                            <  2:1356513/150‑1 (MQ=255)
                                                                                                           tACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGc                                          >  1:1623962/1‑150 (MQ=255)
                                                                                                               gggAGCGCTTTTGTAGGCCGGATAATGCGTTTACGCCGAATCCGGCACCGGCTCAATTCTACAAAACCTATGACACGAGCAGCTGTGTATTCGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCt                                      <  2:695401/150‑1 (MQ=255)
                                                                                                                  aGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCt                                   <  2:166668/150‑1 (MQ=255)
                                                                                                                  aGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCt                                   <  1:1467126/150‑1 (MQ=255)
                                                                                                                   gcgcATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCtt                                  >  2:970886/1‑150 (MQ=255)
                                                                                                                   gcgcATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCtt                                  <  2:496337/150‑1 (MQ=255)
                                                                                                                    cgcATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCCCTCTCCACCACTTGCCCCTGCTtc                                 >  2:478544/1‑150 (MQ=255)
                                                                                                                    cgcATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCTCTGCTtc                                 <  2:1467456/150‑1 (MQ=255)
                                                                                                                    cgcATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTtc                                 >  1:746745/1‑150 (MQ=255)
                                                                                                                    cgcATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACAGATGACACCAGCAGCTGTGTATACTGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTtc                                 <  2:124962/150‑1 (MQ=255)
                                                                                                                        tttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTtcatca                             <  1:671168/150‑1 (MQ=255)
                                                                                                                             aGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGc                        >  1:1666898/1‑150 (MQ=255)
                                                                                                                             aGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGc                        >  1:7181/1‑150 (MQ=255)
                                                                                                                              tgCCGGATGCGGCGTAAACGCCTTATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCTCCAGCa                       >  1:324987/2‑150 (MQ=255)
                                                                                                                               gCCGGATAAGGCGTTTACGCCGCATCCGGCAGCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGTGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCaa                      >  1:1165341/1‑150 (MQ=255)
                                                                                                                               gCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCaa                      >  2:529278/1‑150 (MQ=255)
                                                                                                                                ccGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCaaa                     >  1:139657/1‑150 (MQ=255)
                                                                                                                                     tAAGGCGTTTACGCCGCCTCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGCGCACGCGGTCGGTTAACCCACTCTCCCCCACTTGCCCCTGCTGCCTCCCCAGCAAACGGTc                >  1:248775/1‑150 (MQ=255)
                                                                                                                                       aGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCCCTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCg              >  2:218702/1‑150 (MQ=255)
                                                                                                                                       aGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCg              <  1:1398899/150‑1 (MQ=255)
                                                                                                                                         gCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGcc            <  2:867219/150‑1 (MQ=255)
                                                                                                                                           gTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCTACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCcag          >  1:1019272/1‑150 (MQ=255)
                                                                                                                                           gTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCcag          <  2:617415/150‑1 (MQ=255)
                                                                                                                                            tttACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCcagc         <  1:954217/150‑1 (MQ=255)
                                                                                                                                                cgccgcATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAggc     <  2:238803/150‑1 (MQ=255)
                                                                                                                                                 gccgcATCCGGCACCGGATCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAggcg    >  1:6325/1‑150 (MQ=255)
                                                                                                                                                   cgcATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAggcggg  >  2:1793961/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGTTTTGCTGGTGCAGGATGAAGGTTTTACTGACGTTTTGTACGTTAATCATTTTTCGCCTCGTGGTTCGGGTGTCGGATGCGACGCTGACGCGTCTTATCCGACCTACGGGGAGCGCATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTGCCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGG  >  NC_000913/4317843‑4318139

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑