Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,461,137 | 0 | C | T | 42.9% | 31.1 / 12.6 | 35 | S82F (TCC→TTC) | paaJ | 3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase |
Reads supporting (aligned to +/- strand): ref base C (9/11); new base T (15/0); total (24/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.50e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGCCAG > NC_000913/1461005‑1461287 | tgGTGCGAAACCCGCGTCTCGATGCGGCGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCCGCCGCAGAATGTTTTC‑GGCACAAACCCTAAcc > 1:78596/1‑150 (MQ=255) ccGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGGGGTTTC‑GGCACAACCATTAACCC‑CTTGTGTgg > 1:345440/1‑149 (MQ=255) ctcGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGCGGCGCCGAGGGGTTTC‑GGGGAAAACAATAAACG‑CTTGTGTGGGTCCcggg > 1:51382/1‑150 (MQ=255) cGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCCGAGTGTTTTC‑GGGGCAACCATTAAACG‑CCTGGGTGGGTCCcgg > 2:318020/1‑148 (MQ=255) gCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGGCGCCGAGGGTTTCCGGGCCCACCCATAAACC‑‑CCTGGGTGGGTCCggg > 2:229868/1‑144 (MQ=255) ggggtgtATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGGGTTTCC‑GGCCCAAACATTAACCG‑CTTGTGTGGTTTCGGGCTgg > 1:186964/4‑146 (MQ=255) gAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGCGCCGCAGAGTGGTTCC‑CGCGCCACCACTAACCG‑CCTGTGTGGTTCCGGGCTGGGCg > 2:513034/1‑148 (MQ=255) aTCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGGGGTTTC‑GGCCCAACCATTAACCC‑CTTGTGGGGGTCCCGGGGGGGCGCACCgggg > 1:1416/1‑150 (MQ=255) cGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGtt < 1:20434/150‑1 (MQ=255) cGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGGGTTTTC‑CGCACAACCCTTAAACG‑CTTGTGTGGGTTCGGGCTGGGCGCCCTGGGGgtt > 1:157488/1‑151 (MQ=255) gatgTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCg < 2:626012/151‑1 (MQ=255) gCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGCGCCGCCGAGGGGTTTC‑GGGGCAACCATTAACCG‑CTTGGTGGGGTCCGGGGCGGGCGCCCGGGGGGTTGCCGCCCgg > 2:91890/1‑143 (MQ=255) gCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGCGCCGCCGAGTGTTTTC‑GGCCCAAACATTAAACG‑CTTGTGTGGGTCCGGGGTGGGCGCAATGGGGGTTTCCGCCCCGGGGGTTAAAACgg > 1:380738/1‑151 (MQ=255) cTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGCGCCGCCGAGTGTTTTC‑TGCCCAACCCTTAACCC‑CGTGGGGGGGTCCCGGGGGGGGGGCCTTGGGGTTGCCCCACGGGTGAGTTAAGCggg > 2:56012/1‑151 (MQ=255) aaTCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGggc < 1:680158/150‑1 (MQ=255) aTCAGGCGGGAGAAGATAACCGTAACGTAG‑CCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGTGTTTTC‑GGGACAACCATTTACCGCCTGGTGGGGTCCCGGG‑TGGGCGCACTGGGGGTTGCCGCACGGGCGGTTAAAACggg > 1:322230/1‑147 (MQ=255) gagaAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGAttt < 2:477464/151‑1 (MQ=255) gaagaTAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGGGTTTTC‑GGCACAAACATTAAACC‑CTTGTGTGGGTTCCGGCTGGGCGCCCCGGGGGTTGCCGCACGGGCGGTTTAAACGGGGGGTGGGGAttt > 1:556657/1‑149 (MQ=255) aaCCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGTGTTTCC‑GGCCCAACCATTAACCG‑CCTGTGGGGGTCCGGGCTGGGCGCACTGGGGGTTGCCGCCCGGGCGGTTAAAGCGGGGGGTGGGGCGTTGCTGAt > 2:18635/1‑149 (MQ=255) aCCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGTGTTTTC‑GGCACCACCATTAACCG‑CCTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGGGATTAAAACGGGGGGTGGCGATTTTCTGAt > 1:777930/1‑148 (MQ=255) ccGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGGTGCCGCAGAGTGTTTTC‑CGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGGTTTCCGCACGGGCGGTTAAAACGGGCGATGGCGGTTTGCTGATCGc > 1:813591/1‑150 (MQ=255) aGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGt < 1:801520/150‑1 (MQ=255) aGCCCGGATGGCGACTTTACTGGCGGGGCCGCCCCCAAGGGGTTCC‑GGCCCAACCCTTAACCG‑CTTGTGTGGTTCCGCGGCGGGCGCACTGGGGGTTGCCCCACCGGGGGTTTAAGCGGGCGGCGGCGACTTGCTGGTccccccgggcggg > 2:43090/1‑142 (MQ=255) cccGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGa < 1:87407/151‑1 (MQ=255) ccGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAg < 2:133194/151‑1 (MQ=255) cGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTTC‑GGCACAACCATTAACCG‑CTTGTGTGGGTCCGGGCTGGGCGCACTGGGGGTTGCCGCACGGGCGATTAAAACGGGCGATGGCGATTTGCTGATCGCCGGGGGCGTGGAg > 1:354049/1‑150 (MQ=255) tGGCGACTTTACTGGCGGGGGTGCCGCAGAGGGTTTTC‑CGGACAAACATTAAACG‑CTTGTGTGGGTTCGGGCGGGACGCACTGGGGGTTGCCGCACGGGCGGTTTAAGGGGGCGGGGGCGATTTGCTGATCGCCGGGGGCGGGGGGGCaa > 2:523331/1‑150 (MQ=255) cGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGtt < 2:699829/151‑1 (MQ=255) ggggCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGttt < 2:408143/151‑1 (MQ=255) ggggCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGttt > 1:657166/1‑151 (MQ=255) gCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGTTGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATg < 1:318078/151‑1 (MQ=255) ccGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATgg < 1:102319/151‑1 (MQ=255) gCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGc > 2:88643/1‑151 (MQ=255) agagTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCa > 1:36577/1‑150 (MQ=255) ttCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGcc > 1:35484/1‑151 (MQ=255) cc‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGCCAg < 2:826112/151‑1 (MQ=255) | TGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGCCAG > NC_000913/1461005‑1461287 |
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: ‑ |