Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,918,299 A→G 6.0% intergenic (‑34/+17) 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,2990AG 6.0% 397.0 / 11.0 135intergenic (‑34/+17)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 A (63/61);  new base G (6/2);  total (70/65)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

Only 100 of 138 total aligned reads displayed.

TGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCCAGCAGA  >  NC_000913/3918151‑3918448
                                                                                                                                                    |                                                                                                                                                     
 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)
          cGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAgccgcc                                                                                                                                            <  1:920460/150‑1 (MQ=255)
            aTTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAg                                                                                                                                          <  1:1307729/150‑1 (MQ=255)
               tACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACg                                                                                                                                       >  1:83356/1‑150 (MQ=255)
               tACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACg                                                                                                                                       >  1:868819/1‑150 (MQ=255)
                aCGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGGCCGCAAggcat                                                                                                                                      <  2:410044/150‑5 (MQ=255)
                    aTTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTAcc                                                                                                                                  <  2:1682383/150‑1 (MQ=255)
                       ttGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCGTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTAGCGGCCTGCCCTAAGGCACGCCGCCAGCCGTTACCaag                                                                                                                               >  2:85253/1‑148 (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)
                               ggcgACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGtt                                                                                                                       >  2:389130/1‑150 (MQ=255)
                                 cgACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGc                                                                                                                     >  2:716705/1‑150 (MQ=255)
                                 cgACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGc                                                                                                                     <  1:568504/150‑1 (MQ=255)
                                        ctcCATCGCTTTAGTGCTCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAg                                                                                                              <  2:651514/150‑1 (MQ=255)
                                        ctcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGTCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAg                                                                                                              >  2:1140151/1‑150 (MQ=255)
                                        ctcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGCGCTTCTCCTCAATGCCTTGCGGACTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGAATGGGTTGCTTTGAAg                                                                                                              >  1:1055319/1‑150 (MQ=255)
                                        ctcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAg                                                                                                              <  1:146341/150‑1 (MQ=255)
                                         tcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCCCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATAGGGTTTCTTTGAAgg                                                                                                             >  2:510011/1‑150 (MQ=255)
                                           cATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCCGCCTGCACTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGATTTGAAGGaa                                                                                                           >  2:1547943/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)
                                                       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)
                                                             ttGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGc                                                                                         <  2:198686/150‑1 (MQ=255)
                                                              tGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCt                                                                                        >  1:1340852/1‑150 (MQ=255)
                                                               gCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCtt                                                                                       >  1:946665/1‑150 (MQ=255)
                                                                cGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGTGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTg                                                                                      <  2:1412072/150‑1 (MQ=255)
                                                                cGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTg                                                                                      <  2:41020/150‑1 (MQ=255)
                                                                  tgtTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGcc                                                                                    <  2:1350280/150‑1 (MQ=255)
                                                                  tgtTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGcc                                                                                    <  1:1464476/150‑1 (MQ=255)
                                                                    ttCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCtt                                                                                  <  1:1206616/150‑1 (MQ=255)
                                                                       tGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTACGGCCAGCCGCCCGCCGGTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCCGCGTGccgtggg                                                                               >  2:1712414/1‑145 (MQ=255)
                                                                       tGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGa                                                                               <  1:243948/150‑1 (MQ=255)
                                                                           cgcgTGCGATCTTACTAAGTATCTCTTTTGCGCCGGCCAAGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtc                                                                           <  1:323657/150‑1 (MQ=255)
                                                                           cGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTACTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtc                                                                           <  2:1499474/150‑1 (MQ=255)
                                                                                   aTCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGt                                                                   <  1:752024/150‑1 (MQ=255)
                                                                                   aTCTTACTACGTATCGCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGt                                                                   <  1:1721024/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)
                                                                                         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)
                                                                                                  tctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCGTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGt                                                    <  1:203232/150‑1 (MQ=255)
                                                                                                   ctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTTGTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGtt                                                   >  1:1557016/1‑150 (MQ=255)
                                                                                                    ttttGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTg                                                  <  2:1605436/150‑1 (MQ=255)
                                                                                                     tttGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGa                                                 >  1:729867/1‑150 (MQ=255)
                                                                                                        gcgcCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGAtct                                              >  2:26719/1‑150 (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:1050484/1‑150 (MQ=255)
                                                                                                                cATGAGCTGCTCCTCATTGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGGTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATc                                      <  2:1147222/150‑1 (MQ=255)
                                                                                                                  tGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCaa                                    >  1:43533/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)
                                                                                                                           cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGTGGATGCCTTTCAGCTTGCCTTCGTTTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTg                           <  2:895310/150‑1 (MQ=255)
                                                                                                                           cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTg                           >  1:368308/1‑150 (MQ=255)
                                                                                                                           cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTg                           >  1:695811/1‑150 (MQ=255)
                                                                                                                           cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCGTTGAAGGAATCGAGGATGCCTATCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGGCTGGTTGATCTCTTGCATCAACGGAGCGTg                           <  2:1237460/150‑1 (MQ=255)
                                                                                                                            ctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGa                          >  2:1259912/1‑150 (MQ=255)
                                                                                                                             tcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGAt                         <  1:971941/150‑1 (MQ=255)
                                                                                                                              cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATc                        <  2:43533/150‑1 (MQ=255)
                                                                                                                              cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATc                        >  2:609878/1‑150 (MQ=255)
                                                                                                                              cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATc                        >  1:416288/1‑150 (MQ=255)
                                                                                                                              cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATc                        <  1:54628/150‑1 (MQ=255)
                                                                                                                              cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCATCGAGTTCGTCGTTGTAGCCCGCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATc                        >  2:1138459/1‑150 (MQ=255)
                                                                                                                               aaTGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGCCGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATTCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCa                       >  2:720559/1‑150 (MQ=255)
                                                                                                                                  gCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACgg                    <  1:478686/150‑1 (MQ=255)
                                                                                                                                   ccTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGt                   <  2:868819/150‑1 (MQ=255)
                                                                                                                                    cTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTc                  <  1:1712414/150‑1 (MQ=255)
                                                                                                                                     ttGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCg                 <  1:1140151/150‑1 (MQ=255)
                                                                                                                                     tggCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGATGATCTCTTGCATCAACGGAGCGTGATCACGGTCg                 >  1:697803/3‑150 (MQ=255)
                                                                                                                                       gCGGCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGAc               >  1:1380955/1‑150 (MQ=255)
                                                                                                                                         ggCTTGCCTTAGGGCAGGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCCCCGGTCTGGTTGATCTCTTGCTTCAACGGAGCGTGATCACGGTCGACGt             >  1:1155342/1‑150 (MQ=255)
                                                                                                                                              gCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGcc        <  1:1050484/150‑1 (MQ=255)
                                                                                                                                              gCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGGCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGcc        <  1:1547943/150‑1 (MQ=255)
                                                                                                                                               cccTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCca       >  2:1000084/1‑150 (MQ=255)
                                                                                                                                               cccTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCca       >  2:397093/1‑150 (MQ=255)
                                                                                                                                                ccTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCcag      >  1:1415787/1‑150 (MQ=255)
                                                                                                                                                 cTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCcagc     <  1:1318604/150‑1 (MQ=255)
                                                                                                                                                 cTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCcagc     <  1:992914/150‑1 (MQ=255)
                                                                                                                                                   aaGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCAGCAACGGCGCGTGATCACGGTCGACGTACGCcaggag   >  2:114588/1‑150 (MQ=255)
                                                                                                                                                    aGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCACCGGAGCGTGATCACGGTCGACGTAAGCCAGCaga  >  2:158149/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCCAGCAGA  >  NC_000913/3918151‑3918448

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: