| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,058,958 | 0 | G | A | 56.7% | 7.5 / 32.0 | 30 | L57L (CTG→TTG) | yihO | MFS transporter |
| Reads supporting (aligned to +/- strand): ref base G (9/4); new base A (0/17); total (9/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.00e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.90e-04 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGC‑ATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA > NZ_CP009273/4058870‑4059038 | aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGC‑ATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa < 2:545545/90‑1 (MQ=255) gCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGC‑ATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGAt > 1:450903/1‑90 (MQ=255) gCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGC‑ATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGc < 2:151360/90‑1 (MQ=255) cAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGC‑ATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCAt > 1:235269/1‑90 (MQ=255) cAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGC‑ATATCGGTAAACGCGGAAAAAAACTTCGCCACCAGAAAGATAATCCCGCCAt < 2:730699/90‑1 (MQ=255) atttttttCCGCGAGTCGAGAAAAAACCCGGTGAGC‑ATATCGGTAAACGCGGTAAAAAACTTCGCCCCCAGAAAGATAATCCCGCCATAg < 2:222055/87‑1 (MQ=39) gcgAGTCGAGTAAAAATCCGGTGAGC‑ATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa > 1:662516/1‑90 (MQ=255) gcgAGTCGAGTAAAAACCGGGTGAGC‑ATATCGGTAAACGCGGAAAAAAACTTCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCa < 2:688036/90‑1 (MQ=255) gcgAGTCGAGAAAAAACCCGGTGAGC‑ATATCGGTAACCGCGGAAAAAACCTCCGCCACAAAAAAGATAATCCCGCCATAGTAAGCAGGCa < 2:450903/90‑1 (MQ=255) gcgAGTCAAGAAAAAACCCGGGGAGA‑ATATCGGTAACCGCGGAAAAAACCTCCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCa < 2:668477/90‑1 (MQ=255) aaaCCCGGTGAGC‑ATATCGGAAACCGCGGAAAAAACCTTCCCCCCAAAAAAGATAACCCCGCCATAGTAAGCAGGCACCCCTAATTCATc < 2:135288/90‑1 (MQ=255) aaaCCCGGTGAGC‑ATACCGGAAAACGCGAAAAAAACCTCCCCCCCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATc < 1:279449/90‑1 (MQ=255) aaaCCCGGGGAGC‑ATATCGGAAAACGCGGTAAAAAACTCCCCCCCCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc < 2:274967/90‑1 (MQ=255) aaaacGGAGAGCAATATCGTTAACCGCGGAAAAAAACTTCCCCACACAAAACATCAACCCGCCATAGTAAGCAGGCATCCCTAATTCATc < 1:606684/86‑1 (MQ=255) aaTCCGGGGAGC‑ATATGGGAAAACGGGAAAAAAACCTCCCCCACAAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCg < 2:704913/90‑1 (MQ=255) cGGTGAGC‑ATATCGGTAAACGCGGTAAAAAACTTCGCCACCAG‑AAGGTAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata > 1:403196/1‑90 (MQ=255) c‑ATATCGGTAAACGCGGTAAAAAACTTCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt < 1:705046/90‑1 (MQ=255) tatGGGAAACCGGGGTAAAAAACTCCCCCCCAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga < 2:528223/90‑1 (MQ=255) tatGGGAAAACGGGAAAAAAACCTCCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga < 1:77915/90‑1 (MQ=255) tatGGGAAAACGGGAAAAAAACCTACCCCCCAAAAAAGACACCCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga < 2:613662/90‑1 (MQ=255) cGGAAAACGCGGAAAAAAACTCCCCCACCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga < 2:248843/90‑1 (MQ=255) gaaaaCGCGGAAAAAAACTTCCCCACCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga < 2:620288/88‑1 (MQ=255) aaaCGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATa > 1:528223/1‑90 (MQ=255) cgcgGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACag > 2:8997/1‑90 (MQ=255) gcgGAAAAAAACTCCGCCACCAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACaga < 2:134230/90‑1 (MQ=255) cgGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACagag > 2:279478/1‑90 (MQ=255) gTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTa > 1:368224/1‑90 (MQ=255) tccccacaaAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga < 1:302363/87‑1 (MQ=255) tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga > 1:429467/1‑90 (MQ=255) cccccacaaAAAATATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga < 2:394662/87‑1 (MQ=255) | AATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGC‑ATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA > NZ_CP009273/4058870‑4059038 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |