Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,012,922 | 0 | T | G | 59.1% | 14.4 / 25.4 | 22 | S135A (TCG→GCG) | ssnA | putative aminohydrolase SsnA |
Reads supporting (aligned to +/- strand): ref base T (1/8); new base G (13/0); total (14/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.81e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACC > NZ_CP009273/3012833‑3012994 | aCAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt < 1:19228/90‑1 (MQ=255) aCAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt < 2:486361/90‑1 (MQ=255) aGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCg < 1:474721/90‑1 (MQ=255) gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt < 1:150958/87‑1 (MQ=255) gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGt > 2:150958/1‑87 (MQ=255) gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGgcgc > 2:552010/1‑90 (MQ=255) gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCg < 2:329975/89‑1 (MQ=255) gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCgg > 1:219747/1‑89 (MQ=255) gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCgg > 1:268491/1‑89 (MQ=255) gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCgg > 1:649381/1‑89 (MQ=255) gCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCg > 1:329975/1‑89 (MQ=255) ggACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTccgtg > 1:502085/1‑86 (MQ=255) ggACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGgtc > 1:403583/1‑90 (MQ=255) ggACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGctc > 2:83732/1‑90 (MQ=255) ggACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGctc > 1:94208/1‑90 (MQ=255) aCTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGCCGTATATCGGCGGGGCGCtctc > 2:652469/1‑90 (MQ=255) gagCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACG‑CATTTTTAAAAGt < 1:637838/90‑1 (MQ=255) ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGGGCGCTCCCCCCTTCGCGGCG‑CCCTTTTTTAAAATGGGCTGggcgcggg > 1:387309/1‑88 (MQ=255) ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCCCATTTCGCGACG‑CCTTTTTTAAAATTTGGCTGcgcgcggg > 2:157527/1‑88 (MQ=255) ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCCCAATTCGCGGCGCCTTTTTTAAAATTTGGCc‑gcgcgcggg > 2:404397/1‑88 (MQ=255) ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCCCTCCCCATTTCGCGGCC‑CCTTTTTTAAAATTTGGCTGcgcgcggg > 2:162835/1‑88 (MQ=255) aTCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGAcc < 2:543817/90‑1 (MQ=255) cGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCCCTCCCCCTTTCGCGCCG‑CATTTTTTAAAATTGGGCTGCGCGCGCGGAcccg > 2:440814/1‑88 (MQ=255) cACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGtttt < 2:387591/90‑1 (MQ=255) cATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGaa < 2:231645/90‑1 (MQ=255) ctctCCGGCGTATATCGGCGGGTCGCTCCCCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAactac < 2:104998/90‑1 (MQ=255) cGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGAcc < 1:399765/90‑1 (MQ=255) | ACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCACATTGCGCGACG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGACC > NZ_CP009273/3012833‑3012994 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |