| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,374,036 | 0 | T | G | 28.6% | 29.2 / 13.5 | 21 | L272R (CTC→CGC) | ycjO | putative ABC transporter membrane subunit YcjO |
| Reads supporting (aligned to +/- strand): ref base T (3/12); new base G (0/6); total (3/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.26e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAA > NC_000913/1373973‑1374103 | atCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTtc < 1:260274/68‑1 (MQ=255) tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGGTCTCGGTGGGGGTCTtca > 1:4141562/1‑68 (MQ=255) tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGGTCTCGGGGGGGCTCTtct > 1:2535576/1‑67 (MQ=255) tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGGTCTtcc > 1:3904685/1‑67 (MQ=255) tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTcca > 1:3608202/1‑68 (MQ=255) tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGGGGTGCTCttta > 1:2595337/1‑66 (MQ=255) tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGGGGTGCCCTtca > 1:1438333/1‑68 (MQ=255) tctcTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGGGGGGCTCTtca > 1:3331519/1‑68 (MQ=255) aCAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGGGGTGCTCTTCATCatt > 1:3573142/1‑68 (MQ=255) cctttgcctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCt < 1:387412/68‑1 (MQ=255) ctttgcctttAATGATTTAGGAAAAGCGGCGGCGAACTCGGGGGTGCTCTTCATCATTATTTTCCCTg > 1:3192424/1‑68 (MQ=255) ttgcgtttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc < 1:81622/68‑1 (MQ=255) ttgcctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTc < 1:3744485/68‑1 (MQ=255) tttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTct < 1:2117345/68‑1 (MQ=255) tttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTct < 1:567478/68‑1 (MQ=255) tttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTct < 1:998109/68‑1 (MQ=255) aTTTAAGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGAcc < 1:4174022/68‑1 (MQ=255) ttAGGAAAAGCGGCGGCGATCTCGGTGGT‑CGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAgg < 1:3473481/68‑1 (MQ=255) ttAGAAAAAGGGGCGTCGATCTCGGTGGT‑CGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAgg < 1:2656267/68‑1 (MQ=255) tAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg < 1:297562/68‑1 (MQ=255) tAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAgg < 1:2640781/68‑1 (MQ=255) aGGCAAGGCGGCGGCGATCTCGGTGGT‑CGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGaa < 1:3380252/68‑1 (MQ=255) gAAAAGCGGCGGCGGTCTCGGGGGTGGTTTTTATCATTATTTTTGCTGTTACTCTTCTTACCACGaaa > 1:1681246/1‑68 (MQ=16) aaaGCGGCGGCGATCTCGGTGGT‑CGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACgg < 1:3168278/68‑1 (MQ=255) aaGCGGCGGCGATCTCGGTGGTGCTCTTCTTCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACgg < 1:1858166/68‑1 (MQ=255) aGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACggg < 1:512505/68‑1 (MQ=255) ggcggcgATCTCGGTGGT‑CGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGaa < 1:1999761/68‑1 (MQ=255) cgcgtcgATCTCGGTGGT‑CGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGaa < 1:4181195/67‑1 (MQ=255) gcgccgATCTGGGGGGTGCTCTTCATCATTATTTTCGC‑GTCATTCTGCTGACCAGGAAACGGGTGAAc < 1:2073098/68‑1 (MQ=255) tgtcgATCTCGGTGGT‑CGCCGTATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAcc < 1:2980227/65‑1 (MQ=255) gcgATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCt < 1:2732091/68‑1 (MQ=255) cGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGc < 1:1678370/68‑1 (MQ=255) gtggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGcaa > 1:2525973/1‑68 (MQ=255) ggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTAAATGCaaaca > 1:3534070/1‑68 (MQ=255) gCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCACCAAAt > 1:3292700/1‑68 (MQ=255) tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATaa > 1:3677682/1‑68 (MQ=255) tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATaa > 1:3446694/1‑68 (MQ=255) tctTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATaa > 1:1522468/1‑68 (MQ=255) | ATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAA > NC_000913/1373973‑1374103 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |