Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,379,045 | 0 | T | G | 55.6% | 5.4 / 14.5 | 18 | G224G (GGT→GGG) | ycjX | YcjX family protein |
Reads supporting (aligned to +/- strand): ref base T (4/4); new base G (10/0); total (14/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.29e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.08e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCCACCACTGTAAAGAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAGGCCGATAAGCATACCAATGCCGG > NZ_CP009273/1378971‑1379134 | ctcCACCACTGTAAAGAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCGCGCTGCaa > 2:151122/1‑90 (MQ=255) acTGTAAAGAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCCCGCCGCAAttttttc > 2:424713/1‑88 (MQ=255) gagCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCg > 2:280256/1‑90 (MQ=255) gagCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGGGCCCGCGCTGCAATTCTTCCCCTGGCCg > 2:585409/1‑90 (MQ=255) gCAGGGGCTGCACTTTATTCCGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGa < 2:339632/90‑1 (MQ=255) gCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCCCGCCGCAATTTTTCCCCTGGCCggg > 2:178154/1‑89 (MQ=255) ggggCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGt > 1:559185/1‑90 (MQ=255) ggggCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGt > 1:114305/1‑90 (MQ=255) ggggCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCCCGCTGCAATTTTTTCCCTGGCCGGgtgt > 2:175449/1‑90 (MQ=255) ggggCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCCCGCTGCAATTCTTCCCGTGGCCGGgtgt > 2:337920/1‑90 (MQ=255) ggCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGGGCCCGCGCTGCAATTCTTCCCGTGGGCGGGGGTCg > 2:358607/1‑90 (MQ=255) gCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGa < 2:559185/90‑1 (MQ=255) ttATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGGGCCCGCGCCGCAATTCTTTCCCTGGGCGGTTGGTGACACCTtgg > 1:397422/1‑90 (MQ=255) gCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGTGCCCGCGCTGGAATTCTTCCCGTGGCCGGCTGTCTAAACCTGGGGCGGGTc > 1:687820/1‑90 (MQ=255) tGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGGGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCaa > 1:342236/1‑90 (MQ=255) gagaTATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAGGCCGATAAGc < 1:333267/90‑1 (MQ=255) gagaTATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAGGCCGAAAAGc > 1:343413/1‑90 (MQ=255) tGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAGGCCGATAAGCATACCa < 1:424713/90‑1 (MQ=255) tGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGGTGTCGATACCTGGGGGGAGTCCAAACTTGCGCAGGCCGATAAGCATACCAATGCCgg > 1:141272/1‑90 (MQ=255) | CTCCACCACTGTAAAGAGCAGGGGCTGCACTTTATTCAGCCTGGGCGCTTTGTCTTGCCGGGAGATATGGCAGGTGCGCCCGCGCTGCAATTCTTCCCGTGGCCGGATGTCGATACCTGGGGCGAGTCCAAACTGGCGCAGGCCGATAAGCATACCAATGCCGG > NZ_CP009273/1378971‑1379134 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |