Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,545,211 A→C 41.7% intergenic (‑303/+574) 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,2110AC41.7% ‑1.0 / 10.5 12intergenic (‑303/+574)rpsF/rlmB30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base C (2/3);  total (5/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

ATTGGTGAATATTTGAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTG  >  NC_002947/5545062‑5545243
                                                                                                                                                     |                                
aTTGGTGAATATTTGAACAGCCTGAAACATATTTACAGGTTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGaa                                  <  1:123518/150‑1 (MQ=255)
     tGAATATTTGAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCaaa                             <  1:117992/150‑1 (MQ=255)
     tGAATATTTGAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCAACAGTTGATACGCCAGTTTCCGGGCTTTGCGAAGCaaa                             >  1:257577/1‑150 (MQ=255)
         tatTTGAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGccc                         <  1:111833/150‑1 (MQ=255)
         tatTTAAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTACACACTCAAGCGGCCACCACTTGATACTCCAGTTTCCGGGCTTTGCTAATCAAAGccc                         <  2:687363/150‑1 (MQ=18)
               aaCAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCt                   >  1:55445/1‑150 (MQ=255)
                             tatTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGccc                         >  1:188506/1‑130 (MQ=255)
                             tatTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGccc                         <  2:188506/130‑1 (MQ=255)
                                ttACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTtggcttctg  <  2:260332/150‑1 (MQ=255)
                                             gCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCg                        >  1:187334/1‑115 (MQ=255)
                                             gCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCg                        <  2:187334/115‑1 (MQ=255)
                                                 tttCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCg                        >  1:132866/1‑111 (MQ=2)
                                                 tttCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCg                        <  2:132866/111‑1 (MQ=255)
                                                                                                                                                     |                                
ATTGGTGAATATTTGAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTG  >  NC_002947/5545062‑5545243

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