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_002947 | 5,545,211 | 0 | A | C | 41.7% | ‑1.0 / 10.5 | 12 | intergenic (‑303/+574) | rpsF/rlmB | 30S 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: ‑ |