Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,058,958 | 0 | G | A | 48.9% | 35.5 / 34.1 | 47 | L57L (CTG→TTG) | yihO | MFS transporter |
Reads supporting (aligned to +/- strand): ref base G (21/3); new base A (0/23); total (21/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.81e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTA > NZ_CP009273/4058871‑4059044 | aTTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGaa < 1:632318/90‑1 (MQ=255) ttGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAg > 1:1317132/1‑90 (MQ=255) ttGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAg > 1:546983/1‑90 (MQ=255) ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATc > 2:752815/1‑90 (MQ=255) ttAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATc > 2:1042011/1‑90 (MQ=255) cAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCAt > 2:518284/1‑90 (MQ=255) ttttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa > 2:504122/1‑68 (MQ=255) ttttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGa < 1:504122/68‑1 (MQ=255) ttttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAg > 1:994420/1‑90 (MQ=255) tttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGc > 1:1012280/1‑90 (MQ=255) tttCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGc > 2:1039017/1‑90 (MQ=255) ccGCGAGTCGAGAAAAAACCCGGTGAGCATATCGGTAACCGCGGAAAAAAACTCCGCCACCAAAAAGATAATCCCGCCATAGTAAGCAgg < 1:264004/90‑1 (MQ=255) cgggagtcgaaaAAAAATCGGGGGAGCATATGGGAAAACGCGGAAAAAACCTTCGCCACAAAAAAGATAACCCCGCCATAGTAAGCAGGc < 2:1204802/87‑1 (MQ=255) ccggagccgagaAAAAACCCGGTGAACATAACGAAAAACCGGGAAAAAACCTCCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGc < 1:1586363/87‑1 (MQ=255) gggagtcgagaAAAAACCCGGGGAGCATATGGGAAAACGGGAAAAAAACCTCCCCCACAAAAAAGATACCCCCGCCATAGTAAGCAGGCa < 2:1496337/88‑1 (MQ=255) gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa > 1:263332/1‑90 (MQ=255) gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa > 2:276307/1‑90 (MQ=255) gcgAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa > 2:370358/1‑90 (MQ=255) gcgAGTCGAGAAAAAACCCGGTGACCATATCGGAAACCGGGGTAAAAACCTCCCCCCCCAAAAAAATAACCCCGCCATAGTAAGCAGGCa < 1:1401400/90‑1 (MQ=255) cggagaAAAAACCGGGGGAGCATATCGTAAAGCGCGAAAAAAACCTTCCCCACCAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCt < 1:1365245/88‑1 (MQ=255) aacccGGTGACAATACCGGAAAACGGGAAAAAAACCTTCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCg < 2:994420/87‑1 (MQ=255) cGGTGAGCATATCGAAAAACGCGAAAAAAACCTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGtat < 2:1290360/90‑1 (MQ=255) aatatCGAAAAACGAGAAAAAAACCTCCCCCACAAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAAttt < 2:835291/89‑1 (MQ=255) tatCGGTAACCGGGGAAAAAACCTCCCCCACAAAAAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga < 1:863762/90‑1 (MQ=255) tatCGAAAAACGCGAAAAAAAACTTCCCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTga < 1:772479/90‑1 (MQ=255) tCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg > 1:454980/1‑90 (MQ=255) gcGGAAAACGCGAAAAAAACCTCCCCCACAAAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagg < 2:1132474/89‑1 (MQ=255) gggaaaaCGCGAAAAAAACCTCCCCCACCAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTgagga < 1:1340549/86‑1 (MQ=255) gTAACCGCGGAAAAAACCTCCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga < 2:1412326/90‑1 (MQ=255) gTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga > 1:11226/1‑90 (MQ=255) gTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga > 1:778834/1‑90 (MQ=255) aaaaaCGCGAAAAAAATCTCCGCCCCAAAAAAAATACCCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGgaga < 1:991876/88‑1 (MQ=255) aaaaCGCGAAAAAAACCTCCCCCCCAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGAt < 1:1327505/89‑1 (MQ=255) aaaCGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATa > 2:176996/1‑90 (MQ=255) aaCGCGAAAAAAAACTCCCCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAc < 2:778834/90‑1 (MQ=255) gcgGAAAAAAACTTCGCCACAAAAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACaga < 1:248443/90‑1 (MQ=255) ggaaaaaacCTTCGCCACAAAAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGt < 1:20529/87‑1 (MQ=255) gTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTa > 1:1496879/1‑90 (MQ=255) acaaacCTCCGCCAAAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAcc < 2:797531/88‑1 (MQ=255) aaaaacCTCCCCCACAAAAAAAATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTAcc < 1:1406082/90‑1 (MQ=255) ctccGCCACAAAAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCa < 2:11226/87‑1 (MQ=255) cTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCa > 2:31635/1‑90 (MQ=255) cTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCa > 1:1223302/1‑90 (MQ=255) tcccccacaaAAAAAAAAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCag < 2:1285459/86‑1 (MQ=255) tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa < 2:1217420/39‑1 (MQ=255) tCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCa > 1:1217420/1‑39 (MQ=255) ccaGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGcctgt > 1:797531/1‑86 (MQ=255) | ATTTGCCTTTAGGCCCAATATTTTTCCGCGAGTCGAGTAAAAATCCGGTGAGCATATCGGTAAACGCGGTAAAAAACTTCGCCACCAGAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTA > NZ_CP009273/4058871‑4059044 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |