Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 499,311 | 0 | A | C | 53.8% | 13.0 / 15.6 | 13 | V281G (GTG→GGG) | fsr | fosmidomycin MFS transporter |
Reads supporting (aligned to +/- strand): ref base A (6/0); new base C (0/7); total (6/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.83e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAG > NZ_CP009273/499225‑499400 | aaaTCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTTATCCCCTACAGGCCCGCCGATCACCg > 1:121523/1‑90 (MQ=255) aaaCGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGc > 2:4427/1‑90 (MQ=255) gagCCCCAAATCACATATTTCCGCCCAA‑TTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAAc > 2:174933/1‑90 (MQ=255) gagCCCCAAATCACATATTTCCGCCCAA‑TTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAAc > 2:251990/1‑90 (MQ=255) ccaaaaCCACAT‑TTTCCCCCCAATTTTTCCCCCCAACGGCCCCGCCGATCCCCGTCCCTGCCGCAACGGCAAACAGGAAGGCAAACAGAt < 2:121523/87‑1 (MQ=255) caTTTTCCCCCCCAT‑TTTTACCCCCTACAGGCCCGCCGATCCCCGTCCCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAg < 1:4427/90‑1 (MQ=255) ttCCCCCCAATTTTTACCCCCAACGGCCCCGCCGATCCCCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCAtt < 1:82433/90‑1 (MQ=255) gCCCAA‑TTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGga > 2:196362/1‑90 (MQ=255) ccccACAGCCCCCCCGATCCCCTTCCCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGaa < 1:174933/90‑1 (MQ=255) ccgccgATCCCCGTCCCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATc < 1:228271/90‑1 (MQ=255) ccgccgATCCCCGTCCCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATc < 2:77948/90‑1 (MQ=255) ccgccgATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATc > 2:277205/1‑90 (MQ=255) ccgccgACCCCCGTCCCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATc < 1:227459/90‑1 (MQ=255) aCCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAg < 1:277205/90‑1 (MQ=255) | AAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAA‑TTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAG > NZ_CP009273/499225‑499400 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |