Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 4,196,219 | A→C | 100% | intergenic (‑120/+113) | thiC ← / ← rsd | phosphomethylpyrimidine synthase/regulator of sigma D |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,196,219 | 0 | A | C | 89.7% | 90.0 / 1.8 | 35 | intergenic (‑120/+113) | thiC/rsd | phosphomethylpyrimidine synthase/regulator of sigma D |
Reads supporting (aligned to +/- strand): ref base A (1/3); new base C (15/16); total (16/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.08e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGC > NC_000913/4196089‑4196351 | gttgCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCcgcg < 2:190180/143‑1 (MQ=255) ttgCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTcggcgg > 2:342613/1‑143 (MQ=255) aCATAGCTCAGTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCg < 1:165274/143‑1 (MQ=255) aTAGCTTATTTTTAAAAGTTTTGTTCGTGTTTTGTGAGACGACGGATCAAGCAAGAGACGAGCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGaa < 1:24652/143‑1 (MQ=255) aaaaaGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGtctcc > 1:17417/1‑142 (MQ=255) aaaaaGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACAGTCTCa > 2:139290/1‑143 (MQ=255) aaaGTTAGCGCCGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGc < 1:180602/143‑1 (MQ=255) aaGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCTATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGcc < 1:81846/143‑1 (MQ=255) aGTGAAGCACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCa < 2:355836/143‑1 (MQ=255) ggAGGTGGCTGGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAg < 2:328686/143‑1 (MQ=255) gACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAgg < 1:497534/143‑1 (MQ=255) aGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGa > 1:466289/1‑143 (MQ=255) gacgacGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTTCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGaa > 2:552979/1‑143 (MQ=255) acgacgGATGAAGCAAGAGACGAGCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGaaa < 1:310914/143‑1 (MQ=255) gCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGCCAAGAAATTAGCCCCGCaa > 2:557923/1‑143 (MQ=255) aaGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCCGTTAAGGCACTCCGACACGAACTTAGCCCCGCAAAt > 1:87256/1‑143 (MQ=255) gaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATgggg < 1:541959/143‑1 (MQ=255) tCGCCGCAGGGGGGATGCGATAGCTGAGTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTg < 2:90380/143‑1 (MQ=255) ccgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAAt < 2:218657/143‑1 (MQ=255) ccgcAGGGGCCATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAACTAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAAt < 1:214095/143‑1 (MQ=255) gcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAAAGGTCTCAGCCAGTTAAGGCACTCCGACAAGAACTTAGACCCGCAAATGCGGCATTGAATgc > 2:190459/1‑142 (MQ=255) gtggCGATGCGATAGCAGATCAACTCTTGTCCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAACCGTCCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTaaa < 1:545961/141‑1 (MQ=255) gtAGCCGAGTACCTCTGGTGCCCTTCGCAGGTATTAGCCGGCTCAGGTCCCGCGGAGCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTaa < 2:477540/142‑1 (MQ=255) aTTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCa > 1:76437/1‑143 (MQ=255) aTTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCa > 2:14623/1‑143 (MQ=255) ccTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAgg > 1:407122/1‑142 (MQ=255) ccTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAgg < 2:407122/142‑1 (MQ=255) ccTTCGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTACGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTAAGCGTTAACAGCGCAGAACGCAAGCAgca > 2:514856/1‑141 (MQ=255) cGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCa > 1:517972/1‑96 (MQ=255) cGCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCa < 2:517972/96‑1 (MQ=255) gCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGttt < 1:549867/143‑1 (MQ=255) gCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGttt > 2:457739/1‑143 (MQ=255) gCAGGTATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGttt > 2:164793/1‑143 (MQ=255) gggaTTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGAc < 1:311173/140‑1 (MQ=255) tATTAGCCTGCTCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGAcgc > 1:473223/1‑143 (MQ=255) | GTTGCAGACATAGCTCATTCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATTGAATGTAAATTACGCGTTAACAGCGCAGAACTCAAGCAGGATGTTTGACGC > NC_000913/4196089‑4196351 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |