Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,173,480 C→G 6.0% intergenic (+24/+24) ygiN → / ← parE quinol monooxygenase/DNA topoisomerase IV, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,173,4800CG 6.0% 428.6 / 14.7 133intergenic (+24/+24)ygiN/parEquinol monooxygenase/DNA topoisomerase IV, subunit B
Reads supporting (aligned to +/- strand):  ref base C (67/58);  new base G (4/4);  total (71/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01

Only 100 of 134 total aligned reads displayed.

GAGCAGTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATC  >  NC_000913/3173331‑3173623
                                                                                                                                                     |                                                                                                                                               
      tGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGAcc                                                                                                                                           <  1:1542270/150‑1 (MQ=255)
      tGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAACGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGAcc                                                                                                                                           <  1:910242/150‑1 (MQ=255)
        ggAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCgg                                                                                                                                         >  1:1023010/1‑150 (MQ=255)
        ggAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCgg                                                                                                                                         <  1:705904/150‑1 (MQ=255)
         gAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGc                                                                                                                                       >  2:1471839/1‑150 (MQ=255)
          aaaGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGTCAGGATGGCCGGGc                                                                                                                                       <  1:465444/150‑1 (MQ=255)
             gCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGTCAGGATGGCCGGGCaaa                                                                                                                                    <  2:441237/150‑2 (MQ=255)
             gCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAAt                                                                                                                                    <  2:31035/150‑1 (MQ=255)
              cATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATg                                                                                                                                   >  1:278353/1‑150 (MQ=255)
               aTCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGTCAGGATGGCCGGGCaaaca                                                                                                                                  <  1:474940/150‑4 (MQ=255)
                tCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGtt                                                                                                                                 <  1:812130/150‑1 (MQ=255)
                 cgcgCATCTTGACGCGCATCTGCAAACCCCGCACAGGAAGGGGGATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATTAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTc                                                                                                                                <  2:270149/150‑1 (MQ=255)
                 cgcgCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTc                                                                                                                                <  1:519448/150‑1 (MQ=255)
                 cgcgCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTc                                                                                                                                >  2:171056/1‑150 (MQ=255)
                  gcgcATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGAAATGTTCt                                                                                                                               >  1:1536487/1‑150 (MQ=255)
                  gcgcATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCGGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCt                                                                                                                               <  1:1638195/150‑1 (MQ=255)
                  gcgcATCTTGAAGCGCATCTGCAAACCCCGAACATGAAGGCGTATAGCGAATCCGTAAAATGTGAAGTGATGGAGATGAATATCCGTATTCTGCAGCCAGGGATGTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCTGGCAATGTTCt                                                                                                                               >  1:1543174/1‑150 (MQ=255)
                  gcgcATCTAGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGACGCCGTAAAAGGTGACGTGCGGGAGAGGAATCTCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCt                                                                                                                               <  1:111353/150‑1 (MQ=255)
                    gcATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCttt                                                                                                                             <  2:1626491/150‑1 (MQ=255)
                    gcATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGAGGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCttt                                                                                                                             <  1:414247/150‑1 (MQ=255)
                     cATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTc                                                                                                                            <  2:1380176/150‑1 (MQ=255)
                       tCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCt                                                                                                                          <  1:989396/150‑1 (MQ=255)
                       ccTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCt                                                                                                                          <  2:792868/149‑1 (MQ=255)
                        cTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCtt                                                                                                                         <  2:936880/150‑1 (MQ=255)
                          tGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTa                                                                                                                       <  1:1160897/150‑1 (MQ=255)
                             aGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAAc                                                                                                                    >  1:1167962/1‑150 (MQ=255)
                             aGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAAc                                                                                                                    <  2:680032/150‑1 (MQ=255)
                              gcgcATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAAcc                                                                                                                   >  2:1170268/1‑150 (MQ=255)
                               cgcATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCt                                                                                                                  >  2:1377547/1‑150 (MQ=255)
                                gcATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTc                                                                                                                 >  1:1649265/1‑150 (MQ=255)
                                  aTCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCaa                                                                                                               >  2:1527661/1‑150 (MQ=255)
                                   tCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAAt                                                                                                              >  1:1342484/1‑150 (MQ=255)
                                       cAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGCCCGGGCAATGTTATGTCCTTTAAACCTCAATCTcc                                                                                                          >  2:701740/1‑150 (MQ=255)
                                         aaCCCCTCATATTAAGGCTTATATCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGc                                                                                                        <  2:1322613/150‑1 (MQ=255)
                                         aaCCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGc                                                                                                        <  1:233038/150‑1 (MQ=255)
                                             ccGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATg                                                                                                    >  2:1036568/1‑150 (MQ=255)
                                              cGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAATCCTCTATCTCCGCCATGt                                                                                                   >  1:632943/1‑150 (MQ=255)
                                               gCACATGAAGGAGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGACTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTc                                                                                                  <  1:1170268/150‑1 (MQ=255)
                                                  caTGAAGGCTTATAGCGAAGCCGTAATAGGTGACGGGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATGTCCGCCATGTCGcc                                                                                               <  2:128902/150‑1 (MQ=255)
                                                   aTGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATTTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCt                                                                                              <  2:1691073/150‑1 (MQ=255)
                                                   aTGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCt                                                                                              <  2:1404027/150‑1 (MQ=255)
                                                        ggCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTtctc                                                                                         <  1:356736/150‑1 (MQ=255)
                                                         gCGTATAGCGAGGCCGTAAAAGGTGACGTGCTGGAGATGAATCTCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTctct                                                                                        <  1:1457268/150‑1 (MQ=255)
                                                         gCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTctct                                                                                        <  1:1448732/150‑1 (MQ=255)
                                                          cGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCtt                                                                                       >  2:750431/1‑150 (MQ=255)
                                                           gTATAGCGAAGCCGTAACAGGTGACGTGCGGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTGGCCTTTCTCTTg                                                                                      <  2:1242098/150‑1 (MQ=255)
                                                           gTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCGTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTg                                                                                      <  1:327541/150‑1 (MQ=255)
                                                           gTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTg                                                                                      >  1:1005458/1‑150 (MQ=255)
                                                           gTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTg                                                                                      <  2:967065/150‑1 (MQ=255)
                                                                gCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAAcc                                                                                 <  1:1534171/150‑1 (MQ=255)
                                                                      ccGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTgc                                                                           <  1:487282/150‑1 (MQ=255)
                                                                        gTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCAGCCATGTCGCCTTTCTCTTGCAACCAGTTgcgg                                                                         >  2:1575943/1‑150 (MQ=255)
                                                                         tAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTgcggc                                                                        >  1:925141/1‑150 (MQ=255)
                                                                         tAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTgcggc                                                                        <  1:943719/150‑1 (MQ=255)
                                                                          aaaaGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTgcggcg                                                                       >  1:208394/1‑150 (MQ=255)
                                                                           aaaGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGa                                                                      >  1:1298884/1‑150 (MQ=255)
                                                                             aGGTGACGTGCTGGAGATGAATATCCGTATGCTGCAGCCAGGGATTTAAACCAGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTGCTTTCCTTTTATCCTCAATCTCCGCCTTGTCGCCTTTCTCTTGCAACCAGTGGCGGCGATc                                                                    <  1:35066/150‑1 (MQ=255)
                                                                               gTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCtt                                                                  >  2:1243222/1‑150 (MQ=255)
                                                                               gTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCtt                                                                  >  1:522971/1‑150 (MQ=255)
                                                                                  aCGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCg                                                               >  2:1124579/1‑150 (MQ=255)
                                                                                     tGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGc                                                            >  1:748131/1‑150 (MQ=255)
                                                                                       cTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGt                                                          >  2:1468469/1‑150 (MQ=255)
                                                                                         ggAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGttt                                                        >  1:597918/1‑150 (MQ=255)
                                                                                         ggAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATACTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGttt                                                        <  1:1535493/150‑1 (MQ=255)
                                                                                              tGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCg                                                   <  2:110800/150‑1 (MQ=255)
                                                                                              tGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCg                                                   >  2:1573491/1‑150 (MQ=255)
                                                                                              tGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCg                                                   >  1:1647465/1‑150 (MQ=255)
                                                                                              tGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCg                                                   <  2:440960/150‑1 (MQ=255)
                                                                                               gAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGAGGCGATCTTCCGAGCGTTTCTTCGc                                                  >  2:870768/1‑150 (MQ=255)
                                                                                                aaTATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGcc                                                 <  1:1124579/150‑1 (MQ=255)
                                                                                                  tatCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCAGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCcag                                               <  1:708254/150‑1 (MQ=255)
                                                                                                   atCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCcagc                                              >  2:1027284/1‑150 (MQ=255)
                                                                                                    tCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGACATGTCTCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCcagcc                                             >  2:496133/1‑149 (MQ=255)
                                                                                                     ccGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCcagcag                                            <  1:780381/150‑1 (MQ=255)
                                                                                                     ccGTATTCTGCAGCCAGGGATTTAATCATGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCAATGTCGCCTTTCTATTGCAACCAGTTTCGGAGATCTTCCTAGCGTTTCTTCGCcagcaa                                            >  1:316204/1‑149 (MQ=255)
                                                                                                       gTATTCTGCAGCCAGGGATTTACTCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCCgcagca                                          >  2:1069440/1‑150 (MQ=255)
                                                                                                          ttCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCatat                                       <  1:1088671/150‑1 (MQ=255)
                                                                                                           tCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATc                                      <  2:164509/150‑1 (MQ=255)
                                                                                                             tGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGTGCGTTTCTTCGCCAGCAGCATATCca                                    <  1:171056/150‑1 (MQ=255)
                                                                                                                  ccAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCatcatc                               >  1:1322510/1‑150 (MQ=255)
                                                                                                                   cAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcg                              >  2:1562150/1‑150 (MQ=255)
                                                                                                                    aGGGATTTTATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcgc                             >  1:557861/1‑150 (MQ=255)
                                                                                                                     gggATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcgcg                            >  2:1640231/1‑150 (MQ=255)
                                                                                                                     gggATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcgcg                            <  1:1438570/150‑1 (MQ=255)
                                                                                                                     gggATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcgcg                            >  2:321432/1‑150 (MQ=255)
                                                                                                                        aTTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCa                         >  1:1013925/1‑150 (MQ=255)
                                                                                                                        aTTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCa                         >  2:566254/1‑150 (MQ=255)
                                                                                                                              tCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGc                   <  2:499213/150‑1 (MQ=255)
                                                                                                                                cTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTg                 >  1:939380/1‑150 (MQ=255)
                                                                                                                                 tGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGa                >  1:1551728/1‑150 (MQ=255)
                                                                                                                                  gCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGAt               >  2:682/1‑150 (MQ=255)
                                                                                                                                    cgtGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGCTCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCg             <  1:76129/148‑1 (MQ=255)
                                                                                                                                      tGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTc           <  2:191828/150‑1 (MQ=255)
                                                                                                                                       gTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCt          <  2:261715/150‑1 (MQ=255)
                                                                                                                                       cattGCCCGGTCAGGATGGCCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCt          >  2:874135/3‑150 (MQ=255)
                                                                                                                                        attGCCCGGTCAGGATGGCCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCCTCGCGTCAGTCCGCTGATCGTCtt         >  2:905239/2‑150 (MQ=255)
                                                                                                                                          tGCCCGGTCAGGATGGCCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTACTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTATCATCGCCAGCAGCATATCCAACATCGCGTCAGTACGCTGATCGTCTtca       >  1:608835/1‑150 (MQ=255)
                                                                                                                                           gCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTtcat      >  2:1138577/1‑150 (MQ=255)
                                                                                                                                           gCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCCtcaa      >  2:125822/1‑149 (MQ=255)
                                                                                                                                            cccGGCCATCCTGACCGGGCAATGTTCTTTGCTTTAACCCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATCTCCATCATCGCGTCAGTACGCTGATCGTCTtcata     >  2:835503/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GAGCAGTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATC  >  NC_000913/3173331‑3173623

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: