Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,506,035 | 0 | G | A | 40.4% | 44.5 / 23.9 | 47 | L16L (CTG→TTG) | ydcO | BenE family inner membrane putative transporter |
Reads supporting (aligned to +/- strand): ref base G (19/9); new base A (0/19); total (19/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.46e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAA > NC_000913/1505889‑1506178 | ccATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑c < 2:733776/151‑1 (MQ=255) tAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑cccc > 2:631292/1‑151 (MQ=255) gagTGAGGGTACTGGCCCCCATGGCCAGCCCCAGCGCCGTCATCCAGCCAGAGTTTTGTGCGGGGGTGGCCCGGGCGACAATCGCTTTTTGCCAAATTTTTGCCCCTAAAGGGGGG‑TAGCAAATA‑AATACCGCGAAAAC‑CCCCGCCAGTa < 1:613435/150‑1 (MQ=255) cggCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCACCCAGCCAGAGTTTTGGGCGTGGGGGGCTCCGGAGAAAACCGCGGTTTGCCAGTTTATTCCCGCTGAACTGGCG‑TAACAAATA‑AATCCCGCCAGAAA‑CCCCGCCAGTAGCgtg < 1:641350/149‑1 (MQ=255) gTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCgtgggtg > 2:1073676/1‑151 (MQ=255) tACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCgtgggtgg > 2:787595/1‑151 (MQ=255) aCTGCCCCCTATTGCCCCCCCCGGCCCCGCCACCAACCAAGAGTTTGGGGGGGTGGTGGCCCGGCCAAAAACCGCTCCTGCCAAGATTATTGCCGCGAAACGGGGG‑TA‑CCAATTAAATCCCGCAAGAAC‑CCCCGCCAGTAGCgtgggtg < 1:761439/149‑1 (MQ=255) tGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTggagg > 1:339541/1‑151 (MQ=255) gACCCCCATGGCCCCCCCCGGCCCCTTCACCCAGCCAGAGATTTGTGCAGGGGTGGCTCCGGCGACACCCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCAAATA‑AATCCCGCAAAAAC‑CCCCGCCAGTAGCGTGGGTggagg < 1:461723/150‑1 (MQ=255) cGCCCATGGCCAGCCCCGGCGCCGCCCCCCAGCGAGAAATTTGTGCAGGGGTGGCTCCGGCAACAATCGCTTCTTGCCAGATTATTCCCGCTGATCGGGCG‑TAGCCATTT‑AATCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTggagg < 2:164003/148‑1 (MQ=255) ccccAGCGCCTTCACCCAGCCAGAGTTTTGTGCGGTGGGGGCCCGGGCAACAACCGCTGCTGCCCAGATTATTCCCGCTGAACGGGGG‑TAGCCATTT‑AATCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGa < 2:975086/148‑1 (MQ=255) cgcCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTc > 2:799404/1‑151 (MQ=255) ccccGCCACCCGGCCAGAGTTTTGTGCGGGGGGGGCCCCGCCAACATTCGCTGTTTCCAAGATTTTGCCCGCTAACGTGGGG‑TAGCAAATT‑AAACCCCCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGt < 2:816961/148‑1 (MQ=255) gTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCt > 1:887339/1‑151 (MQ=255) ccagcaagAGATTTGTGCGGTGGGGGCTCCGGCGAAAATCGCTGCTTGCAAGATTATTGCCGCTGAACTGGCG‑TACCAAATA‑AATCCCCCCAAAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCg < 2:366732/148‑1 (MQ=255) gccagAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGtt > 2:1039516/1‑151 (MQ=255) ccagAGTTTTGTGCAGGGGTGGCCCCGGCGACAATCCTTGCTTCCCATATTATTGCCGCTGAACTGGCG‑TAGCAAATT‑AACCCCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTta < 2:69684/151‑1 (MQ=255) caggagtTTTGTGCGGTGGTGGCCCGGGAGAAATCCGCTGTTTGCCAGATTATTGCCGCTGACCTGGCG‑TACCAAATA‑AATCCCCCAAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGt < 2:409491/146‑1 (MQ=255) gagtTTTGGGCGGTGGTGCCCCCGGCGACACTCGCTGCTTCCAAGTTTATTGCCG‑TGGACGGGGGTTACCAAATA‑AATCCCGCAAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGtt < 1:841901/147‑1 (MQ=255) gagtTTTGGGCAGTGGTGGCTCCGGCGACAATCGCTGTTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCAAATT‑AATCCCGCAAAAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat < 2:475361/149‑1 (MQ=255) gagaTTTGTGCGGGGGTGGCTCGGGCAACAACCGCTGCTGGCCAGATTATTGCCGCTGAACTGGCG‑TAGCAAATA‑AACCCCGCAAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat < 1:1019169/149‑1 (MQ=255) gaTTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg > 1:837251/1‑151 (MQ=255) gaTTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg > 1:286070/1‑151 (MQ=255) tttGTGCAGTGGTGGCCCCGGCGACAATCGCTGCTTGCCAGATTATTCCCGTTGAACGGGCG‑TAGCAAATA‑AATCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGtt < 2:335558/143‑1 (MQ=255) gtgCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCAGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTAc > 2:161613/1‑149 (MQ=255) ggtggCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGt > 2:1009542/1‑151 (MQ=255) ggtggCTCCGGCGAAAATCGCTGTTTGCCAGATTATTGCGGCTGAACGGGCG‑TAGCAAATA‑AATCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAcc < 1:359496/149‑1 (MQ=255) tCCGGCAACAATCGCTGTTTGCCAGTTTTTTGCCGCTGAACGGGTG‑TCGCCAATT‑AATCCCGCAAAACA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAg < 2:194514/147‑1 (MQ=255) cGGCGAAAATCGCTGTTTGCTATTTTTTTGCCGCATAACGGGCG‑TAACAAATA‑AATCCCGCAAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGc < 1:640615/146‑1 (MQ=255) ggCGACATTCGCTGCTTGCCAGTTTATTGCCGCTAACCTGGCG‑TAGCAATTT‑AATCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGt < 2:1099319/151‑1 (MQ=255) ggCAACAATCGCTGCTTGCCAGATTATTCCCGCGGAACGGGGG‑TGGCAAATA‑AACACCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGt < 1:486331/151‑1 (MQ=255) aCAACCGCGGCTTGCAAGTTTATTGCGGGTGAGCGGGCG‑TACCCAATA‑AATCCCGCCAGAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTcc < 1:813974/149‑1 (MQ=255) tgctTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGc > 1:149056/1‑151 (MQ=255) gctTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGc > 1:658790/1‑150 (MQ=255) ccAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGt > 2:651538/1‑151 (MQ=255) caaattttttGCCGTTAAACTGGGT‑TATCAATTT‑AAAACCGCAAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAg < 2:498763/142‑1 (MQ=255) caaaatttttCCCGCTGAGCGGGCG‑TACCCAATT‑AACCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcaca < 2:770481/144‑1 (MQ=255) cAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTg > 1:932917/1‑151 (MQ=255) agacgggtggtgcaaaTT‑AA‑ACCCCCAAAAACCCCCGCCAGTTGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtttt < 1:718448/137‑1 (MQ=255) aCTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGc > 2:1011177/1‑151 (MQ=255) tGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCg > 2:888143/1‑151 (MQ=255) ggCG‑TACCAAATA‑AATCCCGCAAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCg < 2:780697/150‑1 (MQ=255) gCG‑TAGCAAATA‑AACCCCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCgga < 2:439759/151‑1 (MQ=255) gCG‑GAGCAAATA‑AACCCCGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtt < 2:649971/143‑1 (MQ=255) tt‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCAt > 2:583997/1‑151 (MQ=255) aaTACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATgg > 2:746646/1‑151 (MQ=255) ccGCCAAAAC‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGaaa < 1:813481/150‑1 (MQ=255) | CCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAA‑CCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGAGGAGTCGCATGGAAA > NC_000913/1505889‑1506178 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |