Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 711,081 | 0 | A | C | 52.9% | 22.5 / 31.6 | 34 | T72P (ACC→CCC) | BW25113_RS03560 | lipoprotein |
Reads supporting (aligned to +/- strand): ref base A (15/1); new base C (0/18); total (15/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.62e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.89e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGAGGATTAACGGGTAGCATTGCGGTAAATATTGCACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTT > NZ_CP009273/710996‑711162 | aagagGATTAACGGGTAGCATGGCGGTAAATATTCCCCTTTTCGGGTGAAAGAAAAATTGCGGTGGTAATACCTTCCCTGATAATCCccg < 2:1749550/89‑1 (MQ=255) agagGATTAACGGGTAGCATTGCGGTAAATATTGCACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACccgc > 1:506937/1‑90 (MQ=255) gagGATTAACGGGTAGCATTGCGGTAAATATTGCACCTTTCCGGTGTAAGACAAATTTCCGTGGTGATAACTTACCTGATAATACccgcc > 2:1428967/1‑90 (MQ=255) gagGATTAACGGGTAGCATTGCGGTAAATATTGCACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACccgcc > 1:14300/1‑90 (MQ=255) gATAAACGGGTAGCATTGGGGTAAATATTCCCCTTTTCCGGTGTAAGAAAAATTGCGGGGGTGATAACTTCCCTGAAAATCCCCGCCGCt < 2:1233377/90‑1 (MQ=255) aaCGGGTAGCATTGCGGTAAATATTGCACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCa > 2:729654/1‑90 (MQ=255) aggggTAGCATTGCGGTAAATTTTGCCCCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTCCCTGATAATCCCCGCCGCTGGCAg < 1:619092/88‑1 (MQ=255) tGCGGTAAATATTGCACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCgg > 2:1653896/1‑90 (MQ=255) gCGGTAAATATTGCACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCggg > 2:1471950/1‑90 (MQ=255) gCGGTAAATATTCCCCCTTTCCGGTGAAAGAAAAATTGCCGGGGTGATAACTTCCCTGATAATCCCCGCCGCTGGCAGTTGTCAGGCggg < 1:691137/90‑1 (MQ=255) ttGCACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAg > 1:144087/1‑90 (MQ=255) cACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAac > 1:384426/1‑90 (MQ=255) tttCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAacttac > 1:735208/1‑90 (MQ=255) tttCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAacttac > 1:243821/1‑90 (MQ=255) tttCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAacttac > 2:557968/1‑90 (MQ=255) ttcccGGTGTAAGAAAAATTGCGGTGGTAATAACTTCCCTGATAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAacttac < 2:384426/87‑1 (MQ=255) ccGGTGAAAAAAAATTTGCGGGGGTAAAAACTTCCCTAAAAATCCCCGCCGCTGGAAGTTGTCAGGCGGGAACGAAAAAGAActtactta < 2:310438/90‑1 (MQ=255) ggTGAAAGAAAAATTGCGGGGGTAAAAACTTACCTAATAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATc < 2:142594/90‑1 (MQ=255) gtAAAAAAATTTCCCGTGGTAATAATTTACCTGATAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTc < 1:288077/90‑1 (MQ=255) gaaaGACAATTTGCCGGGGTGATACCTTCCCTAATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTc < 1:523018/88‑1 (MQ=255) gaaaGAAAATTTCCCGTGGTGATAATTTCCCTAAAAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTc < 1:25661/88‑1 (MQ=255) gaaaGAAAAATTGCCGTGGTGATAATTTACCTGATAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTc < 1:934311/88‑1 (MQ=255) aaGAAAATTTGCGGTGGTAATAACTTACCTAATAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCat < 2:477261/90‑1 (MQ=255) aaaaaTTGCGGTGGTAATAACTTACCTAATAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCataat < 1:1471950/88‑1 (MQ=255) aCAATTTGCGGTGGTAATAACTTACCTAATAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCataat < 2:97306/90‑1 (MQ=255) aaaTTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAatat > 2:1374545/1‑90 (MQ=255) aaTTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATg > 1:97306/1‑90 (MQ=255) aaTTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATg > 1:477261/1‑90 (MQ=255) tGCGGTGGTGATAATTTACCTGAAAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCa < 1:1574017/90‑1 (MQ=255) tggtgATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGcaca > 2:1225224/1‑90 (MQ=255) aaaaCTTACCTAAAAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGaaa < 2:243821/88‑1 (MQ=255) ttACCTGATAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTc < 2:1758186/90‑1 (MQ=255) tcccTAATAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCg < 1:1243654/88‑1 (MQ=255) cTGATAATCCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGtt < 2:144087/90‑1 (MQ=255) | CAGAGGATTAACGGGTAGCATTGCGGTAAATATTGCACCTTTCCGGTGTAAGACAAATTGCCGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTT > NZ_CP009273/710996‑711162 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |