Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,371,314 | 0 | G | A | 44.6% | 40.1 / 32.3 | 56 | F475F (TTC→TTT) | menD | 2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylic‑acid synthase |
Reads supporting (aligned to +/- strand): ref base G (25/6); new base A (0/25); total (25/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.57e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAACATCGCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATC > NZ_CP009273/2371227‑2371397 | gAACATCTCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGaa > 2:963644/1‑90 (MQ=255) gAACATCGCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGaa > 1:339109/1‑90 (MQ=255) tCGCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAAttt > 2:171905/1‑90 (MQ=255) ggcggcgTGCTAAAAAGGGACTTTTGGCGGCATCAGATAGAACCGCTCACGCTCGTTTTGCGGGGTTGGCAACAGCGAGAAATTTTGccc < 2:1618058/90‑1 (MQ=255) ggcggcgTGCAAAAAATGGAGTTTTTGGGGAATAAGAGGAAAACGCGCACGCTGGCTTTGGGGTGTTGGGAACAGAGAGAAAATTTGccc < 1:1459632/90‑1 (MQ=255) cggcgTGCAAAAAATGGACGTTTTGCGGCATAAGATAGAAACCCCCCCGCTCGTTTTGGGGCGTGGGCAACAGCGAGAAAATTTGCccgc < 2:928363/90‑1 (MQ=255) ggcgTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCccgcc < 1:3315459/90‑1 (MQ=255) cgTGCTAAAAAGGGACGTTTGGCGGAATAAAAAAGAACCCCCCCCGCTCGTTTTGCGGCGTTGGCAACAGCGAAAAAATTTCCCCGCCgt < 2:2034316/90‑1 (MQ=255) gTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCgtt > 2:1522464/1‑90 (MQ=255) gTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCgtt > 2:3315873/1‑90 (MQ=255) aaaTGGTCTTTTTGCGGCATAAGAAAAAACCCCTCGCGCTCGTTTTGGGGCGTTGGCAACGGGGAAAAAATTTCCCCGCCGTTGTTGTTc < 2:1451211/90‑1 (MQ=255) aaaTGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTc < 1:1257223/90‑1 (MQ=255) aaaTGGACGTTTTGCGGCATCAGAAAGAACCGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTc < 1:2256374/90‑1 (MQ=255) tGGACTTTTTGGGGCATCAGATAAAACCGCCCCCCCTCGTTTGGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcacc < 1:1522208/90‑1 (MQ=255) tGGACGTTTTGGGGCATCATAAAAAAACGCTCACGCTCGCCTTGCGGAGTTGGCAACAGCGTGAAAATTTACCCGCCGTTGCCGTTgacc > 2:1616087/1‑90 (MQ=255) tGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcacc > 1:3402443/1‑90 (MQ=255) tGGACGTTTGGCGACATAAGAAAAAACCCCCCCCGCTCGTTTTGCGGCGTTGGAAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcacc < 2:2838274/90‑1 (MQ=255) gggACGTTTGGGGGCATCAGATAAAACCCCCCACGCTCGTTTTGGGGCGTGGGCAACAGAGAAAAATTTTCCCCGCCGTTGTTGTTcacc < 2:2384435/89‑1 (MQ=255) gggggTTTTGCGGAATAAGATAAAACCCCCCCCGCTCGTTTTGGGGCGTTGGAAACAGGGAAAAAATTTGCCCGCCGTTGTTGTTcacca < 2:3281235/86‑1 (MQ=255) gccGTTTGGCGGCATCAGATAAAACCGCCCCCGCTCGTTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcaccac < 2:1890583/88‑1 (MQ=255) gatgTTTTGCGGCATCAAAAAGAACCGCTCCCCCTCGTTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcaccac < 2:496422/87‑1 (MQ=255) gACTTTTGGCGGCACCAGATAAAAACCCTCCCGCTCGCTTGGCGGCGTTGGCAACAGCGAAAAATTTTGCCCGCCGTTGTTGTTcaccac < 1:2841724/90‑1 (MQ=255) gACGTTTTGCGGCATCAGAAAAAAACGCTCCCGCTCGCTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcaccac < 1:2966977/90‑1 (MQ=255) gACGTTTGGCGGCATCAGATAGAACCGCCCCCGCTCGTTTTGCGGCGTTGGCAACAGCGAAAAATTTTGCCCGCCGTTGTTGTTcaccac < 2:356107/90‑1 (MQ=255) aCGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACcaca > 2:3001615/1‑90 (MQ=255) aCGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACcaca > 1:2235392/1‑90 (MQ=255) tttGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTa > 2:497167/1‑90 (MQ=255) ttggCGGCATAAGAAAAAACCCCCCCCGCTCGTTTGGGGGGGTTGGCAACAGGGAAAAATTTTGCCCGCCGTTGTTGTTCACCACAATTa < 2:3402443/87‑1 (MQ=255) aTCAAATAAAAACCCTCCCCCTCGTTTTGCGGAGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAgc < 1:781198/90‑1 (MQ=255) tCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAgcg > 1:1547110/1‑90 (MQ=255) aaCGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGaaa > 1:2417546/1‑90 (MQ=255) cgctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAAcc > 1:3241490/1‑90 (MQ=255) gctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcacca < 1:2848577/63‑1 (MQ=255) gctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcacca > 2:2848577/1‑63 (MQ=255) gctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCt > 1:1707018/1‑90 (MQ=255) gctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCt > 1:2781257/1‑90 (MQ=255) ctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTg > 1:1445136/1‑90 (MQ=255) ctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTg > 1:911523/1‑90 (MQ=255) tcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGa > 2:85193/1‑90 (MQ=255) cccgctcgttTTGGGGCGTGGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGAc < 1:3076586/88‑1 (MQ=255) cccgctcgttTGGCGGCGTTGAAAACAGAGAAAAATTTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGAc < 2:1707018/88‑1 (MQ=255) cacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGAc > 2:2777608/1‑90 (MQ=255) ccgctcgttTTGCGGGGTGGGAAACAGGAAAAAAATTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACg < 2:598749/89‑1 (MQ=255) cgctcgctTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGc > 1:1700322/1‑90 (MQ=255) ttGCGGCTTGGGAAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGc < 1:2599828/90‑1 (MQ=255) gcggcggcgGCAACAGAGAAAATCTTGCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCa < 2:212462/90‑1 (MQ=255) gcggcgTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCa > 2:2525019/1‑90 (MQ=255) ggcgTTGCCAACAGCGAAAAATTTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAgc < 2:2514168/90‑1 (MQ=255) cgTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAgcgc > 1:3345877/1‑90 (MQ=255) ttGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGt > 2:1337218/1‑90 (MQ=255) ttGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGt > 2:2711263/1‑90 (MQ=255) gaaaCGGCGAAAAATTTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTga < 2:784734/88‑1 (MQ=255) gCAACAGCGAAAAATTTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTga < 2:48544/90‑1 (MQ=255) gCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTga < 1:2711263/90‑1 (MQ=255) gCAACACAGAAAAATTTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTga < 2:2933896/90‑1 (MQ=255) aaGCGAAAAATTTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATc < 2:339109/89‑1 (MQ=255) | GAACATCGCGGCGGCGTGCTCAAAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATC > NZ_CP009273/2371227‑2371397 |
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: ‑ |