Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 298,960 | 0 | A | C | 45.0% | 50.9 / 12.5 | 40 | Y659D (TAT→GAT) | paoC | PaoABC aldehyde oxidoreductase, Moco‑containing subunit |
Reads supporting (aligned to +/- strand): ref base A (17/5); new base C (0/18); total (17/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.35e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCAC > NC_000913/298814‑299103 | ataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATAcccc < 1:704310/151‑1 (MQ=255) acCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCa < 1:177404/151‑1 (MQ=255) gCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGAt > 1:820211/1‑151 (MQ=255) cgcgCTCGCGCCGCACAGCCCCAGCTCCCCGACCCCTTTGGCCTTCACCGGGGAGGATATGGGGTCGGTATCATCCAGAAAAACCCCCTCCTGTTTTGGA‑ATATCCGCATGAACCGGCCCCTCATACCCCGCCATATCGTGATTAAc < 1:532029/147‑1 (MQ=255) cTCACGCCGCACAGGCCCAGCTCCCCGACCCCTTTGGCCTTCATCGGGGAGGTTAGGGGGTCGGTATCACCCAGAAAAATCCCCTCCTTTTTTGGG‑ATACCCGCAGAAACCGGCCCCTCATCCCCCGCCATATCGTGATTAAc < 1:84991/143‑1 (MQ=255) cTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGc > 1:696522/1‑151 (MQ=255) cgccgcACAGGCCCAGCTCACCGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGAAAAATCCCCTCCTGTTTTGGG‑ATACCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCaa < 1:612540/151‑1 (MQ=255) ccgcACAGGCCCAGCCCCCCGCCCCCTTGGGCCTCAGTCGGGGAGGATAGGGGGTCGTTATCACCCAGAAAAACCCCCTCCTGTTTTGGG‑ATACCCACAGGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCaaa < 2:715232/150‑1 (MQ=255) cgcACAGGCCCACCTCCCCCACACCTTTGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGCCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACg < 2:634824/151‑1 (MQ=255) gCCCAGCTCACCGACAACTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATcc > 2:393400/1‑151 (MQ=255) ccAGCTCACCGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCAc < 1:645669/151‑1 (MQ=255) ccAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCAc > 2:304426/1‑151 (MQ=255) ttcgtcttacagcGGGAAGGATATGGGTGCGTTCCTATCCATGAAATCCCCCCCCTTTTTTGGC‑CTTTCCACGGGAACCGACCCTCCATACCCCGCCATATCGTGATTAACGAATTAGCCCAACCGGTCATCCAACGCCAGCccctcc < 2:286143/136‑1 (MQ=255) ggCCTTCATCGGGGAGAATATGGGGTCGGTATCACCCAGAAAAATCCCCTCCTGTTTTGGG‑ATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGt < 1:66962/151‑1 (MQ=255) ttCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcg > 1:354999/1‑151 (MQ=255) cATGGGGGAGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTTTTTGGGA‑ATACCCGCATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgc < 1:734734/150‑1 (MQ=255) ggggAGGATATGGGGTCGGTACCACCCAGAAAAACCCCCCCCTGTTTTGGG‑ATACCCGCATGAACCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccca < 2:65963/151‑1 (MQ=255) tatGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat > 2:369004/1‑150 (MQ=255) tatGGGGGCGTAATCACCCA‑GAAAAACCCCCCCTGGTTTGGGAATATCCGCAGGACCGGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgccc < 2:422850/148‑1 (MQ=255) tGGGGTGGGTATCAACCAGGAAAACCCCCTCTTTTTTTGGG‑ATAGCCTCAGGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat < 2:573578/148‑1 (MQ=255) ggTCGGTATCATCCAGAAAAATCCCCCCCTGTTTGGGG‑ATATCCGCATGACCGGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCAt < 1:553054/151‑1 (MQ=255) gttATCATCCAGAAAAATCCCCTCCTGTTTTGGG‑ATCTCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcg < 2:814243/149‑1 (MQ=255) ggTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcg > 2:536402/1‑151 (MQ=255) atcCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCACGCCCATAGTCATTGCGCCAAt > 2:514445/1‑150 (MQ=255) gaaaaaTCCCCTCTTGTTTGGGG‑ATGTCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg < 1:550016/144‑1 (MQ=255) cACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgc > 1:538215/1‑151 (MQ=255) cctcttTTTTGGGA‑ATACCCGAATAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcg < 1:732517/150‑1 (MQ=255) tGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt > 1:798878/1‑151 (MQ=255) ggtttGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATCACCTGGCCGCGCGCAGattt > 2:192843/3‑151 (MQ=255) gTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt > 2:690565/1‑151 (MQ=255) atCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg > 1:755650/1‑151 (MQ=255) aTGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGt > 2:586223/1‑151 (MQ=255) tGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTc > 2:516927/1‑151 (MQ=255) gAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTcc > 2:309060/1‑151 (MQ=255) ggCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgca < 1:210959/150‑1 (MQ=255) gcccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgca < 1:236137/146‑1 (MQ=255) ccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt < 1:85023/143‑1 (MQ=255) ccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgcagca < 1:754660/149‑1 (MQ=255) ccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgcagca < 1:2357/149‑1 (MQ=255) aCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac > 1:267214/1‑151 (MQ=255) | ATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCAC > NC_000913/298814‑299103 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |