Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,571,940 | 0 | A | C | 54.5% | 7.6 / 10.7 | 11 | V309G (GTG→GGG) | hsdM | type I restriction‑modification system methyltransferase |
Reads supporting (aligned to +/- strand): ref base A (4/1); new base C (0/6); total (4/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.52e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGACCGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAA > NZ_CP009273/4571856‑4572024 | aTGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCaccaccacc < 2:82752/90‑1 (MQ=255) tGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGAccg > 2:474977/1‑90 (MQ=255) gACACTTACCCATAAGGTCACGCCGAATGTCGGTGCCTTGGCCGCCTTCAAACAGCACGTTATCCGGCCCCCCCACCGCCGCACGAccgc < 1:534889/90‑1 (MQ=255) cacTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCCCCACCGCCGCACGAccgccg < 1:512190/90‑1 (MQ=255) aGGTCACGTCGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGACCGCCGGGATGCAGCGt > 1:401795/1‑89 (MQ=255) aGGTCACGTCGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGACCGCCGGGATGCAGCGt < 2:401795/89‑1 (MQ=255) cgacgaATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGACCGCCGGGATGCAGCGTTTCGata > 2:482895/1‑90 (MQ=255) ggTCCCTTTGCCGCCTTCAAACAGCACTTTACCCGGCCCCCCCCCCGCCGCACGACCGCCGGGATGCAGCGTTTCGATAATATGCTGCAt < 1:307206/90‑1 (MQ=255) ttGCCGCCTTCAAACAGCACTTTATCCGGCACCCCCACCGCCGCACGACCGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGcac < 1:326778/90‑1 (MQ=255) ccgctTTAAAACAGCACTTTACCCGGCCCCCCCACCGCCGCCCGACCGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAAc < 1:474977/90‑1 (MQ=255) ttCAAACAGCACGTTATCCGGCCCCCCCACCGCCGCACGACCGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACtgttt < 1:259392/90‑1 (MQ=255) gCACCACCACCGCCGCACGACCGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACaa > 2:396559/1‑90 (MQ=255) | ATGGTGTGCAGATGACACTTATCCATCAGGTCACGACGAATGTCGGTGCCTTTGCCGCCTTCAAACAGCACGTTATCCGGCACCACCACCGCCGCACGACCGCCGGGATGCAGCGTTTCGATAATATGCTGCATAAAGCACAACTGTTTGTTGCTGGTCGGGTGAACAA > NZ_CP009273/4571856‑4572024 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |