| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,006,268 | 0 | T | A | 72.2% | 24.7 / 10.8 | 19 | intergenic (‑92/‑145) | metR/metE | HTH‑type transcriptional regulator MetR/5‑methyltetrahydropteroyltriglutamate‑‑ homocysteine S‑methyltransferase |
| Reads supporting (aligned to +/- strand): ref base T (3/2); new base A (5/8); total (8/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.08e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGA > NZ_CP009273/4006186‑4006344 | tCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGt > 2:112751/1‑90 (MQ=255) tCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCtttt > 2:123690/1‑90 (MQ=255) ataatTTGCGTTTGAGGAATATACAGTACCCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATg < 2:317213/90‑1 (MQ=255) ttGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAAACCTTCGTCTTTTAAATTTATGGTGCGTTGGCt < 2:360444/90‑1 (MQ=255) ggAATAAACAGAAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAAACCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGt < 2:344092/90‑1 (MQ=255) atatACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAAACCTTCGTCTTTTAAATTTATGGTGCGTCGGCTGCGTTtc > 1:297370/1‑90 (MQ=255) ataCAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTtctc > 2:430016/1‑90 (MQ=255) gTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAAACCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCAccc > 1:144375/1‑90 (MQ=255) aCCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAAACCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCgg > 1:489129/1‑90 (MQ=255) ccGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAAACCTTCGTCTTTTAAATTTATGGTGCGTCGGCTGCGTTTCTCCACCCCGGt < 2:297370/90‑1 (MQ=255) aaTTATGGATGTGTAAACATCTGGACGGCTAAAACCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCactt < 2:378489/90‑1 (MQ=255) aTTATGGATGTGTAAACATCTGGACGGCTAAAAACCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCactta < 1:291042/90‑1 (MQ=255) tGGATGTGTAAACATCTGGACGGCTAAAAACCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTc > 2:354177/1‑90 (MQ=255) tgtAAACATCTGGACGGCTAAAAACCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAg < 1:354177/90‑1 (MQ=255) cTGGACGGCTAAAAACCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCact > 2:123602/1‑69 (MQ=255) cTGGACGGCTAAAAACCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCact < 1:123602/69‑1 (MQ=255) cTGGACGGCTAAAAACCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCgggg < 2:173716/90‑1 (MQ=255) tGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGa < 2:11939/90‑1 (MQ=255) tGGACGGCTAAAAACCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGa < 2:266509/90‑1 (MQ=255) | TCACTTCGGCATGAATAATTTGCGCTTGAGGAATATACAGTAACCGCCAATTATGGATGTGTAAACATCTGGACGGCTAAAATCCTTCGTCTTTTAAATTTATGGTGCGTTGGCTGCGTTTCTCCACCCCGGTCACTTACTTCAGTAAGCTCCCGGGGA > NZ_CP009273/4006186‑4006344 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |