Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 49,644 | T→C | 40.2% | K155K (AAA→AAG) | ksgA ← | S‑adenosylmethionine‑6‑N',N'‑adenosyl (rRNA) dimethyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 49,644 | 0 | T | C | 40.2% | 66.4 / 112.0 | 154 | K155K (AAA→AAG) | ksgA | S‑adenosylmethionine‑6‑N',N'‑adenosyl (rRNA) dimethyltransferase |
Reads supporting (aligned to +/- strand): ref base T (69/23); new base C (18/44); total (87/67) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.38e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01 |
Only 100 of 154 total aligned reads displayed.
GGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGC > minE/49594‑49709 | gggATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATTCGCTTTGCTGTTCGGt < 1:413213/63‑1 (MQ=255) gggATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGt < 1:2124983/63‑1 (MQ=255) gggATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGt < 1:2198403/63‑1 (MQ=255) gggATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGt < 1:231356/63‑1 (MQ=255) gggATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGt < 1:3211533/63‑1 (MQ=255) gggATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGt < 1:2890955/63‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:3090692/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:230960/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:2418390/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:2487608/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:2543669/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:2730536/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:3004064/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:1039260/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:3262452/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:3338927/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:3409274/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:544495/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:572534/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:705668/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:757952/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:976342/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:144890/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:1094665/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:1201445/69‑1 (MQ=255) tCACGTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:1553924/69‑1 (MQ=255) cacgTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:108978/68‑1 (MQ=255) cacgTTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCCTTGCTGTTCGGTCCTGCAACCa < 1:359533/68‑1 (MQ=255) gCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGAtt > 1:2991190/1‑69 (MQ=255) gCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGAtt > 1:1345286/1‑69 (MQ=255) aGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcacca > 1:3042113/1‑69 (MQ=255) aGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcacca > 1:1739034/1‑69 (MQ=255) aGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcacca > 1:653210/1‑69 (MQ=255) aGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcacc > 1:1008835/1‑68 (MQ=255) gTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcaccac < 1:2294546/69‑1 (MQ=255) gTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcaccac < 1:3302852/69‑1 (MQ=255) gTATTGCGCCATGACGCTTAATAGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTcaccac > 1:2194128/1‑69 (MQ=255) gcgcCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGa < 1:1927108/56‑1 (MQ=255) gcgcCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGa < 1:274593/56‑1 (MQ=255) cATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCttttt > 1:1557931/1‑68 (MQ=255) cATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCttttt > 1:967326/1‑68 (MQ=255) cATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCttttt > 1:553018/1‑68 (MQ=255) cATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCttttt > 1:440101/1‑68 (MQ=255) aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTg > 1:2143825/1‑68 (MQ=255) aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc > 1:2433937/1‑69 (MQ=255) aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc > 1:39457/1‑69 (MQ=255) aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc > 1:2535788/1‑69 (MQ=255) aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc > 1:1682682/1‑69 (MQ=255) aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc > 1:1883802/1‑69 (MQ=255) aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc > 1:482479/1‑69 (MQ=255) aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc > 1:1140362/1‑69 (MQ=255) aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc > 1:3326096/1‑69 (MQ=255) aTGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGc > 1:1536620/1‑69 (MQ=255) aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGc > 1:1403197/1‑38 (MQ=255) aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACctc > 1:2200310/1‑59 (MQ=255) aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa > 1:3153286/1‑68 (MQ=255) aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa > 1:3363671/1‑68 (MQ=255) aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa > 1:3428436/1‑68 (MQ=255) aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa > 1:3053471/1‑68 (MQ=255) aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa > 1:2968023/1‑68 (MQ=255) aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa > 1:1835420/1‑68 (MQ=255) aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa > 1:2014897/1‑68 (MQ=255) aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa > 1:2068774/1‑68 (MQ=255) aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa > 1:2180163/1‑68 (MQ=255) aCGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCaa > 1:150572/1‑68 (MQ=255) cTTAATCGACCATACGCTTTGCTTTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATa < 1:2967161/69‑1 (MQ=255) cTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATa < 1:3242411/69‑1 (MQ=255) cGACCATACGCTTTGCTGTTCGGTCCTGCACCCAGACGATTCACCACCTCTTTTTGCAACATAAAGTg < 1:1306263/68‑1 (MQ=255) cGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTg < 1:2013358/68‑1 (MQ=255) cGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTg < 1:901566/68‑1 (MQ=255) cGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTg < 1:1486618/68‑1 (MQ=255) cGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTg < 1:609941/68‑1 (MQ=255) aCCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATa < 1:1218013/61‑1 (MQ=255) ccATACGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACttc > 1:3144255/1‑47 (MQ=39) ccATACGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACttc > 1:2506161/1‑47 (MQ=39) ccATACGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACttc > 1:2640117/1‑47 (MQ=39) ccATACGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACttc > 1:3134878/1‑47 (MQ=39) ccATACGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACttc > 1:1090053/1‑47 (MQ=39) ccATACGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACttc > 1:2806389/1‑47 (MQ=39) ccATACGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACttc > 1:2842489/1‑47 (MQ=39) cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtt > 1:3413462/1‑44 (MQ=255) cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtt > 1:745180/1‑44 (MQ=255) cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtt > 1:2438474/1‑44 (MQ=255) cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc > 1:2416337/1‑69 (MQ=255) cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc > 1:130487/1‑69 (MQ=255) cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc < 1:506250/69‑1 (MQ=255) cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc > 1:2578115/1‑69 (MQ=255) cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc > 1:2763902/1‑69 (MQ=255) cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc > 1:1230853/1‑69 (MQ=255) cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc > 1:1234797/1‑69 (MQ=255) cGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGc > 1:3262207/1‑69 (MQ=255) cGCTTTGCTGTTCGGTCCTGAAACCAGACGATTCACCACCTCtt > 1:370581/1‑44 (MQ=255) cGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACTTCTTTTTGCAACATAAAGTGCATGTCGGc < 1:313708/69‑1 (MQ=255) cGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACTTCTTTTTGCAACATAAAGTGCATGTCGGc < 1:2156663/69‑1 (MQ=255) cGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACTTCTTTTTGCAACATAAAGTGCATGTCGGc < 1:591138/69‑1 (MQ=255) cGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACTTCTTTTTGCAACATAAAGTGCATGTCGGc < 1:3307524/69‑1 (MQ=255) cGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACTTCTTTTTGCAACATAAAGTGCATGTCGGc < 1:886154/69‑1 (MQ=255) cGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACTTCTTTTTGCAACATAAAGTGCATGTCGGc < 1:948451/69‑1 (MQ=255) cGCCTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACTTCTTTTTGCAACATAAAGTGCATGTCGGc < 1:1391572/69‑1 (MQ=255) agCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCtt > 1:2925812/2‑44 (MQ=255) | GGGATCACATTGCAATAGTATTGCGCCATGACGCTTAATCGACCATACGCTTTGCTGTTCGGTCCTGCAACCAGACGATTCACCACCTCTTTTTGCAACATAAAGTGCATGTCGGC > minE/49594‑49709 |
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: ‑ |