Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 471,538 T→A 5.4% N467I (AAC→ATC) ‡ rpoD ← RNA polymerase sigma 70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947471,5380TA 5.4% 1044.6 / 33.0 312N467I (AAC→ATC) ‡rpoDRNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base T (150/145);  new base A (6/11);  total (156/156)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 322 total aligned reads displayed.

CGTCACCGATCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCT  >  NC_002947/471389‑471687
                                                                                                                                                     |                                                                                                                                                     
    aCCGATCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAg                                                                                                                                                   >  1:5829614/1‑150 (MQ=255)
    aCCGATCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAg                                                                                                                                                   <  1:6809526/150‑1 (MQ=255)
    aCCGATCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAg                                                                                                                                                   >  1:7471945/1‑150 (MQ=255)
    aCCGATCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAg                                                                                                                                                   >  1:2490022/1‑150 (MQ=255)
    aCCGATCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCAGGCGGGAAATACGGTTGAg                                                                                                                                                   >  2:423089/1‑150 (MQ=255)
        aTCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCttg                                                                                                                                               >  2:3510425/1‑150 (MQ=255)
           ggggTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCttgttg                                                                                                                                            <  1:1483230/150‑1 (MQ=255)
              gTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATc                                                                                                                                         <  1:2848775/150‑1 (MQ=255)
              gTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATc                                                                                                                                         >  2:8231085/1‑150 (MQ=255)
               tttCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCg                                                                                                                                        >  2:7059174/1‑150 (MQ=255)
                   cATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGtctc                                                                                                                                    <  2:6718588/150‑1 (MQ=255)
                   cATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGtctc                                                                                                                                    <  1:1500781/150‑1 (MQ=255)
                   cATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGtctc                                                                                                                                    <  1:1213441/150‑1 (MQ=255)
                    aTGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGAAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCg                                                                                                                                   >  1:7035662/1‑150 (MQ=255)
                     tGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGa                                                                                                                                  <  2:1061382/150‑1 (MQ=255)
                     tGGAGAGCGGCGCTTTGGCGATCTTCAATACCTTGCGGCTTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGa                                                                                                                                  >  2:5948496/1‑150 (MQ=255)
                       gagaTCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATc                                                                                                                                >  2:5344768/1‑150 (MQ=255)
                        agaTCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCa                                                                                                                               <  2:1242179/150‑1 (MQ=255)
                         gaTCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCAt                                                                                                                              <  2:6255926/150‑1 (MQ=255)
                             ggCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGc                                                                                                                          >  1:1569315/1‑150 (MQ=255)
                              gCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCa                                                                                                                         >  1:6302968/1‑150 (MQ=255)
                               ctctTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCAc                                                                                                                        <  1:8801244/150‑1 (MQ=255)
                                       cGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGc                                                                                                                <  1:6199397/150‑1 (MQ=255)
                                        gATCTTCAATACCTTGCTGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCg                                                                                                               <  1:7639192/150‑1 (MQ=255)
                                        gATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCg                                                                                                               <  1:2098771/150‑1 (MQ=255)
                                        gATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCg                                                                                                               >  1:4651840/1‑150 (MQ=255)
                                          tCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGa                                                                                                             <  1:3409453/150‑1 (MQ=255)
                                          tCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGa                                                                                                             <  1:4336892/150‑1 (MQ=255)
                                           cTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGAt                                                                                                            <  2:4484314/150‑1 (MQ=255)
                                             tCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATgg                                                                                                          >  1:8108191/1‑150 (MQ=255)
                                                  aCCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGgtgcgt                                                                                                     <  2:9137831/150‑1 (MQ=255)
                                                     ttGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGcc                                                                                                  <  2:7161417/150‑1 (MQ=255)
                                                     ttGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGcc                                                                                                  >  2:3722815/1‑150 (MQ=255)
                                                     ttGCGGATTTTGTCCTCAGACATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGcc                                                                                                  >  2:7415880/1‑150 (MQ=255)
                                                      tGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCt                                                                                                 <  2:2380902/150‑1 (MQ=255)
                                                         ggATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGt                                                                                              <  1:1912667/150‑1 (MQ=255)
                                                          gATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTc                                                                                             >  2:8446056/1‑150 (MQ=255)
                                                          gATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTc                                                                                             <  2:3250560/150‑1 (MQ=255)
                                                               tGTCCACAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGa                                                                                        <  1:8325586/150‑1 (MQ=255)
                                                                    tCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGaa                                                                                   >  2:1619448/1‑150 (MQ=255)
                                                                    tCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTACTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGaa                                                                                   >  1:1068016/1‑150 (MQ=255)
                                                                     cAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAAc                                                                                  <  2:4651926/150‑1 (MQ=255)
                                                                     cAGGCATTTCCAGGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCATCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAAc                                                                                  <  1:8710634/150‑1 (MQ=255)
                                                                      aGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGc                                                                                          >  2:8636003/1‑141 (MQ=255)
                                                                       ggCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACgg                                                                                <  2:138759/150‑1 (MQ=255)
                                                                        gCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACggg                                                                               <  2:7625074/150‑1 (MQ=255)
                                                                         cATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGt                                                                              <  1:8561250/150‑1 (MQ=255)
                                                                         cATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGt                                                                              <  2:8684431/150‑1 (MQ=255)
                                                                             tCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATc                                                                          <  2:7999891/150‑1 (MQ=255)
                                                                                aTGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGcc                                                                       <  1:1858864/150‑1 (MQ=255)
                                                                                aTGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGcc                                                                       >  1:7239583/1‑150 (MQ=255)
                                                                                 tGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCt                                                                      <  1:4792432/150‑1 (MQ=255)
                                                                                    gcTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGAc                                                                   <  2:2153870/150‑1 (MQ=255)
                                                                                      tCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGa                                                                    >  1:3022812/1‑147 (MQ=255)
                                                                                      tCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGa                                                                    <  2:3022812/147‑1 (MQ=255)
                                                                                      tCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACgg                                                                 >  1:1375024/1‑150 (MQ=255)
                                                                                      tCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACgg                                                                 >  2:4969332/1‑150 (MQ=255)
                                                                                          ccAGTTCTTCCGGCGTTGGTTCGCGACCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGa                                                                                                             <  1:8394851/102‑1 (MQ=255)
                                                                                            agttctTCCGGCGTTGGTTCGCGACCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATccac                                                           <  2:1965704/147‑1 (MQ=255)
                                                                                            aGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATccac                                                           <  2:3334597/150‑1 (MQ=255)
                                                                                              ctctTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCCTCTGGCGGGAACTACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATccacca                                                         <  1:5340217/150‑1 (MQ=255)
                                                                                                 ttCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGt                                                      <  1:3757260/150‑1 (MQ=255)
                                                                                                 ttCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGt                                                      <  2:8884731/150‑1 (MQ=255)
                                                                                                 ttCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGt                                                      >  1:2198320/1‑150 (MQ=255)
                                                                                                   ccGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTgg                                                    >  2:2585927/1‑150 (MQ=255)
                                                                                                   ccGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTgg                                                    <  2:6945446/150‑1 (MQ=255)
                                                                                                   ccGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTgg                                                    <  1:2572393/150‑1 (MQ=255)
                                                                                                     ggAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCa                                                  <  1:3939625/150‑1 (MQ=255)
                                                                                                           ggTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTc                                            <  2:2918139/150‑1 (MQ=255)
                                                                                                           ggTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTc                                            <  2:8012225/150‑1 (MQ=255)
                                                                                                           ggTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTc                                            <  1:8200623/150‑1 (MQ=255)
                                                                                                             ttCGCGACCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCga                                          <  1:5512133/150‑1 (MQ=255)
                                                                                                               cgcgACCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCgaga                                        >  2:3880965/1‑150 (MQ=255)
                                                                                                               cgcgACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCgaga                                        <  2:576711/150‑1 (MQ=255)
                                                                                                                gcgACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGaa                                       <  1:9037674/150‑1 (MQ=255)
                                                                                                                gcgACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGaa                                       <  2:6651357/150‑1 (MQ=255)
                                                                                                                 cgACCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCGTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAAc                                      <  2:8182277/150‑1 (MQ=255)
                                                                                                                 cgACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAAc                                      <  2:5831731/150‑1 (MQ=255)
                                                                                                                 cgACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAAc                                      <  1:3521895/150‑1 (MQ=255)
                                                                                                                   aCCCATTTCCTGCTGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACtt                                    <  1:5203405/150‑1 (MQ=255)
                                                                                                                     ccATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGt                                  >  1:5604441/1‑150 (MQ=255)
                                                                                                                      cATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTa                                 <  1:2765663/150‑1 (MQ=255)
                                                                                                                            atGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACcacg                           >  2:3061703/2‑150 (MQ=255)
                                                                                                                             tGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGa                          >  2:6408312/1‑150 (MQ=255)
                                                                                                                             tGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATCCACCATGTCGCATAGGTCGAGAACTTGTAACCACGa                          >  1:3728342/1‑150 (MQ=255)
                                                                                                                               cagcaTCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACg                        >  2:2341752/1‑150 (MQ=255)
                                                                                                                                agcaTCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACgg                       <  1:6505306/150‑1 (MQ=255)
                                                                                                                                 gcaTCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGt                      >  1:3068770/1‑150 (MQ=255)
                                                                                                                                  caTCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTa                     <  2:4397315/150‑1 (MQ=255)
                                                                                                                                   aTCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTAt                    <  2:7349692/150‑1 (MQ=255)
                                                                                                                                   aTCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTAt                    <  2:1546727/150‑1 (MQ=255)
                                                                                                                                    tCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGAc                         <  2:3814288/144‑1 (MQ=255)
                                                                                                                                    tCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTAtt                   >  2:1834601/1‑150 (MQ=255)
                                                                                                                                     cTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTc                  <  2:7670378/150‑1 (MQ=255)
                                                                                                                                      gggCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCg                 <  1:2272221/149‑1 (MQ=255)
                                                                                                                                         cGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAAc              <  1:3234057/150‑1 (MQ=255)
                                                                                                                                           ggAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACtt            <  2:2808054/150‑1 (MQ=255)
                                                                                                                                              aaTACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTc         <  2:6171593/150‑1 (MQ=255)
                                                                                                                                                   ggTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCg    >  1:1704136/1‑150 (MQ=255)
                                                                                                                                                    gTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTGTTCGAACTTGTCCACCGc   >  1:4964969/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTCACCGATCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTCCACCGCT  >  NC_002947/471389‑471687

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