Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,173,484 A→G 6.1% intergenic (+28/+20) 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,4840AG 6.1% 366.6 / 11.4 132intergenic (+28/+20)ygiN/parEquinol monooxygenase/DNA topoisomerase IV, subunit B
Reads supporting (aligned to +/- strand):  ref base A (65/58);  new base G (4/4);  total (70/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 = 2.27e-01

Only 100 of 134 total aligned reads displayed.

GTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAAC  >  NC_000913/3173336‑3173633
                                                                                                                                                    |                                                                                                                                                     
gTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGAc                                                                                                                                                      >  2:1088671/1‑150 (MQ=255)
 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)
                  ccTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCt                                                                                                                                    <  2:792868/149‑1 (MQ=255)
                   cTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCtt                                                                                                                                   <  2:936880/150‑1 (MQ=255)
                     tGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTa                                                                                                                                 <  1:1160897/150‑1 (MQ=255)
                        aGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAAc                                                                                                                              <  2:680032/150‑1 (MQ=255)
                        aGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAAc                                                                                                                              >  1:1167962/1‑150 (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                                                                            >  1:522971/1‑150 (MQ=255)
                                                                          gTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCtt                                                                            >  2:1243222/1‑150 (MQ=255)
                                                                             aCGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCg                                                                         >  2:1124579/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:440960/150‑1 (MQ=255)
                                                                                         tGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCg                                                             >  2:1573491/1‑150 (MQ=255)
                                                                                         tGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCg                                                             >  1:1647465/1‑150 (MQ=255)
                                                                                         tGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCg                                                             <  2:110800/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)
                                                                                                     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                                      <  1:1438570/150‑1 (MQ=255)
                                                                                                                gggATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcgcg                                      >  2:1640231/1‑150 (MQ=255)
                                                                                                                gggATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATcgcg                                      >  2:321432/1‑150 (MQ=255)
                                                                                                                   aTTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCa                                   >  2:566254/1‑150 (MQ=255)
                                                                                                                   aTTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCa                                   >  1:1013925/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)
                                                                                                                                     tGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTtca                 <  2:1023010/150‑1 (MQ=255)
                                                                                                                                      gCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCCtcaa                >  2:125822/1‑149 (MQ=255)
                                                                                                                                       cccGGCCATCCTGACCGGGCAATGTTCTTTGCTTTAACCCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATCTCCATCATCGCGTCAGTACGCTGATCGTCTtcata               >  2:835503/1‑149 (MQ=255)
                                                                                                                                                    aCCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAAc  <  1:984496/150‑1 (MQ=255)
                                                                                                                                                    aCCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAAc  <  1:742980/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGGGAAAGCATCGCGCATCTTGAAGCGCATCTGCAAACCCCGCACATGAAGGCGTATAGCGAAGCCGTAAAAGGTGACGTGCTGGAGATGAATATCCGTATTCTGCAGCCAGGGATTTAATCCTGCCTTGTTTGCCCGGCCATCCTGACCGGGCAATGTTCTTTCCTTTAAACCTCAATCTCCGCCATGTCGCCTTTCTCTTGCAACCAGTTGCGGCGATCTTCCGAGCGTTTCTTCGCCAGCAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAAC  >  NC_000913/3173336‑3173633

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: