Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,906,769 | 0 | A | C | 25.5% | 96.9 / 10.6 | 47 | Y98D (TAT→GAT) | rlmA | 23S rRNA m(1)G745 methyltransferase |
Reads supporting (aligned to +/- strand): ref base A (28/7); new base C (0/12); total (28/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.64e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTG > NC_000913/1906630‑1906895 | aaCGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATa < 2:123870/140‑1 (MQ=255) aCGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCTGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATAc > 1:1016426/1‑140 (MQ=255) gTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTcaccac < 2:1190848/140‑1 (MQ=255) gCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGc < 2:475017/140‑1 (MQ=255) gCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGc > 2:241741/1‑140 (MQ=255) acacAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAg > 1:939035/1‑140 (MQ=255) cacaAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGc > 1:376396/1‑140 (MQ=255) cacaAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGc > 1:399340/1‑140 (MQ=255) ccaAAAGGTGCCCGGCGGATGGCGTTCCGCCGCGGCTTTTAGCGCTCCCTTCGAAACATCCAGACAAACCGTGGTGATTTCGGGCATCGCACCGGCAAAG‑GCGTGTGTGTAATCCCCTTCCCCGCAGCCAATATCCAGcg < 2:1271289/139‑2 (MQ=255) acaAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCa < 1:997886/140‑1 (MQ=255) aaGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCAccgcc < 2:1372351/140‑1 (MQ=255) gACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTgg > 1:1271289/1‑140 (MQ=255) cGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGTTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCAAATATCCAGCa < 1:979980/125‑1 (MQ=255) cGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTat > 1:262829/1‑140 (MQ=255) cGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCAAATATCCAGCa > 2:979980/1‑125 (MQ=255) aTAGCTTTTCGCCGCCGCTTTTATCGCTCCCTTCGAAACACCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATCCCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTatcat < 2:1035667/140‑1 (MQ=255) aGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTAtcatca > 1:881077/1‑140 (MQ=255) tttCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGcc < 2:1016426/140‑1 (MQ=255) ccgcTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTa < 2:373763/122‑1 (MQ=255) ccgcTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTa > 1:373763/1‑122 (MQ=255) aTCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATCCCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGc < 2:1377126/140‑1 (MQ=255) tACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTTTAATCCCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACaa < 2:347268/140‑1 (MQ=255) tACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTc > 2:146481/1‑127 (MQ=255) tACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTc < 1:146481/127‑1 (MQ=255) aCCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAAt > 2:553350/1‑140 (MQ=255) cTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTg > 2:418502/1‑140 (MQ=255) cTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTg > 1:1072058/1‑140 (MQ=255) cTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGCGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTg > 1:378002/1‑140 (MQ=255) tCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCa > 2:843679/1‑140 (MQ=255) gacaggACCAAACGTGGTGATTCCGGGCAA‑GCCACGGCAAATGGGGTTTGGATAATCCCCTCCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCa < 2:318914/135‑1 (MQ=255) aTCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCa > 2:488531/1‑140 (MQ=255) tCCAGACCAAACGTGGTGATTTCGGGTAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAg > 2:983334/1‑140 (MQ=255) tCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAg > 2:1329981/1‑140 (MQ=255) tCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAg > 1:990051/1‑140 (MQ=255) tCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAg > 2:338274/1‑140 (MQ=255) ccAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGc < 2:376396/140‑1 (MQ=255) ccAGACAAAACGTGGTGATTTCGGGCACCGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGc < 2:381691/140‑1 (MQ=255) aGACCAAACGTGGTGTTTTGGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATCCCCTTCACCACAGCCAAAATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCgg < 2:65679/140‑1 (MQ=255) gccaaaaCGTGGTGATTTCGGGCAACGCATCGGAAAAT‑GCGTGTGTGTAATCCCCTTCCCCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGc < 2:185265/136‑1 (MQ=255) gACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGc < 2:602460/140‑1 (MQ=255) aCGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATaa > 1:597339/1‑140 (MQ=255) gtgATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGTTaactgtc > 1:318914/1‑135 (MQ=255) gtgATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTcc > 2:413104/1‑140 (MQ=255) tgATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCg < 2:597339/140‑1 (MQ=255) tgATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCg < 1:553350/140‑1 (MQ=255) tgATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCg < 2:262829/140‑1 (MQ=255) ttCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTc > 1:779785/1‑134 (MQ=255) ttCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTc < 2:779785/134‑1 (MQ=255) ggCAACGCATCGGCAAAT‑GCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTaa < 1:971654/140‑1 (MQ=255) aCGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAAt > 1:1037839/1‑140 (MQ=255) aaaT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTg < 2:990051/140‑1 (MQ=255) aaT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGc > 2:478314/1‑140 (MQ=255) aT‑GCGTGTGTGTAATCCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCt < 1:241741/140‑1 (MQ=255) agccGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCt < 1:269105/137‑1 (MQ=255) gtgtgtAATCCCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTg < 2:949425/136‑1 (MQ=255) gtgtgtAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTg > 1:949425/1‑136 (MQ=255) | AACGGCAAACGGTGGCTGGAAGCGACACAAAAAGTGACCTGCGGATAGCGTTTCGCCGCCGCTTTTATCGCTACCTTCGAAACATCCAGACCAAACGTGGTGATTTCGGGCAACGCATCGGCAAAT‑GCGTGTGTGTAATACCCTTCACCACAGCCAATATCCAGCACCGCCGTGGCCTTATCATCAAGCCGTTCCCTCAGTTGGGCGACAATTGCATCACGCAGCGGCTGATAATGTCCGGCATCTAAGAATGCGCGGCGTGCTTG > NC_000913/1906630‑1906895 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |