Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,918,296 C→T 7.7% intergenic (‑31/+20) 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,2960CT 7.7% 362.0 / 11.4 118intergenic (‑31/+20)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 (58/50);  new base T (3/6);  total (61/57)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 126 total aligned reads displayed.

GCCATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCCAGC  >  NC_000913/3918147‑3918445
                                                                                                                                                     |                                                                                                                                                     
gccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGccc                                                                                                                                                       <  1:1043348/150‑1 (MQ=255)
gccATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGccc                                                                                                                                                       <  2:809447/150‑1 (MQ=255)
  cATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTa                                                                                                                                                     >  1:717989/1‑150 (MQ=255)
  cATGCGATCCTGCGATTTACGCAGTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTGCTGGACGCTGGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTa                                                                                                                                                     <  1:303113/150‑1 (MQ=255)
    tGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAg                                                                                                                                                   >  1:1589992/1‑150 (MQ=255)
      cGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGc                                                                                                                                                 <  2:223568/150‑1 (MQ=255)
      cGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGc                                                                                                                                                 <  2:1394324/150‑1 (MQ=255)
      cGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGc                                                                                                                                                 >  1:209685/1‑150 (MQ=255)
        aTCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTAAGTGATCTTTTGCGGGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTAGGGCAGGCCGCAAGGCat                                                                                                                                               <  1:1259258/150‑2 (MQ=255)
            tGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGgccg                                                                                                                                           <  2:928/150‑1 (MQ=255)
                aTTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAg                                                                                                                                       <  1:1236971/150‑1 (MQ=255)
                 tttACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAgg                                                                                                                                      <  2:503443/150‑2 (MQ=255)
                        aTTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTAcc                                                                                                                               >  2:285540/1‑150 (MQ=255)
                         ttttGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCa                                                                                                                              >  2:789098/1‑150 (MQ=255)
                            tGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGa                                                                                                                           >  2:1462423/1‑150 (MQ=255)
                              gAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGAtt                                                                                                                         >  2:1566653/1‑150 (MQ=255)
                                  cggcgACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGt                                                                                                                     <  2:877859/150‑1 (MQ=255)
                                   ggcgACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGtt                                                                                                                    >  2:200560/1‑150 (MQ=255)
                                        ccATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGAAGTTACCAGGTTTGGGTTGCttt                                                                                                               >  2:1591174/1‑150 (MQ=255)
                                          aTCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGAGCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGa                                                                                                             <  1:1159243/150‑1 (MQ=255)
                                            ctcCATCGCTTTAGTGATCTTTTGCGGGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAg                                                                                                           <  2:782703/150‑1 (MQ=255)
                                             tcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAgg                                                                                                          >  1:1534166/1‑150 (MQ=255)
                                             tcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAgg                                                                                                          <  1:94629/150‑1 (MQ=255)
                                              ccATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGa                                                                                                         >  2:1104933/1‑150 (MQ=255)
                                                 tCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATc                                                                                                      <  1:672534/150‑1 (MQ=255)
                                                 tCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATc                                                                                                      <  1:1441318/150‑1 (MQ=255)
                                                 tCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATc                                                                                                      <  2:54655/150‑1 (MQ=255)
                                                   gCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGa                                                                                                    >  1:1275886/1‑150 (MQ=255)
                                                   gCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGa                                                                                                    >  2:589860/1‑150 (MQ=255)
                                                    cTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAg                                                                                                   >  2:1546939/1‑150 (MQ=255)
                                                       tAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGAt                                                                                                >  1:799122/1‑150 (MQ=255)
                                                        aGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGACAGACGTTACCAGGAATGTGTTGCTTTTAATGAATCGatgatg                                                                                               >  1:1299291/1‑150 (MQ=255)
                                                          tgctCTTTTGCGTGTTCTGGACGCTTGCGATATTACTACGTATCTATATTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGTAATCGAGGATGcc                                                                                             <  2:1609204/147‑1 (MQ=255)
                                                          tGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGcc                                                                                             >  1:757864/1‑150 (MQ=255)
                                                           gATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCt                                                                                            >  2:486934/1‑150 (MQ=255)
                                                              cTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAACCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAACGAATCGAGCAAGCCttta                                                                                         >  1:1426857/1‑149 (MQ=255)
                                                                tttGCGTGTGCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAg                                                                                       <  2:790764/150‑1 (MQ=255)
                                                                 ttGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGCGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGc                                                                                      >  1:999616/1‑150 (MQ=255)
                                                                 ttGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGc                                                                                      <  1:758442/150‑1 (MQ=255)
                                                                       gtTCTGGACGCTTGCGATCTTACTTCGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCt                                                                                <  2:302035/150‑1 (MQ=255)
                                                                       gtTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAAGCGAGGATGCCTTTCAGCTTGCCt                                                                                >  2:679903/1‑150 (MQ=255)
                                                                        ttCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCtt                                                                               <  1:1177773/150‑1 (MQ=255)
                                                                         gcTGGACGCTTGCGATCTTACTAAGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTc                                                                              <  1:251879/149‑1 (MQ=255)
                                                                          cTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCg                                                                             <  2:1389237/150‑1 (MQ=255)
                                                                          cTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCg                                                                             <  2:1357052/150‑1 (MQ=255)
                                                                            ggACGCTTGCGATCTTACTACGTATCTCTTATGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCAGCCCAAAGGCACGCCCCCACACGTAACCAGGAATGGGTTGCTTTTAACGAATCGACGATGCCTTTCATCTAGCCTTCGAt                                                                           >  1:1492365/1‑150 (MQ=255)
                                                                                 cTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcgt                                                                      >  1:398431/1‑150 (MQ=255)
                                                                                   tGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcgtcg                                                                    <  2:718181/150‑1 (MQ=255)
                                                                                    gCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTcgtcgt                                                                   >  1:1234649/1‑150 (MQ=255)
                                                                                     cGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCCATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGCCGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGtt                                                                  >  1:819896/1‑150 (MQ=255)
                                                                                         cTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCTGGATTGGGTTGCTTTGAAGGAATCGAGGATGTCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAg                                                              <  1:225933/150‑1 (MQ=255)
                                                                                           tactacGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGcc                                                            >  2:762154/1‑150 (MQ=255)
                                                                                             ctacGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATTCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGTTTTAGTCGTTGTAGcccc                                                          >  2:326966/1‑148 (MQ=255)
                                                                                              tacGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCAcc                                                         <  2:519024/150‑1 (MQ=255)
                                                                                                  tATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGAt                                                                           >  1:1220849/1‑128 (MQ=255)
                                                                                                  tATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGAt                                                                           <  2:1220849/128‑1 (MQ=255)
                                                                                                     ctctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGTTCTgg                                                  >  2:220866/1‑150 (MQ=255)
                                                                                                     ctctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTgg                                                  >  2:1612533/1‑150 (MQ=255)
                                                                                                      tctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGt                                                 <  1:1082210/150‑1 (MQ=255)
                                                                                                      tctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCCGCATGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCGTTAAGCGTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCGGGt                                                 >  2:1086232/1‑150 (MQ=255)
                                                                                                        ttttGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCAcc                                                         <  2:1483434/140‑1 (MQ=255)
                                                                                                           tGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGAtc                                            <  2:430251/150‑1 (MQ=255)
                                                                                                               ccGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCtt                                        <  2:1087048/150‑1 (MQ=255)
                                                                                                                cGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTg                                       <  2:1193838/150‑1 (MQ=255)
                                                                                                                cGGCCATGAGCATCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTGGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTg                                       <  1:1055464/150‑1 (MQ=255)
                                                                                                                 ggCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGt                                                 <  1:512679/139‑1 (MQ=255)
                                                                                                                 ggCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGt                                                 >  2:512679/1‑139 (MQ=255)
                                                                                                                 ggCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCAAAGGCAAGCCGCCAGACGTTCCCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGGCTGGTTGATCTCTTGc                                      <  2:903346/150‑1 (MQ=255)
                                                                                                                  gCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAATGAATCGAGGATGCCTTTCAGCTTGACTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCa                                     >  2:1020300/1‑150 (MQ=255)
                                                                                                                  gCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCa                                     <  1:195555/150‑1 (MQ=255)
                                                                                                                      tGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCaa                                 >  1:789045/1‑150 (MQ=255)
                                                                                                                      tGAGCTTCTCCTCAATGCCTGGCGGCCGGCCCTAAGGCAAGCCGCCCGACGGTACCAGGATTGGGTTGCTTTGAAGGATTCGAGGATGCCTTTCATCTTGCCTTCGATGTCGTCGTTGTAGCCCCCGGTAGTGTTGATCTCTTgaagcca                                 >  2:1311841/1‑144 (MQ=255)
                                                                                                                        aGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACg                               >  2:355815/1‑150 (MQ=255)
                                                                                                                         gCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTGGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACgg                              >  2:1632169/1‑150 (MQ=255)
                                                                                                                           ttCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAg                            >  2:1192344/1‑150 (MQ=255)
                                                                                                                              tcctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGAGGTTACCAGGATTGGGTTGCTTTGAAGGAAGCGAGGATGCCTTTCCGCTTGCCTTCGATTTCGTCGTTGTCGCCTCCGGTCTGGTTGATCTCTTGCATCAACGGAGCGt                         >  1:251953/1‑150 (MQ=255)
                                                                                                                               cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGTTTGGGTTGCTTTGAAGGAATCGAGGATGCCATTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTg                        >  1:505123/1‑150 (MQ=255)
                                                                                                                                ctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGa                       >  2:156395/1‑150 (MQ=255)
                                                                                                                                ctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGa                       <  1:1034364/150‑1 (MQ=255)
                                                                                                                                ctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCGTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGa                       >  1:702941/1‑150 (MQ=255)
                                                                                                                                 tcAATGCGTTGCGGCCTGCGCTAAGGCAAGCCGCCAGACGTTACGTGGTATGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCAGCGGGCTGGTTGATCTCTTGCATCAACGGAGCGTGAt                      <  2:465507/150‑1 (MQ=255)
                                                                                                                                  cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATc                     <  2:154277/150‑1 (MQ=255)
                                                                                                                                  cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATc                     <  1:1462423/150‑1 (MQ=255)
                                                                                                                                  cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATc                     >  2:650215/1‑150 (MQ=255)
                                                                                                                                  cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGGCTGGTTGATCTCTTGCATCAACGGAGCGTGATc                     <  2:1018081/150‑1 (MQ=255)
                                                                                                                                   aaTGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTTGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCa                    >  2:143279/1‑150 (MQ=255)
                                                                                                                                   aaTGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCa                    <  2:1275886/150‑1 (MQ=255)
                                                                                                                                    aTGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCAc                   <  1:200560/150‑1 (MQ=255)
                                                                                                                                       ccTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAAGCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGt                <  1:1432671/150‑1 (MQ=255)
                                                                                                                                        ctggCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTc               >  1:556596/4‑150 (MQ=255)
                                                                                                                                        cTTGGGGCCTGACCGAAGGAAATACGCCAGAAGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGACTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACAGGTCTGGTTGATCTCTTGCATAAACGGAGCGTGATCACGGTc               <  2:1492365/150‑1 (MQ=255)
                                                                                                                                          ggCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTCTAGCCACCGGTCTGGTTGCTCTCTTCCATCAACGGAGCGTGATCAGGGTCGa             >  2:1077038/2‑150 (MQ=255)
                                                                                                                                            cGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACg           >  1:1218476/1‑150 (MQ=255)
                                                                                                                                              gCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTa         >  1:1059959/1‑150 (MQ=255)
                                                                                                                                                cTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAg       <  1:326966/150‑1 (MQ=255)
                                                                                                                                                   cccTAAGGCAACCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTACGATTTCCTCGTTGAAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATAACGGTCGACGTAAGCca    >  1:1264762/1‑150 (MQ=255)
                                                                                                                                                    ccTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCcag   >  1:1174327/1‑150 (MQ=255)
                                                                                                                                                     cTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCcagc  <  1:1546939/150‑1 (MQ=255)
                                                                                                                                                     cTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCcagc  >  1:366732/1‑150 (MQ=255)
                                                                                                                                                     cTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCcagc  <  1:455317/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCCAGC  >  NC_000913/3918147‑3918445

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: