Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 898,598 | 0 | T | A | 23.8% | 100.9 / 13.7 | 42 | A27A (GCT→GCA) | rlmC | 23S rRNA m(5)U747 methyltransferase |
Reads supporting (aligned to +/- strand): ref base T (25/7); new base A (0/10); total (25/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.32e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCGTTCCTGTCAGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCTAAAACCGCCGATCTTAAAAATCTGCTCGCCGACTTTCCGGTTGAGGAA > NC_000913/898550‑898646 | tCGTTCCTGTCAGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCt < 1:5045060/49‑1 (MQ=255) tCGTTCCTGTCAGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCt < 1:3649597/49‑1 (MQ=255) gTTCCTGTCAGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCTaa > 1:1035994/1‑49 (MQ=255) gTTCCTGTCAGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCTaa < 1:6686155/49‑1 (MQ=255) ttCCTGTCAGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCTaaaa > 1:658504/1‑50 (MQ=255) tCCTGTCAGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCTAAAAc > 1:4187309/1‑50 (MQ=255) tCCTGTCAGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCTAAAAc > 1:2909480/1‑50 (MQ=255) cTGTCAGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCTAAAAccg < 1:5185211/50‑1 (MQ=255) cTGTCAGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCTAAAAccg < 1:3370270/50‑1 (MQ=255) cTGTCAGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCTAAAAcc < 1:409587/49‑1 (MQ=255) tGTCAGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCTAAAAccgc > 1:2150596/1‑50 (MQ=255) tCAGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCTAAAAccgc < 1:3056774/48‑1 (MQ=255) cAGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCTAAAACCGCCGa > 1:1204905/1‑50 (MQ=255) aGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCTAAAACCGCCGAt > 1:637701/1‑50 (MQ=255) tGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCTAAAACCGCCGATCt > 1:3803118/1‑50 (MQ=255) gtaTAATGCAGCCGATTCCAGAGCAACTCTCCGCTAAAACCGCCGATCt > 1:2543544/3‑49 (MQ=255) tGCAGCCGATTCCAGAGCAACTCTCCGCTAAAACCGCCGATCTTAAAAAt > 1:4948272/1‑50 (MQ=255) tGCAGCCGATCCCAGACAAACCCTCCGCAAAAACCGCCGATCTTAAAAAt < 1:5715667/50‑1 (MQ=25) gCAGCCGATTCCAGAGCAACTCTCCGCTAAAACCGCCGATCTTAAAAATc > 1:3469080/1‑50 (MQ=255) cAGCCGATTCCAGAGCAACTCTCCGCTAAAACCGCCGATCTTAAAAATCt > 1:2492737/1‑50 (MQ=255) aGCCGATTCCAGAGCAACTCTCCGCTAAAACCGCCGATCTTAAAAATCTg > 1:242519/1‑50 (MQ=255) aGCCGATTCCAGAGCAACTCTCCGCAAAAACCGCCGATCTTAAAAATCTg < 1:130319/50‑1 (MQ=255) cggATTCGAGAGCACCCCCCCGCAAAAACCGCCGATCTTAAAAATCTGCt < 1:3188630/48‑1 (MQ=255) ttttCCAGCGCACCTCCCCGAAAAAGCCGCCGATCTTAAAAATCTGCTCg < 1:3386954/48‑1 (MQ=255) gATTCCAGAGCAACTCTCCGCTAAAACCGCCGATCTTAAAAATCTGCTCg > 1:4963638/1‑50 (MQ=255) ttCCAGAGCAACTCTCCGCTAAAACCGCCGATCTTAAAAATCTGCTCGcc > 1:2779007/1‑50 (MQ=255) ccAGAGCAACTCTCCGCTAAAACCGCCGATCTTAAAAATCTGCTCGCCGa > 1:2352537/1‑50 (MQ=255) cAGAGCACCCCTCCGCAAAAACCGCCGATCTTAAAAATCTGCTCGCCGa < 1:5127054/49‑1 (MQ=38) cAGAGCAACTCTCCGCTAAAACCGCCGATCTTAAAAATCTGCTCGCCGAc > 1:2459476/1‑50 (MQ=255) cAGAGCAACTCCCCGTTAAAACCGCCGATCTTAAAAATCTGCTCGCCGAc < 1:6222039/50‑1 (MQ=255) cAGACCACCTCCCCACAAAAACCGCCGATCTTAAAAATCTGCTCGCCGAc < 1:3667157/50‑1 (MQ=255) gagCAACTCCCCGCAAAAACCGCCGATCTTAAAAATCTGCTCGCCGACtt < 1:2166125/50‑1 (MQ=255) agCACCTCCCCGCAAAAACCGCCGATCTTAAAAATCTGCTCGCCGACttt < 1:2612098/50‑1 (MQ=255) gCAACTCTCCGCTAAAACCGCCGATCTTAAAAATCTGCTCGCCGACTTTc > 1:5103461/1‑50 (MQ=255) acctccCCGCTAAAACCGCCGATCTTAAAAATCTGCTCGCCGACTTTCCg < 1:6658211/48‑1 (MQ=255) acctccCCGCAAAAACCGCCGATCTTAAAAATCTGCTCGCCGACTTTCCg < 1:1871377/48‑1 (MQ=255) aaCTCTCCGCTAAAACCGCCGATCTTAAAAATCTGCTCGCCGACTTTCCg > 1:5527577/1‑50 (MQ=255) aCTCTCCGCTAAAACCGCCGATCTTAAAAATCTGCTCGCCGACTTTCCgg > 1:3659331/1‑50 (MQ=255) tctcCGCTAAAACCGCCGATCTTAAAAATCTGCTCGCCGACTTTCCGGtt > 1:6405608/1‑50 (MQ=255) ctcCGCTAAAACCGCCGATCTTAAAAATCTGCTCGCCGACTTTCCGGtt > 1:4441923/1‑49 (MQ=255) tcCGCAAAAACCGCCGATCTTAAAAATCTGCTCGCCGACTTTCCGGtt < 1:4061346/48‑1 (MQ=255) ccGCTAAAACCGCCGATCTTAAAAATCTGCTCGCCGACTTTCCGGTTGAg > 1:5275/1‑50 (MQ=255) cGCTAAAACCGCCGATCTTAAAAATCTGCTCGCCGACTTTCCGGTTGAg > 1:6681783/1‑49 (MQ=255) cTAAAACCGCCGATCTTAAAAATCTGCTCGCCGACTTTCCGGTTGAGGaa > 1:1554878/1‑50 (MQ=255) | TCGTTCCTGTCAGTGGATAATGCAGCCGATTCCAGAGCAACTCTCCGCTAAAACCGCCGATCTTAAAAATCTGCTCGCCGACTTTCCGGTTGAGGAA > NC_000913/898550‑898646 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |