Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,810,641 | C→G | 11.4% | C335W (TGC→TGG) | thrS → | threonine‑‑tRNA ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,810,641 | 0 | C | G | 11.4% | 103.5 / 4.0 | 35 | C335W (TGC→TGG) | thrS | threonine‑‑tRNA ligase |
Reads supporting (aligned to +/- strand): ref base C (13/18); new base G (2/2); total (15/20) | |||||||||||
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.93e-01 |
CTACCAGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGACGATGC > NC_002947/2810500‑2810784 | cTACCAGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCa < 1:42381/150‑1 (MQ=255) cactaccaaaTCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAg < 2:258038/143‑1 (MQ=255) ccAGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGgt < 2:732195/150‑1 (MQ=255) ccAGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGgt > 1:60437/1‑150 (MQ=255) aGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGgtgt < 1:93442/150‑1 (MQ=255) tCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACTGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTAcc > 2:740552/1‑150 (MQ=255) gCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGAcc > 1:736072/1‑150 (MQ=255) aTCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTg > 2:1067843/1‑150 (MQ=255) ccTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGcc < 2:1017228/150‑1 (MQ=255) gggAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGc < 1:398737/150‑1 (MQ=255) ggCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAg > 1:102586/1‑150 (MQ=255) gCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGt < 2:904734/150‑1 (MQ=255) cACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAAATGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTc > 1:145943/1‑150 (MQ=255) cTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTAc < 2:861248/118‑1 (MQ=255) cTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTAc > 1:861248/1‑118 (MQ=255) ccAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGcc > 2:147359/1‑124 (MQ=255) ccAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGAATACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGcc < 1:147359/124‑1 (MQ=255) tATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCAc > 1:611918/1‑150 (MQ=255) tCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCg < 1:513310/150‑1 (MQ=255) ttCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCt > 1:445672/1‑150 (MQ=255) gtcgGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCa > 2:1176076/1‑150 (MQ=255) ccGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGcatcat < 1:162366/150‑1 (MQ=255) aaaCCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTc < 1:1176076/150‑1 (MQ=255) aaaCCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTc < 2:736072/150‑1 (MQ=255) cGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTAc > 1:1134989/1‑150 (MQ=255) gATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACa < 2:256982/150‑1 (MQ=255) aTGAACTGCCCGGGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAg < 2:145943/150‑1 (MQ=255) gAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGa < 2:1134989/150‑1 (MQ=255) accTGCACGTGGCACGGGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGgcgc > 1:491938/3‑110 (MQ=255) accTGCACGTGGCACGGGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGgcgc < 2:491938/108‑1 (MQ=255) accTGCACGTGGCACGGGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGAc > 2:458485/3‑150 (MQ=255) tGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGACGAt < 2:611918/150‑1 (MQ=255) cccGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGAc > 2:658149/1‑49 (MQ=255) cccGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGAc < 1:658149/49‑1 (MQ=255) cacGTGGCACGGGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGACGATgc < 1:458485/148‑1 (MQ=255) | CTACCAGGAAATCAAGACCCCGCAGGTTGTCGACCGCATCCTGTGGGAGCGTTCCGGCCACTGGTCCAACTATGCCGAAAACATGTTCACCACTTCGTCGGAAAGCCGTGACTACGCGGTAAAACCGATGAACTGCCCGTGCCACGTGCAGGTGTTCAACCAGGGCCTGAAAAGCTACCGCGACCTGCCGCTGCGTCTCGCCGAGTTCGGTGCCTGCCACCGTAACGAGCCGTCCGGCGCCCTGCACGGCATCATGCGCGTGCGTGGCTTCGTACAGGACGATGC > NC_002947/2810500‑2810784 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |