Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,628,625 | 0 | A | C | 22.4% | 97.6 / 10.7 | 49 | V105G (GTC→GGC) | pdeF | cyclic di‑GMP phosphodiesterase PdeF |
Reads supporting (aligned to +/- strand): ref base A (31/7); new base C (0/11); total (31/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.09e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAA > NC_000913/2628489‑2628761 | cGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGAccc < 1:761591/140‑1 (MQ=255) gCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCt > 1:146395/1‑140 (MQ=255) cAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGAcc > 1:610925/1‑136 (MQ=255) cAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGAcc < 2:610925/136‑1 (MQ=255) cagtCCTAAAAATGCAGAAAACGGAAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCTGGAAAATTAAAC‑GGGTATCCCCATGTGAAA‑CGTTGTTCCGCGGGGGAGTAAAGCCCCTG‑TAAc < 2:349467/137‑1 (MQ=255) aaaTCCTAAAAATGCAGCAAACTGAAAAAGTATAAGAAACAGCGTTGCAGGAAACAAAATCTGCCAGAATATACGCGGGGAAATTAAAC‑GGGTATCACCATGTAAAA‑CGTTGTTACGCGGGGGAGTAAAGCCCCTG‑TAAc < 2:670228/139‑1 (MQ=255) tCCTAAAAATGCAGCAAACTGAAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCTGGGAAATAAAAC‑GGGTATCCCCATGTAAAC‑CGTTGTTCCGCGGGGGAGAAAAGCCCCTG‑TAACCGc < 2:701543/140‑1 (MQ=255) aaaaTGCAGAAAACGGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCAAGAATATCCGCGGGAAAATTAAAC‑GGTTATCCCCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGCCCCTG‑TAACCGCCCCaa < 2:430916/139‑1 (MQ=255) aaTGCAGCAAACCGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACa < 1:130524/140‑1 (MQ=255) aaTGCAGAAAACGGAAAAAGTATCAAAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATCCGCTGGAAAATTAACC‑GGGTACCCCCATGTGAAA‑CGTTGTTACGCCGGGGAGAAAAGCCCCTG‑TAACCGCCCCAACa < 2:1083773/140‑1 (MQ=255) cagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCCCCATGTAAAA‑CGTTGTTACGCGGGGGAGTAAAGCCCCTG‑TAACCGCCCCAACAgaga < 2:429480/140‑1 (MQ=255) cagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACagctg > 1:349467/1‑137 (MQ=255) cagcaAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGAGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAgaga > 1:779157/1‑140 (MQ=255) cagaaAACTGAAAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATCCGCGGGAAAATAAAAC‑GGGTATCACCATGAAAAA‑CTTTGTTACGCGGGGGAGTAAAGCCCCTG‑TAACCGCCCCAACAgaga < 2:779157/140‑1 (MQ=255) cagaaAACTGAAAAAGTATAAAAAACAGCGTTGCAGGAAACAAAATCTGCAAAAATATACGCGGGAAAATAAAAC‑GGGTATCACCATGTAAAA‑GGTTGTTACGCCGGGGAGTAAAGCCCCTG‑TAACCGCCCCAACAgaga < 2:402199/140‑1 (MQ=255) aaaGTACCAGAACCAGGTTTGCAGAAAACAAAGCCTCCAAAAATACGCGCTGGAAAATAAAGC‑GGTTACCA‑CAGGTAAAAGTGTTGTTGGGCGGGGGAG‑AAAAACCCGGATACCCGCCCCAACAGAGAatgataatgata < 2:1296163/140‑1 (MQ=255) tCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATaaaa > 1:1186404/1‑140 (MQ=255) aGAAACAGCGTTGAAGAAAACACATTCTGCCAAAATATGCGCGGGGAAATTAAAC‑GGGTATCCCCATGTAAAA‑CTTTGTTACGCGGGGGAAAAAAGCCCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAAtt < 2:1122679/140‑1 (MQ=255) aaaCAGGTTTGCAGAAAACAAAATCGCCCAAAATATACGCTGGAAAATAAAAC‑GGGTATCCCCTTGTAAAA‑CTTTGTTGCGCGGGGGAGAAAAGCCCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg < 2:1055680/140‑1 (MQ=255) aCAGCGTCGCAGGAAACAAAATCTCCAAGAATATACGCGGGGAATTAAAAC‑GGGTATCACCATGTAAAA‑CGTTGTTCCGCCGGCGAATAAAGCCCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGcc < 2:864940/140‑1 (MQ=255) aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt > 1:864940/1‑140 (MQ=255) aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt > 2:1090095/1‑140 (MQ=255) aGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGt > 1:1055680/1‑140 (MQ=255) cGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTa > 1:573563/1‑140 (MQ=255) ttGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAAt > 1:1076623/1‑140 (MQ=255) aaCACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGt > 1:190973/1‑140 (MQ=255) aTCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTa > 1:237118/1‑140 (MQ=255) tCTGCAAGAATATACG‑GGGGAAAATAAACGGGGTATCCCCATGTAAAA‑CGTTTTTACGCGGGGGAGTAAAGCCCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaa < 2:1010770/140‑1 (MQ=255) cTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaa > 1:109163/1‑140 (MQ=255) ccAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAc > 1:1149976/1‑140 (MQ=255) cAGAATATACGCGGGAAAATTAAAC‑GGGTATCACCATGTAAAA‑CTTTGTTACGCGGGGGAGAAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAcc < 1:546183/140‑1 (MQ=255) aaTATTCGCTGGAAAATAAAAC‑GGTTATCACCATGTAAAA‑CGTTGTTACGCCGGGGAAAAAAGCCCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg < 2:1155596/140‑1 (MQ=255) aaTATCCGCGGGAAAATAAACC‑GGGTATCCCCATGTAAAA‑TGTTGTTACGCGGGGGAAAAAAGCCCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg < 2:1149976/140‑1 (MQ=255) aaTATACGCGGGGAAATTAAAG‑GGGTATCACCTTGTAAAA‑CGTTGTTACGCGGGGGAGAAAAGCCCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg < 2:190973/140‑1 (MQ=255) aaTATACGCGGGAAAATTAAAC‑GGGTATCCCCATGTAAAA‑CGTTGTTACGCGGGGGAGAAAAGCCCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg < 2:146395/140‑1 (MQ=255) tataCGCGGGGAAATAAACC‑GGTTATCACCATGTAAAC‑CTTTGTTCCGCGGGGGAGAAAAGCCCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct < 2:237118/140‑1 (MQ=255) gggAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat > 2:1228489/1‑140 (MQ=255) ggAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt > 1:200717/1‑140 (MQ=255) ggAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAtt > 1:1063899/1‑140 (MQ=255) gAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAttt > 2:237364/1‑140 (MQ=255) gAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAttt > 2:1242082/1‑140 (MQ=255) gAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAttt > 2:488019/1‑140 (MQ=255) gAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAG‑GAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTg > 2:912515/1‑140 (MQ=255) c‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa > 1:783724/1‑140 (MQ=255) ggggTATCCCAATGTGAAC‑CTTTTTTCCGCGGGGGAGAAAAGCCCCTT‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaa < 2:1057520/139‑1 (MQ=255) aTGTAAAA‑CGTTTTTCCGCGGGGGAGAAAAGCCCCTT‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg < 2:1279400/140‑1 (MQ=255) tgtgAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGc > 1:1057520/1‑140 (MQ=255) gtgAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt > 2:1091521/1‑134 (MQ=255) gtgAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt < 1:1091521/134‑1 (MQ=255) gtgAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg > 1:947154/1‑140 (MQ=255) tgAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGa > 2:2265/1‑140 (MQ=255) tgAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGa > 1:463910/1‑140 (MQ=255) tgAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGa > 1:180949/1‑140 (MQ=255) tgAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGa > 1:1279400/1‑140 (MQ=255) tgAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGa > 2:921447/1‑140 (MQ=255) gggAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc > 2:1169930/1‑140 (MQ=255) gggAGAAAAGCCCTTG‑AAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGccc < 2:1186404/140‑1 (MQ=255) aaGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaa < 1:488019/140‑1 (MQ=255) aaGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGaa < 1:237364/140‑1 (MQ=255) | CGCCAGTAATCCTACAAATGCAGCAAACTGGAAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAAC‑GGGTATCACCATGTGAAA‑CGTTGTTACGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAATGCCGCCCAGCTGAA > NC_000913/2628489‑2628761 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |