Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 277,985 | 0 | T | G | 40.0% | 56.8 / 109.9 | 75 | intergenic (‑291/‑4) | yagA/yagE | IS481‑like element ISEc19 family transposase/2‑keto‑3‑deoxygluconate aldolase |
Reads supporting (aligned to +/- strand): ref base T (31/14); new base G (19/11); total (50/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.27e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACC > NZ_CP009273/277899‑278075 | aCCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCa < 2:644884/91‑1 (MQ=255) aCCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCa < 2:681864/91‑1 (MQ=255) tcCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCg < 1:559343/91‑1 (MQ=255) gTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTc < 2:425873/91‑1 (MQ=255) tGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCcgcg > 2:792446/1‑91 (MQ=255) tGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCcgcg > 2:715074/1‑91 (MQ=255) tGCGTCATACTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCACGGAGAGCTCACGCCCAAGTCcgcg > 2:406403/1‑91 (MQ=255) gCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgt > 2:772408/1‑91 (MQ=255) gCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgt < 1:530515/91‑1 (MQ=255) gCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCgt > 2:771998/1‑91 (MQ=255) gCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCgt > 1:388921/1‑91 (MQ=255) cGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTCTAATGTAGAACAAAATGCTTACGCAAGGAGATCTCATGCCGCAGTCCGCgtt > 1:788675/1‑91 (MQ=255) aTAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTc < 1:339389/91‑1 (MQ=255) gTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑g < 2:13992/91‑1 (MQ=255) gTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑g < 2:568317/91‑1 (MQ=255) gTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑g < 1:472907/91‑1 (MQ=255) gATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGaa < 2:424620/91‑1 (MQ=255) aTCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGCTCAC‑GGGAAt > 2:450206/1‑91 (MQ=255) tCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATc > 2:203375/1‑91 (MQ=255) tCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATc < 2:342219/91‑1 (MQ=255) cGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCg > 1:723284/1‑66 (MQ=255) cGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCg < 2:723284/66‑1 (MQ=255) gTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAACTGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATGAt > 2:314470/1‑91 (MQ=255) gTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCAt > 1:289476/1‑91 (MQ=255) taaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgt > 2:572506/1‑73 (MQ=255) taaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCgt < 1:572506/73‑1 (MQ=255) taaataaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTc > 2:390838/1‑91 (MQ=255) taaaCGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTccccc < 1:429820/91‑1 (MQ=255) aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAgct < 1:777335/50‑1 (MQ=255) aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCACGGGGAATCATTCCCCCt > 1:369436/1‑91 (MQ=255) aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGaa < 2:253971/79‑1 (MQ=255) aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGaa > 1:253971/1‑79 (MQ=255) aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTg > 2:674836/1‑91 (MQ=255) aacgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTg > 1:786906/1‑91 (MQ=255) acgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑gg > 1:213476/1‑75 (MQ=255) acgaacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑gg < 2:213476/75‑1 (MQ=255) aacgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGtctc < 1:287741/91‑1 (MQ=255) acgCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTcc > 1:325863/1‑91 (MQ=255) gTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCAt > 2:478257/1‑91 (MQ=255) aTGTAGAACAAAATGATTCAGCGAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTAccg > 1:682201/1‑91 (MQ=255) aTGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTAccg > 1:682295/1‑91 (MQ=255) aTGTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTAccg > 2:49570/1‑91 (MQ=255) gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc > 2:846158/1‑91 (MQ=255) gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc > 2:43714/1‑91 (MQ=255) gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc > 2:762257/1‑91 (MQ=255) gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc > 2:371224/1‑91 (MQ=255) gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc > 1:746723/1‑91 (MQ=255) gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc > 2:205331/1‑91 (MQ=255) gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc > 1:800702/1‑91 (MQ=255) gTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACAATTTTTAccgcc > 1:759685/1‑91 (MQ=255) gTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCgt > 1:728020/1‑47 (MQ=255) gTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCgt < 2:728020/47‑1 (MQ=255) gTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATc < 2:806221/61‑1 (MQ=255) gTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATc > 1:806221/1‑61 (MQ=255) gTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc > 2:44781/1‑91 (MQ=255) gTAGAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTAccgcc > 2:574500/1‑91 (MQ=255) tAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTAccgccg < 2:786906/91‑1 (MQ=255) gAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGAc > 2:729640/1‑91 (MQ=255) gAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGAc > 2:521299/1‑91 (MQ=255) gAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGAc > 1:380009/1‑91 (MQ=255) gAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGAc > 2:516175/1‑91 (MQ=255) gAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTc < 2:681551/39‑1 (MQ=39) gAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTc > 1:681551/1‑39 (MQ=39) gAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGAc > 2:56380/1‑91 (MQ=255) gAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGAc > 2:195356/1‑91 (MQ=255) gAACAAAATGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGAc > 1:855802/1‑91 (MQ=255) aaaTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCa > 1:159906/1‑91 (MQ=255) aaaTGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTcca > 1:621232/1‑69 (MQ=255) aaaTGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTcca < 2:621232/69‑1 (MQ=255) aaaTGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCa > 1:863815/1‑91 (MQ=255) aaaTGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCa > 1:718868/1‑91 (MQ=255) aaaTGATTCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCa > 2:647950/1‑91 (MQ=255) aTGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCGCCATTTTTACCGCCGACGGCCAGc > 1:320230/1‑91 (MQ=255) aTTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCg > 1:757476/1‑91 (MQ=255) ttCAGCAAGGAGAGCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGa < 1:478257/91‑1 (MQ=255) cAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATa < 1:792446/91‑1 (MQ=255) cAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATa < 1:583448/91‑1 (MQ=255) aGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCg > 1:648494/1‑91 (MQ=255) aTCTAATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCAc < 2:453255/91‑1 (MQ=255) tctcATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCAcc > 1:246941/1‑91 (MQ=255) | ACCTCTCCTTGTTCTGCGTCATAGTATGATCGTTAAATAAACGAACGCTGTTCTATAATGTAGAACAAAATGATTCAGCAAGGAGATCTCATGCCGCAGTCCGCGTTGTTCAC‑GGGAATCATTCCCCCTGTCTCCACCATTTTTACCGCCGACGGCCAGCTCGATAAGCCGGGCACC > NZ_CP009273/277899‑278075 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |