Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 786,614 | T→G | 10.3% | C328G (TGC→GGC) | aroG → | 3‑deoxy‑7‑phosphoheptulonate synthase, Phe‑sensitive |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 786,614 | 0 | T | G | 10.3% | 80.7 / 2.5 | 35 | C328G (TGC→GGC) | aroG | 3‑deoxy‑7‑phosphoheptulonate synthase, Phe‑sensitive |
Reads supporting (aligned to +/- strand): ref base T (18/13); new base G (2/2); total (20/15) | |||||||||||
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.93e-01 |
TGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCA > NC_000913/786472‑786755 | tgTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCt < 1:666245/143‑1 (MQ=255) gTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTg > 1:114674/1‑143 (MQ=255) gTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCCTGATGGTGGAAAGCCATAGGGCGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGGGCATCACCGATGCCTg > 2:169753/1‑143 (MQ=255) gACCTTTGCCATCATATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTg < 1:130024/143‑1 (MQ=255) aCGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTgg < 1:530347/143‑1 (MQ=255) cGTGTGCCAGCATATTGCCGGTGGCGAAAAGGCCATCATTGGCGTGATTTTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTggg < 1:701668/143‑1 (MQ=255) gcagATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCCTGGAAGATACCg > 1:608212/1‑143 (MQ=255) aTTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAGGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGc < 1:719820/143‑1 (MQ=255) cGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGt > 2:707283/1‑143 (MQ=255) cGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACAGTCACAGAATCCGTGATGCCTGTATCGGCTGGGAAGATACCGATGCTCTGt < 1:749181/143‑1 (MQ=255) cGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCAGCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGt > 2:81537/1‑143 (MQ=255) ggTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGCGAGCGGGGAGCCGCTGGCCTACGGTGAGAGCCTCACCGATGCCTGCATCGGCTGGGGAGATACCGATGctgcgtt > 1:328757/1‑138 (MQ=255) gTGGCGAAAAGGCCATTATTGGGGTGAGGGGGGAAAGCCATCTGGTGGAAGGCAAGCAGAGCCTCGAGAGCGGGGAGCCGCTGGGCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTa < 1:240262/143‑1 (MQ=255) cGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCAGGCATCggtgat > 2:191566/1‑114 (MQ=255) cGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCAGGCATCggtgat < 1:191566/118‑5 (MQ=255) gAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCCCCGATGCCTGCATCGGCTGGCGAGATACCGATGCTCTGTTACGTCa > 2:404810/1‑143 (MQ=255) cATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGa > 1:224705/1‑143 (MQ=255) tggAGAGCCATCGGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCCCCGATGCCGGCATCGGCTGGGAAGATACCGATGCTTTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGt < 2:492004/143‑1 (MQ=255) gAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAACGCGTcg > 2:543222/1‑143 (MQ=255) gAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCAGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTcg < 1:214949/143‑1 (MQ=255) aGCCAGCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCgg < 1:694163/143‑1 (MQ=255) cTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTCATTGTCGGATGCGCCGt > 1:565900/1‑143 (MQ=255) agagCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGCAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCaga < 2:403257/143‑1 (MQ=255) cGGGGACCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGACACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTgg < 1:254089/143‑1 (MQ=255) gcggAGCCGCTCGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGCCACCCATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGc < 2:348269/141‑1 (MQ=255) ccGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATcc > 1:681645/1‑143 (MQ=255) gCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGAAGCTCTGATACGTCAACTGGCGACTGCAGTAAAAGCGCGTCGCGGGTAAGGACTAATTGTCGGCTGCGCCGGAAGAGAGGCGtcagcga > 2:102592/1‑137 (MQ=255) cTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATc > 1:322152/1‑143 (MQ=255) aCGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTGAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATcac > 1:384882/1‑143 (MQ=255) agagCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATcaccaaccg > 1:566742/1‑139 (MQ=255) cATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGGAAGGGTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGCATCACCACAGGCCTg > 1:146146/1‑143 (MQ=255) agccGATGCAGGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATa > 2:712939/3‑143 (MQ=255) ccGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAg < 2:224705/143‑1 (MQ=255) aTGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGCCGAATGCAGTAAAAGCGCGTCGCGGGTAAGCATTAATTGTCGGATGCGCCGGCAGCGTTGCGTATCCGATGAATCACCACAGGCCTGATAAGAcg > 2:770633/1‑143 (MQ=255) tcTGCTTCGGCTGGGAAGAGACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGCGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTCGCGTATTCGATTAATCACCACAGGCCTGATAAGTCgcgc < 2:367319/142‑1 (MQ=255) cTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAACGCGCGTCGCGGGTACGGTTTAATTGTCGGATGCGCCGTCACAGTGGCGCATCCGATCAATCACCCCCGGCCTGACACGCCGCGCa > 2:265588/1‑143 (MQ=255) | TGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGTAAGAGCATCACCGATGCCTGCATCGGCTGGGAAGATACCGATGCTCTGTTACGTCAACTGGCGAATGCAGTAAAAGCGCGTCGCGGGTAAGGTTTAATTGTCGGATGCGCCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCA > NC_000913/786472‑786755 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |