Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 4,333,478 | G→A | 54.5% | R555H (CGC→CAC) | yjdA → | conserved hypothetical protein with nucleoside triphosphate hydrolase domain |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 4,333,478 | 0 | G | A | 54.5% | ‑1.6 / 26.8 | 22 | R555H (CGC→CAC) | yjdA | conserved hypothetical protein with nucleoside triphosphate hydrolase domain |
Reads supporting (aligned to +/- strand): ref base G (5/5); new base A (5/7); total (10/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01 |
CGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTT > W3110S.gb/4333411‑4333541 | cgcgACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGccc > 1:14430/1‑71 (MQ=255) gcgACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGcccc < 1:839035/71‑1 (MQ=255) gcgACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGcccc < 1:1080741/71‑1 (MQ=255) cgACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCACCCCa < 1:438992/71‑1 (MQ=255) aCTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCAtt < 1:1250456/71‑1 (MQ=255) tACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAg > 1:92868/1‑70 (MQ=255) aCGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCACCCCATTGAGc < 1:1704940/70‑1 (MQ=255) gTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGt < 1:430833/70‑1 (MQ=255) tttGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCa > 1:895765/1‑70 (MQ=255) cGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATgg < 1:1735596/71‑1 (MQ=255) cTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCACCCCATTGAGCAACAAGTCAAAATgg < 1:187682/67‑1 (MQ=255) cTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCACCCCATTGAGCAACAAGTCAAAATgg < 1:804204/67‑1 (MQ=255) ttGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCACCCCATTGAGCAACAAGTCAAAATGGAAtt > 1:749674/1‑70 (MQ=255) tGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAAtt > 1:1707172/1‑69 (MQ=255) gAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCACCCCATTGAGCAACAAGTCAAAATGGAATTGAg > 1:1203924/1‑71 (MQ=255) tttGCCGGGTGTTAACCGATGTTATTCACCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAggg > 1:1381501/1‑71 (MQ=255) ggTGTTAACCGATGTTATTCACCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTc < 1:66303/71‑1 (MQ=255) tgtTAACCGATGTTATTCACCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCt < 1:118859/71‑1 (MQ=255) gtTAACCGATGTTATTCACCCCATTGAGCAACAAGTCAAAATGGAAtt > 1:186559/1‑48 (MQ=255) ttAACCGATGTTATTCACCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTgg > 1:593517/1‑71 (MQ=255) ttATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCAttt > 1:196756/1‑70 (MQ=255) ttATTCACCCCATTGAGCAACAAGTCAAAATGGaa < 1:1083614/35‑1 (MQ=255) | CGCGACTTTACGCTGCGTTTTGCCGAGCTTGAATCCACCCTTTGCCGGGTGTTAACCGATGTTATTCGCCCCATTGAGCAACAAGTCAAAATGGAATTGAGCGAGTCAGGGTTTCGTCCTGGGTTTCATTT > W3110S.gb/4333411‑4333541 |
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: ‑ |