Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,545,683 T→C 20.5% 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,6830TC20.5% 72.5 / 14.3 39intergenic (‑775/+102)rpsF/rlmB30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Reads supporting (aligned to +/- strand):  ref base T (17/14);  new base C (4/4);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

ATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAACC  >  NC_002947/5545580‑5545830
                                                                                                       |                                                                                                                                                   
aTAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAAc                                                                                                       >  2:38564/1‑150 (MQ=38)
         tATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCAcc                                                                                              >  1:971829/1‑150 (MQ=255)
           ttCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCt                                                                                            >  2:978228/1‑150 (MQ=255)
              ggcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAggg                                                                                                                               >  1:137847‑M2/81‑112 (MQ=255)
              ggcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAggg                                                                                                                               <  2:137847‑M2/32‑1 (MQ=255)
               gcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGaa                                                                                                                                  >  2:532413‑M2/80‑108 (MQ=255)
               gcggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGaa                                                                                                                                  <  1:532413‑M2/29‑1 (MQ=255)
               gcggcTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAg                                                                                        <  2:1065330/150‑1 (MQ=255)
                 ggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGAt                                                                                                                      <  2:789652‑M2/41‑1 (MQ=255)
                 ggttggcccggtagggccagccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTGAt                                                                                                                      >  1:789652‑M2/78‑118 (MQ=255)
                 ggttggcccggtagggccaaccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTg                                                                                                                        >  2:884982‑M2/78‑116 (MQ=255)
                 ggttggcccggtagggccaaccgccttttctttttcaaaatgtcagcgtcgctgcggaactgaagccgtacacctcaTGTTTGCTGCGCGACAGCGCAGCCCGGAAGGGCTGGGTg                                                                                                                        <  1:884982‑M2/39‑1 (MQ=255)
                    tGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGAGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCg                                                                                   >  1:727595/1‑150 (MQ=255)
                        cTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGtt                                                                               >  2:433494/1‑150 (MQ=255)
                                  cAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGaa                                                                     >  2:419507/1‑150 (MQ=255)
                                    gTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAaga                                                                   <  2:971829/150‑1 (MQ=255)
                                            aTGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCATGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGaa                                                           <  2:61494/150‑1 (MQ=255)
                                               cTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCa                                                                                                <  1:856578/110‑1 (MQ=255)
                                               cTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCa                                                                                                >  2:856578/1‑110 (MQ=255)
                                                 tttCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGc                                                      <  1:419507/150‑1 (MQ=255)
                                                      aGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAgg                                                       <  1:34754/144‑1 (MQ=255)
                                                      aGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAgg                                                       >  2:34754/1‑144 (MQ=255)
                                                          ctcGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCAGCGCTGAGGCCGGAACAGGCTGAGAGTCa                                             <  1:163842/150‑1 (MQ=255)
                                                                          aTCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCCGGTTGCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTCGACCTGCCGc                             >  1:979807/1‑150 (MQ=25)
                                                                           tCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCt                            >  2:291939/1‑150 (MQ=255)
                                                                            cTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTg                           >  2:438617/1‑150 (MQ=255)
                                                                              gACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGAc                         >  1:920904/1‑150 (MQ=255)
                                                                                  tCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTTATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCAc                     >  1:662910/1‑150 (MQ=255)
                                                                                    aTGCTCATCATGTTTGTTTTGCGACAGCGCAGCCGGGAAGGGCTGGGTGATCTGTTACAGGGGAACGCAATCACGCTGTAGGAGCGGGTTCACCGGCGAAGAGTCCGGAACAGGCTGAGAGTCTGCGCTTGACCTGCCGCTGACggtcgg                   <  2:283637/150‑6 (MQ=255)
                                                                                      gcgcATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCt                 >  1:58048/1‑150 (MQ=255)
                                                                                       cgcATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCtt                <  1:438617/150‑1 (MQ=255)
                                                                                        gcATCCTGTTTGCTGTGAGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTc               <  2:727595/150‑1 (MQ=255)
                                                                                          aTCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCGTTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGa             >  1:84426/1‑150 (MQ=255)
                                                                                           tCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGaa            <  1:1120571/150‑1 (MQ=255)
                                                                                           tCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGaa            >  2:398677/1‑150 (MQ=255)
                                                                                             cTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCCGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACa          <  1:978228/150‑1 (MQ=255)
                                                                                             cTGTTTGCTGTGAGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACa          >  2:1013571/1‑150 (MQ=255)
                                                                                              tGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAg         <  1:1133110/150‑1 (MQ=255)
                                                                                               gTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACACAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGa        <  2:526501/150‑1 (MQ=255)
                                                                                                     tgtgCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCg                              <  2:38241/122‑1 (MQ=255)
                                                                                                     tgtgCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCg                              >  1:38241/1‑122 (MQ=255)
                                                                                                     tgtgCGACAGCGCAGCCCGGAAGGGCTGAGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAAcc  >  2:261173/1‑150 (MQ=255)
                                                                                                      gtgCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAgg                                                       <  2:1032517/96‑1 (MQ=255)
                                                                                                      gtgCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAgg                                                       >  1:1032517/1‑96 (MQ=255)
                                                                                                       |                                                                                                                                                   
ATAGCAAAGTATTCGGCGGCTGACCTACCCGGATCAGTCGGCGAATGCTTTTCTAGATCTCGATGTCAGTGCAGATCTGACGTCATGCGCATCCTGTTTGCTGTGCGACAGCGCAGCCCGGAAGGGCTGGGTGATCTCCTACAGGGGAACGCAATCACCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCGGAACAGGCTGAGAGTCAGCGCTTGACCTGCCGCTGACGCACGGCTTCGAACAGACAAACC  >  NC_002947/5545580‑5545830

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: