Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,162,884 | 0 | T | G | 72.7% | 46.9 / 12.8 | 22 | I231R (ATA→AGA) | yegS | lipid kinase YegS |
Reads supporting (aligned to +/- strand): ref base T (2/4); new base G (7/9); total (9/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTCAGGCCGGTGGCGGTCAGCAATTGTGTCCGAACGCGTTAATTAACGATGGCTTGCTGCAACTGCGCATTTTTACCGGCGATGAAATACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGAATATTATCGAAGGCGCTTCGTCGTGGTTTGATATTCAGGCAC > NZ_CP009273/2162797‑2162973 | gTCAGGCCGGTGGCGGTCAGCAATTGTGTCCGAACGCGTTAATTAACGATGGCTTGCTGCAACTGCGCATTTTTACCGGCGATGAAAGAc > 1:817623/1‑90 (MQ=255) ggTCAGCAATTGTGTCCGAACGCGTTAATTAACGATGGCTTGCTGCAACTGCGCATTTTTACCGGCGATGAAAGACTTCCGGCTCTCGTa > 2:26876/1‑90 (MQ=255) tCAGCAATTGTGTCCGAACGCGTTAATTAACGATGGCTTGCTGCAACTGCGCATTTTTACCGGCGATGAAATACTTCCGGCTCTCGTATc < 2:226646/90‑1 (MQ=255) aGCAATTGTGTCCGAACGCGTTAATTAACGATGGCTTGCTGCAACTGCGCATTTTTACCGGCGATGAAAGACTTCCGGCTCTCGTATCaa > 1:708330/1‑90 (MQ=255) gTTAATTAACGATGGCTTGCTGCAACTGCGCATTTTTACCGGCGATGAAATACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGa > 1:95266/1‑90 (MQ=255) taattaaCGATGGCTTGCTGCAACTGCGCATTTTTACCGGCGATGAAAGACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATa < 1:867143/90‑1 (MQ=255) aaCGATGGCTTGCTGCAACTGCGCATTTTTACCGGCGATGAAAGACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCg < 2:217268/90‑1 (MQ=255) cGATGGCTTGCTGCAACTGCGCATTTTTACCGGCGATGAAATACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGaa < 2:926628/90‑1 (MQ=255) cGATGGCTTGCTGCAACTGCGCATTTTTACCGGCGATGAAAGACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGaa < 2:189777/90‑1 (MQ=255) cGATGGCTTGCTGCAACTGCGCATTTTTACCGGCGATGAAAGACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGaa > 2:722916/1‑90 (MQ=255) gATGGCTTGCTGCAACTGCGCATTTTTACCGGCGATGAAAGACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGAat < 2:870786/90‑1 (MQ=255) ctgcAACTGCGCATTTCTACCGGCGATGAAAGACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGAATATTATCGaa < 2:1039478/90‑1 (MQ=255) aCTGCGCATTTTTACCGGCGATGAAAGACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGAATATTATCGAAGgcgc < 1:722916/90‑1 (MQ=255) tGCGCATTTTTACCGGCGATGAAATACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATATCCCGAATATTATCGAAGGCGCtt < 2:779354/90‑1 (MQ=255) gcgcATTTTTACCGGCGATGAAATACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGAATATTATCGAAGGCGCTtc > 2:431358/1‑90 (MQ=255) gcgcATTTTTACCGGCGATGAAAGACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGAATATTATCGAAGGCGCTtc > 2:531235/1‑90 (MQ=255) gcgcATTTTTACCGGCGATGAAAGACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGAATATTATCGAAGGCGCTtc > 1:869866/1‑90 (MQ=255) cATTTTTACCGGCGATGAAAGACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGAATATTATCGAAGGCGCTtcgtc < 2:708330/90‑1 (MQ=255) ttACCGGCGATGAAAGACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGAATATTATCGAAGGCGCTTCGTCGTGGt > 2:267665/1‑90 (MQ=255) ggCGATGAAAGACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGAATATTATCGAAGGCGCTTCGTCGTGGTTTGat < 2:417975/90‑1 (MQ=255) ggCGATGAAAGACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGAATATTATCGAAGGCGCTTCGTCGTGGTTTGat < 2:868541/90‑1 (MQ=255) aaTACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGAATATTATCGAAGGCGCTTCGTCGTGGTTTGATATTCAGGc < 2:603307/90‑1 (MQ=255) tACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGAATATTATCGAAGGCGCTTCGTCGTGGTTTGATATTCAGGcac < 2:1090921/90‑1 (MQ=255) | GTCAGGCCGGTGGCGGTCAGCAATTGTGTCCGAACGCGTTAATTAACGATGGCTTGCTGCAACTGCGCATTTTTACCGGCGATGAAATACTTCCGGCTCTCGTATCAACCTTAAAATCTGACGAAGATAACCCGAATATTATCGAAGGCGCTTCGTCGTGGTTTGATATTCAGGCAC > NZ_CP009273/2162797‑2162973 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |