Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,387,038 | 0 | A | T | 55.6% | 1.3 / 18.8 | 18 | S860R (AGT→AGA) | yhdP | AsmA2 domain‑containing protein YhdP |
Reads supporting (aligned to +/- strand): ref base A (1/7); new base T (10/0); total (11/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.51e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.11e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCAC > NZ_CP009273/3386956‑3387117 | aCAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGGGGTTTtctgtctg > 1:202688/1‑90 (MQ=255) aGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGGGGTTTTCTTTCTgcc > 1:334084/1‑90 (MQ=255) gggCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGGGGTTTTCTTTCTgccc > 1:211725/1‑89 (MQ=255) gggCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGCGGTTTTCTTTCTtccc > 1:698890/1‑89 (MQ=255) cagccaCTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGCGGTTTTCTTTCTGccccccc > 1:51474/1‑89 (MQ=255) gccaCTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGGGGGTTTCTGTCTGcccccccca > 2:564449/1‑87 (MQ=255) aTTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGGGTTTTTCTTTCTGCCGCCCCAATAACACGATCGAg > 1:562517/1‑90 (MQ=255) aTCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTc < 1:184131/90‑1 (MQ=255) gAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGGGGGTTTTTCTTTCTGCCGCCCCAAAAGCACGATCGAGCGGCAGCTTTTTACCCAg > 1:659373/1‑90 (MQ=255) aGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGc < 2:541117/90‑1 (MQ=255) ttCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGGTTTTCTGTCTTCCGCCCCAAAAACACCATCGAGCGGCAGCTTTTTTCCGAGCaa > 2:106858/1‑90 (MQ=255) cACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGc < 1:42893/90‑1 (MQ=255) gTTCCGGCAACGGCGGGGGCGGTTTTCTGTTTTCCCCCCCAAAAACACGATCGAGCGGCCGCCTTTTTCCGCGCAAACACCGGCCATTGa > 2:397284/1‑90 (MQ=255) ttCCGGCAACTGCGGGGGGGGTTTTCTGTTTGCCGCCCCAAAAGCACCAACCGGCGGCAGCGTTTTTCCGAGCAACCAGCGGCTATTTaa > 2:617710/1‑90 (MQ=255) gCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGAt < 2:82149/90‑1 (MQ=255) aaCGGCGGGAGCGTTTTACTGTCTGCCGCCCAAAAAGCACGATCGAGCGTCAGCTTTTTACCGAGCAACCAGCGGCTATTGAAATGATTa > 1:39617/1‑90 (MQ=255) cggcggGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATc < 1:307345/90‑1 (MQ=255) cggcggGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATc < 1:21030/90‑1 (MQ=255) ggAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCAc < 1:397284/90‑1 (MQ=255) | ACAGGGCCAGCCACTCGGCACCATTCATCGGCGGCATATTGAGTTCAACGCCACTTTGTTCCGGCAACGGCGGGAGCGTTTTACTGTCTGCCGCCCAAATAGCACGATCGAGCGTCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCAC > NZ_CP009273/3386956‑3387117 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |