Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 341,039 | 0 | C | A | 29.8% | 94.6 / 17.3 | 47 | intergenic (+337/‑76) | BW25113_RS01680/BW25113_RS01685 | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (29/4); new base A (0/14); total (29/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.96e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGT‑GTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCA > NZ_CP009273/340919‑341173 | aaTGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGT‑GTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCt > 2:217950/1‑138 (MQ=255) ttggttgACATTCATATAAAAAAAATCATAATCCCATCATTTTTGT‑GTGGGTTTTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTGACAAAAAATAGGCTGAAA‑GGCTAATATCAAAAACCGCCATTAGCTCATCAgg < 1:187687/139‑1 (MQ=255) tggttgACATTCATATGAAAAAAATCATAATTCCATCATGTTTGT‑GTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAgg > 1:219477/1‑138 (MQ=255) gttgACATTCATATAAAAAAAACAAAAATCCCATCATGTTTGT‑GGGGGTTTTTTTC‑TGTACCTTACGCACCGCACCCACCCCAGACAGAAAATATGCTGAAG‑GCCTCATATCAAAAACCGCCATTAGCTCATCAGGaa < 2:150720/138‑1 (MQ=255) gACATTCATATGAAAAAAATCATAATTACATCATGTTTGT‑GTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAagag > 2:160037/1‑138 (MQ=255) aCATTCATATAAAAAAAAACATAATTCCATCATGTTTGT‑GTGGGGTCTTTTC‑TGTATCTTCCGCACCCCACACAAGCCTGACAGAAAATATGCGGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGc < 1:141103/138‑1 (MQ=255) cATTCATATGAAAAAAATCATAATTCCATCATGTTTGT‑GTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAagag > 2:54144/1‑136 (MQ=255) catagaaaaaaaaTCAAAATCCCATC‑TGTTTGTGGTGGGTTTTTTTC‑TTTATCTTCACCAACCCCCCAAAGCCTGACAAAAAATAGGCTGAAA‑GGCAAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacg < 2:219477/132‑1 (MQ=255) tatGAAAAAAATCATAATTCCATCATGTTTGT‑GTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacgac > 1:255524/1‑138 (MQ=255) aaaaaagCCATAATCCCAGTAGGTTTGG‑GTGGGG‑CTGTTCTTGTCTGTTCCACACCGCACTCAAGCTTGCCAAAAAATAGGC‑GTAAGGGCCCATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCa < 1:145438/139‑1 (MQ=255) aaTCCCATCATTTTTGT‑GTGGGTTTTTTTT‑TTTATCTACCCCACCCCACCCACGCCTGAAAAAAAAGGTGTTGAAA‑GCCACATAAAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTa < 2:270751/139‑1 (MQ=255) aTTCCATTATTTTTGT‑GTGGGTTTTTTTT‑TTTACATTACGCCTCGCACTCAAGCCGGACAAAAAATTTGCTTAAG‑GGAGTATATAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTa < 1:160037/138‑1 (MQ=255) atcatGTTTGT‑GTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGATCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTAc > 1:231255/1‑139 (MQ=255) tcatGTTTGT‑GTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACg > 1:168684/1‑139 (MQ=255) tcatGTTTGT‑GTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCCCATCAGGAAGAGCAGACGACTACCAGGATTGTTGTATggcagt > 1:119113/1‑135 (MQ=255) tggggggTTTTTTC‑TGTATCTTACGCACCGCCCTCAACCCTGAAAGAAAATAGGCTGAAA‑GGCTCATAAAAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGa < 1:217950/135‑1 (MQ=255) aTCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTcc > 2:206700/1‑138 (MQ=255) tACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATCTCAAAAACCGCCACTAGCTCATCAGGAAGAGCAGACGACACCCAGGATTCTTGCATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtt > 2:112843/1‑139 (MQ=255) tACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtt > 1:252757/1‑139 (MQ=255) tACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtt > 1:187864/1‑139 (MQ=255) tACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtt > 1:36668/1‑139 (MQ=255) tACGCATCGCACTCAAGCCTGAAAGAAAATATGCTGTAA‑GGCTCATCTCACAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACta > 2:279461/1‑138 (MQ=255) cGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACtttt > 1:206102/1‑139 (MQ=255) aTCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCa > 1:9123/1‑138 (MQ=255) cGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaaa > 2:161714/1‑139 (MQ=255) cGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaaa > 1:65984/1‑139 (MQ=255) cGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaaa > 2:103675/1‑139 (MQ=255) cGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaaa > 1:59475/1‑139 (MQ=255) cGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATAGCGGTCCAACTTTTCaaaa > 2:50487/1‑139 (MQ=255) cGCACTCAAGCCTGACAGAAAATATGATGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCaaaa > 2:93309/1‑139 (MQ=255) cACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAAct > 1:129650/1‑139 (MQ=255) ccttaaCAGAAAAGTGGCTGTAG‑GGCACAAACAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGa < 2:45942/134‑1 (MQ=255) ccTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGaa > 1:93654/1‑139 (MQ=255) aCAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAAt > 1:228261/1‑136 (MQ=255) aCAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAAt > 1:186389/1‑136 (MQ=255) caaaaaaTAGGCGGTAA‑GGCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGa < 1:49741/136‑1 (MQ=255) cAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAAt > 2:162346/1‑135 (MQ=255) aGAAAATATGCTGTAA‑GGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGAGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGAt < 1:218430/139‑1 (MQ=255) atgatGCTGTAA‑GTCTAATATAAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcccc < 1:112643/136‑1 (MQ=255) aatgaacTGTAG‑GCCCCATACAAAAAGCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGAGcccc < 2:254141/133‑1 (MQ=255) aagaagCTGTAC‑GCCTCAAACAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcccc < 2:49511/134‑1 (MQ=255) aaTATGCTGTAA‑GCCTCATAAAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTTTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcccc < 2:59475/139‑1 (MQ=255) aaTAAGCAGTAA‑GCCACAAACAAAAATCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATcccc < 2:65984/139‑1 (MQ=255) cTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCt > 2:221193/1‑139 (MQ=255) cTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCt > 1:186119/1‑139 (MQ=255) atatAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCa < 2:252757/139‑1 (MQ=255) ataaCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCa < 1:4635/139‑1 (MQ=255) | AATGTGTTGGTTGACATTCATATGAAAAAAATCATAATTCCATCATGTTTGT‑GTGGGGTCTTTTC‑TGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAA‑GGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCAACTTTTCAAAACTCTTCTGAATCAGATCCCCTTATTATCATCCCTGCA > NZ_CP009273/340919‑341173 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |