Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 890,381 | 0 | A | C | 52.4% | 21.5 / 25.4 | 21 | intergenic (+29/‑66) | potF/potG | spermidine/putrescine ABC transporter substrate‑binding protein PotF/putrescine ABC transporter ATP‑binding subunit PotG |
Reads supporting (aligned to +/- strand): ref base A (9/1); new base C (0/11); total (9/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.40e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTC > NZ_CP009273/890311‑890465 | gtgCGCCCCCGCGCGGGGACCAAAGTGAAGAGCGAAAAATAACCCGCAGTCGTAGATGCCGGGGGGGCGCCCCACACCCGCCGGCAATTc < 1:395194/90‑1 (MQ=255) gcgcCCCCGCGCGTGGACCAAAGTAAAGAGCGGAAAATAACCCGCAGTCGTAGATGCGGGGGGGGCGCCCCACCCCCGCCGGCAATTCGc < 2:113970/90‑1 (MQ=255) gcgcCCCCGCGCGGGGACAAAAGTGAAGAGCGGAAAATAACCCGCAGTCGTAGAGGCGGGGGGGGCGCCCCACCCCCGCCGGCAATTCGc < 2:206207/90‑1 (MQ=255) cccccGCGCGTGGACCAAAGTTAAGAGCGAAAAATAACCCCCAGTCGAAGAGGCGGGAGGGGCGCCCCACCCCCGCCGGCAATTCGcccc < 2:591977/88‑1 (MQ=255) gcgcgTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTAt > 2:563366/1‑90 (MQ=255) gcgcgTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTAt > 2:520333/1‑90 (MQ=255) gcgcgTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTAt > 2:25038/1‑90 (MQ=255) cgTGGACAAAAGTGAAGAGCGGAAAAAAACCCGCAGTCGTAGATGCCGGGGGGGCGCCCCACCCCCGCCGGCAATTCGCACCATTATGGt < 2:508273/90‑1 (MQ=255) gTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTg > 2:116794/1‑90 (MQ=255) tGGACAAAAGGGAAGAGGGAAAAATAACCCGCAGTCGTAGATGCGGGGGGGGCCCCCCACCCCCGCCGGCAATTCGCACCATTATGGTgc < 2:294626/90‑1 (MQ=255) gggACAAAAGTGAAGAGCGAAAAATAACCCGCAGTCGTAGAGGCGGGGGGGGCCCCCCACCCCCGCCGGCAATTCGCACCATTATGGTgc < 2:608719/89‑1 (MQ=255) gACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTgcgc > 1:655183/1‑90 (MQ=255) aaGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACAt > 1:563229/1‑90 (MQ=255) gggCGAAAAAAAATCCCCAGTCGAAGAGGCGGGGGGGCCGCCCCCCCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTc < 2:655183/88‑1 (MQ=255) gagCGGAAAATAACCCGCAGTCGTAGATGCGGGGGGGGCGCCCCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTc < 2:146416/90‑1 (MQ=255) cgaaaaaTAACCCGCAGCCGTAGAGGCGGGGGGGGCCCCCCCCACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAAt < 2:650955/87‑1 (MQ=255) aataacCCGCAGTGGTAGAGGCGGGGGGGGCCCACCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCgg < 2:573007/90‑1 (MQ=255) atCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGagag > 1:425167/1‑90 (MQ=255) tCCGCAGTCGTAGATGCGGGGGGGGCGCCCCCCCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCATTGCCGGAGagc < 1:373200/90‑1 (MQ=255) tGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACg > 1:701468/1‑90 (MQ=255) ggCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTc > 1:628281/1‑90 (MQ=255) | GTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTC > NZ_CP009273/890311‑890465 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |