Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,160,437 | 0 | A | C | 72.2% | 31.8 / 10.7 | 58 | S161R (AGC→CGC) | PP_3661 | membrane protein |
Reads supporting (aligned to +/- strand): ref base A (1/14); new base C (26/13); total (31/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.29e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTG‑CGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCT > NC_002947/4160302‑4160564 | aCCAGTACCTTGGTTTTCGGCCTGTGGTGTGGCTTGAATGTTTTAATGCTCGACCTCGATACCGCCTTGGTGACCGCGGCGGTCAGCGCCATCTGCGGCGCCGCCGCGGTGCTGTCATTCTAATCGGCGGTG‑CGCa < 2:2487918/136‑1 (MQ=255) cccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGACTCGGGGCGG‑GGGCGCGcccc > 1:1467328/1‑136 (MQ=255) cccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGG‑GGCCGcccccc > 2:2267980/1‑136 (MQ=255) ccTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGG‑CGCCGCGccccc > 2:1805384/1‑136 (MQ=255) tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGACTCGGGGCGG‑GGGCGCGccccccccaa > 1:2601096/1‑132 (MQ=255) tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGG‑GGGCGCGcccccccccc > 1:160073/1‑132 (MQ=255) tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGG‑GGGCGCGcccccccaaa > 1:1031295/1‑132 (MQ=255) tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGG‑CGCccccccccccccca > 1:2028212/1‑132 (MQ=255) tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGG‑CGCCGCGgcccccccaa > 1:1611927/1‑132 (MQ=255) tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCCCCGCCGTGCTGGCATTCGAATCGGCGCCG‑CGCCCCGccccccccac > 1:2664565/1‑132 (MQ=255) tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCAGCGGCGCCGCCCCCGTGCTGGCAATCGAATCGGCGGGG‑GGCCGCGccccgccaaa > 1:2586195/1‑130 (MQ=255) tgATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCACCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCCG‑CGCCGccccccccccaa > 1:1066810/1‑132 (MQ=255) aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCTGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGG‑GGCCGCGCCCCCCCCaaaa > 1:2929839/1‑136 (MQ=255) aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCCATCGAATCGGGGCAG‑CGGCGCGGCCCCCaccaaa > 1:2693648/1‑131 (MQ=255) aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGG‑GGCGGGGCCCCCCCCaaaa > 1:1945905/1‑136 (MQ=255) aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGG‑GGGCGCGcgccccccccaa > 1:1261502/1‑130 (MQ=255) aTCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCCCCGTGCTGGCAATCGAATCGGGGCGG‑GGCGCGGCCCCCCCCaaaa > 1:1187601/1‑136 (MQ=255) ggCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGTCCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCAG‑CGCCGCCCCCCCCCAaaaaacc > 2:2884607/1‑133 (MQ=255) ggCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCTTCTGCGGCGCCGCCGCCCTGCTGGCAATCGAATCGGGGCTG‑CAGCGCTCCCCCCACaaaaacc > 2:654726/1‑133 (MQ=255) ggCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGCG‑CGCCGCCCCCCCCCAaaaaacc > 1:2885154/1‑133 (MQ=255) ggCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGCCTG‑CGCCGCCCCCCCCCTAAAAgcc > 2:1199227/1‑135 (MQ=255) ccTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAG‑CCCCCCGCCCCCCAAAAAAgaccc > 2:278948/1‑132 (MQ=255) ccTGTGGTGTGGCATGAAAATGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCCCAG‑CCCCCCCCCCCCCACaaaacacac > 2:2347426/1‑131 (MQ=255) cTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGG‑CGGGGCGCCCCCCCCaaaaaaaaca > 1:1836739/1‑130 (MQ=255) tgtggtgtggCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGG‑CGCCGCGCCCCCCCAAAAAAAGCCaa > 1:821267/1‑135 (MQ=255) gtggCATTAAAGTGTGAAAGCTGCACCGCGAAACCGCCTTGCTGACGGCCGCGGGCAGCGCCTTCGGGGGCGCCGCCGCGGTGGGGGATTTCTAAGCGGGGCTG‑CGCAGCGCCCCCCATAAAAGCGCCATGGCAGt < 2:2101509/136‑1 (MQ=255) cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGGCCGCGCCCCCCCCCAAAAAACCCCCATGGC‑GTCGGa > 2:1152553/1‑136 (MQ=255) cATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTG‑CGCCCCGCCCCCCCCAAAAACGCCCAGGCAGTCggg > 2:654854/1‑135 (MQ=255) ttcaaaGTGTTCACTCTCGCCCCCCAAACCGCATGGCGGCCCGCCCCCGCCTGCGCCAGCGCCGGCGCCGCCGCCGTGCTGGAATTCCAATCGGCGCTG‑CGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGaa < 2:1642011/133‑1 (MQ=255) aTGAATGTGTTCACGGACGACCGCTAACCCGCCTTGCTGCCCGCCGCCGGCTGCGCCATCGGCGGCGCCGCCGCGGTGCTGGAATTCGACGCGGCGCTG‑CGCCGCCCCCCCCATAAAAGCGCCATGGCAGTCGGaa < 2:2953936/136‑1 (MQ=255) aaaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGCGG‑GGGCGCGCCCCCCCAAAAAACGCCCAGGGAGTCGGGAGCg > 1:1351824/1‑136 (MQ=255) aaGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGGGG‑CGGCGCGCCCCCCCCAAAAGACCCCAGGGAGGCGGAAGCgg > 1:1575132/1‑135 (MQ=255) tgtgtTCATGCGCGACCGCTAACCCGCTTTGCTGCCCGCCGCCGGCAGCGCCATGGGGGGGGCCGCCGCCGTGCGGGAATGCGATGCGGCGCTG‑CGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa < 2:1227617/135‑1 (MQ=255) tgtgtTCAAGCTCGACCGCTATCCCGCGTTGCTCCCCGCCGCGGGCAGCGCAATCTGGGGCGCCGCCGCGGTGGTGCCTTTCTAATGGGCGCGG‑CGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa < 2:1499447/135‑1 (MQ=255) tgtgtTCAAGCTCCCCCGCTATCCCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCGGGCATTCTAAGCGGCGCGG‑CGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTa < 2:289569/135‑1 (MQ=255) gtgtTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGGGGGG‑GGCCGCGCCCCCCACAAAAAGCCCAGGGGAGgcggaggcggag > 1:1933704/1‑136 (MQ=255) tgtTCAAGCTCGACCGCTATCCCGCCTTGCTGCCCGCCGCCGGCTGCGCCATGTGCGGCGCCGCGGCCGTGCTGGCATTCGATTCGGCGCTG‑CGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGt < 2:1662964/136‑1 (MQ=255) tgtTCAAGCTCGACCCCGATACCGCCTTGCTGACCGCCCCCGGCAGCGCCATCTGCGGCGCCGCCGCGGTGGTGGCATTCGAATCGGCGCTG‑CGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGt < 2:1208585/136‑1 (MQ=255) aaGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCGG‑CGCCGCGCCCCCCCAAAAAACCCCAAGGCCGTAGGAAGCGGAGTGcatg > 1:2798286/1‑133 (MQ=255) cccccGCGATCCCGCGTTGCTCCCCGCCGCCGGCAGCGCCATCGGCGGCGCCGCCGCCGTGCTGGAATGCGAAGCGGCGGGG‑CGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGt < 2:2032885/133‑1 (MQ=255) taTACCGCTTTGCTGCCCGCCGCGGCCCGCGCTTTCGGCGGCGCCGCCGCGGTGTGGGCATTCGACGCGGCGCTG‑CGCCGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGt < 2:2609630/135‑1 (MQ=255) tACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCAG‑CGCCGCGCCCCCCAAAAAAACGCCAAGGCAGGGGGAAGCGGAGTGGGGTTGGGGAccgtttgg > 2:1934857/1‑130 (MQ=255) ccccGCCTGGCCCACCCCCGCCGCCGGCGCCAGCTGCGCCCCGGCGGCCGTGCTTGCATTCGATGCGCCGCTG‑CGCCCCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCg < 2:2367425/134‑1 (MQ=255) ccgGCCGCGCCATCGGGGGCGCCGCCGCGGTGCTGCCATCCCAAGCGGCGCGG‑CGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGc < 2:1128995/136‑1 (MQ=255) ccaccggcggcgCCGCCGCCGGGCTGCCATCCGCCCCGGCCCGG‑CGCACCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGc < 1:2346761/130‑1 (MQ=255) gggggcgCCGCGGCCGTGCTGCCTTTCTCCGCGGCGCGG‑CGCACCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggct < 1:988375/133‑1 (MQ=255) gggggcgCCGCCGCGGTGCTGGAATTCGTAGCGGCGCGG‑CGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggct < 2:524254/133‑1 (MQ=255) ggcggcgCCGCCGCGGTGCGGCAATTCGAATCGGCGCGG‑CGCCGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggct < 1:1825461/135‑1 (MQ=255) ggcggcgCCCCCGCCGTGCTGGCATTCGACGCGGCGCCG‑CGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggct < 2:2410358/135‑1 (MQ=255) gggcgCCGCGGGGGTGCTGCCATTCTCCGGGGCGCGG‑CGCCGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctgg < 2:3026804/135‑1 (MQ=255) gggcgCCGCCGCCGGGCTGGCATTCGAATCGGCGCTG‑CGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctgg < 2:2929839/135‑1 (MQ=255) ccgccgccggCGTGCTGGCATTCGAAGCGGCGCTG‑CGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCggctggct < 1:1771179/135‑1 (MQ=255) ggccgccgccGTGCTGGCTTTCGAATCGGCGCTG‑CGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGACGTTCCTTTACCCACTGGCGATCAATGCCGgctggctg < 2:3051598/135‑1 (MQ=255) cgccggcgcGGGGGTGGCATTCGAAGCGGCGCTG‑CGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGgctggctg < 2:2063938/136‑1 (MQ=255) ggcggTGCTGTCTTTCGACGGGGCGCCG‑CGCCGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTg < 1:1577423/132‑1 (MQ=255) tacgcggcgcgg‑cgCCGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCg < 2:1836739/132‑1 (MQ=255) catgcGGCTCTG‑CGCAGCCCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCg < 2:43538/132‑1 (MQ=255) aTCGGCGCAG‑CGCAGCGCCCCCCAAAAAAACGCCAATGCAGTCGGAAGCGTAGTGGTGTTAGGTAACTTGTCGATGTTGCTTTAACCACTGGCGATAAAATCCGGCTGGCTTCATATGGATACCATGGGCGCCggg > 2:1216713/1‑136 (MQ=255) ggcgcgg‑cgcAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctg < 2:2664565/136‑1 (MQ=255) cgcTG‑CCCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGgctgct < 1:1209134/136‑1 (MQ=255) | ACCAGTACCCTGCTGATCGGCCTGTGGTGTGGCATGAAAGTGTTCAAGCTCGACCGCGATACCGCCTTGCTGACCGCCGCCGGCAGCGCCATCTGCGGCGCCGCCGCCGTGCTGGCATTCGAATCGGCGCTG‑CGCAGCGCCCCCCATAAAAGCGCCATGGCAGTCGGAAGCGTAGTGCTGTTCGGTACCTTGTCGATGTTCCTTTACCCACTGGCGATCAATGCCGGCTGGCTGCATCTGGATACCATGGGCGCCGGGCTGCT > NC_002947/4160302‑4160564 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |