Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 300,711 | 0 | G | C | 44.7% | 43.5 / 20.3 | 38 | P75R (CCG→CGG) | paoC | PaoABC aldehyde oxidoreductase, Moco‑containing subunit |
Reads supporting (aligned to +/- strand): ref base G (5/16); new base C (17/0); total (22/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.04e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCG > NC_000913/300574‑300857 | cGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTgcgc < 1:6715/151‑1 (MQ=255) gAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATTAGAGCCCAGACGCCCCGCGCctttt > 1:699316/1‑146 (MQ=255) gTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGGATTACAACCCGTACCCCCCGCGCCTTTTTTgcggcg > 2:562269/1‑150 (MQ=255) gTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGAAAGCCCGTACGCCCCGCGCcttttttgggg > 2:570799/1‑142 (MQ=255) gCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTAt < 2:11373/151‑1 (MQ=255) cACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCa < 2:402938/151‑1 (MQ=255) cgcAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTg < 2:693893/151‑1 (MQ=255) aaTGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGGGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGGATTAAAACCCCGACCCCCCGCGCCTTTTTCGCGGCG‑GCCCGATCCA‑GGGCggggggg > 2:343389/1‑146 (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:737866/149‑1 (MQ=255) ggCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑ggcgg < 1:660525/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) | CGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCG‑TCCGTATCAA‑GGGCGGTGAGGCGTCCTTTGGCAA‑TGGCGGAACCGACGATATAGCCATA‑GGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAGTTTTCAGCGGTCCG > NC_000913/300574‑300857 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |