Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,243,741 | G→A | 13.7% | intergenic (+30/+69) | tsaA → / ← rnt | peroxiredoxin/ribonuclease T |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,243,741 | 0 | G | A | 13.7% | 125.6 / 10.2 | 51 | intergenic (+30/+69) | tsaA/rnt | peroxiredoxin/ribonuclease T |
Reads supporting (aligned to +/- strand): ref base G (24/20); new base A (4/3); total (28/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01 |
CACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCG > NC_002947/1243598‑1243888 | cACGAAGAGCACGGCTAAGTGTGCCCTGCCTACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGc < 3:98507/148‑1 (MQ=255) cGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCgg > 8:304054/1‑149 (MQ=255) aCGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGttttttca < 4:256357/149‑3 (MQ=39) gCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGc > 8:197194/1‑149 (MQ=255) cAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGcc > 3:298294/1‑148 (MQ=255) gAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGcc < 8:35924/138‑1 (MQ=255) aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt < 6:264822/123‑1 (MQ=255) aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTAt > 5:264822/1‑123 (MQ=255) aaaaaGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCgg < 3:17802/149‑1 (MQ=255) gACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGccc > 4:285251/1‑149 (MQ=255) ggCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAAt > 1:158431/1‑149 (MQ=255) gCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTATTCGCGGGCGCGCCCGCGaa > 6:255544/1‑147 (MQ=255) cATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTATTCGCGGGCGCGCCCGCGaagag < 5:255544/149‑4 (MQ=255) tGAACGGTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTTTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCt < 7:299941/149‑1 (MQ=255) gAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCtt > 2:50103/1‑149 (MQ=255) aaCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCtt < 7:304054/148‑1 (MQ=255) aCGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCtt > 8:62287/1‑147 (MQ=255) cGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg < 4:298294/149‑1 (MQ=255) gCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg > 1:276112/1‑148 (MQ=255) ttCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTTCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAg > 5:246538/1‑146 (MQ=255) tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCa > 6:219384/1‑142 (MQ=255) tGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCa < 5:219384/142‑1 (MQ=255) aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATc < 3:285251/148‑1 (MQ=255) aaGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCg < 7:197194/149‑1 (MQ=255) tgctgcCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACGGGCCTTCGTGGCCGTTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGaaa < 1:50103/148‑1 (MQ=255) cTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgcc > 8:129113/1‑149 (MQ=255) tGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGACCGCGAATGCCTTTAGCAGATCAATCGTTGAAATAGCGCCAgccg > 7:285056/1‑149 (MQ=255) gAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccgc > 2:297274/1‑149 (MQ=255) cGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCAt > 1:60717/1‑68 (MQ=17) cGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCAt < 2:60717/68‑1 (MQ=17) gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAAt > 7:253533/1‑124 (MQ=255) gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAAt < 8:253533/124‑1 (MQ=255) gagaACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCATAAGCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAgccgcc > 1:108291/1‑148 (MQ=255) cGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctt > 8:75362/1‑149 (MQ=255) cTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctttc < 2:276112/149‑1 (MQ=255) tAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa < 7:62287/148‑1 (MQ=255) aaTTGCTGCGTGCATGAAAAAAGCGGCCTTCGTGGCGGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGGAGCTCAATGGTTGAAATCGCGGCCGCCGCCCATTTGTTTCCAGCGGTTGAc > 7:53452/1‑148 (MQ=255) ttGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCGGATCAATCGtt > 1:27697/1‑107 (MQ=255) ttGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGtt < 2:27697/107‑1 (MQ=255) gcGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCg > 2:215914/1‑148 (MQ=255) cGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGc < 2:158431/148‑1 (MQ=255) cGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCa < 5:56570/149‑1 (MQ=255) cGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCCAATGCCATTAGCAGATAAATAGTTTAAATACCGCCAGCCGCCCATTTCATTACAGCGGTTGACGATGccccc > 6:56570/1‑146 (MQ=255) gCATAAAAAAACCGGCCACGAAGGCCGGTTTTTTCATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa > 3:157922/1‑149 (MQ=38) ataaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa > 2:67314/4‑147 (MQ=255) gAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCAtttctttc > 4:195099/1‑122 (MQ=255) aaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa < 3:324684/144‑1 (MQ=255) aaaaaaaCCGGCCACGAAGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGaa > 4:324684/2‑145 (MQ=255) aaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa < 4:10366/130‑1 (MQ=255) aaaaaCCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGa > 3:10366/1‑130 (MQ=255) cGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTtctc < 7:75362/148‑1 (MQ=255) gTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCg < 7:129113/148‑1 (MQ=255) | CACGAAGAGCACGGCGAAGTGTGCCCTGCCAACTGGAAAAAAGGCGACAAGGGCATGAACGCTTCGCCTGAAGGCGTTGCTGCCTACCTGAGCGAGAACGCTGGCAAGCTGTAATTGCTGCGTGCATGAAAAAACCGGCCTTCGTGGCCGGTTTTTTTATGCCT‑‑TCCTGTGCCGGCCTCTTCGCGGGCGCGCCCGCGAATGCCTTTAGCAGATCAATCGTTGAAATCGCGCCAGCCGCCCATTTCTTTCCAGCGGTTGACGATGCCGCAGAACAGCTCAGCGGTTTTCTCG > NC_002947/1243598‑1243888 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |