| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 300,745 | 0 | T | G | 47.2% | 48.4 / 13.3 | 36 | T64P (ACC→CCC) | paoC | PaoABC aldehyde oxidoreductase, Moco‑containing subunit |
| Reads supporting (aligned to +/- strand): ref base T (3/16); new base G (17/0); total (20/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.35e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATGCGGTCATGGGGACGACCGACAACCTTCAGC > NC_000913/300596‑300893 | cgcAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTg < 2:693893/151‑1 (MQ=255) gatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTgg < 2:236999/150‑1 (MQ=255) gatAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCCGGACCCCCCGCGCCTTTTTCGCGGGG‑TCCGGATCAA‑GGGCGGGGAGGCgggcctt > 1:266290/1‑142 (MQ=255) gCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAAATACAGCCCGTACGCCCCGCGCCTTTTTCGCGGGG‑GCCGGATTAAGGGGGGGGGGGGGGCCCTTTGGCa > 1:275426/1‑147 (MQ=255) cTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAACCCGGACGCCCCGCGCCTTTTTCGCGGGG‑GCCGGATCAA‑GGGCGGGGAGGCGGCCCTTTGCCA‑Tgggg > 1:90153/1‑149 (MQ=255) aTAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGAAAGCCAGTACGCCCCGCGCTTTTTTGGCGGGG‑TCCGGATCAA‑GGGGGGGGAGGGGGCCTTTTGCAA‑TGGCgg > 1:525631/1‑148 (MQ=255) aGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATTACAGCCCGTAAGCCCCGCGCTTTTTTTGCGGGG‑GCCGTATCAA‑GGGGGGGGGGGCGGCCCTTTGCCA‑TGGCGGGAcc > 1:452090/1‑150 (MQ=255) gTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAAcc < 2:548170/149‑1 (MQ=255) ccgccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAAATAAAACCAGTAAGCCCGGCGCTTTTTTTGCGGGG‑GCCGGATCAA‑GGGGGGGGAGGGGGCCCTTTGCAA‑AGGCgg > 1:660422/1‑139 (MQ=255) aCAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAAAGGCAGCCAGGACTCCCCGCGCCTTTTGCGGGGCGCGCCGCACAAG‑GGGGGGGGAGGGGGCCTTTGGCAA‑‑GGGGGAGCCCACGacattagccg > 2:85607/1‑144 (MQ=255) ccTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑ggcg < 1:96835/151‑1 (MQ=255) tGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAACCCGTACGCCCCGCGCCTTTTTCGCGGCG‑GCCGTATCAA‑AGGGGGGGAGGCGGCCCTTGGCAA‑TTGCGGGAACGACGGTATAACCATaagggg > 1:265632/1‑145 (MQ=255) ggCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑ggcgg < 1:660525/149‑1 (MQ=255) ggCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑ggcgg < 1:737866/149‑1 (MQ=255) cGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTTCGCGGCG‑GCCCTTTCAA‑GGGCGGGGGGGCGTCCTTTTGCAA‑TGGCGGGACCGACGATAAAGCCATA‑Agcggcg > 1:598622/1‑149 (MQ=255) cGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGtt < 2:284020/151‑1 (MQ=255) gtgtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGGACGCCCCGCGCCTTTTTCGCGGCG‑GCCGGATCAA‑GGGCGGGGGGGCGTCCTTTGGGCA‑TTGGGGAACCGGCGATATAGCCATA‑AGCGGCGGTgg > 1:526495/1‑151 (MQ=255) tgtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATTACAACCCGTACGCCCCGCGCTTTTTGCGCGGCG‑GCCGGATCCA‑GGGGGGGGGGGGGTCCTTTGGCAA‑TGGCGGGAACCACGGTATAACCATAGGGGGGCGTTgg > 2:692018/1‑151 (MQ=255) gtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTggg < 2:525616/150‑1 (MQ=255) tttttGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAACCAGGAAGCCCGGCGCTTTTTGCGGGGGG‑GCCGGATCAA‑GGGGGGGGGGGGGGCCTTTTGCAA‑AGGGGGAAACGGCGAAAAAGACATA‑GGGGGGGGTggggg > 1:548269/1‑151 (MQ=255) tttGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAGCCAGTACGCCCGGCGCTTTTTTCGCGGCG‑TCCGTATCAA‑GGGCGGGGAGGCGTCCTTTTGCAA‑TGGCGGGACCGACGATATAGCCATA‑GGCGGCGGTggggg > 1:330433/1‑149 (MQ=255) tGCCGAGTGCCCCGGCGTTACTGGCGGTAAAGACAGACAGCACCCCCCGCGCTTTTTTAGGGGCG‑TGCGTATCCA‑GGGCGGGGGGGCGGCCCTTTGCAATTGGGGGAACCG‑CCAACTAAACCTA‑GGGGGGGGTGGGGGCTTTTTCa > 2:684771/1‑146 (MQ=255) gCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCa < 2:499585/150‑1 (MQ=255) ccGAGTGCCCCGGCGTTACTGGCGGTAATGACAACCCGGAAGCCCCGGGCTTTTTGCGCGGCG‑GCCGGATCAA‑GGGCGGGGAGGCGGCCTTTGGCAA‑TGGCGGGACCGACGATAAAGCCATA‑GGCGGGGTTGGGGGGTTCTTTCTTCCAtt > 1:396604/1‑151 (MQ=255) gAGTGCCCCGGCGTTACTGGCGGTAATGACAACCCGTACGCCCCGGGCCTTTTGCGCGGCG‑GCCGGATCAA‑GGGGGGGGAGGGGGCCCTTGGCAA‑TGGCGGGACCGACGATATAGACATA‑GGCGGCGGTTGGGGGTTCTTTATGCCCtt > 2:531584/1‑149 (MQ=255) gTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTa < 1:492257/151‑1 (MQ=255) ccGGCGTTACTGGCGGTAATTAAAAGCAGTAAGCCCCGCGCTTTTTGCGCGGCG‑GCCGGATTAA‑AGGCGGGGAGGGGTCCCTTGGCAA‑AGGCGGGACCGACGCTATAAACAAA‑GGCGGGGGTGGGGGGGTTTTTATGCCCtt > 2:266433/1‑142 (MQ=255) cGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTa < 2:542660/151‑1 (MQ=255) gCCAGTACGCCCGGCGCTTTTTGCGCGGCT‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAGCCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTc < 2:710545/150‑1 (MQ=255) ccAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTc > 2:330875/1‑149 (MQ=255) gCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCg < 2:96586/151‑1 (MQ=255) ttGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATGCGGTc < 2:599678/150‑1 (MQ=255) gcg‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATGCGGTCATGGGgac < 2:263411/151‑1 (MQ=255) gTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTCCTTCATCCCATTCGTAGGGGTAGCGTGCCGGGCCGGTATTTTTACGGGGCCCGCCGCTGGGGTATTGGGgccgac > 1:22051/1‑148 (MQ=255) aTCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATGCGGTCATGGGGAcgaccgac > 1:671228/1‑150 (MQ=255) ggTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATGCGGTCATGGGGACGACCGACAACCTTCAGc < 1:729605/151‑1 (MQ=255) | CGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCGTCGATGCGGTCATGGGGACGACCGACAACCTTCAGC > NC_000913/300596‑300893 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |