Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,033,631 | C→T | 11.7% | intergenic (+18/+26) | idi → / ← lysS | isopentenyl diphosphate isomerase/lysine tRNA synthetase, constitutive |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,033,631 | 0 | C | T | 11.7% | 143.4 / 9.1 | 51 | intergenic (+18/+26) | idi/lysS | isopentenyl diphosphate isomerase/lysine tRNA synthetase, constitutive |
Reads supporting (aligned to +/- strand): ref base C (21/24); new base T (4/2); total (25/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.19e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01 |
TGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAA > NC_000913/3033488‑3033770 | tgtgATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCgggg > 4:169784/1‑149 (MQ=255) gtgATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGt > 4:27328/1‑149 (MQ=255) aTTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTtgt < 2:670007/149‑1 (MQ=255) aTTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTtgt < 3:663959/149‑1 (MQ=255) tttAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGGCAAATGTCGGGGttttt < 4:14782/149‑4 (MQ=255) aGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGGCAAATGTCGGGGttt > 2:602022/1‑140 (MQ=255) aGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGGCAAATGTCGGGGttt < 1:602022/141‑2 (MQ=255) gATGTATTACACGGTATTGATGGCACGCCGTGGGCGTTCAGTCCGTGGATGTTGATGCATGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACTCAGCTTAGATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATaa > 4:581962/1‑149 (MQ=255) gTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGt < 1:651123/149‑1 (MQ=255) cGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTa < 2:295923/149‑1 (MQ=255) cGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTa < 4:661155/149‑1 (MQ=255) gCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGa < 4:359655/149‑1 (MQ=255) ggCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCgg < 4:361208/149‑1 (MQ=255) ggCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTACCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCgg < 2:635238/149‑1 (MQ=255) gCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCggg > 4:442070/1‑149 (MQ=255) tCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAAt > 1:66395/1‑149 (MQ=255) cGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATaa < 3:169784/149‑1 (MQ=255) gCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGt > 1:569639/1‑149 (MQ=255) gCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGt > 1:327933/1‑149 (MQ=255) ggCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATg > 1:14348/1‑149 (MQ=255) ggCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATg < 4:421129/149‑1 (MQ=255) ggCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATg > 2:232273/1‑149 (MQ=255) aaaTCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGtt < 2:14348/149‑1 (MQ=255) aaTCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTg < 3:442070/149‑1 (MQ=255) cgcgAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTg < 2:327933/149‑1 (MQ=255) aGCCAGAAAACGATTATCTGCATTTACCCAGCATAAATAAAAAAACCCCGACATTTGCCGGAGTTGTGAGCATAACGTAATGCTTCTTTTACCGGACGCATCGCTGGGAACACAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAg < 3:581962/149‑1 (MQ=255) aGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAtt > 2:229816/1‑149 (MQ=255) aGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCAtt > 2:322800/1‑149 (MQ=255) gAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTg > 3:352769/1‑136 (MQ=255) gAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTg < 4:352769/136‑1 (MQ=255) gAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTa < 1:229816/149‑1 (MQ=255) gATTATCTGCATTTACCCAGCTTAAATAAAAAAACCTCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATg < 3:528285/118‑1 (MQ=255) gATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATg > 4:528285/1‑118 (MQ=255) ttATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACgg < 4:526694/149‑1 (MQ=255) tATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGt < 1:322800/149‑1 (MQ=255) tACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTAACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAAt < 3:607493/149‑1 (MQ=255) aaataaaAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACct < 4:251949/149‑1 (MQ=255) ataaaAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACctgc > 1:630770/1‑149 (MQ=255) ataaaAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACct > 3:288488/1‑147 (MQ=255) ataaaAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACct < 4:288488/147‑1 (MQ=255) aaaaaaaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATa > 4:557195/1‑121 (MQ=255) aaaaaaaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATa < 3:557195/121‑1 (MQ=255) tcacaaCCCCGGCAAATGTCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGt > 2:627704/5‑149 (MQ=255) aaaaaaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAc > 1:306203/1‑142 (MQ=255) aaaaaaCCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGAc < 2:306203/142‑1 (MQ=255) caaCCCCGGCAAATGTCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGTTCGATACCAATTCCCAGACCTGCTGTagg > 1:27278/2‑149 (MQ=255) aaCCCCGGCAAATGTCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTcggc > 4:155827/1‑149 (MQ=255) cGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTaa > 4:66901/1‑149 (MQ=255) cGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTaa > 2:603504/1‑149 (MQ=255) cGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTaa > 2:474461/1‑149 (MQ=255) cGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTaa > 1:182813/1‑149 (MQ=255) | TGTGATTTAGCAGATGTATTACACGGTATTGATGCCACGCCGTGGGCGTTCAGTCCGTGGATGGTGATGCAGGCGACAAATCGCGAAGCCAGAAAACGATTATCTGCATTTACCCAGCTTAAATAAAAAAACCCCGACATTTGCCGGGGTTGTGAGCATAACGTAATGCTTATTTTACCGGACGCATCGCCGGGAACAGAATAACGTCGCGGATGGTATGGCTGTTGGTGAACAGCATTACCATACGGTCGATACCAATTCCCAGACCTGCTGTCGGCGGTAA > NC_000913/3033488‑3033770 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |