| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,248,382 | 0 | C | G | 27.6% | 42.0 / 10.3 | 33 | A51P (GCA→CCA) | lysP | lysine transporter |
| Reads supporting (aligned to +/- strand): ref base C (15/6); new base G (8/0); total (27/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.48e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACC > NC_000913/2248320‑2248448 | ttCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCcccgc < 2:411449/68‑1 (MQ=255)ttCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCcccgc > 2:2635/1‑68 (MQ=255) tCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGGcccgcc > 2:429336/1‑68 (MQ=255) ccGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCgg > 2:646291/1‑68 (MQ=255) cGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGa > 1:565591/1‑68 (MQ=255) gagaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGAc > 2:458896/1‑68 (MQ=255) gagaCTGGTCATCAGGAAGTAAACAATCAGGCCAAACAGCATAAACGAGAGCAATGACACGACCGGAc > 1:577465/1‑68 (MQ=255) tcatcaGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGGCCCGCCCGGACCTGCCTga > 2:529791/1‑68 (MQ=255) ggAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGACCCGCCCGGACCTGCCTGAGAAATc > 2:778245/1‑68 (MQ=255) gTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGGCCCGCCCGGACCTGCCTGAGAAATCGTTg > 2:315708/1‑68 (MQ=255) aaCCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTgcgc > 1:56415/1‑68 (MQ=255) aGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGc < 2:294452/68‑1 (MQ=255) gCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGcaa > 2:611846/1‑68 (MQ=255) gCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGcaa > 1:857803/1‑68 (MQ=255) aaTCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGcaacaa < 1:242360/68‑1 (MQ=255) cAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACaaaaa > 1:182897/1‑68 (MQ=255) acgCAAATACGAGAGCAATGACACGCCAGGACCTGAATGAGAACTCGTTCAGCAAGAGGCACaaaaaa > 2:573199/3‑68 (MQ=255) gCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGa < 2:952015/68‑1 (MQ=255) atatACGAGAGCAATGGCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGAcc > 2:857156/1‑68 (MQ=255) atatACGAGAGCAATGGCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGAcc > 2:689737/1‑68 (MQ=255) ataCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGcaaca < 1:550005/58‑1 (MQ=255) ataCGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGcaaca > 2:550008/1‑58 (MQ=255) cGAGAGCACTGCCCCGCCCGGACCTGCCTGAGAACTCGTTGCGCCAGAGGCAACAAAAAGCCCTGTAc > 1:715469/1‑68 (MQ=255) cGAGAGCAATGGCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTAc > 2:791696/1‑68 (MQ=255) cGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGtgc > 1:321150/1‑66 (MQ=255) cGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTAc > 2:738094/1‑68 (MQ=255) agCAATGGCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGAt > 1:34736/1‑68 (MQ=255) agCAATGGCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGAt > 2:363621/1‑68 (MQ=255) gCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATg > 1:436552/1‑68 (MQ=255) gCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATg > 1:102869/1‑68 (MQ=255) gCAATGACCCGCCCGGACCTGCCTGAGAAATCGTAGCGACTTAGTCAACAAAAAGCCCTGTACCGATg > 1:482686/1‑68 (MQ=255) aaTGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGa < 2:946044/68‑1 (MQ=255) gCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAAcc > 1:680681/1‑68 (MQ=255) | TTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCTGAGAAATCGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACC > NC_000913/2248320‑2248448 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |