Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,667,126 | 0 | A | C | 48.1% | 22.9 / 12.4 | 27 | Y241D (TAT→GAT) | hcaT | putative 3‑phenylpropionate transporter |
Reads supporting (aligned to +/- strand): ref base A (11/3); new base C (0/13); total (11/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.39e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAA > NC_000913/2666998‑2667255 | aaTGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCAccgccg > 1:631475/1‑140 (MQ=255) cgcAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGctgcc < 1:42563/140‑1 (MQ=255) aTCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCa > 2:201809/1‑140 (MQ=255) cGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGt > 2:659926/1‑140 (MQ=255) aaCAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACCCAGCGCAAAGATAATGCCTCCCGCCACCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCt < 2:184718/140‑1 (MQ=255) aaCAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACCCAGCGAAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCt < 2:974084/140‑1 (MQ=255) cAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTaa < 1:201809/140‑1 (MQ=255) tatCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCg > 1:179905/1‑140 (MQ=255) cgcgTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTaa > 2:639382/1‑140 (MQ=255) gCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGc < 2:179905/140‑1 (MQ=255) aCTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCg < 1:920577/140‑1 (MQ=255) aCTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCg < 1:639382/140‑1 (MQ=255) aCTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCg > 2:508258/1‑140 (MQ=255) ggcggAAAAGTTTATTCCTCAGCGCAAAGATAATGACTTCCGCCCCCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGcc < 1:351304/140‑1 (MQ=255) ggcggAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGcc < 1:644140/140‑1 (MQ=255) ggcggAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGcc < 1:61848/140‑1 (MQ=255) cggAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGcccc < 2:808213/140‑1 (MQ=255) cggAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGcccc > 1:742409/1‑140 (MQ=255) aaGTTTATTACCCAGCGCAAAGATAAGGACTCCCGCCACCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCa < 1:125826/140‑1 (MQ=255) aGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCa > 1:441422/1‑140 (MQ=255) tgccTTCCGCCCCCCGCCCCAGCGACCAAAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAgcgc < 2:441422/137‑1 (MQ=255) gCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTg > 2:584332/1‑140 (MQ=255) cGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCa < 1:584332/140‑1 (MQ=255) cGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCACCAGt < 2:910222/122‑1 (MQ=255) cGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCACCAGt > 1:910222/1‑122 (MQ=255) acccACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCaa < 2:742409/138‑1 (MQ=255) aCCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAAc < 2:863987/130‑1 (MQ=255) aCCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAAc > 1:863987/1‑130 (MQ=255) | AATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAA > NC_000913/2666998‑2667255 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |