Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,545,683 T→C 14.8% intergenic (‑775/+102) rpsF ← / ← rlmB 30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,545,6830TC14.8% 181.3 / 23.2 82intergenic (‑775/+102)rpsF/rlmB30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base C (7/5);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 7.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

GGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAACCCC  >  NC_002947/5545573‑5545832
                                                                                                               |                                                                                                                                                     
gggAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACa                                                                                                                <  1:1791190/150‑1 (MQ=32)
 ggAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAg                                                                                                               >  1:1417169/1‑150 (MQ=33)
  gAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAgg                                                                                                              <  2:1417169/150‑1 (MQ=34)
     ccATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGa                                                                                                           <  1:299789/150‑1 (MQ=37)
      cATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGaa                                                                                                          >  2:1634735/1‑150 (MQ=38)
          gCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAgg                                                                                                              >  2:1407130/1‑142 (MQ=34)
          gCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAgg                                                                                                              <  1:1407130/142‑1 (MQ=34)
            aaaaagaaaaggcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTAcgggcgagggcgc                                                                                                     >  1:1374367‑M2/91‑144 (MQ=255)
               aagaaaaggcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACagtggt                                                                                                           <  2:41195‑M2/54‑2 (MQ=255)
               aagaaaaggcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACagtggt                                                                                                           >  1:41195‑M2/88‑140 (MQ=255)
                tATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCAcc                                                                                                >  2:1439172/1‑150 (MQ=255)
                 aTTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCAccc                                                                                               >  2:1186319/1‑150 (MQ=255)
                  ttCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGCCGTGGATAATGTAGGAGGCGCATGCTCATCCTGTTTGCTTGGCGACAGCGCAGACCAGTAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCt                                                                                              <  1:532179/150‑1 (MQ=9)
                   tCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTg                                                                                             <  1:2183669/150‑1 (MQ=255)
                    aaggcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGaa                                                                                                                                    >  1:1277628‑M2/83‑111 (MQ=255)
                    aaggcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGaa                                                                                                                                    <  2:1277628‑M2/29‑1 (MQ=255)
                     ggcggcTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTa                                                                                           <  1:1186319/150‑1 (MQ=255)
                     aggcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTa                                                                                                                  >  1:2104730‑M2/82‑128 (MQ=255)
                      gcggcTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAg                                                                                          >  2:2128242/1‑150 (MQ=255)
                      gcggcTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAg                                                                                          <  1:1634735/150‑1 (MQ=255)
                         ggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGt                                                                                                                           >  2:1237287‑M2/78‑115 (MQ=255)
                         ggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGt                                                                                                                           <  1:1237287‑M2/38‑1 (MQ=255)
                         ggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGAtctc                                                                                                                     >  2:1768019‑M2/78‑121 (MQ=255)
                         ggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGAtctc                                                                                                                     <  1:1768019‑M2/44‑1 (MQ=255)
                           ttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGt                                                                                                                           <  2:2172679‑M2/38‑1 (MQ=255)
                           ttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGt                                                                                                                           >  1:2172679‑M2/76‑113 (MQ=255)
                             aCCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCggg                                                                                   >  1:822466/1‑150 (MQ=255)
                             aCCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCggg                                                                                   <  1:1439172/150‑1 (MQ=255)
                                tACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAggg                                                                                                             <  1:2046171/121‑1 (MQ=35)
                                tACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAggg                                                                                                             >  2:2046171/1‑121 (MQ=35)
                                tACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTc                                                                                <  1:1622954/150‑1 (MQ=255)
                                 aCCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCa                                                                               >  2:1155696/1‑150 (MQ=255)
                                  cccGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCAc                                                                              <  1:70083/150‑1 (MQ=255)
                                         cAGTTGGCGAACTCTTTTCTATATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCAGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGaa                                                                       <  2:605169/150‑1 (MQ=18)
                                         cAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGaa                                                                       <  2:1319159/150‑1 (MQ=255)
                                         cAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGaa                                                                       <  1:707227/150‑1 (MQ=255)
                                         cAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGaa                                                                       <  1:1155696/150‑1 (MQ=255)
                                            tCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAcgag                                                                    >  2:711909/1‑150 (MQ=255)
                                            tCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAagag                                                                    >  1:246085/1‑150 (MQ=255)
                                             cGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAgg                                                                   <  2:1091069/150‑1 (MQ=255)
                                              ggCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCCGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGc                                                                  >  1:266711/1‑150 (MQ=255)
                                               gCGAATGCTTTTCTAGATCTGGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTa                                                                                                                  <  1:1897264/101‑1 (MQ=14)
                                               gCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTa                                                                                                                  >  2:1897264/1‑101 (MQ=31)
                                                 gAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACCGGGGACCGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCgg                                                               <  1:537818/150‑1 (MQ=255)
                                                  aaTGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGa                                                              >  1:720240/1‑150 (MQ=255)
                                                   aTGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGaa                                                             >  1:1366028/1‑150 (MQ=255)
                                                     gCTTTTCTAGATCTCGATGTCTGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACa                                                           >  1:518362/1‑150 (MQ=255)
                                                       ttttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAgg                                                         <  2:878118/150‑1 (MQ=255)
                                                           cTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCt                                                                                              <  1:1319687/109‑1 (MQ=255)
                                                           cTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCt                                                                                              >  2:1319687/1‑109 (MQ=255)
                                                           cTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTga                                                     <  1:1803293/150‑1 (MQ=255)
                                                              gATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGagag                                                  >  1:970326/1‑150 (MQ=255)
                                                               aTCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGt                                                 >  1:1416462/1‑150 (MQ=255)
                                                                tctcGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGTAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTc                                                >  1:248123/1‑150 (MQ=255)
                                                                  tcGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAg                                              >  1:1915043/1‑150 (MQ=255)
                                                                    gATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAAc                                                            <  1:969926/134‑1 (MQ=255)
                                                                    gATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAAc                                                            >  2:969926/1‑134 (MQ=255)
                                                                     aTGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTga                                                     <  2:1032043/140‑1 (MQ=255)
                                                                     aTGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTga                                                     >  1:1032043/1‑140 (MQ=255)
                                                                        tCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAgg                                                         >  2:454256/1‑133 (MQ=255)
                                                                        tCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAgg                                                         <  1:454256/133‑1 (MQ=255)
                                                                         cAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGa                                       >  1:1517204/1‑150 (MQ=255)
                                                                         cAGTGCAGATCTGACGTCATCTTCAGCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGCGTGATCTCCTACAGGGGCACCTGATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGa                                       <  1:711909/150‑1 (MQ=25)
                                                                          aGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCg                                                                <  1:1630505/124‑1 (MQ=255)
                                                                          aGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCg                                                                >  2:1630505/1‑124 (MQ=255)
                                                                           gTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGAcc                                     <  1:1340198/150‑1 (MQ=255)
                                                                             gCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTg                                   >  2:1008838/1‑150 (MQ=255)
                                                                              cAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGc                                  >  2:929991/1‑150 (MQ=255)
                                                                                 aTCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGc                               >  1:1143135/1‑150 (MQ=255)
                                                                                    tGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGa                            <  2:246548/150‑1 (MQ=255)
                                                                                            aTGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACgg                     >  2:701680/1‑150 (MQ=255)
                                                                                               ggcATCCTGTTTAATGTGAGACCGAGCAGCCCGGACGGGATGGGTTCAATCCCACAGGCGCACGAAATCACCCTGTAGGTGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCtt                  <  2:1653564/149‑1 (MQ=25)
                                                                                               cgcATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCtt                  <  1:701680/150‑1 (MQ=255)
                                                                                               cgcATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCtt                  <  2:1733218/150‑1 (MQ=255)
                                                                                                    ccTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAAc             <  2:276560/150‑1 (MQ=255)
                                                                                                      tGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAg           <  2:822466/150‑1 (MQ=255)
                                                                                                      tGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAg           >  1:2273554/1‑150 (MQ=255)
                                                                                                      tGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAg           >  1:1817680/1‑150 (MQ=255)
                                                                                                        tttGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGAc         <  2:2273554/150‑1 (MQ=255)
                                                                                                        tttGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGAc         >  1:240508/1‑150 (MQ=255)
                                                                                                            cTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAAc     >  2:2274987/1‑150 (MQ=255)
                                                                                                            cTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAAc     <  1:2274987/150‑1 (MQ=255)
                                                                                                             tgtgCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACa            <  2:150297/142‑1 (MQ=255)
                                                                                                             tgtgCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACa            >  1:150297/1‑142 (MQ=255)
                                                                                                             tgtgCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAAcc    >  2:1484508/1‑150 (MQ=255)
                                                                                                               tgCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAAcccc  >  2:379666/1‑150 (MQ=255)
                                                                                                               tgCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAAcccc  >  1:312599/1‑150 (MQ=255)
                                                                                                               |                                                                                                                                                     
GGGAACCATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAACCCC  >  NC_002947/5545573‑5545832

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