| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 300,721 | 0 | G | T | 38.2% | 42.4 / 13.1 | 34 | Q72K (CAA→AAA) | paoC | PaoABC aldehyde oxidoreductase, Moco‑containing subunit |
| Reads supporting (aligned to +/- strand): ref base G (7/14); new base T (13/0); total (20/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.60e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGC‑GCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATGCGGTCAT > NC_000913/300579‑300870 | gTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCTGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGC‑GCTTTTTGCGCGGCGt < 1:356905/151‑1 (MQ=255) cGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGC‑GCTTTTTGCGCGGCGTcc < 2:197126/149‑1 (MQ=255) ccACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGC‑GCTTTTTGCGCGGCGTCCGTATc < 1:444918/151‑1 (MQ=255) aCTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAAAGACAGCCAGTACGCCCCGC‑GCCTTTTTTGCGGCGGCCGGATc > 2:499003/1‑149 (MQ=255) cAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTGATTACAGCCCCGACGCCCCGC‑GCCTTTTTCGCGGGGGCCGGATC‑AAggg > 1:257664/1‑149 (MQ=255) aGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGGAACCCCCGC‑GCCTTTTTCGCGGCGGCCGGATC‑CAAGGCgg > 2:503397/1‑151 (MQ=255) tGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATTAAAGCCAATACGCCCCGC‑GCTTTTTTTGCGGGGGCCGTATC‑AAGGGGGGGGAgggg > 2:566240/1‑149 (MQ=255) gataatAGTGGTCAATAGAGGGGCCGCCTAACAGCATGGCGGTTTTTTTGTCGACTTTTCAGAGTGACCCTGCTATACTGTTGTTAATGACAGCCAGTACGCCCGGC‑GGTTTTTGCGCGGCGTCCGCATC‑AAGTGCGTTGAGGCGACCtt < 1:540108/150‑1 (MQ=255) cTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGTATTACAGACAAGAAGGCCGGCGGCTTTTTTTCCGG‑GGCCGGATT‑AAGGGCGGGGGGGGGGCCTTTGGcc > 1:350193/1‑144 (MQ=255) gggCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTGCTGGCGGTAATGACAGCCAGTACGCCCGGC‑GCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACcgacga < 1:100789/151‑1 (MQ=255) gccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAACCAGTACGCCCCGC‑GCTTTTTTCGCGGCGGCCGGATCAAAGGGGGGGGGGGGGCCCTTGGGCAA‑GGC‑GGGAAcgacga > 2:539591/1‑149 (MQ=255) cgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATTGAAGCCCGCACGCCCCGC‑GGCTTTTTCGCGGCCGCCCTATC‑CAGGGCGGGGGGGCGGCCCTTTGGGATGGC‑GGGACccacgaaa > 1:297312/1‑147 (MQ=255) ccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGC‑GCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGc < 1:477716/151‑1 (MQ=255) aCAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATTTAAGCCAGTACCCCCCGC‑GCCTTTTTCGCGGCGTCCGGATC‑CAGGGGGGGGAGGGGGCCCTTTGCCATGGCGGGAACcgaggtataagca > 1:183925/1‑139 (MQ=255) cAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGC‑GCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAg < 2:634777/151‑1 (MQ=255) ggCGGTGTTTTTGTCTCCTGTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGC‑GCTTTTTGCGCGGCGTCGGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAGgcggcg < 2:243424/151‑1 (MQ=255) ggCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGC‑GCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAggcggc < 2:309933/150‑1 (MQ=255) cGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAACCAGTACGCCCCGC‑GCTTTTTGCGCGGCGGCCGGATC‑AAGGGCGGGGAGGCGGCCTTTGGCCATGGC‑GGGAACGAAGATATAGCCAAAGgcggcggt > 2:80299/1‑149 (MQ=255) tgtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATTACAGCCCGTACGCCCCGC‑GCCTTTTTCGCGGCGTCCCTATC‑AAGGGGGGGGGGGCGTCCCTTGGCAATGGC‑GGGACCGGCGATATAGCCACAGGCGGCGGTggg > 1:533868/1‑151 (MQ=255) tttttGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGC‑GCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAGGCGGCGTTgggg < 2:334984/150‑1 (MQ=255) tttttGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGC‑GCTTTTTTCGCGGCGGCCGGATC‑AAGGGGGGGGAGGCGGCCCTTGGCAATGGC‑GGGAACGGCGAGAAAGCCAAAGGGGGGGGTggggg > 1:679108/1‑151 (MQ=255) tttttGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGGAAGCCCCGC‑GCTTTTTTCGCGGGGGCCGGATC‑AAGGGGGGGGAGGGGTCCTTTTGCAATGGG‑GGGACCGACGAAATAGCCATAGGCGGCGGTggggg > 2:562726/1‑151 (MQ=255) cGAGTGCCCCGGCGTTACTGGCGGTAAAGACAGCCCGTACGCCCCGC‑GCCTTTTTCGCGGCGTCCGCATC‑CAGGGCGGGGAGGCGGCCTTTGGCAATGGC‑GGGACCGCCCATATAGCCATAGGCGGCGGTTGGGGGGTCTTCATGCCCtt > 1:158138/1‑150 (MQ=255) ggCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGC‑GCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAg < 1:262112/151‑1 (MQ=255) gCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGC‑GCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGc < 1:322670/151‑1 (MQ=255) ggcaaTGACAGCCCGTACCCGCGGC‑GCTTTTTTCGCGGCGTCCGTCTC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGCACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCCTGCCATTCGTAGGCGt > 1:205106/4‑137 (MQ=255) gTAATGACAGCCAGTACGCCCGGC‑GCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCgg > 2:223895/1‑151 (MQ=255) aaTGACAGCCAGTACGCCCGGC‑GCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTa < 2:326457/151‑1 (MQ=255) aGTACGCCCGGC‑GCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCg > 1:286036/1‑151 (MQ=255) tACGCCCGGC‑GCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGt > 2:447586/1‑151 (MQ=255) gCCCGGC‑GCTTTTTGCGCGGGGGCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTct < 1:366969/150‑2 (MQ=255) gcTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATGc > 2:331762/1‑151 (MQ=255) ttttGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATGCGGt > 2:518862/1‑151 (MQ=255) ttGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATGCGGTc > 1:564325/1‑150 (MQ=255) gcgcGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATGCGGTCag < 1:472050/150‑2 (MQ=255) | GTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGC‑GCTTTTTGCGCGGCGTCCGTATC‑AAGGGCGGTGAGGCGTCCTTTGGCAATGGC‑GGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATGCGGTCAT > NC_000913/300579‑300870 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |