| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,789,346 | Δ1 bp | 34.3% | coding (483/1269 nt) | ygaF → | hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,789,346 | 0 | A | . | 34.3% | 40.7 / 43.4 | 35 | coding (483/1269 nt) | ygaF | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base A (10/10); new base . (9/3); total (21/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.67e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01 | |||||||||||
TGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACGGCGTGGTGAT > minE/1789277‑1789411 | tgtgCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCaga < 1:371015/70‑1 (MQ=255) gtgCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCaga < 1:1034482/69‑1 (MQ=255) tCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGa < 1:1129242/71‑1 (MQ=255) tCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGa < 1:1801158/71‑1 (MQ=255) tCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGa < 1:2071490/71‑1 (MQ=255) cGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTa > 1:80641/1‑71 (MQ=255) tCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGGGGCGGTGa < 1:385981/58‑1 (MQ=255) aGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCaga > 1:552525/1‑48 (MQ=255) tACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAG‑GGCGGTGAAATTATCTATAACGCCg > 1:271398/1‑69 (MQ=255) ccGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGaa < 1:492360/70‑1 (MQ=255) cgcgACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGt < 1:371674/71‑1 (MQ=255) gTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATc < 1:1716663/51‑1 (MQ=255) aCGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCt > 1:561283/1‑71 (MQ=255) aCGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTAAGCGGGCt > 1:1654342/1‑71 (MQ=255) aCGGCGGCGATGGCAAAAATCTTCCAGTCCAG‑GGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGc < 1:2084327/69‑1 (MQ=255) aCGGCGGCGATGGCAAAAATCTTCCAGTCCAG‑GGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCt > 1:1433298/1‑70 (MQ=255) ggcgATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAAt < 1:1684480/70‑1 (MQ=255) cgATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAg < 1:1700885/71‑1 (MQ=255) gATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAg > 1:1270814/1‑70 (MQ=255) tGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACg > 1:1049873/1‑45 (MQ=255) ggCAAAAATCTTCCAGTCCAG‑GGCGGTGAAATTATCTATAACGc > 1:222964/1‑44 (MQ=38) ggCAAAAATCTTCCAGTCCAG‑GGCGGTGAAATTATCTATAACGc > 1:1037241/1‑44 (MQ=38) cAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGt > 1:1779342/1‑49 (MQ=255) cAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGt > 1:2073218/1‑49 (MQ=255) aaaaTCTTCCAGTCCAG‑GGCGGTGAAATTATCTATAACGCCGAAGt < 1:1541904/46‑1 (MQ=39) aaaaTCTTCCAGTCCAG‑GGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAAc > 1:1344907/1‑71 (MQ=255) aaTCTTCCAGTCCAGGGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGa > 1:663271/1‑60 (MQ=255) aaTCTTCCAGTCCAGGGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGa > 1:1982771/1‑60 (MQ=255) aaTCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACg > 1:1578923/1‑71 (MQ=255) aaTCTTCCAGTCCAG‑GGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAg > 1:909029/1‑60 (MQ=255) aaTCTTCCAGTCCAG‑GGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACgg > 1:2143823/1‑71 (MQ=255) aaTCTTCCAGTCCAG‑GGCGGTGAAATTATCTATAACGCCGAAGTCAGAGGGCTTAATGAg > 1:1707771/1‑60 (MQ=255) aTCTTCCAGTCCAG‑GGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACGGc > 1:1024169/1‑71 (MQ=255) ccagtccagAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACGGCgtg < 1:1169036/70‑1 (MQ=255) ccagtccagAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACGGCgt < 1:724581/69‑1 (MQ=255) gtccag‑GGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACGGCGTGGTGAt < 1:1473257/71‑1 (MQ=255) cagAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCg > 1:233937/1‑38 (MQ=255) cagAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGa > 1:987787/1‑48 (MQ=255) | TGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACGGCGTGGTGAT > minE/1789277‑1789411 |
| 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: ‑ |