Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 204,556 | A→G | 100% | intergenic (+156/‑28) | alaV → / → rrlH | tRNA‑Ala/23S ribosomal RNA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 204,556 | 0 | A | G | 87.1% | 133.4 / 7.5 | 76 | intergenic (+156/‑28) | alaV/rrlH | tRNA‑Ala/23S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base A (10/0); new base G (29/37); total (39/37) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.03e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
CTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCT > minE/204503‑204617 | cTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACTCGATGATGGATCGCAAGAAAc > 1:720765/1‑66 (MQ=11) cTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGa < 1:266285/63‑1 (MQ=255) cTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATc > 1:2962856/1‑69 (MQ=255) cTGAACAATGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATc > 1:2256100/1‑69 (MQ=255) aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGtt > 1:496312/1‑67 (MQ=25) aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGtt > 1:388263/1‑67 (MQ=25) aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGtt > 1:2397509/1‑67 (MQ=25) aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGtt > 1:541844/1‑67 (MQ=25) aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGtt > 1:1215931/1‑67 (MQ=25) aaaGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGtt > 1:263186/1‑67 (MQ=25) aGTTGTTTGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTtgtg < 1:1857957/68‑1 (MQ=255) aGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTtgtg < 1:3009114/68‑1 (MQ=25) aGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTtgtg < 1:2310571/68‑1 (MQ=25) aGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTtgtg < 1:415897/68‑1 (MQ=25) tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc < 1:2736356/68‑1 (MQ=25) tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc < 1:737290/68‑1 (MQ=25) tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc < 1:3063870/68‑1 (MQ=25) tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc < 1:3253540/68‑1 (MQ=25) tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc < 1:2324132/68‑1 (MQ=25) tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc < 1:3374924/68‑1 (MQ=25) tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc < 1:1558341/68‑1 (MQ=25) tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc < 1:2435209/68‑1 (MQ=25) tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc < 1:1123118/68‑1 (MQ=25) tGAGTCTCTCAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGc < 1:1346909/68‑1 (MQ=25) tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAg > 1:1059686/1‑69 (MQ=25) tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAg > 1:1909679/1‑69 (MQ=25) tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAg > 1:3354438/1‑69 (MQ=25) tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAg > 1:1038151/1‑69 (MQ=25) tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAg > 1:3229825/1‑69 (MQ=25) tctcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAg > 1:2521328/1‑69 (MQ=25) tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGc > 1:2817025/1‑68 (MQ=25) tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg > 1:3266250/1‑69 (MQ=25) tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg > 1:3334054/1‑69 (MQ=25) tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg > 1:2023376/1‑69 (MQ=25) tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg > 1:1483747/1‑69 (MQ=25) tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg > 1:605410/1‑69 (MQ=25) tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg > 1:2716096/1‑69 (MQ=25) tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg > 1:1055969/1‑69 (MQ=25) tcAAATTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg > 1:2837402/1‑69 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:2901164/67‑1 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:388577/67‑1 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:506467/67‑1 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:597652/67‑1 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:1012858/67‑1 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:2779144/67‑1 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:1206128/67‑1 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:1508595/67‑1 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:1803927/67‑1 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:1027351/67‑1 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:198503/67‑1 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:2029654/67‑1 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:2092389/67‑1 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:2109007/67‑1 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:2323388/67‑1 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa > 1:3203029/1‑69 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa > 1:2128705/1‑69 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa > 1:2955887/1‑69 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa > 1:1877849/1‑69 (MQ=25) aaaTTTTCGCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTa > 1:1917337/1‑69 (MQ=25) ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg > 1:3178112/1‑69 (MQ=35) ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg > 1:114611/1‑69 (MQ=35) ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg > 1:458026/1‑69 (MQ=35) ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg > 1:2958314/1‑69 (MQ=35) ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg > 1:1010527/1‑69 (MQ=35) ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg > 1:1997439/1‑69 (MQ=35) ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg > 1:2130271/1‑69 (MQ=35) ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg > 1:3321832/1‑69 (MQ=35) ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg > 1:938742/1‑69 (MQ=35) ttCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACgg > 1:2083933/1‑69 (MQ=35) gCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg < 1:947607/68‑1 (MQ=25) gCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg < 1:844293/68‑1 (MQ=25) gCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg < 1:981/68‑1 (MQ=25) gCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg < 1:449054/68‑1 (MQ=25) gCAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg < 1:2137363/68‑1 (MQ=25) cAACACGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg < 1:2066169/67‑1 (MQ=25) acacGATGATGGATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTg < 1:2843214/65‑1 (MQ=25) tgatggATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 1:3290557/50‑1 (MQ=255) tgatgaATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCt > 1:1373339/1‑67 (MQ=35) tgatgaATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCt > 1:357248/1‑67 (MQ=35) | CTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGCAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCT > minE/204503‑204617 |
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: ‑ |