Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,918,291 C→T 5.2% intergenic (‑26/+25) atpG ← / ← atpA F1 sector of membrane‑bound ATP synthase, gamma subunit/F1 sector of membrane‑bound ATP synthase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,918,2910CT 5.2% 440.9 / 7.5 135intergenic (‑26/+25)atpG/atpAF1 sector of membrane‑bound ATP synthase, gamma subunit/F1 sector of membrane‑bound ATP synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (63/65);  new base T (5/2);  total (68/67)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

GGCCGCCATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTA  >  NC_000913/3918143‑3918438
                                                                                                                                                    |                                                                                                                                                   
ggCCGCCATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCt                                                                                                                                                    >  1:1718919/1‑150 (MQ=255)
 gccgccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTg                                                                                                                                                   <  1:1150813/150‑1 (MQ=255)
  ccgccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGc                                                                                                                                                  >  1:204861/1‑150 (MQ=255)
  ccgccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGc                                                                                                                                                  <  1:1391891/150‑1 (MQ=255)
  ccgccATGCGATCCTGCGATTTACCCATTTTGGAAGCGGCGTCCATCTCCAGCGCTTTCGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTTTCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGc                                                                                                                                                  <  2:231114/150‑1 (MQ=255)
  ccgccATGCGATCCTGCGATACACGCATTTGGGAAGCGGCGACCATCTCCATCGCTATATGGATCTTTTGCCTGGTCTGGACCCTTGCGATCTTACTACGTATCTCTTATGCGCCGGCCATGAGCTTCACCTCAATGCCTTGCGGCCTGc                                                                                                                                                  <  2:264679/150‑1 (MQ=255)
    gccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGccc                                                                                                                                                <  2:1140936/150‑1 (MQ=255)
    gccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGccc                                                                                                                                                <  2:1013391/150‑1 (MQ=255)
     ccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCt                                                                                                                                               >  1:1351688/1‑150 (MQ=255)
        tGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAg                                                                                                                                            <  1:1307536/150‑1 (MQ=255)
         gCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCGTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAgg                                                                                                                                           >  1:1681208/1‑150 (MQ=255)
          cGATCCTGCGATTTACGCTTTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGTTCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCGGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTATGGc                                                                                                                                          <  1:615470/150‑1 (MQ=255)
          cGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGc                                                                                                                                          <  1:1036446/150‑1 (MQ=255)
                tGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAgccg                                                                                                                                    >  1:1483256/1‑150 (MQ=255)
                       tACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACg                                                                                                                             >  1:83356/1‑150 (MQ=255)
                        aCGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAggcat                                                                                                                            <  2:410044/150‑5 (MQ=255)
                         cGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCAGCCCTAAAGCAAGCCGCCCGACGtt                                                                                                                           >  1:1462271/1‑150 (MQ=255)
                            aTTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTAcc                                                                                                                        <  2:1682383/150‑1 (MQ=255)
                               ttGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCGTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTAGCGGCCTGCCCTAAGGCACGCCGCCAGCCGTTACCaag                                                                                                                     >  2:85253/1‑148 (MQ=255)
                                 ggAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGAt                                                                                                                   >  2:1482035/1‑150 (MQ=255)
                                 ggAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGACGCCAGACGTTACCAGGAt                                                                                                                   >  2:568504/1‑150 (MQ=255)
                                 ggAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTAGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGAt                                                                                                                   >  1:895310/1‑150 (MQ=255)
                                  gAAGCGGCGACCATCTGCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGGTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGACCTTAGGCAAGCCGCCAGACGTTACCAGGAtt                                                                                                                  <  1:1053893/150‑1 (MQ=255)
                                    aGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTgg                                                                                                                >  2:355452/1‑150 (MQ=255)
                                     gcggcgACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTTCCAGGATTggg                                                                                                               <  1:941349/150‑1 (MQ=255)
                                      cggcgACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGt                                                                                                              >  1:207134/1‑150 (MQ=255)
                                         cgACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGc                                                                                                           >  2:716705/1‑150 (MQ=255)
                                            ccATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCttt                                                                                                        >  1:1486729/1‑150 (MQ=255)
                                                ctcCATCGCTTTAGTGCTCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAg                                                                                                    <  2:651514/150‑1 (MQ=255)
                                                ctcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGTCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAg                                                                                                    >  2:1140151/1‑150 (MQ=255)
                                                ctcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGCGCTTCTCCTCAATGCCTTGCGGACTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGAATGGGTTGCTTTGAAg                                                                                                    >  1:1055319/1‑150 (MQ=255)
                                                 tcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCCCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATAGGGTTTCTTTGAAgg                                                                                                   >  2:510011/1‑150 (MQ=255)
                                                        cTTTAGGGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAg                                                                                            <  1:1039976/150‑1 (MQ=255)
                                                         tttAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAgg                                                                                           <  1:1222968/150‑1 (MQ=255)
                                                              tgttCTTTTGCGTGTTCTGGAGGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGcc                                                                                      <  2:831324/147‑1 (MQ=255)
                                                              tGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGcc                                                                                      <  1:1596459/150‑1 (MQ=255)
                                                               gATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCt                                                                                     >  2:1318604/1‑150 (MQ=255)
                                                               gATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGCCAAGCCGCAAGACCTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCt                                                                                     >  2:622746/1‑150 (MQ=255)
                                                                aTCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTTCCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCtt                                                                                    >  2:754751/1‑150 (MQ=255)
                                                                  cTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTc                                                                                  <  2:586548/150‑1 (MQ=255)
                                                                  cTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTc                                                                                  <  2:384879/150‑1 (MQ=255)
                                                                   ttttGCGTGGGCTGGACTCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCa                                                                                 <  2:118677/150‑1 (MQ=255)
                                                                      tGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCt                                                                              >  1:1340852/1‑150 (MQ=255)
                                                                       gCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCtt                                                                             >  1:946665/1‑150 (MQ=255)
                                                                        cGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTg                                                                            <  2:41020/150‑1 (MQ=255)
                                                                          tgtTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGcc                                                                          <  1:1464476/150‑1 (MQ=255)
                                                                          tgtTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGcc                                                                          <  2:1350280/150‑1 (MQ=255)
                                                                             tCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTc                                                                       <  1:884105/150‑1 (MQ=255)
                                                                               tGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTACGGCCAGCCGCCCGCCGGTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCCGCGTGccgtggg                                                                     >  2:1712414/1‑145 (MQ=255)
                                                                               tGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGa                                                                     <  1:243948/150‑1 (MQ=255)
                                                                                   cGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTACTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtc                                                                 <  2:1499474/150‑1 (MQ=255)
                                                                                       tGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcgtcg                                                             <  2:1564283/150‑1 (MQ=255)
                                                                                            tCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTa                                                        <  1:885016/150‑1 (MQ=255)
                                                                                             cTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAg                                                       <  2:929643/150‑1 (MQ=255)
                                                                                              ttACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGc                                                      >  1:757579/1‑150 (MQ=255)
                                                                                                 ctacGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGCCGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCCGCTTGCCTTCGATTTCGTCGTTGTAGCCAc                                                   >  1:536440/1‑150 (MQ=255)
                                                                                                 ctacGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTAGTTGTAGCCAc                                                   >  1:1561835/1‑150 (MQ=255)
                                                                                                  tacGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCAcc                                                  <  2:1126217/150‑1 (MQ=255)
                                                                                                      tATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCTCCGGTc                                              >  2:430886/1‑150 (MQ=255)
                                                                                                       aTCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCt                                             >  2:1619691/1‑150 (MQ=255)
                                                                                                       aTCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCt                                             >  2:1151851/1‑150 (MQ=255)
                                                                                                       aTCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCAGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTTGGTTGCTTTTAAGGAATCGAGGATGCCTTTGAGCTTGCCGTAGTTTTAGTCGTTGTACCCAACGGTCt                                             >  2:1516373/1‑150 (MQ=255)
                                                                                                         ctctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTgg                                           <  1:1125138/150‑1 (MQ=255)
                                                                                                          tctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCGTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGt                                          <  1:203232/150‑1 (MQ=255)
                                                                                                           ctTTTGCGCCGGCCATGAGCTTCCCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGtt                                         <  1:727845/150‑1 (MQ=255)
                                                                                                            ttttGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTg                                        <  2:1605436/150‑1 (MQ=255)
                                                                                                             tttGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGa                                       >  1:729867/1‑150 (MQ=255)
                                                                                                               tGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGAtc                                     <  1:788623/150‑1 (MQ=255)
                                                                                                               tGCGCCGGCCATGAGCGTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTTCCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGAtc                                     <  2:1140518/150‑1 (MQ=255)
                                                                                                                 cgcCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGAtctc                                   <  2:604820/150‑1 (MQ=255)
                                                                                                                   ccGGCCATGAGCTTCTCCTCATTGCGTTGCGGCCTGCCCTATGGCAAGCCGCCAGACGTTACCAGGCTTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCtt                                 <  2:1055319/150‑1 (MQ=255)
                                                                                                                    cGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTg                                <  2:946665/150‑1 (MQ=255)
                                                                                                                       ccATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTATGAAGGAATCGAGGATGCCTATCAGCTAGCCTTCGATTTCGTCGTTGTAGACACCGGTCTGGTTGCTCTCTTGCAt                             >  1:264300/1‑150 (MQ=255)
                                                                                                                        cATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATc                            >  1:1037293/1‑150 (MQ=255)
                                                                                                                        cATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATc                            >  2:395190/1‑150 (MQ=255)
                                                                                                                         aTGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGAAGTTACCAGGATTGGGTTGCTTTGAAGGAATAGAGGATGCCTTTCATATTGCCTTCGATTTAGTCGTTGTAGCCACCGGTCAGGTTGATCTCTTtcatca                           >  1:1712754/1‑150 (MQ=255)
                                                                                                                          tGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCGTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTACGATTTCGTCGTTGTAGCCAGCGGTCTGGTTGATCTCTTGGATaaa                          >  1:315932/1‑150 (MQ=255)
                                                                                                                               ttCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAg                     <  2:808464/150‑1 (MQ=255)
                                                                                                                                tctcCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGc                    >  1:123883/1‑150 (MQ=255)
                                                                                                                                   cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGACTCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCCGGTTGATCTCTTGCATCAACGGAGCGTg                 <  1:412098/150‑1 (MQ=255)
                                                                                                                                   cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGTGGATGCCTTTCAGCTTGCCTTCGTTTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTg                 <  2:895310/150‑1 (MQ=255)
                                                                                                                                   cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTg                 >  1:368308/1‑150 (MQ=255)
                                                                                                                                   cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCGTTGAAGGAATCGAGGATGCCTATCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGGCTGGTTGATCTCTTGCATCAACGGAGCGTg                 <  2:1237460/150‑1 (MQ=255)
                                                                                                                                     tcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGAt               <  1:971941/150‑1 (MQ=255)
                                                                                                                                      cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGCGGATGCCTTTCAGCTTGCCTGCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGAGCTCTTGAATCAACGGAGCGTGATc              >  2:1132639/1‑150 (MQ=255)
                                                                                                                                      cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATc              <  2:43533/150‑1 (MQ=255)
                                                                                                                                      cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATc              <  1:54628/150‑1 (MQ=255)
                                                                                                                                      cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATc              >  1:416288/1‑150 (MQ=255)
                                                                                                                                       aaTGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGCCGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATTCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCa             >  2:720559/1‑150 (MQ=255)
                                                                                                                                       aaTGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCa             <  2:1049397/150‑1 (MQ=255)
                                                                                                                                          gCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACgg          >  2:370264/1‑150 (MQ=255)
                                                                                                                                          gCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACgg          <  1:478686/150‑1 (MQ=255)
                                                                                                                                           ccTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGt         <  2:83356/150‑1 (MQ=255)
                                                                                                                                           ccTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGt         <  2:868819/150‑1 (MQ=255)
                                                                                                                                             tggCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGATGATCTCTTGCATCAACGGAGCGTGATCACGGTCg       >  1:697803/3‑150 (MQ=255)
                                                                                                                                              ggCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCgc      >  1:614309/2‑149 (MQ=255)
                                                                                                                                               gCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGAc     >  1:1380955/1‑150 (MQ=255)
                                                                                                                                                 ggCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCCCCGGTCTGGTTGATCTCTTGCTTCAACGGAGCGTGATCACGGTCGACGt   >  1:1155342/1‑150 (MQ=255)
                                                                                                                                                 ggCTTGCCTTAGGGCAGGCCGCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTa  >  2:193665/1‑150 (MQ=255)
                                                                                                                                                  gCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTa  >  1:1179304/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCCGCCATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTA  >  NC_000913/3918143‑3918438

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