Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,825,061 | 0 | G | A | 23.8% | 158.2 / 26.2 | 63 | intergenic (‑124/+79) | ves/spy | HutD family protein Ves/ATP‑independent periplasmic chaperone |
Reads supporting (aligned to +/- strand): ref base G (45/3); new base A (0/15); total (45/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.68e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAC > NC_000913/1825013‑1825110 | aCGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAG‑AAGagaaaag > 1:3072753/1‑49 (MQ=255) aCGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAG‑AAGagaaaag < 1:1744217/49‑1 (MQ=255) aCGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAG‑AAGagaaaag > 1:1407855/1‑49 (MQ=255) cGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAG‑AAGagaaaagaa < 1:3602821/50‑1 (MQ=255) cGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAG‑AAGagaaaagaa < 1:1425126/50‑1 (MQ=255) gTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAG‑AAGagaaaagaaa < 1:2488155/50‑1 (MQ=255) tAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAAc > 1:5971626/1‑50 (MQ=255) aGTGGTGATCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAAcc > 1:5849178/1‑50 (MQ=255) aGTGGTGATCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAAcc > 1:5861514/1‑50 (MQ=255) gtggtgATCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAAccg < 1:2034605/50‑1 (MQ=255) gtggtgATCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAAccg > 1:467920/1‑50 (MQ=255) gtggtgATCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAAccg > 1:1343034/1‑50 (MQ=255) gtggtgATCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAAcc < 1:5716920/49‑1 (MQ=255) tggtgATCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAAccgc > 1:4235376/1‑50 (MQ=255) ggtgATCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAAccgcc > 1:5541596/1‑50 (MQ=255) gtgATCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAAccgcc < 1:6351689/49‑1 (MQ=255) tgATCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGa > 1:1512809/1‑50 (MQ=255) tgATCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGa > 1:3954652/1‑50 (MQ=255) gATCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGAt > 1:4048117/1‑50 (MQ=255) tCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATcc > 1:3769726/1‑50 (MQ=255) cAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCt > 1:5702479/1‑50 (MQ=255) aGGTTTTTTTTACAATTGTTGGAG‑AAAAAAAAAAAAAACCGCCGATCCTg < 1:5019992/50‑1 (MQ=255) tttttttACAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCa > 1:590540/1‑50 (MQ=255) tttttttACAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCa > 1:364020/1‑50 (MQ=255) tttttttACAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCa > 1:14243/1‑50 (MQ=255) tttttACAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCAcc > 1:6389997/1‑50 (MQ=255) ttttACAGTTTTTTGAG‑AAGAAAAAAAAAAACCGCCGATCCTGTCCAcc < 1:1669895/49‑1 (MQ=25) ttttACAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCg > 1:3005405/1‑50 (MQ=255) tttACAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCg > 1:3471061/1‑49 (MQ=255) ttACAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGCTCCTGTCCACCGCa > 1:2824989/1‑50 (MQ=255) ttACAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCa > 1:6439696/1‑50 (MQ=255) aCAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAt > 1:6685537/1‑49 (MQ=255) aCAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAt > 1:675447/1‑49 (MQ=255) cAACTTTTTGAA‑AAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTa < 1:4099105/50‑1 (MQ=255) cAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa > 1:700546/1‑50 (MQ=255) cAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa > 1:6224327/1‑50 (MQ=255) aaCTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTAc > 1:855914/1‑50 (MQ=255) aaCTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTAc > 1:5211610/1‑50 (MQ=255) aatGTTTGAAAAAGAAAAAAAAAAGCCGCCGATCCTGTCCACCGCATTAc < 1:1430559/48‑1 (MQ=19) aCTGTTTGAT‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACt > 1:2490215/1‑50 (MQ=255) aCTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACt > 1:390962/1‑50 (MQ=255) aCTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACt > 1:2271376/1‑50 (MQ=255) aCTGTTTGAA‑AAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACt < 1:543669/50‑1 (MQ=255) cTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTg > 1:1010989/1‑50 (MQ=255) tGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTg > 1:990839/1‑49 (MQ=255) gTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc > 1:691270/1‑49 (MQ=255) gTTGGAG‑AAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGc < 1:5850960/49‑1 (MQ=255) tttGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc > 1:3553155/1‑48 (MQ=255) tttGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc > 1:4341414/1‑48 (MQ=255) ttGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAg > 1:2533830/1‑50 (MQ=255) ttGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAg > 1:4270596/1‑50 (MQ=255) tGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg > 1:6385311/1‑50 (MQ=255) tGAG‑AAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAg < 1:1685187/49‑1 (MQ=39) tGAA‑AAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAg < 1:205456/49‑1 (MQ=255) ggag‑aaggaaaAAAAAAGCCGCCGATCCTGTCCACCGCATTACTGCAAg < 1:5897295/48‑1 (MQ=37) gag‑aagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg > 1:1920393/1‑49 (MQ=255) ag‑aagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg > 1:1404587/1‑48 (MQ=255) ag‑aagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg > 1:342930/1‑48 (MQ=255) ag‑aagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAgg > 1:1208279/1‑48 (MQ=255) ag‑aaaaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa < 1:1310316/50‑1 (MQ=255) g‑aagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg > 1:3984487/1‑50 (MQ=255) aagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGt > 1:822402/1‑50 (MQ=255) aagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGt < 1:5200685/50‑1 (MQ=255) agagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg > 1:2118458/1‑50 (MQ=255) agagaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg < 1:4827191/50‑1 (MQ=255) agaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg < 1:5396616/50‑1 (MQ=255) agaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg < 1:4992822/50‑1 (MQ=255) agaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg < 1:5739688/50‑1 (MQ=255) gagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTg > 1:296856/1‑49 (MQ=255) agaaaagaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGa > 1:1117988/1‑50 (MQ=255) agaaaagaaaACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGa > 1:3091248/1‑50 (MQ=255) gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGa > 1:2663453/1‑49 (MQ=255) gaaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGa > 1:4608021/1‑49 (MQ=255) aaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc > 1:4708100/1‑49 (MQ=255) aaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa < 1:6253982/50‑1 (MQ=255) aaaaaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa < 1:1181549/50‑1 (MQ=255) aagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc > 1:2486069/1‑47 (MQ=255) aagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGAGGAc > 1:4477138/1‑47 (MQ=255) agaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGa > 1:3443928/1‑50 (MQ=255) gaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAc > 1:630698/1‑50 (MQ=255) gaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAc > 1:2296420/1‑50 (MQ=255) | ACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAG‑AAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGAC > NC_000913/1825013‑1825110 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |