Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,834,068 | G→A | 22.2% | L56L (CTG→CTA) | ygbM → | conserved hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,834,068 | 0 | G | A | 22.2% | 66.9 / 15.8 | 36 | L56L (CTG→CTA) | ygbM | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (22/6); new base A (5/3); total (27/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.84e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCTGTTTCCCTATAACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGAC > minE/1834002‑1834135 | tCTGTTTCCCTATAACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTAgcg < 1:362209/70‑1 (MQ=255) tttCCCTATAACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGtt < 1:2493620/71‑1 (MQ=255) cccTATAACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCaa > 1:3361306/1‑71 (MQ=255) ccTATAACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTAGCGCTGTTCAAt < 1:76116/71‑1 (MQ=255) tataACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTg > 1:2962447/1‑57 (MQ=255) tataACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATAc > 1:3040470/1‑71 (MQ=255) tataACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTAGCGCTGTTCaa < 1:3357274/68‑1 (MQ=255) taACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTAGCGCTGTTCAATa > 1:135315/1‑68 (MQ=255) tactCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTcc > 1:3701500/1‑71 (MQ=255) tactCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTcc > 1:3544837/1‑71 (MQ=255) tCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTc < 1:2264087/57‑1 (MQ=255) cACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCaa > 1:101734/1‑57 (MQ=255) cACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCg > 1:3272279/1‑68 (MQ=255) cTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGatat > 1:1187577/1‑71 (MQ=255) cAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGgc > 1:877730/1‑41 (MQ=255) aaaTCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAAt > 1:2414328/1‑71 (MQ=255) aaaTCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAAt > 1:3314554/1‑71 (MQ=255) aTCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGc > 1:1974493/1‑71 (MQ=255) ccAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTa > 1:968820/1‑65 (MQ=255) ccAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTa > 1:1993020/1‑65 (MQ=255) aaGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTaa < 1:466075/62‑1 (MQ=255) aaGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTaa < 1:3028301/62‑1 (MQ=255) aaGCAACTGGAGCAAAACCATCTGACACTAGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCgggg > 1:3360936/1‑69 (MQ=255) aGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGaa > 1:3848287/1‑71 (MQ=255) gAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGAt > 1:504514/1‑70 (MQ=255) gCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATAtt < 1:423225/49‑1 (MQ=255) cAAAACCATCTGACACTAGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATc > 1:1122930/1‑71 (MQ=255) aaaaCCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATcc < 1:3130588/71‑1 (MQ=255) aaaCCATCTGACACTAGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCg > 1:288238/1‑71 (MQ=255) aaaCCATCTGACACTAGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCg > 1:2863498/1‑71 (MQ=255) aCCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCg > 1:2885087/1‑50 (MQ=255) aTCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCtt > 1:2183903/1‑71 (MQ=255) cTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTcc > 1:3901482/1‑71 (MQ=255) cTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTcc > 1:989089/1‑71 (MQ=255) cacTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGa > 1:2211168/1‑71 (MQ=255) acTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGAc > 1:3270243/1‑71 (MQ=255) cTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGc > 1:299797/1‑59 (MQ=255) | TCTGTTTCCCTATAACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGAC > minE/1834002‑1834135 |
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: ‑ |