| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,852,319 | 0 | G | C | 40.6% | 56.4 / 82.9 | 69 | Q285E (CAG→GAG) | rsgA | ribosome small subunit‑dependent GTPase A |
| Reads supporting (aligned to +/- strand): ref base G (18/23); new base C (13/15); total (31/38) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCGTGTTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCAC > minE/2852258‑2852382 | tCGTGTTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTg > 1:4886460/1‑71 (MQ=255)tCGTGTTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTg > 1:1473326/1‑71 (MQ=255)tCGTGTTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTg > 1:1971755/1‑71 (MQ=255)tCGTGTTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTg > 1:1813510/1‑71 (MQ=255) gtgtTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGtt > 1:1250706/1‑71 (MQ=255) tttGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCg > 1:1174973/1‑71 (MQ=255) tttGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCg > 1:4412269/1‑71 (MQ=255) tGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc > 1:845769/1‑71 (MQ=255) tGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc > 1:4317865/1‑71 (MQ=255) tGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc > 1:2349188/1‑71 (MQ=255) tGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc > 1:2195896/1‑71 (MQ=255) gCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc < 1:3059882/70‑1 (MQ=255) gCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc < 1:3732542/70‑1 (MQ=255) gCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc < 1:2325207/70‑1 (MQ=255) gCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGc < 1:1600312/70‑1 (MQ=255) aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg < 1:4293240/61‑1 (MQ=37) aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg < 1:3048949/61‑1 (MQ=37) aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg < 1:296000/61‑1 (MQ=37) aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg < 1:2948145/61‑1 (MQ=37) aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg < 1:2855677/61‑1 (MQ=37) aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg < 1:1811397/61‑1 (MQ=37) aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg < 1:3420681/61‑1 (MQ=37) aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg < 1:3645799/61‑1 (MQ=37) aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg < 1:4082575/61‑1 (MQ=37) aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg < 1:2140367/61‑1 (MQ=37) aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg < 1:4601851/61‑1 (MQ=37) aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg < 1:1904762/61‑1 (MQ=37) aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg < 1:185387/61‑1 (MQ=37) aaTCGCGATATTTACACAGACCTAAATAGTCATGGAATACGACAAATCACTCAGGGCTTTg < 1:3391353/61‑1 (MQ=25) aTCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTg < 1:4467940/60‑1 (MQ=25) tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc > 1:3184812/1‑60 (MQ=25) tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc > 1:4371975/1‑60 (MQ=25) tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc > 1:3579138/1‑60 (MQ=25) tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc > 1:4498201/1‑60 (MQ=25) tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc > 1:4780622/1‑60 (MQ=25) tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc > 1:1173593/1‑60 (MQ=25) tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc > 1:1790294/1‑60 (MQ=25) tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc > 1:1989108/1‑60 (MQ=25) tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc > 1:2762747/1‑60 (MQ=25) tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc > 1:2761603/1‑60 (MQ=25) tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc > 1:2742086/1‑60 (MQ=25) tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc > 1:2894957/1‑60 (MQ=25) tCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAATCACTCAGGGCTTTGtc > 1:2644703/1‑60 (MQ=25) atatTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGcc < 1:1665753/71‑1 (MQ=255) tACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCagag < 1:1606548/70‑1 (MQ=255) gACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAAc < 1:4195010/71‑1 (MQ=255) gACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAAc < 1:4121507/71‑1 (MQ=255) gACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAAc < 1:4012701/71‑1 (MQ=255) aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGAt > 1:3109736/1‑35 (MQ=255) aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACg > 1:2535318/1‑71 (MQ=255) aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACg > 1:411856/1‑71 (MQ=255) aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACg > 1:414741/1‑71 (MQ=255) aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACg > 1:3149772/1‑71 (MQ=255) aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACg > 1:1324699/1‑71 (MQ=255) aaaTAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCAGGCTCCAGGTGCCAGAGGCCGAACTCAc > 1:4782047/1‑70 (MQ=255) aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg < 1:3373312/71‑1 (MQ=255) aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg < 1:510481/71‑1 (MQ=255) aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg < 1:4613349/71‑1 (MQ=255) aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg < 1:4538645/71‑1 (MQ=255) aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg < 1:2265221/71‑1 (MQ=255) aTGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg < 1:2578753/71‑1 (MQ=255) tGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg < 1:1514603/70‑1 (MQ=255) gAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTgggg < 1:72508/71‑1 (MQ=255) ttCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGcc < 1:12341/41‑1 (MQ=255) aCAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAat < 1:2741247/71‑1 (MQ=255) aaGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAatcac < 1:2174942/71‑1 (MQ=255) aaGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAatcac < 1:4408771/71‑1 (MQ=255) aaGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAatcac < 1:2471259/71‑1 (MQ=255) cccTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTg < 1:3075735/55‑1 (MQ=255) | TCGTGTTTGCAATCGCGATATTTACACAGACCTAAATAGTCATGGAATTCGACAAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCAC > minE/2852258‑2852382 |
| 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: ‑ |