Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,001,918 | 0 | A | C | 31.4% | 49.6 / 11.3 | 35 | intergenic (‑129/+192) | fliA/fliC | RNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein |
Reads supporting (aligned to +/- strand): ref base A (22/2); new base C (0/11); total (22/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.87e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGC > NC_000913/2001783‑2002057 | aaTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑cccc > 2:921912/1‑140 (MQ=255) aaTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑cccc > 1:752254/1‑140 (MQ=255) tCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCccgcc > 1:757158/1‑140 (MQ=255) tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCAc > 1:498033/1‑140 (MQ=255) gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAA‑CCCGGCCGTGGGGGGGGAAGCACTTTGCTGaaaa > 2:688164/1‑140 (MQ=255) gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa > 1:994707/1‑140 (MQ=255) gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa > 1:664300/1‑140 (MQ=255) aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTTCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg > 2:577061/1‑140 (MQ=255) aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg > 1:297694/1‑140 (MQ=255) aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg > 1:581542/1‑140 (MQ=255) ggaTGTTTTTCCTTTTCTGGAGGTTTTTTATGCGGACGATTGGTGGGTAAAAGGGGGGTTTATTGGGGGGTAAAAGGAAAATTCCAGGAAA‑AAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt < 2:664300/138‑1 (MQ=255) aTTATCTTTTTCTGCAGGTTTTTTAATGGGACGATTAGTGGGTAAAAGGAGGGGTTATTTGGGGGTTACAGGTAAATCCCAGGCAA‑AAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATg < 2:556264/140‑1 (MQ=255) tCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg > 1:637502/1‑140 (MQ=255) tCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg > 2:211600/1‑140 (MQ=255) cTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc > 1:923640/1‑140 (MQ=255) gggttttttAATCGGACGATTAGGGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAA‑AAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg < 2:444494/136‑1 (MQ=255) ggaggtttAGGGGGTAAAATGGGGGTTTATTTGGGGTTAAAAGGTAAATCCAAGAAAA‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc < 2:757158/134‑1 (MQ=37) ggACGTTTAGGGGGTGAAATGGGGGGTTTTTGGGGGGTTACAGGTAAATTCCAGGAAA‑AAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc < 2:637502/140‑1 (MQ=38) gggttttAGGGGGGAAAATGGGGGTTTGTTTGGGGGTTAAAGGTAAATTCCGGGAAA‑AAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc < 1:688164/135‑1 (MQ=37) gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc > 2:149893/1‑140 (MQ=255) gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc > 1:162328/1‑140 (MQ=255) gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc > 2:457/1‑140 (MQ=255) aCGATTAGTGGGTGAAATGGGGGGTTATTGGGGGTTTACAGGTAAATCCCAGGAAA‑AAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg < 1:446215/140‑1 (MQ=38) cGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg > 1:574560/1‑140 (MQ=255) gggggTGAAATGGGGGTTTTTTGGGGGGTTACAGTAAAATCCAAGGAAA‑AAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa < 2:752254/138‑1 (MQ=25) ggTAAAAGGGGGGGTTATTGGGGGGTTACAGGAAAATTCCAGGCAA‑AAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa < 1:703965/140‑1 (MQ=37) gtaaaaTGAGGGGTTTTTTGGGGGTTACAGGTAAATCCCAGGAAA‑AAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa < 2:923640/137‑1 (MQ=38) gAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg > 2:359240/1‑140 (MQ=255) tttggggggTTACAGGTAAATTCCAGGCAA‑AAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACAAAAATGTGCAAATTCAATa < 1:921912/136‑1 (MQ=38) aTTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa > 2:11875/1‑140 (MQ=255) tGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGTAAATTCAATAAAt < 1:149893/140‑1 (MQ=255) gTTACAGGAAAATTTCAGGAAA‑AAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCaa < 1:577061/140‑1 (MQ=37) cAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa > 2:1001171/1‑140 (MQ=255) ggCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg > 2:235062/1‑140 (MQ=255) ggCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg > 2:911801/1‑140 (MQ=255) a‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagc > 2:289144/1‑140 (MQ=255) | AATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAG‑AAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGC > NC_000913/2001783‑2002057 |
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: ‑ |