| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 67,887 | 2 | . | T | 38.5% | 18.3 / 65.6 | 39 | coding (451/1503 nt) | araA | L‑arabinose isomerase |
| Reads supporting (aligned to +/- strand): ref base . (12/12); new base T (6/9); total (18/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.42e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTG > NC_000913/67748‑68028 | gCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCCTGGGCTTGCTTGGCGTGCCA‑‑GT‑‑g > 2:2284/1‑143 (MQ=255) gCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAAc < 1:573598/143‑1 (MQ=255) aTCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCAc < 1:659212/143‑1 (MQ=255) tCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCaagaa > 2:688029/1‑139 (MQ=255) tGACCGCTACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGt < 2:250326/143‑1 (MQ=255) gACCGCCACGTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGtt < 1:629122/143‑1 (MQ=255) cgtcccGCCTGTTATCGCCAAATCGGCAGACTTTCAGATGACGTGTA‑CCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACg < 2:411847/138‑1 (MQ=255) catggaTGTCATCGCTCAATCGGAAGACTTTGAGATGACGGGTACCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCa < 1:121614/138‑1 (MQ=255) ttATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGCTACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCCTGTTGCCGACGCATAcgcgcg > 2:153660/1‑143 (MQ=255) cGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCACTGT‑‑GAAAGACAAACATGGCATGTTGCTGACGCATACGCGCGCCaa < 1:593580/143‑1 (MQ=255) ccAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCCCCTGGGCTTGTTTATCCTGCCC‑‑GTGGGACCGGGAACCCCGGCCAGTTGCTGGCGTATACGTGCGCCGATg > 1:512486/1‑143 (MQ=255) ggCAGACTTTCAGATGACGGGTATCCTGTTGAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGa < 2:581269/143‑1 (MQ=255) gaaGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGaa < 1:804759/141‑1 (MQ=255) gCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGATGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACt > 1:325721/1‑143 (MQ=255) gACTTTCAGAGGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTcgc < 1:411788/143‑1 (MQ=255) tttCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCACGCCAATGAAACAGAACCcgcc > 2:560944/1‑142 (MQ=255) cAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCAGCCAGGAGCGGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGAcc < 1:103642/143‑1 (MQ=255) gACGGGCATCGTGTTTAGTGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGAACGCCAtgt < 2:37650/143‑1 (MQ=255) cGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGc < 1:165251/143‑1 (MQ=255) gTATCCTGTTTAGAGACCGCCTGACGCAGCCAGGAGCCGATCCGCTCATGCGCTAGTTTATCCTGCCA‑‑GT‑‑GACCGGTCACAACCGGCTGTTGCTGACGCATACGCGCGCCAATGAAGCCGACATCGCGCCCGCCATCGGCAGt > 2:12589/1‑143 (MQ=255) tATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑AT‑‑GACCGTTAACCACGGCATGTTGCGGCCGCCTAGACACGCCTATGCAGTCGAACTCGAGACAGCCAGGGGCCGTc > 2:189233/1‑143 (MQ=255) gTTTAGAGACCGCCTGTCGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACCCGCGACCGCCATGTGCAGTCAggt > 1:797292/1‑143 (MQ=255) gTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttc > 1:156515/1‑143 (MQ=255) tAGGGACCGCCTGACGCATCCAGGAGCCGCTACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttcagg < 2:48172/143‑1 (MQ=255) tGACGCGTCCACGAGCCGATACGCTCTTGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCCGTCACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTc < 2:658274/143‑1 (MQ=255) cgcgTCCCGGGGCCGATACGCTCATGGGCTTGTTTTGCCTGCCA‑‑GT‑‑GGCGGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCat < 2:479264/143‑1 (MQ=255) cAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACAGCCCTGTGCAGTCTGGTGGAGGGTCATAAAGTCAATATCg > 1:158096/1‑143 (MQ=255) aGGAGCCGATACGCTCATGGGCTAGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGAAGCATACGCGCGCCAGTGAAGCCGGACTGGCGACAGCGATGTCCACTGTGGTTCCGCTCCATAACGTCCGAATCGa > 2:402741/1‑143 (MQ=255) caGCCGATACGCTCTTGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATa < 1:270342/142‑1 (MQ=255) ggctgATTCGCGCATGGGCTTGTTTCTCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGATGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTg < 1:376402/139‑1 (MQ=255) cGATACGCTCCTGGGCTTGTTTATCCTTCCA‑‑GT‑‑GCCCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTcc < 1:542117/143‑1 (MQ=255) aTACGCTCATGGGCTTGTTTATCCTGCCAGTGT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACACGTGCAGAGGAATCGGACCTAGGGTCAGCGAAGGTCAGTCTGGTTCAGGTTCATAAAGTCCATCTCGATACTGTcc > 2:169112/1‑143 (MQ=255) gggCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATCAAGTCCATATCGATCCTGTCCCACGGCCGcgcc > 1:28735/1‑143 (MQ=255) cTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCcgcg > 1:340864/1‑143 (MQ=255) ccTGCCG‑‑GT‑‑GCCCGGTACCCACGGCATGTTCCTGCCGCATACGCGGGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTgg < 2:334940/143‑1 (MQ=255) ccTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCGTAGTACCGGTACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTgg < 1:584693/143‑1 (MQ=255) cTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCCCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTggg > 1:504742/1‑143 (MQ=255) tGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCGCGGCAGCGCCGCGGTGAGCTGGgt > 1:701017/1‑143 (MQ=255) cA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCACGGCAATGAACCCGCACTCGCGACGGCCATGAGCCTTCGGGTTCAGCTTCATAAAGTCGAGATCGGTAGTGTAGCACGGCAGCCCCGAGCTGTACttcgtggt > 2:310116/1‑136 (MQ=255) | GCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCA‑‑GT‑‑GACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTG > NC_000913/67748‑68028 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |