Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 2,530,398 | G→A | 11.1% | A26T (GCA→ACA) | mutM → | formamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,530,398 | 0 | G | A | 11.1% | 98.2 / 4.3 | 36 | A26T (GCA→ACA) | mutM | formamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase |
Reads supporting (aligned to +/- strand): ref base G (23/9); new base A (2/2); total (25/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.70e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01 |
ACCCGAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGACC > minE/2530334‑2530461 | aCCCGAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAgtg < 1:1715374/70‑1 (MQ=255) cccGAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATg > 1:2180103/1‑64 (MQ=255) gAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTgcg > 1:835482/1‑71 (MQ=255) gTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCaa > 1:388565/1‑71 (MQ=255) gTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCaa > 1:2007260/1‑71 (MQ=255) ttGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATg < 1:1430973/57‑1 (MQ=255) gAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAAc > 1:1652547/1‑69 (MQ=255) ccAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGc > 1:1068827/1‑71 (MQ=255) aGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAgtggtg > 1:171172/1‑57 (MQ=255) aGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCtt > 1:1144615/1‑71 (MQ=255) ccgccgCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTgc < 1:1975675/71‑1 (MQ=255) cgccgcGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCtt > 1:687685/1‑68 (MQ=255) cgccgcGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTgcg < 1:230877/71‑1 (MQ=255) ccgcGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATg < 1:2011383/44‑1 (MQ=255) ccgcGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATg < 1:498137/44‑1 (MQ=255) ccgcGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATg < 1:82248/44‑1 (MQ=255) gCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTgcgc < 1:1015314/65‑1 (MQ=255) gAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAgtgg > 1:1946299/1‑40 (MQ=255) aCCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCag < 1:1078744/71‑1 (MQ=255) ccGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCag > 1:1949386/1‑70 (MQ=255) cGCATCTCGTTGGTGCAACCATTCTTCATACAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaa < 1:1997424/71‑1 (MQ=255) cATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaaga < 1:1759835/71‑1 (MQ=255) cATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaaga < 1:544075/71‑1 (MQ=255) tcGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaa > 1:2050360/1‑65 (MQ=255) tcGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaa > 1:400398/1‑65 (MQ=255) gTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTg < 1:495068/50‑1 (MQ=255) ggTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCg > 1:394937/1‑71 (MQ=255) gCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGtt > 1:2107962/1‑70 (MQ=255) gCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGtt > 1:356215/1‑70 (MQ=255) aCCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAg > 1:1901802/1‑71 (MQ=255) ccATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGc > 1:541717/1‑71 (MQ=255) ccATTCTTCATACAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGc < 1:1407520/71‑1 (MQ=255) cATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAgaga > 1:1169710/1‑56 (MQ=255) cATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAgaga > 1:1702270/1‑56 (MQ=255) cATTCTTCATACAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGc > 1:870683/1‑70 (MQ=255) cATTCTTCATACAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGc > 1:2172257/1‑70 (MQ=255) ttcttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAagag < 1:1777135/53‑1 (MQ=255) ttcttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGac > 1:329472/1‑71 (MQ=255) ttcttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGac > 1:1811304/1‑71 (MQ=255) cttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCagaa > 1:1553621/1‑48 (MQ=255) cttcATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGacc > 1:444454/1‑70 (MQ=255) | ACCCGAAGTTGAAACCAGCCGCCGCGGCATAGAACCGCATCTCGTTGGTGCAACCATTCTTCATGCAGTGGTGCGCAACGGACGCTTGCGCTGGCCGGTTTCAGAAGAGATCTACCGTTTAAGCGACC > minE/2530334‑2530461 |
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: ‑ |