Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,540,683 | 0 | A | C | 52.9% | 11.1 / 17.0 | 17 | G248G (GGT→GGG) | yddG | aromatic amino acid efflux DMT transporter YddG |
Reads supporting (aligned to +/- strand): ref base A (5/3); new base C (0/9); total (5/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.05e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACATTCCATGCAGCATAAGCAA > NZ_CP009273/1540600‑1540744 | gcgcCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTAcc > 2:764250/1‑90 (MQ=255) gcCTTGCCAGAGCGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTGCAGGCGTAAAATACGAACCTACCGc > 2:773137/1‑90 (MQ=255) ttGCCAGAAGGAGAACGACAGGGGGGCGCGGAGCAGCCCGGCTCCAAGCGTTGGGGGAAGTACAGGCGTAAAATACGACCCTACCGCCAt < 2:426693/90‑1 (MQ=255) gaacgagaacgaCAGGGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGt < 1:576843/90‑1 (MQ=255) gaCAGCGGGGCGCTGAGCAGCCCGGCTGCAAGCGCGGAGGAAAGTACAGGCGTAAAATACGACCCTACCGCCATAATGGTGACATTGCCa < 2:549025/90‑1 (MQ=255) gaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGACCCTACCGCCATAATGGTGACATTGCCa < 1:213259/90‑1 (MQ=255) gaCAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCa < 1:330818/90‑1 (MQ=255) aaCAGCGGGCCCCTGAGCACCACGGCGGCAAGCGCTGAGGAAAGAACAGGCTTAAAATCCGACCCTACCGCCATAATGGTGACATTGCCa < 2:6150/89‑1 (MQ=255) gCGGGGCGCTGAGCAGCACGGTTGCAGGCGGTGAGGAAAGGACAGGCGTAAAATACAACCCTACCGCCATAATGGTGACATTGCCATGCa < 1:656152/90‑1 (MQ=255) gCGGGGCGCTGAGCAGCACGGCTGCAAGCGCGGAGGAAAGTACAGGCGAAAAATCCGACCCTACCGCCATAATGGTGACATTGCCATGCa < 2:406027/90‑1 (MQ=255) gcgcTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAAtata > 2:56997/1‑90 (MQ=255) cTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCaa > 2:481061/1‑83 (MQ=255) cTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCaa < 1:481061/83‑1 (MQ=255) gcagcaCGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACAt > 1:560036/1‑90 (MQ=255) gCAAGCGTTGAGGAAAGTACAGGCTTAAAATACGACCCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACATTCCATgcagca < 1:566898/90‑1 (MQ=255) gcggagaAAAGTACAGGCGAAAAATACAACCCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACATTCCATGCAGCATAAGCa < 2:796307/87‑1 (MQ=255) cTGAGGAAAGTACAGGCGTAAAATACGACCCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACATTCCATGCAGCATAAGCaa < 2:781205/90‑1 (MQ=255) | GCGCCTTGCCAGAACGAGAACGACAGCGGGGCGCTGAGCAGCACGGCTGCAAGCGCTGAGGAAAGTACAGGCGTAAAATACGAACCTACCGCCATAATGGTGACATTGCCATGCAATATACCGACATTCCATGCAGCATAAGCAA > NZ_CP009273/1540600‑1540744 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |