Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,461,121 | 0 | C | G | 45.7% | 34.7 / 11.1 | 35 | L77V (CTG→GTG) | paaJ | 3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase |
Reads supporting (aligned to +/- strand): ref base C (6/13); new base G (16/0); total (22/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.88e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTT > NC_000913/1460971‑1461267 | gatgatCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGc < 1:495917/151‑1 (MQ=255) gctgctATCCCTTTGCGGGAACTGTTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCaga < 1:368904/151‑1 (MQ=255) gctATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAgt < 1:247773/150‑1 (MQ=255) gctATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAgt < 1:430642/150‑1 (MQ=255) aTCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGtt < 2:818479/150‑1 (MQ=255) tCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCGGGCGCcgaggggt > 1:309161/1‑145 (MQ=255) 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) | GATGATCTGGCTGCTATCCCTTTGCGGGAACTGCTGGTGCGAAACCCGCGTCTCGATGCGGAGTGTATCGATGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCC‑GGCACAACCATTAACCG‑CTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTT > NC_000913/1460971‑1461267 |
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: ‑ |