| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,591,634 | T→C | 30.7% | intergenic (+93/+96) | yfeS → / ← cysM | conserved hypothetical protein/cysteine synthase B |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,591,634 | 0 | T | C | 30.7% | 95.8 / 76.6 | 101 | intergenic (+93/+96) | yfeS/cysM | conserved hypothetical protein/cysteine synthase B |
| Reads supporting (aligned to +/- strand): ref base T (39/31); new base C (21/10); total (60/41) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.81e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01 | |||||||||||
Only 100 of 103 total aligned reads displayed.
GTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGT > minE/1591569‑1591704 | ttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTaa > 1:2207788‑M1/27‑71 (MQ=255)ttaagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTaa > 1:4081563‑M1/27‑71 (MQ=255)gTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTaa > 1:3292707/1‑71 (MQ=255) ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc > 1:235094/1‑67 (MQ=39) ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTaa < 1:2770226/70‑1 (MQ=255) ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAAt > 1:4599416/1‑71 (MQ=255) caagcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCtaaa > 1:209120‑M1/26‑70 (MQ=255) caagcgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCtaaa < 1:4873076‑M1/46‑2 (MQ=255) tAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAAtt > 1:2591310/1‑71 (MQ=40) tAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAAtt > 1:2474244/1‑71 (MQ=40) aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa > 1:2446982‑M1/25‑69 (MQ=255) aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa > 1:3968611‑M1/25‑69 (MQ=255) aagcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa > 1:2920001‑M1/25‑69 (MQ=255) agcgtcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCtaaa < 1:2228188‑M1/46‑2 (MQ=255) agcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACgctaaatc > 1:4238166‑M1/24‑68 (MQ=255) agcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaat < 1:2089850‑M1/47‑3 (MQ=255) agcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaat < 1:2082957‑M1/47‑3 (MQ=255) gcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATCa > 1:655331‑M1/23‑71 (MQ=255) gcgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATCa > 1:4811518‑M1/23‑71 (MQ=255) aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTc < 1:1177903/70‑1 (MQ=39) aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCa < 1:3609408/71‑1 (MQ=38) aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCa < 1:2427778/71‑1 (MQ=38) cgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag > 1:187267‑M1/22‑70 (MQ=255) cgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag > 1:3585021‑M1/22‑70 (MQ=255) cgtcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag > 1:5546817‑M1/22‑70 (MQ=255) cgtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag < 1:519117‑M1/50‑2 (MQ=255) gtcgcatcaggcattgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAAtcagg < 1:1299606‑M1/51‑3 (MQ=255) tcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat > 1:3384337‑M1/20‑64 (MQ=255) tcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGctaaat > 1:3391076‑M1/20‑64 (MQ=255) tcgcatcaggcgttgatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAatcaggc > 1:3504489‑M1/20‑68 (MQ=255) tcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACGTACAAAATCACGCTAAatcaggc < 1:3454996‑M1/52‑4 (MQ=255) tcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCtaaa > 1:462860‑M1/20‑64 (MQ=255) tcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAatcaggc > 1:5118810‑M1/20‑68 (MQ=255) tcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAatcaggc > 1:4696563‑M1/20‑68 (MQ=255) tcgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAatcagga > 1:1080715‑M1/20‑68 (MQ=255) cgcatcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGctaaat > 1:754138‑M1/19‑63 (MQ=255) gTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAatata < 1:3988206/71‑1 (MQ=39) tGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAatatat < 1:2293409/71‑1 (MQ=39) catcaggcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATCa < 1:3876066‑M1/49‑1 (MQ=255) gcgcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATAtt > 1:1995596/1‑71 (MQ=39) cTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGc < 1:4304830/70‑1 (MQ=40) gcgttgatgtCGGATGCGGCATAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag > 1:3826488‑M1/11‑59 (MQ=255) gTTTATGCCGGATTCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTc < 1:3353321/56‑1 (MQ=255) gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaat > 1:2029427/1‑70 (MQ=11) gTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTc < 1:2608126/56‑1 (MQ=35) tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAc > 1:2781268/1‑69 (MQ=255) tATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCt > 1:2820729/1‑71 (MQ=255) tATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaaatca > 1:305977/1‑68 (MQ=11) tATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaaatc < 1:4377490/70‑3 (MQ=11) gatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag < 1:4317859‑M1/50‑2 (MQ=255) gatgtCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAATCACGCTAAATcag < 1:228643‑M1/50‑2 (MQ=255) gATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGc > 1:2859712/1‑71 (MQ=11) gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCt < 1:321210/41‑1 (MQ=32) gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAgg > 1:161687/1‑70 (MQ=255) gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGc < 1:487477/71‑1 (MQ=255) gATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGc > 1:3033291/1‑71 (MQ=255) aTGCGGCGTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaa < 1:4827526/59‑1 (MQ=11) aTGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGc > 1:1538705/1‑55 (MQ=255) tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCt < 1:4258487/71‑1 (MQ=255) gcggcgTAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTa > 1:4840718/1‑65 (MQ=11) cggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAc > 1:5689638/1‑58 (MQ=255) cttAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATaa < 1:1974287/69‑1 (MQ=255) gTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAg > 1:1172367/1‑71 (MQ=255) gTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAg < 1:291041/71‑1 (MQ=255) tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaaat < 1:5390949/53‑2 (MQ=11) tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGaaat < 1:3676233/53‑2 (MQ=11) tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAg < 1:1589502/60‑1 (MQ=11) tAAACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCt < 1:3975788/64‑1 (MQ=11) tAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCa < 1:750435/38‑1 (MQ=40) tAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAACCCTTGt < 1:3206/58‑1 (MQ=255) tAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAg < 1:579092/70‑1 (MQ=255) aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCaa > 1:3982589/1‑37 (MQ=40) aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTg > 1:2682872/1‑44 (MQ=38) aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCg > 1:3221703/1‑70 (MQ=255) aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCg > 1:3680747/1‑70 (MQ=255) aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg > 1:2013041/1‑47 (MQ=11) aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg > 1:2702139/1‑47 (MQ=11) aCGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTa < 1:2269605/71‑1 (MQ=11) gCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTAgc > 1:2944531/1‑71 (MQ=11) gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaatt < 1:4338094/48‑3 (MQ=11) cTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAgcg > 1:234498/1‑70 (MQ=255) attATCCGGCCTACAAAATCGTGCTAATTCAATATATTGGAGAAACCTTg < 1:4279856/49‑1 (MQ=255) tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTa > 1:3572312/1‑50 (MQ=255) tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt > 1:1725128/1‑71 (MQ=255) tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt > 1:1372159/1‑71 (MQ=255) tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt > 1:1676644/1‑71 (MQ=255) ttCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGaa < 1:329634/40‑1 (MQ=11) aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCtt < 1:136632/46‑1 (MQ=255) aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt > 1:1384668/1‑70 (MQ=255) aTCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCAt > 1:596644/1‑70 (MQ=255) tACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTg < 1:559899/47‑1 (MQ=255) tACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAgcgc < 1:5312983/60‑1 (MQ=255) cAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATaa > 1:5266805/1‑49 (MQ=255) cAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCa > 1:2144260/1‑59 (MQ=255) aaaTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtt > 1:4591587/1‑67 (MQ=255) aaTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCg > 1:5019737/1‑49 (MQ=255) aaTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCg > 1:5153723/1‑71 (MQ=255) aTCGTGCTAATTCAATATATTGCAGAAACCTTGTa > 1:2087793/1‑35 (MQ=255) tGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGt < 1:5566025/71‑1 (MQ=255) tGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGt < 1:61967/71‑1 (MQ=255) | GTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGT > minE/1591569‑1591704 |
| 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: ‑ |