Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,058,958 | 0 | G | A | 53.8% | 10.9 / 30.3 | 26 | L57L (CTG→TTG) | yihO | MFS transporter |
Reads supporting (aligned to +/- strand): ref base G (6/6); new base A (0/14); total (6/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.01e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAG > NZ_CP009273/4058869‑4059039 | aaaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑Ag < 1:1677/90‑1 (MQ=255) aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGa < 2:466047/90‑1 (MQ=255) aaTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGa < 1:97982/90‑1 (MQ=255) tGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGa < 2:549765/90‑1 (MQ=255) ccTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATa < 2:902415/90‑1 (MQ=255) gCCAAATTTTTTTCGGCGAGCCGAGAAAAAATCCGGTGAGAATATCGGAAAACGCGGAAAAAACCTCCCCCACA‑AAAAAGATACCCCCGc < 2:890406/90‑1 (MQ=255) ccAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATAATCCCGCCa > 1:564428/1‑90 (MQ=255) tttcCCGGGAGCCGAGAAAAAACCCGGTGAGAATATCGGAAAACGCGGAAAAAACCTCCCCCACA‑AAAAAAAAAACCCCGCCATAGTAAg < 1:91523/90‑1 (MQ=255) tttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGaa < 2:868079/68‑1 (MQ=255) tttCCGCGAGTCGAGTAAAAATCCGGCGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGaa > 1:868079/1‑68 (MQ=255) ttcccGCGAGCCGAAAAAAAATCCGGTGAAAATATCGAAAAACGCGGAAAAAACCTCCCCCAAA‑AAAAAAATAACCCCGCCATAGTAAGc < 2:575813/87‑1 (MQ=255) cgcgAGTCGAGAAAAAACCCGGTGAGCATATGGGAAAACGCGGAAAAAAACTTCCCCACC‑AAAAAGATAACCCCGCCATAGTAAGCAGGc < 1:837101/90‑1 (MQ=255) gcgGGTCGAAAAAAAATCGGGTGAGAATAGGGGAAAACGGGAAAAAAACCTCCCCCACA‑AAAAAAATAATCCCGCCATAGTAAGCAGGCa < 1:443219/90‑1 (MQ=255) gcgAGTGGAGAAAAAATCCGGTGAGCATATCGGAAAACGCGGAAAAAAACTCCCCCACA‑AAAAAGATAACCCCGCCATAGTAAGCAGGCa < 1:54504/90‑1 (MQ=255) gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCa > 2:151197/1‑90 (MQ=255) gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCa > 2:837101/1‑90 (MQ=255) gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCa > 2:659376/1‑90 (MQ=255) gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGAAAAAAACTTCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCa < 1:552918/90‑1 (MQ=255) aaacTCGGGTGAGCATATGGGAAAACGCGGTAAAAACCTTCCCCCCC‑AAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCAt < 1:293905/90‑1 (MQ=255) aaacTCCGGTGAAAATA‑CGGAAAACCGGGAAAAAAACCTCCCCCCCAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCAt < 2:188773/90‑1 (MQ=255) aaaaTCGGGTGACCATACCGTAAAACGCGGAAAAAAACTCCCCCCCA‑AAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCAt < 1:155920/90‑1 (MQ=255) ggTGAGCATATAGAAAAACGCGAAAAAAACCTTCCCCACA‑AAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata < 2:879104/90‑1 (MQ=255) tatCGAAAAACGGGGAAAAAACCTCCCCCACA‑AAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga < 2:465674/90‑1 (MQ=255) cgtaaaaCGCGAAAAAAACCCCCGCCACA‑AAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga < 2:326637/86‑1 (MQ=255) gaaaaaaacCTCCGCCACA‑AAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt < 2:537137/87‑1 (MQ=255) cTTCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCa > 1:541441/1‑90 (MQ=255) ccccaca‑aAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCagag < 2:541441/88‑1 (MQ=255) | AAATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAG > NZ_CP009273/4058869‑4059039 |
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: ‑ |