Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,012,922 | 0 | T | G | 52.9% | 6.5 / 13.5 | 17 | S135A (TCG→GCG) | ssnA | putative aminohydrolase SsnA |
Reads supporting (aligned to +/- strand): ref base T (3/5); new base G (9/0); total (12/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.05e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAA > NZ_CP009273/3012833‑3012998 | aCAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt < 1:593610/90‑1 (MQ=255) aGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCg < 1:208420/90‑1 (MQ=255) gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt > 2:58352/1‑87 (MQ=255) gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt < 1:58352/87‑1 (MQ=255) gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGc > 2:467438/1‑90 (MQ=255) gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCgg > 2:690561/1‑89 (MQ=255) ggACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGgcgccc > 2:532289/1‑88 (MQ=255) ggACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGctc > 2:130691/1‑90 (MQ=255) aCTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCtctc < 1:235251/90‑1 (MQ=255) cTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGGGCCCTCCCCACTTcgcggc > 1:184095/1‑88 (MQ=255) cTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCCCCCTGCCCGAc > 1:699173/1‑90 (MQ=255) tCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACG‑CCTTTTTTAAAATTGGGCTGcgcgcg > 2:493134/1‑90 (MQ=255) ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGGGCCCTCCCCACTTCGCGGCGCCCTTTTTAAAATTTGG‑CTgcgcggggt > 2:663371/1‑86 (MQ=255) ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCCCTCCCCACTGCCCGACG‑CCCTTTTTAAAAATGGCCTGcgcgcggt > 1:277936/1‑88 (MQ=255) aTCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGAcc < 1:586677/90‑1 (MQ=255) ccGGCGTATATCGGCGGGGGGCTCTCCCCCTTTCGCGGCC‑CCTTTTTTTAAATTTGCCTGCGCGCGGTGGCCCGTTTTTTAACTACTTAc > 1:248501/1‑90 (MQ=255) ccGGCGTATATCGGCGGGGCGCTCTCCCCACTGCGCGGCGCCTTTTTTAAAAATTGGCC‑GCGCGCGGGGGCCCGTTTTTTAAATACTTAc > 2:408536/1‑90 (MQ=255) cGGCGTATATCGGCGGGTCGCTCTCCCCATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGAcc < 1:444875/90‑1 (MQ=255) cGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGAcc < 2:69908/90‑1 (MQ=255) gCGTATATCGGCGGGGGGCTCTCCCCCTTTCGGGGCGCCTTTTTTAAAAATTGG‑CTGCGCGGGGTGACCCGTTTTGTAAATAATGACCgg > 2:274825/1‑89 (MQ=255) gTATATCGGCGGGTCGCCCTCCCCCTTTCGCCGCG‑CCTTTTTTAAAAATGGCCTGCGCGCGGTGACCCGTTTTTTAAATACTGACCCTaa > 2:426807/1‑90 (MQ=255) | ACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAA > NZ_CP009273/3012833‑3012998 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |