Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,920,796 T→C 6.5% A31A (GCA→GCG atpF ← ATP synthase Fo complex subunit b

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,920,7960TC 6.5% 372.9 / 10.2 122A31A (GCA→GCGatpFATP synthase Fo complex subunit b
Reads supporting (aligned to +/- strand):  ref base T (49/65);  new base C (3/5);  total (52/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

TTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCA  >  NC_000913/3920647‑3920942
                                                                                                                                                     |                                                                                                                                                  
    tCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCa                                                                                                                                                <  1:1531847/150‑1 (MQ=255)
       ccTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTa                                                                                                                                             >  2:1440165/1‑150 (MQ=255)
        cTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTaa                                                                                                                                            >  2:1635595/1‑150 (MQ=255)
        cTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTaa                                                                                                                                            <  2:122771/150‑1 (MQ=255)
         tGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAAt                                                                                                                                           <  1:1093860/150‑1 (MQ=255)
            tCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATggc                                                                                                                                        >  1:1176634/1‑150 (MQ=255)
               atgatTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATggcggc                                                                                                                                     >  2:589009/1‑150 (MQ=255)
                tgatTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGcc                                                                                                                                    <  2:1001151/150‑1 (MQ=255)
                   ttACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATGGCAGCCAtc                                                                                                                                 <  2:1158172/150‑2 (MQ=255)
                      ccTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGt                                                                                                                              <  2:244339/150‑1 (MQ=255)
                       cTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTa                                                                                                                             >  2:694685/1‑150 (MQ=255)
                        tGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTAc                                                                                                                            >  1:903076/1‑150 (MQ=255)
                         gggCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGGAAGGTCAATGTCCTTATGTGCTCGTTCTGCGGAATCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGTCTGCCATTAATGGCGGCCATACGTACt                                                                                                                           <  2:1599214/150‑1 (MQ=255)
                          ggCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACtt                                                                                                                          >  2:384948/1‑150 (MQ=255)
                          ggCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACtt                                                                                                                          >  1:1451227/1‑150 (MQ=255)
                          ggCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACtt                                                                                                                          <  1:1364060/150‑1 (MQ=255)
                            cTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCa                                                                                                                        <  2:1156306/150‑1 (MQ=255)
                            cTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCa                                                                                                                        <  1:585203/150‑1 (MQ=255)
                             ttCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCAt                                                                                                                       <  1:1335345/150‑1 (MQ=255)
                               ccGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGcc                                                                                                                                    <  2:1573733/135‑1 (MQ=255)
                               ccGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGcc                                                                                                                                    >  1:1573733/1‑135 (MQ=255)
                                cgctttcgctttTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGca                                                                                                                    <  2:1451227/150‑1 (MQ=255)
                                cgctttcgctttTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGca                                                                                                                    <  1:368870/150‑1 (MQ=255)
                                 gctttcgctttTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcag                                                                                                                   <  2:687480/150‑1 (MQ=255)
                                  ctttcgctttTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcaga                                                                                                                  <  1:694685/150‑1 (MQ=255)
                                   tttcgctttTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcagaa                                                                                                                 >  1:881849/1‑150 (MQ=255)
                                    ttcgctttTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcagaac                                                                                                                >  2:1336008/1‑150 (MQ=255)
                                     tcgctttTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcagaaca                                                                                                               >  1:1489892/1‑150 (MQ=255)
                                      cgctttTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATGGCAGCCa                                                                                                                                   >  1:341804/1‑129 (MQ=255)
                                       gctttTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcaga                                                                                                                  >  2:559739/1‑145 (MQ=255)
                                       gctttTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGcaga                                                                                                                  <  1:559739/145‑1 (MQ=255)
                                       gctttTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATGGCAGCCa                                                                                                                                   >  2:42502/1‑128 (MQ=255)
                                       gctttTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCCGCCATTAATGGCAGCCa                                                                                                                                   <  1:42502/128‑1 (MQ=255)
                                          tttttCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacg                                                                                                          >  2:157280/1‑150 (MQ=255)
                                           ttttCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacga                                                                                                         >  2:279352/1‑150 (MQ=255)
                                             ttCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAgaacgaac                                                                                                       <  2:126648/150‑1 (MQ=255)
                                              tCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACa                                                                                                      <  2:219422/150‑1 (MQ=255)
                                                  cTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGAc                                                                                                  >  1:304102/1‑150 (MQ=255)
                                                   tGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACa                                                                                                 <  2:700502/150‑1 (MQ=255)
                                                     gtcggtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCagaacagaac                                                                                                           >  2:343991/1‑138 (MQ=255)
                                                      tcggtcgCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAAc                                                                                              <  2:6459/150‑1 (MQ=255)
                                                            gcgcTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtg                                                                                        >  1:1221113/1‑150 (MQ=255)
                                                             cgcTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtgg                                                                                       >  2:433397/1‑150 (MQ=255)
                                                             cgcTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtgg                                                                                       <  1:97719/150‑1 (MQ=255)
                                                               cTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcc                                                                                     <  1:4358/150‑1 (MQ=255)
                                                                   ccTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcctggc                                                                                 >  1:951765/1‑150 (MQ=255)
                                                                   ccTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcctggc                                                                                 >  2:889627/1‑150 (MQ=255)
                                                                    cTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGAtggcctggcc                                                                                <  1:1498286/150‑1 (MQ=255)
                                                                       tGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAg                                                                             <  2:1220862/150‑1 (MQ=255)
                                                                       tGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAg                                                                             <  1:832227/150‑1 (MQ=255)
                                                                        gCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAgg                                                                            >  1:1487900/1‑150 (MQ=255)
                                                                         caaggtcaaggtCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGa                                                                           >  1:215996/1‑150 (MQ=255)
                                                                            ggtcaaggtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttg                                                                        >  1:406548/1‑150 (MQ=255)
                                                                              tcaaggtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGAttgtt                                                                      >  2:917538/1‑150 (MQ=255)
                                                                                 aggtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCg                                                                   <  1:1635595/150‑1 (MQ=255)
                                                                                 aggtcCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCg                                                                   <  1:1477841/150‑1 (MQ=255)
                                                                                      cTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAg                                                              <  2:1558481/150‑1 (MQ=255)
                                                                                       ttATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGa                                                             <  1:298628/150‑1 (MQ=255)
                                                                                         aTGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGAtt                                                           <  1:889627/150‑1 (MQ=255)
                                                                                          tgtgCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTc                                                          <  1:814501/150‑1 (MQ=255)
                                                                                           gtgCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTca                                                         <  2:951765/150‑1 (MQ=255)
                                                                                           gtgCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTca                                                         >  2:1463825/1‑150 (MQ=255)
                                                                                            tgCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTcac                                                        <  2:1433982/150‑1 (MQ=255)
                                                                                            tgCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTcac                                                        <  2:1253125/150‑1 (MQ=255)
                                                                                             gCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTcaca                                                       >  2:479405/1‑150 (MQ=255)
                                                                                               tCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGc                                                     <  1:917538/150‑1 (MQ=255)
                                                                                                   tCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACaa                                                 >  1:955828/1‑150 (MQ=255)
                                                                                                   tCTGCGGAAGCAAGGCCGTCACCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGACTGCCCGAGGATTGTTGCGTTAAGATTCACAGCACaa                                                 >  2:232003/1‑150 (MQ=255)
                                                                                                    cTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAAt                                                <  1:1304374/150‑1 (MQ=255)
                                                                                                     tGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATg                                               >  1:1259313/1‑150 (MQ=255)
                                                                                                     tGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATg                                               >  1:586240/1‑150 (MQ=255)
                                                                                                      gCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGc                                              >  2:1336448/1‑150 (MQ=255)
                                                                                                      gCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGc                                              >  2:35413/1‑150 (MQ=255)
                                                                                                        ggAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCct                                            >  1:243601/1‑150 (MQ=255)
                                                                                                        ggAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCct                                            >  1:244906/1‑150 (MQ=255)
                                                                                                        ggAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCct                                            >  2:142497/1‑150 (MQ=255)
                                                                                                         gAAGCAAGGCCGTCAGCAATTTCTTTTTGACGGTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCctc                                           <  2:685465/150‑1 (MQ=255)
                                                                                                          aaGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCctct                                          <  2:748522/150‑1 (MQ=255)
                                                                                                           aGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTa                                         <  2:304102/150‑1 (MQ=255)
                                                                                                              aaGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTAttt                                      <  2:406548/150‑1 (MQ=255)
                                                                                                               aGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTa                                     <  1:279352/150‑1 (MQ=255)
                                                                                                                ggCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAg                                    >  1:887372/1‑150 (MQ=255)
                                                                                                                   cGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTa                                 <  1:1084742/150‑1 (MQ=255)
                                                                                                                     tCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAAc                               <  1:157280/150‑1 (MQ=255)
                                                                                                                     tCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAAc                               >  2:1098800/1‑150 (MQ=255)
                                                                                                                      cAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACg                              <  1:125328/150‑1 (MQ=255)
                                                                                                                          aaTTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCt                          <  2:881849/150‑1 (MQ=255)
                                                                                                                           aTTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTg                         <  2:1462289/150‑1 (MQ=255)
                                                                                                                             ttCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGat                       <  2:244906/150‑1 (MQ=255)
                                                                                                                                   ttGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGct                 >  1:20667/1‑150 (MQ=255)
                                                                                                                                    tGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGctc                <  2:1497782/150‑1 (MQ=255)
                                                                                                                                     gACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGctct               >  1:10093/1‑150 (MQ=255)
                                                                                                                                       cGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCttt             <  1:1392888/150‑1 (MQ=255)
                                                                                                                                          ttttCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTtaaa          >  2:274099/1‑150 (MQ=255)
                                                                                                                                             tCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTtaaataa       >  2:181722/1‑150 (MQ=255)
                                                                                                                                              cGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTtaaataaa      <  2:1555448/150‑1 (MQ=255)
                                                                                                                                                 tGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGc   <  1:142497/150‑1 (MQ=255)
                                                                                                                                                 tGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGc   <  2:721758/150‑1 (MQ=255)
                                                                                                                                                 tGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGc   <  1:1344780/150‑1 (MQ=255)
                                                                                                                                                  ggCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCa  <  1:589009/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTGTTCGCCTGCTCGATGATTACCTGGGCTTCCGCTTTCGCTTTTTTCAGCTGGTCGGTCGCGCTGGCCTTTGCAAGGTCAAGGTCCTTATGTGCTCGTTCTGCGGAAGCAAGGCCGTCAGCAATTTCTTTTTGACGTTTTTCGATGGCTGCCATTAATGGCGGCCATACGTACTTCATGCAGAACAGAACGAACAGGACAAACGCGATGGCCTGGCCGAGGATTGTTGCGTTAAGATTCACAGCACAATGCCTCTATTTAGTTAACGTTCTGATATTGCTCTTTAAATAAAAGCA  >  NC_000913/3920647‑3920942

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