Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 890,386 | 0 | A | C | 51.8% | 29.7 / 36.0 | 27 | intergenic (+34/‑61) | 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 (10/3); new base C (0/14); total (10/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.39e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCA > NZ_CP009273/890297‑890472 | cGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACcaca > 1:1093091/1‑90 (MQ=255) cGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACcaca > 2:840133/1‑90 (MQ=255) cGAAAATCGACCGGGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACcaca > 1:901099/1‑90 (MQ=255) gcgcCCCCGGGGGGGGACAAAAGTGAAGAGGGAAAAATAATCCGCAGTCGTAGATGCGGGGGGGGCGCCCCCCCCCCGCCGGCAATTCGc < 2:1171746/90‑1 (MQ=255) gcgcCCCCGGGCGGGGACCAAAGTGAAGGGCGGAAAATAATCCGCAGTCGTAGATGCGGGGGGGGCGCCCCACCCCCGCCGGCAATTCGc < 2:306389/90‑1 (MQ=255) gcgcCCCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAACCCGCAGTCGTAGATGCCGGGGGGGCGCACCACCCCCGCCGGCAATTCGc < 2:151324/90‑1 (MQ=255) gcgcCCCCGCGCGTGGACCAAAGTGAAGAGCGAAAAATAATCCGCAGTCGTAGATGCGGGAGGGGCGCACCACCCCCGCCGGCAATTCGc < 1:1198252/90‑1 (MQ=255) gcgcCCCCGCGCGGGGACCAAAGTGAAGAGCGAAAAATAACCCGCAGTCGTAGATGCCGGAGGGGCGCCCCACCCCCGCCGGCAATTCGc < 1:713442/90‑1 (MQ=255) gcgcACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGc > 2:447246/1‑90 (MQ=255) gcccccGCGGGTGGACAAAAGGGAAGGGCGAAAAATAACCCGCAGTCGTAGATGCGGGGGGGGCGCACCACCCCCGCCGCCAATTCGCAc < 2:1201529/87‑1 (MQ=255) gcccccGCGCGTGGACCAAAGTGAAGAGCGGAAAATAACCCGCAGTCGTAGATGCGGGAGGGGCGCACCACACCCGCCGGCAATTCGCAc < 1:954767/87‑1 (MQ=255) ccGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCAtt > 1:823831/1‑90 (MQ=255) gcgcgTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTAt > 1:1044272/1‑90 (MQ=255) acaaaaGTAAAGAGCGAAAAATAACCCGCAGTCGTAGATGCGGGAGGGGCGCACCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCt < 1:605782/87‑1 (MQ=255) ccAAAGCAAAGAGAGAAAAATACCCCGCGTGGGTAGATGCCGGGGGGGCCCACAACACCCGCCGGCAATTCGCACCATTATGGTGCGCtt < 2:901099/90‑1 (MQ=255) caaaaGTGAAGCGCGAAAAAAAACCCGCAGTCGTAGAGGCGGGGGGGGCGCCCCCCCCCCGCCGGCAATTCGCACCATTATGGTGCGCtt < 2:1093091/88‑1 (MQ=255) caaaaGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGGGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCtt < 1:1109708/88‑1 (MQ=255) aaGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACAt > 2:654542/1‑90 (MQ=255) gCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCaa > 1:572136/1‑90 (MQ=255) cccGCAGTCGAAGGTGCGGGGGGGGCCCCCCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGagc < 2:135369/89‑1 (MQ=255) gggTCGTAGAGGCGGGGGGGGCCCCCCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGcc < 2:1246882/88‑1 (MQ=255) tAGATGCGGGGGGGGCCCCCCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAAt < 1:499665/90‑1 (MQ=255) tAGATGCCGGGGGGGCGCCCCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAAt < 1:1325085/90‑1 (MQ=255) aGATGCGGGAGGGGCCCCCCACCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATg < 1:840133/90‑1 (MQ=255) ccGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCt > 2:1011787/1‑90 (MQ=255) ccGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCt > 2:1206904/1‑90 (MQ=255) ccGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCt > 1:1246882/1‑90 (MQ=255) ggggggggCCCCCCCCCCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTa < 2:49324/86‑1 (MQ=255) gcgcACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTGTctctta > 1:814640/1‑88 (MQ=255) gcACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGcc > 2:1114377/1‑90 (MQ=255) caccCCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCa < 1:162001/90‑1 (MQ=255) | CGAAAATCGACCGTGTGCGCACCCGCGCGTGGACCAAAGTGAAGAGCGGAAAATAATCCGCAGTCGTAGATGCCGGAGGGGCGCACCACACCCGCCGGCAATTCGCACCATTATGGTGCGCTTGCACACATTCAATGCCGGAGAGCAGCCGTGAATGACGCTATCCCTCGCCCGCA > NZ_CP009273/890297‑890472 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |