| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,591,415 | G→T | 42.3% | K203N (AAG→AAT) | yfeS → | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,591,415 | 0 | G | T | 42.3% | 74.4 / 51.0 | 59 | K203N (AAG→AAT) | yfeS | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (17/17); new base T (13/12); total (30/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01 | |||||||||||
ATGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGACGGTCAG > minE/1591347‑1591481 | atGATTACCTACTCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGAt > 1:470692/1‑71 (MQ=255)atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGAt > 1:1119067/1‑71 (MQ=255) tGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGATg > 1:309567/1‑71 (MQ=255) gATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGATgg > 1:1136816/1‑71 (MQ=255) aTTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATgg < 1:2052542/70‑1 (MQ=255) tACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGATGGGAc > 1:482481/1‑71 (MQ=255) cTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCa < 1:74592/71‑1 (MQ=255) cGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGATGGGACTCaa < 1:213200/69‑1 (MQ=255) cGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGATGGGACTCaa < 1:1471517/69‑1 (MQ=255) cTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGATgg < 1:1142775/57‑1 (MQ=255) tGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAg < 1:1498125/52‑1 (MQ=255) cGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGa > 1:1442889/1‑47 (MQ=255) ggTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGtc < 1:611384/70‑1 (MQ=255) ggTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGtc < 1:1966913/70‑1 (MQ=255) cAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCaa > 1:477101/1‑53 (MQ=255) aTTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTc < 1:549053/48‑1 (MQ=255) aTTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTc < 1:910301/48‑1 (MQ=255) aTTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCt < 1:1714281/71‑1 (MQ=255) tAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCaa > 1:2183971/1‑48 (MQ=255) tAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGtcatc > 1:1504859/1‑65 (MQ=255) tAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTa > 1:251198/1‑71 (MQ=255) tAAAGTCATGGCTAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTa > 1:360606/1‑71 (MQ=255) aaaGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTa < 1:1676499/70‑1 (MQ=255) tCATGGGTAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCa < 1:1774891/47‑1 (MQ=255) tCATGGGTAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGc > 1:694011/1‑70 (MQ=255) aTGGGTAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaaa < 1:975107/71‑1 (MQ=255) aTGGGTAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaaa < 1:1497865/71‑1 (MQ=255) tGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTc < 1:1121750/38‑1 (MQ=255) tGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTc < 1:878681/38‑1 (MQ=255) ggTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAt > 1:452551/1‑71 (MQ=255) ggTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGAGCGTCATCAGCTTACCGCAAAAAt > 1:123180/1‑71 (MQ=255) tAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTa < 1:867639/59‑1 (MQ=255) tAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttc < 1:1605674/71‑1 (MQ=255) aaaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaaaaa > 1:1658790/1‑66 (MQ=255) aaaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaaaaa > 1:791672/1‑66 (MQ=255) aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc < 1:328123/71‑1 (MQ=255) aTCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttctt > 1:131068/1‑69 (MQ=255) tctcCCATAAACTTAAAAAGATGGGACTCAATGCACTAGcgcg < 1:1545954/43‑1 (MQ=255) tctcCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCaa > 1:2116009/1‑71 (MQ=255) tctcCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCaa > 1:1613736/1‑71 (MQ=255) tcCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGtc > 1:519039/1‑43 (MQ=255) ccATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCaa < 1:1915764/67‑1 (MQ=255) tAAACTTAAAAATATGGGACTCAATGCACTAGCGCGtcatc < 1:1836385/41‑1 (MQ=255) tAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCaa > 1:2157995/1‑64 (MQ=255) aaaCTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGa > 1:65811/1‑71 (MQ=255) aaCTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAg > 1:390082/1‑41 (MQ=255) aaCTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAg > 1:877783/1‑41 (MQ=255) aaCTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaaaaa > 1:1078432/1‑54 (MQ=255) aaCTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaaaaa > 1:1500553/1‑54 (MQ=255) aaCTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACAGCaaaaa > 1:712653/1‑54 (MQ=255) aaCTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGAc > 1:1297928/1‑71 (MQ=255) aCTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGACg > 1:1519220/1‑71 (MQ=255) aaaaTATGGGACACAATGCACTAGCGCGTCATCAGCTTAc > 1:1271470/1‑40 (MQ=255) aaaaGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGACGGTCAg < 1:60154/71‑1 (MQ=255) aaaaGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGACGGTCAg < 1:2143044/71‑1 (MQ=255) aaaaGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGACGGTCAg < 1:1677700/71‑1 (MQ=255) aaaaGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCCCAAAAATTCTTCAATGGGGTGACGGTCAg < 1:1170750/71‑1 (MQ=255) aaaGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATg < 1:47033/57‑1 (MQ=255) aaGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGACGGTCAg < 1:29832/69‑1 (MQ=255) aGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaaa < 1:1954705/43‑1 (MQ=255) | ATGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGACGGTCAG > minE/1591347‑1591481 |
| 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: ‑ |