Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,570,447 | 0 | A | T | 40.0% | 85.3 / 19.2 | 55 | L278Q (CTA→CAA) | hsdS | type I restriction‑modification system specificity subunit |
Reads supporting (aligned to +/- strand): ref base A (2/31); new base T (22/0); total (24/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.12e-13 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCA > NZ_CP009273/4570374‑4570536 | aaCACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACa < 2:1016626/90‑1 (MQ=255) ccAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGGTTTTTTTACTTTTTTAAAAttat > 1:776245/1‑88 (MQ=255) ccAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGGTTTTTTTACTTTTTTAAAATTCt > 1:1318809/1‑90 (MQ=255) aacaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTTAAAATTCTaa > 1:15592/1‑89 (MQ=255) acaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGa < 2:300185/90‑1 (MQ=255) acaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGa < 2:923134/90‑1 (MQ=255) aaaTTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGaaa < 2:477568/90‑1 (MQ=255) aaaCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTTTTTACTTTTTTAAAATTTTAGAAAACGAata > 2:890215/1‑90 (MQ=255) aaCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTTAAAAATCTTAAAAACGAattt > 1:991898/1‑88 (MQ=255) ttCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGt < 2:1129635/90‑1 (MQ=255) ccGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTTAAAATTCTAGAAACCCAATATTGttt > 1:197971/1‑90 (MQ=255) ccGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGttt < 1:323687/90‑1 (MQ=255) ccGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGttt < 1:70323/90‑1 (MQ=255) ttATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGa < 2:218951/90‑1 (MQ=255) tataGCGAGTAAATAAAAGATTTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGAt < 1:1355872/90‑1 (MQ=255) tataGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGAt < 2:96726/90‑1 (MQ=255) tataGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGAt < 2:16211/90‑1 (MQ=255) cGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACa < 2:1415963/90‑1 (MQ=255) gAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTTAACATTTTTGAAAACGAATATCGTTTTTATTTACAt > 1:1350353/1‑90 (MQ=255) gAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTTAAAATTTTAGAAAACCAAAATCGTTTTTATTTACAt > 2:980556/1‑90 (MQ=255) gAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGGTTTGTTTACTTTTTTAAAATTTTAGAAAACGAATATTGTTTTTTTTTACAt > 1:839189/1‑90 (MQ=255) gTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTTTTTACTTTTTTAAAATTCTAGAAAACGAAAATCGTTTTTATTTTCATgg > 2:674540/1‑90 (MQ=255) aaaaGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCa < 1:126556/90‑1 (MQ=255) aaaaGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCa < 1:1259022/90‑1 (MQ=255) aaaGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTTTTAACATTTTAGAAAACGAAAATCGTTTTTATTTACATGGCCAGCaa > 2:860249/1‑89 (MQ=255) aaGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTCACTTTTTTAAAATTTTTGAAAACCAAAATCGTTTTTTTTTACATGGCCAGCAcc > 2:1004310/1‑89 (MQ=255) aGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTGAACATTCTAGAAAACAAATATCGTTTTTATTTACATGGGCAGCACGt > 2:1320283/1‑90 (MQ=255) gATCTCCATCTTGTAATTTGTGGCGGTTTTTTTTACTTTTTTTAAATTTTTTAAAACCAATAATGTTTTTTTTTTCATTGGCAGCACGTa > 1:358297/1‑90 (MQ=255) gATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTTAAAATTCTAGAAACCGAATATCGTTTTTTTTTACATTGGCAGCACGTa > 1:152231/1‑90 (MQ=255) aTCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTAc < 2:839189/90‑1 (MQ=255) aTCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTAc < 1:969186/90‑1 (MQ=255) tctcCATCTTGTAATTTGTGGCGGTTTTGTTCACCTTTTTAACAATTTAAAAAACGAAAATCGTTTTTATTTTCATGGCCAGGACGTACa > 2:671132/1‑90 (MQ=255) cATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAAc < 2:977786/90‑1 (MQ=255) aTCTTGTAATTTGTGGCGGTTTTGTTCACTTTTTTAACAATCTAAAAACCCAAAATCGTTTTTATTTACATTGCCAGCACGGACAGaaaa > 2:504155/1‑88 (MQ=255) cTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTaa < 1:890215/90‑1 (MQ=255) ttGTAATTTGTGGCGGTTTTGTTTTCTTTTTTTAAATTTTTTAAAACGAATATCGTTTTTTTCTTCATTGGCACCACGGACAAGACTaaa > 1:1114083/1‑89 (MQ=255) ttGTAATTTGTGGCGGTTTTGTTTACTTTTTTAAAATTTTAGAAAACGAATATCGTTTTTTTTTACATTGCCAGCACCTACAGAaataat > 2:1301245/1‑90 (MQ=255) ttGTAATTTGTGGCGGTTTTGTTCACTTTCTGAACATTCTAGAAAACCAAAATCGTTTTTATCTACATTGCCAGCACCTACAGAACTaaa > 2:618384/1‑89 (MQ=255) tGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATg < 2:1340580/90‑1 (MQ=255) gTAATTTGTGGCGGTTTTGTTTACTTTTTTAAAATTCTTGAAAACCAATATTGTTTTTTTTTTCATTGGCAGCCCGGACAAAAATAATGc > 2:794872/1‑90 (MQ=255) gTAATTTGTGGCGGGTTTTTTTACTTTCTTGATATTTTTTAAAACGAATATTGTTTTTTTTTTCATTGCCAGCACGGACAAAAATAATGc > 1:1186274/1‑90 (MQ=255) tAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGTCAGCACGTACAGAACTAATGCg < 2:1484043/90‑1 (MQ=255) aaTTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCgt < 2:79979/90‑1 (MQ=255) aaTTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCgt < 1:1366301/90‑1 (MQ=255) aaTTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCgt < 1:368364/90‑1 (MQ=255) aaTTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCgt < 1:674540/90‑1 (MQ=255) aaTTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCgt < 2:628647/90‑1 (MQ=255) aaTTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCgt < 2:358297/90‑1 (MQ=255) aTTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCgta < 2:46792/90‑1 (MQ=255) ttGTGGCGGTTTTGTTCACTTTCTTAACATTCTTGAAAACGAATATCGTTTTTTTTTACATGGCCAGCACGTACAGAACTAATGCgtagt > 2:713928/1‑90 (MQ=255) gtgGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTAt < 2:1114083/90‑1 (MQ=255) gtgGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTAt < 2:1191700/90‑1 (MQ=255) gCGGTTTAGTTCACTTTCTTAAAATTTTTGAAAACCAATATCGTTTTGATTTACATTGCCAGCACCTACAGAAAAAAAGCGGAGGATTgt > 1:369888/1‑89 (MQ=255) cGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGa < 2:15592/90‑1 (MQ=255) ggTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGAt < 2:565774/90‑1 (MQ=255) aGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCa > 2:956770/1‑90 (MQ=255) | AACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCA > NZ_CP009273/4570374‑4570536 |
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: ‑ |