Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,442,524 | C→T | A254V (GCG→GTG) | mnmC → | fused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,442,524 | 0 | C | T | 97.2% | 133.7 / ‑4.0 | 36 | A254V (GCG→GTG) | mnmC | fused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme |
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (18/17); total (19/17) | |||||||||||
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.86e-01 |
ACGCTGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGGG > NC_000913/2442379‑2442668 | aCGCTGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGt > 1:520013/1‑151 (MQ=255) cGCTGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTccgtgccgtg < 1:146887/151‑1 (MQ=255) gCTGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTgg > 2:324863/1‑151 (MQ=255) tGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGtt > 1:52680/1‑151 (MQ=255) cGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTaa < 2:58938/151‑1 (MQ=255) cGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACggg < 2:297727/151‑1 (MQ=255) ccGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGgcagcag < 2:539190/151‑1 (MQ=255) ccGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGgcagcag < 2:100371/151‑1 (MQ=255) tttGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAAc > 1:101816/1‑151 (MQ=255) gTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACggg > 1:572831/1‑151 (MQ=255) cgccgcGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGa > 1:17598/1‑149 (MQ=255) gccgcGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAgc < 1:41782/151‑1 (MQ=255) cgcgGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAgcg < 2:300636/150‑1 (MQ=255) tttGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGAt < 1:429378/150‑1 (MQ=255) cAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATcg < 2:173067/151‑1 (MQ=255) aCGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATcggcgg > 2:233646/1‑151 (MQ=255) cGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATcggcgg > 1:581717/1‑150 (MQ=255) gCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATcggcgg > 2:171118/1‑149 (MQ=255) aTTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTAtt > 1:530735/1‑151 (MQ=255) aTTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTAtt < 1:115813/151‑1 (MQ=255) tCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGAtt > 1:299658/1‑149 (MQ=255) taCGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATcggcggagttg > 2:424911/2‑141 (MQ=255) aaaaaCGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGCGGGATTTGCcgcgcggt > 2:522206/1‑149 (MQ=255) gggCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGCGGGATTGGCAGCGCGTtgttgtt > 2:176963/1‑146 (MQ=255) gcgcAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCATCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTa < 1:569574/151‑1 (MQ=255) cgcAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTAc < 1:321598/151‑1 (MQ=255) gAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCt < 2:350944/151‑1 (MQ=255) aTGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGc < 1:419793/151‑1 (MQ=255) gCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCa < 1:552047/150‑1 (MQ=255) gCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAg < 1:213861/151‑1 (MQ=255) ttGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAgg < 2:508851/149‑1 (MQ=255) ggggTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCCGCGCGGTGTTGTCGCTGGCGCTATTTCGGCGCGGGTGGGAGGGaa > 2:613724/1‑148 (MQ=255) gTGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGCGCGGTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGGAGGGAAAGCttt > 1:10148/1‑151 (MQ=255) tGATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGTCCGCGCGTTGTTTTCGCTGGCGCTATTAACGCGCGGGTGGGAGGGAACGCtt > 1:579219/1‑149 (MQ=255) gATGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTAt < 1:195960/151‑1 (MQ=255) aTGGAACAGACATTACCGCTCCCCTGCTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTTTCAGCGCGTTGTTGTCGCTGGGGCTATTAAGGCGCGGCCTGCCGGGAACGCCTtttt > 1:196465/1‑151 (MQ=255) cTCCGTGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCCGCGCGGTGTTGTTGCCGGGGCTATTACGGCGCGGGCGGGCGGGAAAGCCTTTTTTCGCGGGGGAGGgcccccccccggg > 2:210905/1‑142 (MQ=255) | ACGCTGGCGACATTTACGTCTGCCGGTTTTGTCCGCCGCGGTTTGCAGGACGCCGGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGAGGCCCCCGCACTGGG > NC_000913/2442379‑2442668 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |