Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 471,531 A→T 5.2% I469I (ATT→ATA) ‡ rpoD ← RNA polymerase sigma 70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947471,5310AT 5.2% 1050.5 / 28.7 320I469I (ATT→ATA) ‡rpoDRNA polymerase sigma 70 factor
Reads supporting (aligned to +/- strand):  ref base A (153/150);  new base T (6/11);  total (159/161)
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.98e-01

Only 100 of 324 total aligned reads displayed.

TCTTCGTCGTCACCGATCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTC  >  NC_002947/471382‑471680
                                                                                                                                                     |                                                                                                                                                      
     gtcgtcACCGATCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACg                                                                                                                                                  <  1:4197460/150‑1 (MQ=255)
      tcgtcACCGATCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACgg                                                                                                                                                 >  1:5465401/1‑150 (MQ=255)
        gtcACCGATCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGtt                                                                                                                                               <  1:6099557/150‑1 (MQ=255)
        gtcACCGATCGGGGGTTCCATGGAGATCGGCTCTTTGGCGATCGTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGtt                                                                                                                                               <  2:6914200/150‑1 (MQ=255)
           aCCGATCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAg                                                                                                                                            <  1:6809526/150‑1 (MQ=255)
           aCCGATCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAg                                                                                                                                            >  1:7471945/1‑150 (MQ=255)
           aCCGATCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAg                                                                                                                                            >  1:2490022/1‑150 (MQ=255)
                tCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCttgt                                                                                                                                       <  2:3734668/150‑1 (MQ=255)
                  ggggTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCttgttg                                                                                                                                     <  1:1483230/150‑1 (MQ=255)
                   gggTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGa                                                                                                                                    >  1:5323314/1‑150 (MQ=255)
                     gTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATc                                                                                                                                  <  1:2848775/150‑1 (MQ=255)
                     gTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATc                                                                                                                                  >  2:8231055/1‑150 (MQ=255)
                          cATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGtctc                                                                                                                             <  1:1213441/150‑1 (MQ=255)
                          cATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGtctc                                                                                                                             <  1:1500781/150‑1 (MQ=255)
                            tGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGa                                                                                                                           >  2:7151791/1‑150 (MQ=255)
                            tGGAGAGCGGCGCTTTGGCGATCTTCAATACCTTGCGGCTTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGa                                                                                                                           >  2:5948496/1‑150 (MQ=255)
                             ggAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGAt                                                                                                                          >  2:3234057/1‑150 (MQ=255)
                               agaTCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCa                                                                                                                        <  2:1242179/150‑1 (MQ=255)
                                gaTCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCAt                                                                                                                       <  1:458330/150‑1 (MQ=255)
                                gaTCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCAt                                                                                                                       <  2:6255134/150‑1 (MQ=255)
                                    ggCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGc                                                                                                                   >  1:1569315/1‑150 (MQ=255)
                                     gCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCa                                                                                                                  >  2:5873266/1‑150 (MQ=255)
                                     gCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCa                                                                                                                  >  1:6302968/1‑150 (MQ=255)
                                        ctTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCACTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCg                                                                                                               <  1:3440707/150‑1 (MQ=255)
                                               gATCTTCAATACCTTGCTGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCg                                                                                                        <  1:7639192/150‑1 (MQ=255)
                                                 tCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGa                                                                                                      <  2:2439478/150‑1 (MQ=255)
                                                  cTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGAt                                                                                                     >  1:8474674/1‑150 (MQ=255)
                                                  cTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGAt                                                                                                     <  2:5768831/150‑1 (MQ=255)
                                                   ttCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATc                                                                                                                                  >  1:3331986/1‑120 (MQ=255)
                                                   ttCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATg                                                                                                    >  1:6976816/1‑150 (MQ=255)
                                                         aCCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGgtgcgt                                                                                              <  2:9137831/150‑1 (MQ=255)
                                                           cTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGc                                                                                                         >  1:8721172/1‑137 (MQ=255)
                                                            ttGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGTAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTAGATCATGTGCACCGGAATGCGGATGGTGCGTGcc                                                                                           >  2:3572355/1‑150 (MQ=255)
                                                            ttGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGcc                                                                                           >  2:7961463/1‑150 (MQ=255)
                                                                 gATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTc                                                                                      <  2:3250560/150‑1 (MQ=255)
                                                                  aTTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCg                                                                                     >  1:4713668/1‑150 (MQ=255)
                                                                   ttttGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCgg                                                                                    <  2:8007686/150‑1 (MQ=255)
                                                                        tCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATc                                                                               <  2:4697004/150‑1 (MQ=255)
                                                                           tCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATc                                                                               >  2:6560172/1‑147 (MQ=255)
                                                                           tCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGaa                                                                            >  2:1619448/1‑150 (MQ=255)
                                                                           tCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTACTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGaa                                                                            >  1:1068016/1‑150 (MQ=255)
                                                                             aGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGc                                                                                   <  1:8636003/141‑1 (MQ=255)
                                                                             aGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGc                                                                                   >  2:8636003/1‑141 (MQ=255)
                                                                              ggCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACgg                                                                         <  2:138759/150‑1 (MQ=255)
                                                                              ggCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACgg                                                                         <  1:3011077/150‑1 (MQ=255)
                                                                               gCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACggg                                                                        <  2:7625071/150‑1 (MQ=255)
                                                                                cATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGt                                                                       <  2:8684431/150‑1 (MQ=255)
                                                                                cATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGt                                                                       <  1:1742884/150‑1 (MQ=255)
                                                                                 aTTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTg                                                                      >  2:4500657/1‑150 (MQ=255)
                                                                                 aTTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTg                                                                      >  1:390498/1‑150 (MQ=255)
                                                                                   ttCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGAt                                                                    >  2:5201875/1‑150 (MQ=255)
                                                                                       aTGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGcc                                                                >  1:7239583/1‑150 (MQ=255)
                                                                                        tGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCt                                                               >  2:5446256/1‑150 (MQ=255)
                                                                                          cgcTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATTACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTg                                                              <  2:705026/150‑1 (MQ=255)
                                                                                           gcTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGAc                                                            <  2:2153870/150‑1 (MQ=255)
                                                                                           gcTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGAc                                                            <  1:2671401/150‑1 (MQ=255)
                                                                                             tCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCGATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACgg                                                          >  2:7371407/1‑150 (MQ=255)
                                                                                             tCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGa                                                             <  2:3023022/147‑1 (MQ=255)
                                                                                             tCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACgg                                                          >  1:7762501/1‑150 (MQ=255)
                                                                                             tCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACgg                                                          >  1:1375024/1‑150 (MQ=255)
                                                                                               gCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGAt                                                        >  1:4960268/1‑150 (MQ=255)
                                                                                               agccAGTTCTTCCGGCGTTGGTTCGCGACCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAAt                                                        >  2:3919178/3‑150 (MQ=255)
                                                                                                gccAGTTCTTCCGGCGTTGGTTCGCGACCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATc                                                       <  2:4303712/149‑1 (MQ=255)
                                                                                                 ccAGTTCTTCCGGCGTTGGTTCGCGACCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGa                                                                                                      >  2:8394851/1‑102 (MQ=255)
                                                                                                 ccAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATcc                                                      <  1:8211467/150‑1 (MQ=255)
                                                                                                   agttctTCCGGCGTTGGTTCGCGACCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGAATccac                                                    <  2:1965704/147‑1 (MQ=255)
                                                                                                   aGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATccac                                                    <  2:3334597/150‑1 (MQ=255)
                                                                                                        ttCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGt                                               >  2:1159698/1‑150 (MQ=255)
                                                                                                        ttCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGt                                               <  2:8884731/150‑1 (MQ=255)
                                                                                                        ttCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGt                                               >  1:2198320/1‑150 (MQ=255)
                                                                                                        ttCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGt                                               <  1:2004749/150‑1 (MQ=255)
                                                                                                        ttCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGt                                               >  1:511382/1‑150 (MQ=255)
                                                                                                        ttCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGt                                               <  1:3757260/150‑1 (MQ=255)
                                                                                                            ggAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCa                                           <  1:5544185/150‑1 (MQ=255)
                                                                                                                 cGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGt                                      <  2:1848673/150‑1 (MQ=255)
                                                                                                                  ggTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTc                                     <  2:2918139/150‑1 (MQ=255)
                                                                                                                   gTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCg                                    >  1:8112404/1‑150 (MQ=255)
                                                                                                                    ttCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCga                                   <  1:1863151/150‑1 (MQ=255)
                                                                                                                    ttCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCga                                   <  2:2597752/150‑1 (MQ=255)
                                                                                                                       gcgACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGaa                                <  1:9037674/150‑1 (MQ=255)
                                                                                                                        cgACCGATTTCATGCAGCATCTGGCGCGATGTACGAACGATCTTGTTGATCGTTTCGATCATGTGCACCGGAATACGGATGGTGCGCGCCTGATCGGCAATCGAACGGGTGATTGCCTGACGACTCCACCATGTCGCATAGGTCGAGAAc                               <  2:8181901/150‑1 (MQ=255)
                                                                                                                          aCCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACtt                             <  1:5809920/150‑1 (MQ=255)
                                                                                                                           cccATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTg                            <  1:8385926/150‑1 (MQ=255)
                                                                                                                             cATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTa                          <  1:2765663/150‑1 (MQ=255)
                                                                                                                             cATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTa                          <  1:2177764/150‑1 (MQ=255)
                                                                                                                                 tCCTGCAGCATATGGCGCGAAATACCGTTGAGCTTGTTTATCGTCTCGGTCATGTGCACCGCAAAGACTATGGTGCGAGCATTGTCTGCGATCGAACGGGTGATCGCCTGAAGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCcg                      >  2:6717361/1‑150 (MQ=255)
                                                                                                                                    tGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGa                   >  2:7670364/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:7313847/150‑1 (MQ=255)
                                                                                                                                          aTCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTAt             <  2:1546727/150‑1 (MQ=255)
                                                                                                                                          aTCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTAt             <  2:3816941/150‑1 (MQ=255)
                                                                                                                                            cTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTc           <  2:6375357/150‑1 (MQ=255)
                                                                                                                                            cTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTc           <  2:4452525/150‑1 (MQ=255)
                                                                                                                                             gggCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCg          <  1:2272221/149‑1 (MQ=255)
                                                                                                                                               gCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGaa        <  1:8881591/150‑1 (MQ=255)
                                                                                                                                                cGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAAc       <  1:3234057/150‑1 (MQ=255)
                                                                                                                                                  ggAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACtt     <  2:2808054/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCTTCGTCGTCACCGATCGGGGTTTCCATGGAGATCGGCTCTTTGGCGATCTTCAATACCTTGCGGATTTTGTCCTCAGGCATTTCCATGCGCTCGCCCAGCTCTTCCGGAGTCGGTTCGCGACCCATTTCCTGCAGCATCTGGCGGGAAATACGGTTGAGCTTGTTGATCGTCTCGATCATGTGCACCGGAATGCGGATGGTGCGTGCCTGGTCGGCGATCGAACGGGTGATCGCCTGACGGATCCACCAGGTGGCATAGGTCGAGAACTTGTAGCCGCGACGGTATTCGAACTTGTC  >  NC_002947/471382‑471680

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: