Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 300,711 | 0 | G | C | 58.3% | 5.1 / 27.9 | 24 | P75R (CCG→CGG) | paoC | PaoABC aldehyde oxidoreductase, Moco‑containing subunit |
Reads supporting (aligned to +/- strand): ref base G (3/7); new base C (14/0); total (17/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.47e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-03 |
CGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAG > NC_000913/300566‑300843 | cgcCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCCGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTACTGGCAGCCCATACCCCCCgcgc > 1:549052/1‑150 (MQ=255) gTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCtttttt > 1:96743/1‑150 (MQ=255) gTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCCGTACCCCCCGGGCcttttt > 1:2055/1‑150 (MQ=255) gTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATTACAGCCCGTACCCCCCGCGCTTTTTTTgcggcgg > 1:116105/1‑150 (MQ=255) ccACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATc < 1:130134/151‑1 (MQ=255) cAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAGCCAGGACGCCCGGCGCCTTTTTCGCGGCGGCCGGATCAAgggggggggg > 2:172853/1‑147 (MQ=255) ggCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGCGTCCGTATCAAGGGCGGGGGGGCgg > 2:112531/1‑150 (MQ=255) ggCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAGCCCGTACGCCCCGCGCTTTTTTCGCGGGGGCCGGATCAAgggggggggggggg > 2:188393/1‑143 (MQ=255) tgatgaTAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGAAAGCCAGTACGCCCGGGGCCTTTTTCGCGGCGTCCGGTTCAAAGGGGGGGGGGCGGcc > 1:126854/1‑149 (MQ=255) gCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGGACGCCCCGCGCTTTTTTCGCGGCGGCCGTATCCAGGGGGGGGAGGGGTCCTTTGGCAATggg > 2:537792/1‑150 (MQ=255) cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGGCCGGATCAAGGGCGGGGAGGCGTCCTTTGGCAATGGCggg > 2:253546/1‑150 (MQ=255) cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGGGAGGCGGCCTTTTGCAATGGCggg > 2:551068/1‑150 (MQ=255) tAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAAcc < 1:240870/151‑1 (MQ=255) gccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGGACGCCCCGCGCCTTTTTCGCGGCGTCCGTATCAAGGGCGGGGGGGCGTCCCTTGGGCATGGCGGGACCGGCGATATAg > 1:431124/1‑151 (MQ=255) gtgtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGGGAGGCGTCCTTTGGCAATGGCGGAAACGACGATATAGCCATAGGCGGCGTTgg > 2:343494/1‑151 (MQ=255) ccTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAACCAGTACGCCCCGCGCCTTTTGCGCGGCGCCCGTATCAAGGGCGGGGAGGCGGCCCTTGGCCATGGCGGGACCGACGATATAGCCATAGGGGGCGTTGGGGGcttcttt > 1:493790/1‑148 (MQ=255) tttGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGGGTCCGTCTCAAGGGCGGGGAGGCGCCCTTTTGCCATGGCGGAACCCGCGCTATAACCATAGGGGGCGGTGGGGGGTTCTTCATg > 1:485430/1‑150 (MQ=255) aGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGCGGCCGTATCAAGGGGGGGGAGGCGGCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGGGTTGGGGGGCTCTTCCTGcccctcgg > 1:364256/1‑145 (MQ=255) gCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAg < 1:253542/150‑1 (MQ=255) gCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAg < 1:141708/150‑1 (MQ=255) aCTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGCGTCCGTATCAAGGGCGGGGAGGCGTCCTTTTGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGGTTCTTTCTGCCATTCGTAGGCGTAGcgtgc > 1:457572/1‑150 (MQ=255) cTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtgccg < 2:2055/151‑1 (MQ=255) aTGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAg < 1:551058/151‑1 (MQ=255) aTGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAg < 2:364264/151‑1 (MQ=255) | CGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAG > NC_000913/300566‑300843 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |