Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,506,020 | 0 | C | A | 39.6% | 72.8 / 21.6 | 48 | G21C (GGC→TGC) | ydcO | BenE family inner membrane putative transporter |
Reads supporting (aligned to +/- strand): ref base C (19/10); new base A (0/19); total (19/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.05e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGAACAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGA > NC_000913/1505872‑1506163 | agaACAGGTACGCGATCCCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCa < 2:471139/150‑1 (MQ=255) aaCAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAAtt > 2:947023/1‑151 (MQ=255) aaCAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAAtt < 2:39059/151‑1 (MQ=255) aCAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑a < 2:1044006/151‑1 (MQ=255) cgcgATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGc > 2:92052/1‑151 (MQ=255) ccATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAAc < 2:733776/151‑1 (MQ=255) tAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAAcccc > 2:631292/1‑151 (MQ=255) gagTGAGGGTACTGGCCCCCATGGCCAGCCCCAGCGCCGTCATCCAGCCAGAGTTTTGTGCGGGGGTGGCCCGGGCGACAATCGCTTTTTGCCAAATTTTTGCCCCTAAAGGGGGG‑TAGCAAATA‑AATACCGCGAAAACCCCCGCCAGTa < 1:613435/150‑1 (MQ=255) cggCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCACCCAGCCAGAGTTTTGGGCGTGGGGGGCTCCGGAGAAAACCGCGGTTTGCCAGTTTATTCCCGCTGAACTGGCG‑TAACAAATA‑AATCCCGCCAGAAACCCCGCCAGTAGCgtg < 1:641350/149‑1 (MQ=255) gTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCgtgggtg > 2:1073676/1‑151 (MQ=255) tACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCgtgggtgg > 2:787595/1‑151 (MQ=255) aCTGCCCCCTATTGCCCCCCCCGGCCCCGCCACCAACCAAGAGTTTGGGGGGGTGGTGGCCCGGCCAAAAACCGCTCCTGCCAAGATTATTGCCGCGAAACGGGGG‑TA‑CCAATTAAATCCCGCAAGAACCCCCGCCAGTAGCgtgggtg < 1:761439/149‑1 (MQ=255) tGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTggagg > 1:339541/1‑151 (MQ=255) gACCCCCATGGCCCCCCCCGGCCCCTTCACCCAGCCAGAGATTTGTGCAGGGGTGGCTCCGGCGACACCCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCAAATA‑AATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTggagg < 1:461723/150‑1 (MQ=255) cGCCCATGGCCAGCCCCGGCGCCGCCCCCCAGCGAGAAATTTGTGCAGGGGTGGCTCCGGCAACAATCGCTTCTTGCCAGATTATTCCCGCTGATCGGGCG‑TAGCCATTT‑AATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTggagg < 2:164003/148‑1 (MQ=255) ccccAGCGCCTTCACCCAGCCAGAGTTTTGTGCGGTGGGGGCCCGGGCAACAACCGCTGCTGCCCAGATTATTCCCGCTGAACGGGGG‑TAGCCATTT‑AATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGa < 2:975086/148‑1 (MQ=255) cgcCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTc > 2:799404/1‑151 (MQ=255) ccccGCCACCCGGCCAGAGTTTTGTGCGGGGGGGGCCCCGCCAACATTCGCTGTTTCCAAGATTTTGCCCGCTAACGTGGGG‑TAGCAAATT‑AAACCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGt < 2:816961/148‑1 (MQ=255) gTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCt > 1:887339/1‑151 (MQ=255) ccagcaagAGATTTGTGCGGTGGGGGCTCCGGCGAAAATCGCTGCTTGCAAGATTATTGCCGCTGAACTGGCG‑TACCAAATA‑AATCCCCCCAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCg < 2:366732/148‑1 (MQ=255) gccagAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGtt > 2:1039516/1‑151 (MQ=255) ccagAGTTTTGTGCAGGGGTGGCCCCGGCGACAATCCTTGCTTCCCATATTATTGCCGCTGAACTGGCG‑TAGCAAATT‑AACCCCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTta < 2:69684/151‑1 (MQ=255) caggagtTTTGTGCGGTGGTGGCCCGGGAGAAATCCGCTGTTTGCCAGATTATTGCCGCTGACCTGGCG‑TACCAAATA‑AATCCCCCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGt < 2:409491/146‑1 (MQ=255) gagtTTTGGGCGGTGGTGCCCCCGGCGACACTCGCTGCTTCCAAGTTTATTGCCG‑TGGACGGGGGTTACCAAATA‑AATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGtt < 1:841901/147‑1 (MQ=255) gagtTTTGGGCAGTGGTGGCTCCGGCGACAATCGCTGTTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCAAATT‑AATCCCGCAAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat < 2:475361/149‑1 (MQ=255) gagaTTTGTGCGGGGGTGGCTCGGGCAACAACCGCTGCTGGCCAGATTATTGCCGCTGAACTGGCG‑TAGCAAATA‑AACCCCGCAAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat < 1:1019169/149‑1 (MQ=255) gaTTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg > 1:286070/1‑151 (MQ=255) gaTTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg > 1:837251/1‑151 (MQ=255) tttGTGCAGTGGTGGCCCCGGCGACAATCGCTGCTTGCCAGATTATTCCCGTTGAACGGGCG‑TAGCAAATA‑AATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGtt < 2:335558/143‑1 (MQ=255) gtgCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCAGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTAc > 2:161613/1‑149 (MQ=255) ggtggCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGt > 2:1009542/1‑151 (MQ=255) ggtggCTCCGGCGAAAATCGCTGTTTGCCAGATTATTGCGGCTGAACGGGCG‑TAGCAAATA‑AATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAcc < 1:359496/149‑1 (MQ=255) tCCGGCAACAATCGCTGTTTGCCAGTTTTTTGCCGCTGAACGGGTG‑TCGCCAATT‑AATCCCGCAAAACACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAg < 2:194514/147‑1 (MQ=255) cGGCGAAAATCGCTGTTTGCTATTTTTTTGCCGCATAACGGGCG‑TAACAAATA‑AATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGc < 1:640615/146‑1 (MQ=255) ggCGACATTCGCTGCTTGCCAGTTTATTGCCGCTAACCTGGCG‑TAGCAATTT‑AATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGt < 2:1099319/151‑1 (MQ=255) ggCAACAATCGCTGCTTGCCAGATTATTCCCGCGGAACGGGGG‑TGGCAAATA‑AACACCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGt < 1:486331/151‑1 (MQ=255) aCAACCGCGGCTTGCAAGTTTATTGCGGGTGAGCGGGCG‑TACCCAATA‑AATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTcc < 1:813974/149‑1 (MQ=255) tgctTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGc > 1:149056/1‑151 (MQ=255) gctTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGc > 1:658790/1‑150 (MQ=255) ccAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGt > 2:651538/1‑151 (MQ=255) caaattttttGCCGTTAAACTGGGT‑TATCAATTT‑AAAACCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAg < 2:498763/142‑1 (MQ=255) caaaatttttCCCGCTGAGCGGGCG‑TACCCAATT‑AACCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcaca < 2:770481/144‑1 (MQ=255) cAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTg > 1:932917/1‑151 (MQ=255) aCTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGc > 2:1011177/1‑151 (MQ=255) tGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCg > 2:888143/1‑151 (MQ=255) ggCG‑TACCAAATA‑AATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCg < 2:780697/150‑1 (MQ=255) gCG‑TAGCAAATA‑AACCCCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCgga < 2:439759/151‑1 (MQ=255) gCG‑GAGCAAATA‑AACCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtt < 2:649971/143‑1 (MQ=255) | AGAACAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGA > NC_000913/1505872‑1506163 |
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: ‑ |