| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 599,576 | 0 | T | G | 24.2% | 179.6 / 12.1 | 63 | V288G (GTG→GGG) | cusA | copper/silver export system RND permease |
| Reads supporting (aligned to +/- strand): ref base T (6/41); new base G (15/0); total (22/41) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.83e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTG > NC_000913/599528‑599625 | gcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtg < 1:6385694/50‑1 (MQ=255) ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtgg < 1:2602133/50‑1 (MQ=255) cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtggt < 1:3131086/50‑1 (MQ=255) cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGgggggt > 1:2324496/1‑50 (MQ=255) cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGgggggt > 1:1625542/1‑50 (MQ=255) ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggg > 1:944939/1‑48 (MQ=255) ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtgggg > 1:7228116/1‑48 (MQ=255) ggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAt > 1:4586747/1‑50 (MQ=255) gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATc > 1:2202150/1‑50 (MQ=255) gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATc > 1:7286775/1‑50 (MQ=255) gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATc > 1:3074467/1‑50 (MQ=255) aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt < 1:3964040/50‑1 (MQ=255) aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCt < 1:307146/50‑1 (MQ=255) aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCt > 1:2481655/1‑50 (MQ=255) ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGTTcccg > 1:2262649/1‑48 (MQ=255) ttGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTg < 1:5933955/50‑1 (MQ=255) ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCg < 1:3423727/49‑1 (MQ=255) ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCg < 1:7245626/49‑1 (MQ=255) ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCg < 1:5108997/49‑1 (MQ=255) ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGa < 1:2266216/50‑1 (MQ=255) ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGa < 1:3985127/50‑1 (MQ=255) ccGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGa > 1:6078592/1‑50 (MQ=255) ccGAAATAAACGGCGAAGGCGAAGTGGCGGGAGGGGTGGTGATCCTGCGa < 1:3443113/50‑1 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATc < 1:5475053/50‑1 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATTCTGCGATc > 1:1242713/1‑50 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGTTc > 1:844129/1‑50 (MQ=255) gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGCGGGGATCCTGCGATc > 1:4676565/1‑50 (MQ=255) aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATcc < 1:4807065/50‑1 (MQ=255) aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGCTTCTGCGGTct > 1:6870974/1‑49 (MQ=255) aCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCg < 1:3117577/50‑1 (MQ=255) cTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCgg < 1:3616846/50‑1 (MQ=255) cTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCgg < 1:1455224/50‑1 (MQ=255) cTAAACGGCGAAGGCGAAGTGGCGGGCCGGGTGGTGCTCCCGCGATCCgg > 1:461988/1‑50 (MQ=255) tAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGGGATCCTGCGATCCgg > 1:3690668/1‑49 (MQ=255) tAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTTCGATCCgg > 1:7418277/1‑49 (MQ=37) aaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGc < 1:4752030/48‑1 (MQ=255) aCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCaaa < 1:1718093/50‑1 (MQ=255) cGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCaaaa < 1:406314/50‑1 (MQ=255) ggcgaaggcgaaGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCaaaaa < 1:3265997/50‑1 (MQ=255) ggcgaaggcgaaGTGGCGGGCGGGGGGGGGGTCCTGCGGTCCGGCaaaaa > 1:5337024/1‑50 (MQ=255) ggcgaaggcgaaGTGGCGGGCGGGGGGGGGATCCCGCGAACCGGCaaaaa > 1:740053/1‑50 (MQ=255) gcgaaggcgaagTGGCGGGCGGGGGGGGGATCCTGCGACTCGGCTAAAAc > 1:4763547/1‑50 (MQ=255) cgaaggcgaagTGGCGGGCGGGGGGGGGATCCTGCGATCCGGCAAAAACg > 1:1631873/1‑50 (MQ=255) gaaggcgaagTGGCGGGCGGGGTGGTGGTCCTGCGATCCGGCAAAAACGc < 1:960448/50‑1 (MQ=255) aggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGccc < 1:1176/50‑1 (MQ=255) ggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCg < 1:224882/50‑1 (MQ=255) cgaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCgag < 1:2836586/50‑1 (MQ=255) gaagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCgaga < 1:4141232/50‑1 (MQ=255) aagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGaa < 1:1021299/50‑1 (MQ=255) aagTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGaa < 1:5840995/50‑1 (MQ=255) agTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAg < 1:6242194/50‑1 (MQ=255) gTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGt < 1:3102023/50‑1 (MQ=255) gTGGCGAGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGt < 1:3973636/50‑1 (MQ=255) tGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTg < 1:5740797/50‑1 (MQ=255) gcgggcggGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGa < 1:1858750/49‑1 (MQ=255) gcgggcggGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGa < 1:239370/49‑1 (MQ=255) gcgggcggGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGAt < 1:2520530/50‑1 (MQ=255) cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc < 1:4319764/50‑1 (MQ=255) cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc < 1:1176972/50‑1 (MQ=255) cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc < 1:7368470/50‑1 (MQ=255) cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc < 1:1963409/50‑1 (MQ=255) cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc < 1:623551/50‑1 (MQ=255) cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc < 1:6212058/50‑1 (MQ=255) cgggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc > 1:4852186/1‑50 (MQ=255) gggcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcg > 1:7475279/1‑50 (MQ=255) ggcgggGTGGTGATGCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgc < 1:3686501/50‑1 (MQ=255) gcgggGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgcc < 1:3758008/50‑1 (MQ=255) ggggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccgc > 1:1671453/1‑50 (MQ=255) gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc < 1:6830620/50‑1 (MQ=255) gggTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc > 1:1616664/1‑50 (MQ=255) gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGt < 1:3899878/50‑1 (MQ=255) gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGt > 1:5874813/1‑50 (MQ=255) gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGt > 1:3677941/1‑50 (MQ=255) gtggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGt < 1:1320651/50‑1 (MQ=255) tggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGt > 1:3581444/1‑49 (MQ=255) tggtgATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTg < 1:3809019/50‑1 (MQ=255) | GCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTG > NC_000913/599528‑599625 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |