| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 69,905 | 0 | A | C | 52.9% | 11.5 / 16.5 | 17 | G368G (GGT→GGG) | thiP | thiamine/thiamine pyrophosphate ABC transporter permease ThiP |
| Reads supporting (aligned to +/- strand): ref base A (4/4); new base C (0/9); total (4/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.94e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-02 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAG‑‑CACTACGCACAATACAC > NZ_CP009273/69835‑69984 | aaaaGAGGCGGGTGGCCAGCACAACCCCCGGCAGGGGGAGGATCAACATGCCCCTCATCCCCAGCCCCTGACCCGCCAGCATTTTCTGcc < 2:119325/90‑1 (MQ=255) gCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCg > 2:111505/1‑90 (MQ=255) ccGGTAGCCAGCACAACCCCCGGCAGGGCGAGGATCAACATGCCGCTCATCTCCAGCCCCTGACCCGCCAGCATTTTCTGCCGCGCCCGc < 1:35846/90‑1 (MQ=255) ccGGTACCCACCACACCCCCCGGCAGGGCGAGGATCAACATGCCGCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGc < 2:129995/90‑1 (MQ=255) gCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGAc > 2:199551/1‑90 (MQ=255) ggAGGATCAACATGCCCCTCATCCCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcatca < 1:68401/89‑1 (MQ=255) ggAGGACCAACATCCCCCTCACCCCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcatca < 2:129777/89‑1 (MQ=255) cGAGGATCAACATGCCCCTCATCTCCAGCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcatca < 1:199995/90‑1 (MQ=255) cGAGGATCAACAGGCCCCTCATCCCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcatca < 2:198578/90‑1 (MQ=255) cGAGGATCAACAGCCCGCTCACCTCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGcatca < 1:154162/90‑1 (MQ=255) aTCAACAGGCCCCTCCCCCCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTc < 2:331215/90‑1 (MQ=255) ccGCTCATCTCCAGCACCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAG‑‑CACTACg < 1:338643/90‑1 (MQ=255) cccTCACCCCCACCCCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAG‑‑CACTACGc < 2:238819/88‑1 (MQ=255) aTCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAG‑‑CACTACGCACAAt < 1:111505/90‑1 (MQ=255) aGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCACACTACGCACAATacac > 1:39906/1‑90 (MQ=255) tGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGAc < 1:323574/40‑1 (MQ=255) tGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGAc > 2:323574/1‑40 (MQ=255) | AAAAGAAGCCGGTAGCCAGCACAATCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAG‑‑CACTACGCACAATACAC > NZ_CP009273/69835‑69984 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |