| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,211,283 | A→T | 9.0% | A23A (GCT→GCA) ‡ | cspC ← | stress protein, member of the CspA‑family |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,211,283 | 0 | A | T | 9.0% | 210.6 / 10.6 | 77 | A23A (GCT→GCA) ‡ | cspC | stress protein, member of the CspA‑family |
| Reads supporting (aligned to +/- strand): ref base A (30/40); new base T (2/5); total (32/45) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.93e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTT > minE/1211215‑1211353 | tCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc > 1:81415/1‑71 (MQ=255)tCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc > 1:3223482/1‑71 (MQ=255) cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc > 1:2272288/1‑71 (MQ=255) cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc > 1:2591574/1‑71 (MQ=255) cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc > 1:2992955/1‑71 (MQ=255) cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc > 1:3071137/1‑71 (MQ=255) cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGcc > 1:3032835/1‑71 (MQ=255) cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc < 1:864049/70‑1 (MQ=255) cAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGc < 1:1122496/70‑1 (MQ=255) agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt < 1:436698/71‑1 (MQ=255) agagTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt < 1:2435086/71‑1 (MQ=255) gagTTTTAAAACCATTAGTCTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGt > 1:3538257/1‑70 (MQ=255) tGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGaa < 1:4144428/70‑1 (MQ=255) tGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAg > 1:1055842/1‑71 (MQ=255) gAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCa < 1:294587/55‑1 (MQ=255) ttACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGa < 1:3547024/61‑1 (MQ=255) aCCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAg > 1:2499169/1‑47 (MQ=255) gtcTGGATTGCAGAGAAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCtt < 1:176250/68‑1 (MQ=255) ctcTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCttt < 1:3281052/69‑1 (MQ=9) ctcTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCttt < 1:3928632/69‑1 (MQ=9) ctcTGGATGGCGGTGAAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCttt < 1:998392/69‑1 (MQ=9) cccTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt < 1:519036/70‑1 (MQ=255) ccTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCTTCAGCCGGAGTAATGAAGCCAAAACCTTTa < 1:1024990/71‑1 (MQ=255) ccTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa < 1:174806/71‑1 (MQ=255) tGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGt < 1:2012236/50‑1 (MQ=255) tGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGa > 1:1306103/1‑71 (MQ=255) tGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGACTAATGAAg > 1:2273317/1‑57 (MQ=255) gATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTaa > 1:3880776/1‑50 (MQ=255) gATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTaa > 1:1361331/1‑50 (MQ=255) gATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGa > 1:2403135/1‑53 (MQ=255) aTAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTc < 1:987740/71‑1 (MQ=255) aTAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTc < 1:1077156/71‑1 (MQ=255) tAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg > 1:2372705/1‑42 (MQ=255) tAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg > 1:1498796/1‑42 (MQ=255) tAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCg > 1:3045156/1‑42 (MQ=255) gCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt < 1:1532013/61‑1 (MQ=255) cagagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCAtta < 1:3513115/69‑2 (MQ=255) cagagaAGTGTACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCAtta < 1:2920047/69‑2 (MQ=255) cGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTg < 1:2085230/71‑1 (MQ=255) cGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTg < 1:3129450/71‑1 (MQ=255) ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaa < 1:1566687/53‑1 (MQ=255) ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaa < 1:920578/53‑1 (MQ=255) ggAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTg < 1:1812119/70‑1 (MQ=255) gtgaAATGGACGAAAACGTCTTTGCTGCCATCATCAGGAGTGATAAAGCCaa < 1:2806829/50‑1 (MQ=9) gagaAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTa < 1:1362761/42‑1 (MQ=255) gagaAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTa < 1:3648418/42‑1 (MQ=255) gagaAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGtt > 1:565356/1‑68 (MQ=255) tgtACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCAtta < 1:3980499/63‑2 (MQ=9) tgtACGAACACGTCTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACtt < 1:1475939/71‑1 (MQ=12) tACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaaa < 1:2253170/45‑1 (MQ=255) tACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt < 1:2967863/50‑1 (MQ=255) tACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt < 1:3732320/50‑1 (MQ=255) tACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTaa > 1:2218304/1‑71 (MQ=255) tACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTaa > 1:4031362/1‑71 (MQ=255) aCGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa < 1:1258982/51‑1 (MQ=255) cGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCttt > 1:3870191/1‑49 (MQ=255) gAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa < 1:272721/49‑1 (MQ=255) aaCACATCTTTGCTGCCATCAGCCGGAGTAATGAAGACAAAACCTTTAGACTCGTTGAACCACTTAACCTg < 1:1069805/71‑1 (MQ=255) acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaaaa > 1:791573/1‑41 (MQ=255) acacATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTg < 1:3472441/70‑1 (MQ=255) cacaTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGAc > 1:1542654/1‑71 (MQ=255) caTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa > 1:2468677/1‑44 (MQ=255) caTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGtt > 1:1353606/1‑52 (MQ=255) aTCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTa < 1:878038/43‑1 (MQ=255) tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCaaa < 1:3233421/35‑1 (MQ=255) tCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCtt < 1:412725/40‑1 (MQ=255) cTTTGCTGCCGTCTTCCGGAGTAATGAAACCGAATCCTTTGGACTCATTAAACCACTTAACGTTACCttt > 1:1573532/1‑70 (MQ=255) tttGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTaa > 1:3207005/1‑71 (MQ=255) tgctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTATAATc > 1:3009437/1‑71 (MQ=255) ctgcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCtt < 1:787853/71‑1 (MQ=255) gcCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTg > 1:576175/1‑71 (MQ=255) cATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAAc > 1:641329/1‑52 (MQ=255) aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTg < 1:847488/54‑1 (MQ=255) aTCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCa > 1:2139001/1‑71 (MQ=255) tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAAc < 1:755129/42‑1 (MQ=255) tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGa < 1:3598664/54‑1 (MQ=255) tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGAc < 1:3901107/55‑1 (MQ=255) tCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAt < 1:3714274/71‑1 (MQ=255) aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCa > 1:1931556/1‑42 (MQ=255) aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCt > 1:1060131/1‑55 (MQ=255) aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt > 1:2768151/1‑71 (MQ=255) aGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCAttt > 1:513095/1‑71 (MQ=255) | TCAGCCAGAGTTTTGAAGCCATTACCCTGGATAGCGGAGAAGTGTACGAACACATCTTTGCTGCCATCAGCCGGAGTAATGAAGCCAAAACCTTTAGACTCGTTGAACCACTTAACCTGACCTTTAATCTTTGCCATTT > minE/1211215‑1211353 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |