Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,545,210 | 0 | A | T | 44.4% | ‑0.5 / 19.5 | 18 | intergenic (‑302/+575) | rpsF/rlmB | 30S ribosomal protein S6/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB |
Reads supporting (aligned to +/- strand): ref base A (5/5); new base T (2/6); total (7/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.67e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATTGGTGAATATTTGAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAACTCCT > NC_002947/5545061‑5545256 | gATTGGTGAATATTTGAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGa > 2:89214/1‑150 (MQ=255) aTTGGTGAATATTTGAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGaa > 2:1247774/1‑150 (MQ=255) tatTTGAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGccc < 2:460685/150‑1 (MQ=255) tatTTGAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGccc < 1:529657/150‑1 (MQ=255) atTTGAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCCAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCCGCCACCAGTTACTACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCg < 1:500683/150‑1 (MQ=18) tttGAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCgg < 2:61522/150‑1 (MQ=255) ttGAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGccc > 2:1099815/1‑147 (MQ=255) ttGAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGccc < 1:1099815/147‑1 (MQ=255) aGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGcc > 1:862882/1‑150 (MQ=255) gAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTAATCCACTACCATACAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCt > 1:184478/1‑150 (MQ=17) aCATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCg < 2:1454333/134‑1 (MQ=255) aCATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCTTCGCAAAGCCCg > 1:1454333/1‑134 (MQ=255) aCATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTtgg < 2:862882/150‑1 (MQ=255) aCATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTtgg < 2:168505/150‑1 (MQ=255) gCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCg > 1:448510/1‑150 (MQ=38) tGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGaa < 2:448510/150‑1 (MQ=35) ggCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAActcc < 1:89214/150‑1 (MQ=32) gCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAActcct > 2:116873/1‑150 (MQ=32) ccTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTtggcttctggct < 2:1261128/138‑1 (MQ=32) ccTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTtggcttctggct > 1:1261128/1‑138 (MQ=32) | GATTGGTGAATATTTGAACAGCCTGAAACATATTTACAGGCTGCGGGCCCTTTCGCGGGTGAACCCGCTCTCATCAAACAAAATGCAACCCATTGTTTTCTAAGCACTCAAGCGGCCACCAGTTGATACCCCAGTTTCCGGGCTTTGCGAAGCAAAGCCCGGACCTGCCGGTCTTGGCTTCTGGCTTCGAACTCCT > NC_002947/5545061‑5545256 |
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: ‑ |