| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,058,951 | 0 | G | C | 55.0% | 19.2 / 28.6 | 20 | A59G (GCG→GGG) | yihO | MFS transporter |
| Reads supporting (aligned to +/- strand): ref base G (6/3); new base C (0/11); total (6/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.17e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-05 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCC‑ACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTA > NZ_CP009273/4058870‑4059028 | aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCC‑ACC‑AGa < 2:763331/90‑1 (MQ=255)aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCC‑ACC‑AGa < 2:451114/90‑1 (MQ=255) ttcccGGGAGTCGAAAAAAAACCGGGTGAGCATATCGGAAACCGGGGAAAAAACCTCCCCC‑CCA‑AAAAAAAAAACCCCGCCATAGTAAGc < 1:700668/87‑1 (MQ=255) ttcccGCGAGTCGAGAAAAAACCCGGGGAGCATATCGGAAAACGGGGAAAAAACCTTCCCC‑ACA‑AAAAAAATAATCCCGCCATAGTAAGc < 1:149668/87‑1 (MQ=255) tCCGCGAGTCAAGAAAAAACCCGGTGAGCATATCGGTAAACGCGGAAAAAACCTCCCCC‑CCC‑AAAAAGATAATCCCGCCATAGTAAGCAg < 2:600780/90‑1 (MQ=255) gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCC‑ACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCa > 2:835852/1‑90 (MQ=255) gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCC‑ACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCa < 2:5546/90‑1 (MQ=255) ccgagccgaaaAAAAACCCGGTGAGCATATCGGAAAACGCGAAAAAAACCTCCCCC‑ACC‑AAAAAAATAACCCCGCCATAGTAAGCAGGCa < 1:345002/89‑1 (MQ=255) ggggAGCATATCGGAAAACGCGGAAAAAAACTCCCCC‑CCC‑AAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata < 2:235797/87‑1 (MQ=255) ggTGAGCATATCGTTAAACGCGGTAAAAACCTCCCCC‑ACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata < 2:235796/90‑1 (MQ=255) ggTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCC‑ACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata > 2:668064/1‑90 (MQ=255) gggTAAACG‑GGAAAAAAACCTCCCC‑CCCAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga < 2:752505/89‑1 (MQ=255) gTAAACGCGGTAAAAAACTTCGCC‑ACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga > 2:416978/1‑90 (MQ=255) aaaaaCGCGAAAAAAACCTCCCCC‑ACA‑AAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga < 1:835852/88‑1 (MQ=255) aaaCGCGGTAAAAAACTTCGCC‑ACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATa > 2:607582/1‑90 (MQ=255) accgcgGAAAAAAACTCCCCC‑CCA‑AAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc < 2:48395/88‑1 (MQ=255) ggggaaaaaaCCTCCCCAACA‑AAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACaga < 1:25333/86‑1 (MQ=255) ggTAAAAAACTTCGCC‑ACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt > 2:586584/1‑90 (MQ=255) gaaaaaaacCTCCCCC‑ACA‑AAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt < 1:632400/87‑1 (MQ=255) gTAAAAAACTTCGCC‑ACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTa > 2:666475/1‑90 (MQ=255) | AATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCC‑ACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTA > NZ_CP009273/4058870‑4059028 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |