Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 880,753 | G→A | 51.1% | intergenic (+72/‑122) | yciM → / → pyrF | conserved hypothetical protein/orotidine‑5'‑phosphate decarboxylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 880,753 | 0 | G | A | 51.1% | ‑6.2 / 66.4 | 45 | intergenic (+72/‑122) | yciM/pyrF | conserved hypothetical protein/orotidine‑5'‑phosphate decarboxylase |
Reads supporting (aligned to +/- strand): ref base G (15/7); new base A (11/12); total (26/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.31e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01 |
TTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCTCGC > minE/880684‑880822 | tttAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGgc > 1:6679/1‑71 (MQ=255) tttAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGgc > 1:720363/1‑71 (MQ=255) tttAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGgc > 1:2532717/1‑71 (MQ=255) tttAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGgc > 1:3992169/1‑71 (MQ=255) aaaTCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGaaaa < 1:2650354/71‑1 (MQ=255) tCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGc < 1:686069/71‑1 (MQ=255) tCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGc < 1:891489/71‑1 (MQ=255) cTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGcc < 1:915358/71‑1 (MQ=255) cTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGcc < 1:137728/71‑1 (MQ=255) tAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAg < 1:1181153/70‑1 (MQ=255) caacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGgaga < 1:1864385/70‑1 (MQ=255) caacaTACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGgaga < 1:408269/70‑1 (MQ=255) aacaTACTAATTATTAATGTTCCGTTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGaaa > 1:4045672/1‑71 (MQ=255) tACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGc < 1:2730455/71‑1 (MQ=255) cTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa > 1:3902355/1‑71 (MQ=255) cTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa > 1:193302/1‑71 (MQ=255) cTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa > 1:1636740/1‑71 (MQ=255) ttattaATGTTCCATTGTGCTCCGGCAACGACGACGCAGaa < 1:547236/41‑1 (MQ=255) ttattaATGTTCCATTGTGCTCCGGCAACGACGACGCAGaa < 1:3135587/41‑1 (MQ=255) ttattaATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGgaga < 1:2655695/59‑1 (MQ=255) ttattaATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGgaga < 1:3866670/59‑1 (MQ=255) ttattaATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGc < 1:3344099/65‑1 (MQ=255) attaATGTTCCATTGTGCTCCGGCAACGACGACGCAGaaaa > 1:914953/1‑41 (MQ=255) attaATGTTCCATTGTGCTCCGGCAACGACGACGCAGaaaa > 1:3653553/1‑41 (MQ=255) attaATGTTCCATTGTGCTCCGGCAACGACGACGCAGaaaa > 1:1929480/1‑41 (MQ=255) attaATGTTCCATTGTGCTCCGGCAACGACGACGCAGaaaa > 1:2060598/1‑41 (MQ=255) taATGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGc < 1:431435/61‑1 (MQ=255) aTGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa > 1:685277/1‑61 (MQ=255) aTGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa > 1:1090347/1‑61 (MQ=255) aTGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa > 1:752537/1‑61 (MQ=255) aTGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa > 1:542723/1‑61 (MQ=255) tGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAAtttt > 1:1178124/1‑71 (MQ=255) tGTTCCATTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGCCA‑GGGAGAAATCGCAACTGTTAAttttt > 1:50430/1‑71 (MQ=255) ttCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa > 1:2101127/1‑58 (MQ=255) ttCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa > 1:3008136/1‑58 (MQ=255) ccATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTAt > 1:2194931/1‑71 (MQ=255) aTTGTGCTCCGGCAACGACGACGCAGAAAAAGCCTGcc < 1:746297/38‑1 (MQ=255) ttGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGc < 1:430226/51‑1 (MQ=255) ttGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGt < 1:1969786/57‑1 (MQ=255) tgtgCTCCGGCAACGACGACGCAGAAAAAGCCTGcc < 1:3320028/36‑1 (MQ=255) gtgCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCa > 1:924570/1‑71 (MQ=255) gtgCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCa > 1:925021/1‑71 (MQ=255) gtgCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCa > 1:4080046/1‑71 (MQ=255) gtgCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCa > 1:3122385/1‑71 (MQ=255) gtgCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCa > 1:1221110/1‑71 (MQ=255) gacgGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCaa > 1:2281718/1‑38 (MQ=255) gacgGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACg > 1:3263852/1‑60 (MQ=255) gacgGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACg > 1:2522280/1‑60 (MQ=255) gacgGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCt > 1:3520444/1‑71 (MQ=255) ggCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCTCGc < 1:619196/71‑1 (MQ=255) | TTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATGTTCCATTGTGCTCCGGCAACGACGGCGCAGAAAAAGCCTGCCAGGGGAGAAATCGCAACTGTTAATTTTTTATTTCCACGGGTAGAATGCTCGC > minE/880684‑880822 |
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: ‑ |