Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,325,522 | T→G | 21.2% | K71T (AAA→ACA) | wcaL ← | predicted glycosyl transferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,325,522 | 0 | T | G | 21.2% | 34.9 / 21.7 | 33 | K71T (AAA→ACA) | wcaL | predicted glycosyl transferase |
Reads supporting (aligned to +/- strand): ref base T (15/11); new base G (2/5); total (17/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.25e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01 |
AGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCAT > minE/1325453‑1325590 | aGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGtt < 1:1899144/70‑1 (MQ=255) aGCGCCTGCCATGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCa < 1:178423/70‑1 (MQ=255) cTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCa < 1:1440669/65‑1 (MQ=255) tGCGAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCa < 1:313614/64‑1 (MQ=255) aGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGTa < 1:914505/71‑1 (MQ=255) aGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGTa < 1:1737922/71‑1 (MQ=255) tATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGt > 1:307799/1‑71 (MQ=255) tATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGTAGGt > 1:488094/1‑71 (MQ=255) tttttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGccc > 1:1524933/1‑63 (MQ=255) tttttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGtt < 1:348701/70‑1 (MQ=255) tttttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTc < 1:2119053/71‑1 (MQ=255) tttttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTc < 1:1587248/71‑1 (MQ=255) tttttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTc < 1:1522916/71‑1 (MQ=255) tttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAg > 1:586379/1‑65 (MQ=255) tttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCg < 1:11566/70‑1 (MQ=255) tttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGt > 1:1809817/1‑71 (MQ=255) tttACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTGTCGCCACTTTGCCCGTAGGTTCGt > 1:1199108/1‑71 (MQ=255) gAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAAc > 1:90659/1‑71 (MQ=255) gCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGt < 1:390323/59‑1 (MQ=255) cgcgCAAGGTCTGGCTGGCTCGGTGGCTCAGTTTCGCCACTTTGc > 1:1008760/1‑45 (MQ=255) cgcgCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTcc > 1:153726/1‑59 (MQ=255) cgcgCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTcc > 1:698547/1‑59 (MQ=255) cgcgCAAGGTCTGGCTGGCTCGGTCGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCgg > 1:332743/1‑71 (MQ=255) aGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTg > 1:1652531/1‑38 (MQ=255) aGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGtt > 1:227277/1‑68 (MQ=255) gTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCg < 1:1163386/62‑1 (MQ=255) ctggctggCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGgcag > 1:750795/1‑71 (MQ=255) tggctCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTg < 1:1890501/62‑1 (MQ=255) ggTGGCGCAGTGTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTa < 1:266286/70‑1 (MQ=255) gcgcAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCg > 1:1887256/1‑71 (MQ=255) gcgcAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCg > 1:407504/1‑71 (MQ=255) cgcAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGt > 1:851000/1‑71 (MQ=255) cAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTcc < 1:582727/71‑1 (MQ=255) gTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCAt < 1:1256637/71‑1 (MQ=255) | AGGTTGAGCGCCTGCCAGGTATTTTTACGATGAATGCCGCGCAAGGTCTGGCTGGCTCGGTGGCGCAGTTTCGCCACTTTGCCCGTAGGTTCGTCCTGTAACCAGCGGGTTCTGGCAGCAAGGTTGTATTTCGTCCAT > minE/1325453‑1325590 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |