Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 341,039 | 0 | C | A | 34.2% | 88.2 / 22.7 | 41 | intergenic (+337/‑76) | BW25113_RS01680/BW25113_RS01685 | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (26/1); new base A (0/14); total (26/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.26e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCA > NZ_CP009273/340959‑341125 | catcatGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGc > 2:397550/1‑90 (MQ=255) atcatGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGcc > 2:430329/1‑90 (MQ=255) tttGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGc > 2:169120/1‑90 (MQ=255) ggTTTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGa < 2:26047/90‑1 (MQ=255) ggTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGa > 1:219579/1‑90 (MQ=255) ttttCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAagag > 2:93182/1‑90 (MQ=255) tATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAgacga > 2:131334/1‑90 (MQ=255) ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc > 1:303496/1‑90 (MQ=255) ttACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAc > 2:414608/1‑90 (MQ=255) tACGCATCGCACTCAAGCCTGACAGAAAATACGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc > 1:462846/1‑90 (MQ=255) aCGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCGGACGACAACCa > 2:52839/1‑90 (MQ=255) aCGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCa > 2:44011/1‑90 (MQ=255) tCGCCCTCAAGCCTGAAAAAAAATATGCTGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt < 2:150698/90‑1 (MQ=255) tCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt > 1:27278/1‑90 (MQ=255) tCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt > 1:381523/1‑90 (MQ=255) tCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt > 1:150698/1‑90 (MQ=255) tCGCACTAAAGCCTGAAAAAAATTTTGCTGTAAGGTTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAtt < 2:333291/90‑1 (MQ=255) gCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt > 1:333291/1‑90 (MQ=255) gCACTCAAGCCTGAAAAAAAATATGTTGAAAGCCTAATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt < 2:381523/90‑1 (MQ=255) ccACTCAAGCCAAAAAAAAAATATGCTGAAAGGCTAAAAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGAttgt < 1:430329/89‑1 (MQ=255) cAAGCCTGAAAGAAAATACGTTGAAAGGCTCAAAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt < 2:462846/90‑1 (MQ=255) cAAGCCTGAAAAAAAATATGCTGAAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt < 2:135592/90‑1 (MQ=255) cAAGCCTAAAAAAAAATATGTTGAAAGGCTAATAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt < 1:163663/90‑1 (MQ=255) aaacccTGACAAAAAATATGCTTTAAGGCCCATATAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTAt < 2:27278/86‑1 (MQ=255) aaGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATg > 2:270567/1‑90 (MQ=255) ctaaCAAAAAATATGCTGAAAGCCTAATAAAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTAc < 1:90770/87‑1 (MQ=255) tGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACg > 2:64916/1‑90 (MQ=255) atatGCTGAAAGCCTCAAAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAg < 1:169120/90‑1 (MQ=255) tatGCTGAAAGCCTCATAAAAAAACCCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAgg < 1:131391/90‑1 (MQ=255) gCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCt > 1:131472/1‑90 (MQ=255) gCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCt > 1:131750/1‑90 (MQ=255) tGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTa > 2:218269/1‑68 (MQ=255) tGTAAGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTa < 1:218269/68‑1 (MQ=255) gaaaGGCTCATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc < 1:408515/53‑1 (MQ=255) gTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAAcc > 2:408515/1‑55 (MQ=255) cATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCg > 1:99374/1‑73 (MQ=255) cATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTc > 1:382706/1‑90 (MQ=255) cATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTc > 2:63309/1‑90 (MQ=255) cATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGAACGGGGTTCGCGGCCTCGCTGGCGGTc > 1:142667/1‑90 (MQ=255) cATATAAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCg < 2:99374/73‑1 (MQ=255) tatCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCa > 1:75441/1‑90 (MQ=255) | CATCATGTTTGTGTGGGGTCTTTTCTGTATCTTACGCATCGCACTCAAGCCTGACAGAAAATATGCTGTAAGGCTCATATCAAAAACCGCCATTAGCTCATCAGGAAGAGCAGACGACAACCAGGATTGTTGTATGGTACGGGGTTCGAGGCCTCGATGGCGGTCCA > NZ_CP009273/340959‑341125 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |