Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,001,918 | 0 | A | C | 33.3% | 40.0 / 15.5 | 33 | intergenic (‑129/+192) | fliA/fliC | RNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein |
Reads supporting (aligned to +/- strand): ref base A (20/2); new base C (0/11); total (20/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.03e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGC > NC_000913/2001783‑2002057 | aaTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑cccc > 2:96707/1‑140 (MQ=255) aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCccg > 1:756619/1‑140 (MQ=255) aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCccg > 1:608503/1‑140 (MQ=255) tCACGATAAACAGCCCTGCGTTATATGAGTTATGGGCATGATTATCCTTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGGTTATTG‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAC‑CCccgcc < 2:690731/140‑1 (MQ=255) aTAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTAGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTgg > 2:101924/1‑140 (MQ=255) aaatatGAGTTATCGGCAGAATTATCCTTTTTGGCGGGTTTTTTAATCGGACGTTTAGTGGGAAAAAGGAGGGGTTATTG‑GGGGGTAAAAGGAAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTg < 2:229801/138‑1 (MQ=255) gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa > 1:619104/1‑140 (MQ=255) aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg > 2:101070/1‑140 (MQ=255) aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg > 1:438830/1‑140 (MQ=255) cGGCATGATTATCCGCTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgct > 1:229801/1‑139 (MQ=255) aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt > 2:431424/1‑140 (MQ=255) cTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc > 1:745562/1‑140 (MQ=255) tttttaAATGGGACGATTAGGGGGTGAAATGAGGGTTTATTT‑GGGGGTTACAGGAAAATTCCAGAAAAAAAAACA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa < 2:727936/139‑1 (MQ=255) tttAACGGGACGATTAGGGGGTGAAATGGGGGTTTATTG‑GGGGGTTACAGGAAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg < 2:438830/140‑1 (MQ=255) gACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc > 1:475220/1‑140 (MQ=255) ggTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa > 2:8760/1‑140 (MQ=255) ggTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa > 2:34017/1‑140 (MQ=255) tGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt > 2:17357/1‑140 (MQ=255) gAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg > 2:529672/1‑140 (MQ=255) aTGGGGGGTTATTG‑GGGGGTTACAGGTAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc < 2:475220/140‑1 (MQ=38) gggTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGTAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt > 1:161290/1‑140 (MQ=255) ggTTTTTTGGGGGGTTACAGATAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt < 2:745562/140‑1 (MQ=37) ggggTTACAGGTAAATTCCAGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg < 1:34017/140‑1 (MQ=255) ttagagTAAAATCCAGAAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAt < 2:161290/135‑1 (MQ=25) aCAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc > 2:264251/1‑140 (MQ=255) aCAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc > 1:285140/1‑140 (MQ=255) aCAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc > 1:273249/1‑140 (MQ=255) aCAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc > 1:201557/1‑140 (MQ=255) cAGGTAAATTCCAGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa < 1:264251/140‑1 (MQ=255) aaGGAAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa < 2:619104/139‑1 (MQ=37) aaaTCCCAGGCAGAAAAACC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg < 1:96707/140‑1 (MQ=255) tcccAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg < 2:756619/138‑1 (MQ=255) ccAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc < 1:610663/140‑1 (MQ=255) a‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagc > 2:330382/1‑140 (MQ=255) | AATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTT‑GGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGC > NC_000913/2001783‑2002057 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |