Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,058,967 | 0 | T | C | 61.5% | 7.6 / 24.8 | 26 | I54V (ATT→GTT) | yihO | MFS transporter |
Reads supporting (aligned to +/- strand): ref base T (8/2); new base C (0/16); total (8/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.88e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA > NZ_CP009273/4058879‑4059038 | ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATAATc > 1:846523/1‑90 (MQ=255) tttCGGCGAGTCGAGAAAAAGCCCGGTGAGCATATCGGAAACCGCGGAAAAAACCTTCCCCCCC‑AAAAAAAAAACCCCGCCATAGTAAGc < 2:294052/90‑1 (MQ=255) tttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGc > 1:735266/1‑90 (MQ=255) ttcccGGGAGTCGAGAAAAAACCCGGTGAGCATATCGGAAACCGCGAAAAAAACCTTCCCCACC‑AAAAAGAAAACCCCGCCATAGTAAGc < 2:41623/87‑1 (MQ=255) ttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCa > 2:808515/1‑90 (MQ=255) ccGCGAGTGGAAAAAAAACCCGGAGAGCATATCGAAAAACGCGAAAAAAAACTCCCCCACA‑AAAAAAATAACCCCGCCATAGTAAGCAgg < 1:160147/90‑1 (MQ=255) gggagtcgaaaAAAAATCGGGGGAGAATACGGGAAAACGGGAAAAAAACCTCCCCCACA‑AAAAAAATAACCCCGCCATAGTAAGCAGGCa < 1:614329/88‑1 (MQ=255) gcgAGTCGAAAAAAAACCGGGTGAGCATATCGGAAACCGCGGAAAAAAACTCCCCCACA‑AAAAAAATAATCCCGCCATAGTAAGCAGGCa < 2:221035/90‑1 (MQ=255) gcgAGTCGAAAAAAAACCGGGTGAAAATCTGGGAAACCGCGAAAAAAACCTCCCCCCCA‑AAAAAAAAAACCCCGCCATAGTAAGCAGGCa < 1:374913/90‑1 (MQ=255) gcgAGTAGAGAAAAAACCGGGTGAGAATATGGGAAAACGGGGAAAAAATCTCCGCCCCC‑AAAAAGATAACCCCGCCATAGTAAGCAGGCa < 2:394396/90‑1 (MQ=255) aaaTCCGGTGCCAATATGGGAAAACGCGGAAAAAAACTTCCCCACA‑AAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc < 2:204149/90‑1 (MQ=255) aaaTCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATc > 2:392822/1‑90 (MQ=255) aaTCCGGTGAGCATATCGGAAAACGAGAAAAAAACCTCCCCCACA‑AAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCg < 1:600366/90‑1 (MQ=255) cGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtat > 1:42791/1‑90 (MQ=255) cGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtat > 2:600366/1‑90 (MQ=255) cATATCGGAAAACGCGGAAAAAACCTCCCCCCCA‑AAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt < 2:294314/90‑1 (MQ=255) tatCGATAAACGCGGAAAAAACCCCCCCCACA‑AAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga < 1:666340/90‑1 (MQ=255) tCGGTAACCGCGTAAAAAAACTCCCCCACC‑AAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg < 2:622782/90‑1 (MQ=255) gTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga > 1:380224/1‑90 (MQ=255) atAAACGCGAAAAAAATCTTCCCCGCA‑AAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga < 1:572723/89‑1 (MQ=255) aaaaCGCGGAAAAAAACTCCCCCCCA‑AAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGAt < 2:42791/89‑1 (MQ=255) ggaaaaaaCCTCCCCCCCAAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt < 2:76621/88‑1 (MQ=255) aaaaacCTCCCCCACA‑AAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAcc < 2:600668/90‑1 (MQ=255) aaaaaaCTCCCCCACC‑AAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAcc < 1:275340/90‑1 (MQ=255) ttCCCCACA‑AAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCag < 2:282122/90‑1 (MQ=255) tCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCaga > 1:670531/1‑90 (MQ=255) | TTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACC‑AGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGA > NZ_CP009273/4058879‑4059038 |
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: ‑ |