| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,984,188 | 0 | C | G | 45.2% | 19.6 / 39.4 | 31 | G294G (GGG→GGC) | araG | L‑arabinose ABC transporter ATPase |
| Reads supporting (aligned to +/- strand): ref base C (1/16); new base G (14/0); total (15/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.66e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTT‑‑‑CACCACTGCGAACCGCCAGACTTAT > NC_000913/1984126‑1984247 | atatGAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgccccgc < 1:25940/68‑1 (MQ=255) ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTTATTCGCTACGGCCCGCCCCTGcc > 2:1112204/1‑68 (MQ=255) ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgcc < 2:849631/57‑1 (MQ=255) ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAcgcc > 1:849634/1‑57 (MQ=255) ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTAGGGCCCGCCCCTAcc > 2:985293/1‑68 (MQ=255) ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTAcc > 2:1106379/1‑68 (MQ=255) ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTAcc > 2:983539/1‑68 (MQ=255) ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTAcc > 2:863991/1‑68 (MQ=255) ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTAcc > 1:1093767/1‑68 (MQ=255) ccTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTAcc > 1:297222/1‑68 (MQ=255) ccTGACCGGCGGTGATTTGCGTCCAGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTAcc > 1:498092/1‑68 (MQ=255) ccTGACCGGCGGTGATTAGCGTCACGCCAAACATGCCTTTCATTAATTCGCTACGGCACGACCCTAcc > 1:553846/1‑68 (MQ=255) tGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAg < 2:261006/68‑1 (MQ=255) ccGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTACCAGAcc > 2:353879/1‑68 (MQ=255) ggTGATTTGCGTCCCGCCATACATGCCTTTCATTATTTCGCTACGCTCCGCCCCTACCAGACCAAACa < 2:85856/68‑1 (MQ=255) ggTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACa < 1:514769/68‑1 (MQ=255) tGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGc < 2:868946/68‑1 (MQ=255) ggATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGc < 2:1074040/67‑1 (MQ=255) gATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGcc < 1:182555/68‑1 (MQ=255) tttGCGTCCCGCCAAACATGCCTTTCGTTAATTCGCTACGGCCCGCCCCTACCAGGCCAAACAGCCca > 2:1140592/1‑68 (MQ=255) tttGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGGCCCGCCCCTACCAGACCAAACAGCCca > 1:580063/1‑68 (MQ=255) cGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCcaacaa < 2:117380/68‑1 (MQ=255) cGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCcaacaa < 2:297222/68‑1 (MQ=255) cGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCcaacaa < 2:378243/68‑1 (MQ=255) aaaCATGCCTTTCATTAATTCGCTACGGCCCGCCCCTACCAGACCAAACAGCCCAACAATTT‑‑‑caccac > 1:364968/1‑68 (MQ=255) aaCATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTT‑‑‑CACCACt < 2:801471/68‑1 (MQ=255) acTTTCATTAATTCGCTACGGCCAGCCAATACAAGAACAAGCAGCAAACAAATTT‑‑‑CAAAACAGCGAAc > 1:780642/2‑68 (MQ=255) cTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTT‑‑‑CACCACTGCGAAcc < 2:371150/68‑1 (MQ=255) tttCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTT‑‑‑CACCACTGCGAACCg < 1:436871/68‑1 (MQ=255) ttCATAAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTT‑‑‑CACCACTGCGAACCGc < 1:251006/68‑1 (MQ=255) tAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTT‑‑‑CACCACTGCGAACCGCCAGAc < 1:380935/68‑1 (MQ=255) tAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTT‑‑‑CACCACTGCGAACCGCCAGAc < 2:586901/68‑1 (MQ=255) ttCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTTCACCACCACTGCGAACCGCCAGAc < 2:1038255/68‑1 (MQ=255) tCGCTACGGCCCGCCCCTACCAGACCAAACAGCCCAACAATTT‑‑‑CACCACTGCGAACCGCCAGACTTAt > 1:589108/1‑68 (MQ=255) | ATATAAACCTGACCGGCGGTGATTTGCGTCCCGCCAAACATGCCTTTCATTAATTCGCTACGCCCCGCCCCTACCAGACCAAACAGCCCAACAATTT‑‑‑CACCACTGCGAACCGCCAGACTTAT > NC_000913/1984126‑1984247 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |