Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 69,900 | 0 | A | C | 58.3% | 7.1 / 10.4 | 12 | V370G (GTG→GGG) | thiP | thiamine/thiamine pyrophosphate ABC transporter permease ThiP |
Reads supporting (aligned to +/- strand): ref base A (4/1); new base C (0/7); total (4/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.01e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCA > NZ_CP009273/69831‑69976 | aGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTc > 1:30798/1‑90 (MQ=255) gCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCg > 1:269247/1‑90 (MQ=255) ccGGTACCCAGCACAACCCCCGGCATGGCGAGGATCAACATGCCGCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGc < 1:41809/90‑1 (MQ=255) cccGGCATGGCGAGGATCAACATGCCGCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCAt < 1:59401/90‑1 (MQ=255) cGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTcgcg < 1:138504/78‑1 (MQ=255) cGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTcgcg > 2:138504/1‑78 (MQ=255) ggCAGGGGGGGGACCAACATGCCCCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATagc < 2:269247/90‑1 (MQ=255) gggggATCAACATGCCCCTCCTCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcatca < 1:144758/87‑1 (MQ=255) cGCTCATCTCCAGCCCCTCCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc < 2:194879/90‑1 (MQ=255) cGCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc < 1:3764/90‑1 (MQ=255) cGCTCACCCCCACCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGc < 2:30798/90‑1 (MQ=255) gCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGca > 1:109431/1‑90 (MQ=255) | AGTAAAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACTACGCA > NZ_CP009273/69831‑69976 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |