| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_007779 | 884,768 | 0 | C | G | 46.9% | 31.5 / 22.5 | 32 | A225G (GCG→GGG) | cmr | multidrug efflux system protein |
| Reads supporting (aligned to +/- strand): ref base C (2/15); new base G (15/0); total (17/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.42e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGC > NC_007779/884670‑884868 | gCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGggcgc < 3:1260299/101‑1 (MQ=255) cTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGggcgctg < 3:628177/101‑1 (MQ=255) tGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGggcgctgg < 1:881258/101‑1 (MQ=255) gCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTggg < 3:1557389/101‑1 (MQ=255) cACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCGGGCGCTggggt > 3:1195975/1‑99 (MQ=255) aCGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGAtt < 1:352296/101‑1 (MQ=255) cgTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGt < 1:348299/101‑1 (MQ=255) aGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTc < 1:608907/101‑1 (MQ=255) gCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCGGGCGCGGGGGGTCGTTTGTCTg > 3:924615/1‑101 (MQ=255) gagaAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCTGGCGCCGGGGTGCGTTTGTTTGcc > 3:1485350/1‑101 (MQ=255) aaaCTGTCACCGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGtt < 1:229714/101‑1 (MQ=255) tGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCGGGGGGCGGGGTTCGTTAGTTTGCCGTTGCTg > 1:726968/1‑101 (MQ=255) tGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTg < 1:1067501/101‑1 (MQ=255) gTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGGCGGGGCGGGGGGTCGTTTGTCTGCCGCTGCTgg > 3:594999/1‑101 (MQ=255) cACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCGGGCGCTGGGGTTGGTTCATGTGCCGTTGCTGGCg > 3:924127/1‑101 (MQ=255) cACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCCGGCGCCGGGGGTTGTTTGTCTGcccgttctgggg > 1:1179422/1‑91 (MQ=255) tGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTgg < 1:953644/101‑1 (MQ=255) aaaGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCTGGCGGTGGGGATAGTTAGTCTGCCGTTGCttgggggggt > 3:670830/1‑92 (MQ=255) aaGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGACCGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATc < 3:639559/101‑1 (MQ=255) aaGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATc < 3:298910/101‑1 (MQ=255) cGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGt < 1:441087/101‑1 (MQ=255) gACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGGTGGCGCTGGGGTTCGTTAGTCTGCCGCTGCTGGCGGGGGTCGccccggccgcccg > 1:325129/1‑91 (MQ=255) tataAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCTGGCGCGGGGGTTCGTTAGTTTGTCGTTTTTGGGGTGGGTCGCCCCCTCGCCggatat > 1:951858/1‑96 (MQ=255) tataAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTat < 1:1012339/101‑1 (MQ=255) tataAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTat < 3:281254/101‑1 (MQ=255) gCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCTGGCGCCGGGGTTCGTTAGTCTGCCGTTTCTGGCGTGGGTCGCCCCGTCCCCGACTCtcatca > 3:1048885/1‑101 (MQ=255) tGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCGGGCGCTGGGGTTCGTTTGTCTGCCGTTGCTGGCGGGGGTCGCCCCGTCGCCGGTTATCATCAtt > 1:52561/1‑101 (MQ=255) ggTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGGGGGCGCCGGGGTTCGGTTGGCTTCCGCTGCTTGCGGGGGTTGCCCACACGCCGCGTATTATCATTa > 3:1286327/1‑101 (MQ=255) tGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCTGGCGCTGGGGTTCGTTTGTCTGCCGTTGCTGGCGGGGGTCGCCCCGTCGCCGCGTATTATCATTAcc > 1:812170/1‑101 (MQ=255) tGAAGAACGGCCGCTTTGTGGCGGGGGGGGTGGGGCGGGGATTCGTTCTTCTGCCGTTCCGGGTCTGGGTCGCCCAGCCGCCCCTTatgatatatacaggc > 3:344130/1‑91 (MQ=255) gcaCGGCCGCTTTGTGGCGGGTGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGc < 1:517542/99‑1 (MQ=255) gCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTg < 1:915364/101‑1 (MQ=255) cGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGag < 1:843109/101‑1 (MQ=255) gggggCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAat > 3:86340/1‑101 (MQ=255) gggCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAatat > 1:679050/1‑101 (MQ=255) cgctggcgctggGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGc > 1:1211003/1‑101 (MQ=255) | GCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGC > NC_007779/884670‑884868 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |