Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 704,276 | T→C | 9.1% | L123L (CTA→CTG) | ymcA ← | conserved hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 704,276 | 0 | T | C | 9.1% | 153.8 / 6.8 | 55 | L123L (CTA→CTG) | ymcA | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (28/22); new base C (3/2); total (31/24) | |||||||||||
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.43e-01 |
GATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGC > minE/704215‑704346 | gATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTaa < 1:3910614/71‑1 (MQ=255) gATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTaa > 1:102284/1‑71 (MQ=255) aTACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAAc > 1:1144146/1‑71 (MQ=255) aCCGTCGAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACcttctt > 1:532811/1‑71 (MQ=255) aCCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCACACGCCACTTGCGGTAGCCAGTa > 1:3084523/1‑63 (MQ=255) aCCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACcttc < 1:3864751/69‑1 (MQ=255) ccGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGc < 1:4122212/57‑1 (MQ=255) gTCAAACAGACCCGTACCAGCAATATAACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTcc > 1:3779271/1‑70 (MQ=255) aaaCAGCCCCGTACCGCCGATATCCCGCGCACCAACCGCCACTTGCGGCAGCCAGTAACTCTCTTCCcac > 1:3502649/1‑70 (MQ=255) aaaCAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCca > 1:1997378/1‑69 (MQ=255) aaaCAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCca > 1:2456672/1‑69 (MQ=255) aaCAGCCCCGTACCGCCGATATCCCGCGCACCAACCGCCACTTGCGGCAGCCAGt < 1:1429573/55‑1 (MQ=255) aCAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTaa > 1:622001/1‑56 (MQ=255) cAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGc > 1:1491974/1‑71 (MQ=255) tACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGt < 1:1116272/45‑1 (MQ=255) tACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTATCCTTCTTCCCACAGGCGAAGTTTa < 1:2250353/71‑1 (MQ=255) ccAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTaa > 1:347774/1‑70 (MQ=255) ccAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTaa > 1:2822752/1‑70 (MQ=255) ccAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTaa > 1:2687091/1‑70 (MQ=255) aGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAg < 1:3394286/63‑1 (MQ=255) gCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATc > 1:3844021/1‑71 (MQ=255) aaTATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCg < 1:1360885/70‑1 (MQ=255) atatCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGaa > 1:558766/1‑71 (MQ=255) tatCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAAt < 1:923859/71‑1 (MQ=255) tCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACctt > 1:441709/1‑40 (MQ=255) tCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCca < 1:3641498/47‑1 (MQ=255) tCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCg < 1:3827924/53‑1 (MQ=255) tCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGt > 1:822092/1‑57 (MQ=255) tCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTaaaa < 1:2337680/63‑1 (MQ=255) tCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATg > 1:2516100/1‑70 (MQ=255) aCGTTTACCAAACGCCACTTGCGGTAGCCAGTAAcc > 1:448107/1‑36 (MQ=255) aCGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAtt > 1:3129330/1‑71 (MQ=255) aCGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAtt > 1:3979769/1‑71 (MQ=255) cGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGAttt < 1:1233849/71‑1 (MQ=255) gTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTg > 1:1693256/1‑71 (MQ=255) gTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTg > 1:491916/1‑71 (MQ=255) tACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCca > 1:3164502/1‑40 (MQ=255) tACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTc < 1:650526/70‑1 (MQ=255) tACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCt > 1:4029543/1‑71 (MQ=255) ccAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCg < 1:3318279/44‑1 (MQ=255) cAACCGCCACTTGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGttt < 1:2605070/49‑1 (MQ=16) cAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGa < 1:2675394/57‑1 (MQ=255) cAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTg < 1:266453/71‑1 (MQ=255) aaaCGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCt < 1:3525145/67‑1 (MQ=255) accGCCACTTGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGaa > 1:611367/3‑56 (MQ=255) accGCCACTTGCGGCAGCCAGTAACTCTCTTCCCACAAACGCAGTTTGAGATCGa > 1:1838457/3‑55 (MQ=255) aaCGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTa > 1:666771/1‑71 (MQ=255) aCGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTAt < 1:574692/71‑1 (MQ=255) aCGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTAt < 1:751886/71‑1 (MQ=255) cACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTAt > 1:995098/1‑67 (MQ=255) aCTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGt < 1:1677144/64‑1 (MQ=255) aCTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTg < 1:194918/71‑1 (MQ=255) aCTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTg < 1:3033067/71‑1 (MQ=255) cTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGa > 1:4055974/1‑71 (MQ=255) gCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTa > 1:1718357/1‑61 (MQ=255) tAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATcgccgc > 1:2069970/1‑71 (MQ=255) | GATACTCACCGTCAAACAGACCCGTACCAGCAATATCACGTTTACCAAACGCCACTTGCGGTAGCCAGTAACCTTCTTCCCACAGGCGAAGTTTAAAATCGAATGATTTGTCTTTGTATGACTGATCGCCGC > minE/704215‑704346 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |