Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,012,931 | 0 | A | C | 54.5% | 2.2 / 11.1 | 11 | T138P (ACA→CCA) | ssnA | putative aminohydrolase SsnA |
Reads supporting (aligned to +/- strand): ref base A (0/5); new base C (6/0); total (6/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.16e-03 | |||||||||||
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. |
TTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCT‑CCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAAC > NZ_CP009273/3012846‑3013002 | ttGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGTGTTCT‑Ccccacc > 1:263452/1‑88 (MQ=255) cTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCT‑CCCCCTTTcgcggc > 2:70315/1‑88 (MQ=255) gAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCT‑CCACATTGCGCGACG‑Ca < 1:70315/90‑1 (MQ=255) ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGGGCCCT‑CCCCCTTTCGCGGCGCCTTTTTTTAAATTTGGCc‑gcgcggggt > 1:302258/1‑86 (MQ=255) ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCCCT‑CCCCCTTGCGCGGCG‑CATTTTTTTAAATTGGCCTGcgcgcggt > 1:275013/1‑88 (MQ=255) aTCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGGTCTCCCCCATTCCCGGCCC‑CCTTTTTAAAAATTGGCC‑GCGCGCGGTGGcc > 1:121364/1‑90 (MQ=255) ctctCCGGCGTATATCGGCGGGTCGCTCT‑CCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAactac < 1:156850/90‑1 (MQ=255) ctcCGGCGTATATCGGCGGGTCGCTCT‑CCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTg < 2:302258/90‑1 (MQ=255) gCGTATATCGGCGGGGCGCCCT‑CCCCATTGCGCGGCG‑CATTTTTTAAAAATTGCCTGCGCCCGGTGACCCGTTTTTTAAAAACTGAccct > 2:232892/1‑88 (MQ=255) atatCGGCGGGTCGCTCT‑CCCCATTGCGCGACG‑CATTTTTTAAAATTGGCCTGCGCGCGGTGACCCGTTTTTTAAATACTGACCGTaaaa > 2:124097/1‑88 (MQ=255) atCGGCGGGTCGCTCT‑CCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTaacaac < 1:124097/90‑1 (MQ=255) atCGGCGGGTCGCTCT‑CCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTaacaac < 1:222324/90‑1 (MQ=255) | TTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCT‑CCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACCGTAACAAC > NZ_CP009273/3012846‑3013002 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |