Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,058,958 | 0 | G | A | 55.2% | 2.6 / 23.8 | 29 | L57L (CTG→TTG) | yihO | MFS transporter |
Reads supporting (aligned to +/- strand): ref base G (8/5); new base A (0/16); total (8/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.00e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTT‑CGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGA > NZ_CP009273/4058878‑4059025 | tttAGCCCCAATTTTTTTCCGCGAGTCGAGAAAAAACCGGGGGAGCATATCGGTAACGGGGGAAAAAACCTC‑CGCCACAAAAAAAATaac < 1:656390/90‑2 (MQ=255) ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTT‑CGCCACCAGAAAGATAATc > 2:293739/1‑90 (MQ=255) ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTT‑CGCCACCAGAAAGATAATc > 2:3151/1‑90 (MQ=255) ttttcCCGCGAGCCGAGAAAAAATCCGGTGAGCATATCGGAAAACGCGGAAAAAAACTT‑CGCCACCAAAAAGAAAATCCCGCCATAGTaa < 2:160746/90‑1 (MQ=255) ttttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTT‑CGCCACCAGAAAGATAATCCCGCCATAgctgt > 2:319057/1‑86 (MQ=255) tttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTT‑CGCCACCAGAAAGATAATCCCGCCATAGTAAGc > 1:32066/1‑90 (MQ=255) ccGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTT‑CGCCACCAGAAAGATAATCCCGCCATAGTAAGCAgg > 1:168229/1‑90 (MQ=255) cgcgAGTCGAGAAAAAACCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTT‑CGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGc < 2:216363/90‑1 (MQ=255) gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTT‑CGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa > 1:620407/1‑90 (MQ=255) gcgAGTCGAGTAAAAACCCGGTGAGCATATCGGTAAACGCGGTAAAAACCTT‑CGCCCCCAAAAAGATAACCCCGCCATAGTAAGCAGGCa < 2:47241/90‑1 (MQ=255) aaaTCCGGTGAGAATACCGGAAAACGCGGAAAAAACCTC‑CGCCCCAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATc < 1:120464/90‑1 (MQ=255) aaaCCCGGTGAGCATATCGGAAAACGCGGAAAAAAACTC‑CCCCACCAAAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATc < 2:620407/90‑1 (MQ=255) ggggAGCATATCGGAAAACGGGAAAAAAACCTC‑CCCCACAAAAAAGAAACCCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata < 2:596111/87‑1 (MQ=255) ggTGAGCATATCGGTAAACGCGGAAAAAAACTT‑CGCCACCAGAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtata < 2:152861/90‑1 (MQ=255) cATATCGGAAACCGCGGAAAAAAACTT‑CGCCCCAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt < 2:584247/90‑1 (MQ=255) tatCGGAAACCGCGGAAAAAAACTC‑CCCCACAAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga < 2:25341/90‑1 (MQ=255) tatCGGAAACCGCGAAAAAAACCTC‑CCCCCCAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga < 1:293739/90‑1 (MQ=255) tatCGGAAAACGCGGTAAAAAACTC‑CCCCACCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga < 2:32066/90‑1 (MQ=255) tCGGAAACCGGGGAAAAAAACTC‑CCCCCCCAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg < 2:63553/90‑1 (MQ=255) tCGGAAACCGCGGAAAAAAACTT‑CCCCACCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg < 2:109954/90‑1 (MQ=255) gaaaaCGAGGAAAAAAACTC‑CCCCAAAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga < 1:457842/88‑1 (MQ=255) gTAAACGCGGTAAAAAACTT‑CGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa < 2:165331/58‑1 (MQ=255) gTAAACGCGGTAAAAAACTT‑CGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa > 1:165331/1‑58 (MQ=255) gTAAACGCGGAAAAAAACTT‑CGCCCCCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga < 1:427079/90‑1 (MQ=255) aCGCGGTAAAAAACTT‑CGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACa > 1:179502/1‑90 (MQ=255) gcgGAAAAAACCTC‑CCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACaga < 1:109904/90‑1 (MQ=255) gcgGAAAAAAACTT‑CGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACaga < 2:302040/90‑1 (MQ=255) gcgGAAAAAAACTC‑CCCCCCAAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACaga < 2:419889/90‑1 (MQ=255) gggTAAAAAATTTCCCCCACGAAAAAGATACCCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACaga < 2:662763/89‑1 (MQ=255) | TTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTT‑CGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGA > NZ_CP009273/4058878‑4059025 |
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: ‑ |