| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_012967 | 3,011,193 | 0 | T | G | 36.4% | 49.0 / 15.4 | 44 | intergenic (+164/‑179) | ECB_RS14865/ECB_RS14870 | DUF4942 domain‑containing protein/IS110 family transposase |
| Reads supporting (aligned to +/- strand): ref base T (20/8); new base G (16/0); total (36/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.65e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TAGAAGCATTCTTCGATTCCGAGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTAAACACTCCCGAAT > NC_012967/3011129‑3011264 | tAGAAGCATTCTTCGATTCCGAGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGgttttttt > 2:4256532/1‑69 (MQ=255) aGAAGCATTCTTCGATTCCGAGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTTCtt > 2:1966843/1‑72 (MQ=255) aaGCATTCTTCGATTCCGAGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTTGtttt > 2:4879705/1‑72 (MQ=255) aaGCATTCTTCGATTCCGAGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTTCtttt > 2:1515533/1‑72 (MQ=255) aaGCATTCTTCGATTCCGAGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTTCtttt > 1:3982099/1‑72 (MQ=255) ttcttcGATTCCGAGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGTTTTTGCTTTTGGATc > 1:3396733/1‑72 (MQ=255) ttcttcGATTCCGAGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTTCTTTTTGATc > 1:5909161/1‑72 (MQ=255) ttttcGTTTCGGAGTCCGGCCCCACTATTTCTAAAAACCCCGCCCACAG‑GG‑G‑GGGTTTTTGCTTTTGGATC‑t < 1:5920491/70‑1 (MQ=255) tcttcGCTTCCGAGTCCGGCACCACTAATTCATCAGAACCCGCCCACAC‑GG‑C‑AGGTATTTGCTTTTGGATC‑t > 1:3887913/1‑72 (MQ=255) cttcGATTCCGAGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GGCC‑GGGTTTTTCTTTTTGGGTC‑t > 1:2948824/1‑72 (MQ=255) attcGATTCCGAGTCCGGCACCNCTAATTCTTAAGAACCCACCCACAA‑GG‑C‑GGGTTTTTGCTTTTGGATC‑tt < 1:2096130/71‑1 (MQ=255) gATTCCGAGTCCGGCACCACTAATTCTTAAGNACCCGCNCCCAA‑AG‑CGGGGTTTTTTTTTTTGGttt‑ttttt > 1:435049/1‑70 (MQ=255) gATTCCGAGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑CGGGGTTTTTCCTTTTGG‑TC‑TTttta > 2:1626245/1‑72 (MQ=255) gATTCCGAGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGTTTTTGCTTTTGGATC‑TTtata > 2:5099835/1‑72 (MQ=255) gATTCCGAGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑AG‑C‑GGGGTTTTTCTTTTGGGTC‑TTttta > 1:2307965/1‑72 (MQ=255) ttCCGAGTCCTGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTGCTTTTGGATC‑TTtata‑t‑a > 1:3442282/1‑72 (MQ=255) cGAGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTGCTTTTTGATC‑TTTATa‑t‑aaaa > 2:1954995/1‑72 (MQ=255) aGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTTCCTTTTGATC‑TTTTTa‑t‑aaaata > 2:211621/1‑72 (MQ=255) aGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTGCTTTTGGATC‑TTTATt‑t‑taaata > 2:5388181/1‑72 (MQ=255) aGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTGCTTTTGGATC‑TTTATa‑t‑caaata > 2:5816730/1‑72 (MQ=255) tCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑AG‑C‑GGGGTTTTTCTTTTTGGTC‑TTTTTATt‑aaaataa > 1:6872548/1‑72 (MQ=255) ccGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTGCTTTTTGGTC‑TTTTTA‑T‑AAAATAAAc > 1:5438858/1‑72 (MQ=255) caccacTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGTTTTTGCTTTTGGGTC‑TTTATA‑T‑AAAATAAAA‑tttg > 2:1906668/1‑70 (MQ=255) accacTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGTTTTTGCTTTTTGATC‑TTTATA‑T‑AAAATAAAC‑TTgtg > 1:1679799/1‑70 (MQ=255) ccacTAATTCAAAACACCCCGCCCACAG‑GG‑G‑GGGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGa < 1:7329068/72‑1 (MQ=255) cacTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTTCTTTTTGGTCTTTTTTT‑T‑AAAAAAAAA‑TTgggg > 2:4520452/1‑71 (MQ=255) cacTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTTCTTTTTGGTC‑TTTTTT‑TAAAAATAAAA‑TTgggg > 2:1408492/1‑71 (MQ=255) cTAATTCTTAAGAACCCGCCCACAA‑GG‑CGGGGTTTTTGTTTTTGG‑TC‑TTTTTA‑T‑TAAAAAAACATTggggga > 2:4963148/1‑69 (MQ=255) cTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTTTTTTTTGATC‑TTTATA‑T‑TAAATAAAC‑CTTGGAGCa > 2:3500199/1‑72 (MQ=255) cTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTTCTTTTTGATC‑TTTATA‑T‑AAAATAAAC‑TTTGGAGCa > 1:6894686/1‑72 (MQ=255) aaTTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGGTTTTGCTTTTGGATC‑TTTATA‑T‑TAAATAAAC‑TTGGGTGCATa > 2:3581645/1‑72 (MQ=255) cTTAAGAACCCGCCCACAA‑AG‑CGGGGGTTTTCTTTTTGGTTT‑TTTTTT‑T‑AAAAAAAAC‑TTGGGGGCATAACg > 2:3460704/1‑72 (MQ=255) ttAAGAACCCGCCCACAA‑GG‑CGGGGGTTTTGTTTTTGGTTT‑CTTTTT‑T‑AAAAAAAAA‑TTGGGGGAATaatac > 2:6874535/1‑69 (MQ=255) aacaCCCCCCCAAAAGGG‑C‑GGGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCAt < 2:6065644/69‑1 (MQ=38) aaCCCGCCCACAA‑GG‑C‑GGGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAg > 1:5707137/1‑72 (MQ=255) aaCCCGCCCACAA‑GG‑C‑GGGGTTTTTCTTTTTGATC‑TTTATA‑T‑AAAATAAAA‑TTGGGAGCATAACGCCATCAg > 2:7027111/1‑72 (MQ=255) cccccccAACAG‑GG‑G‑GGGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAgg < 1:3633747/71‑1 (MQ=255) aCCCGCCCACAA‑GG‑C‑GGGGTTTTGCTTTTTGATC‑TTTATA‑T‑AAAATAAAA‑TTGGGAGCATAACGCCCTCacg > 1:2005056/1‑70 (MQ=255) ccGCCAACAG‑GG‑C‑GAGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTc < 2:1655114/72‑1 (MQ=255) cccccAGAG‑GG‑G‑GGGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCt < 1:5633063/70‑1 (MQ=255) ccACAA‑GG‑C‑GGGTTTTTTCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTaa > 2:5991945/1‑72 (MQ=255) ccACAA‑GG‑C‑GGGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTaa > 2:1403651/1‑72 (MQ=255) cacaA‑GG‑C‑GGGTTTTTGCTTTTTGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTaaa > 2:3324883/1‑72 (MQ=255) cacaA‑GG‑C‑GGGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTaaa > 2:6991062/1‑72 (MQ=255) acaA‑GG‑C‑GGGTTTTTGCTTTTTGAAC‑TTTATA‑T‑AAAATAAAC‑TTTGGAGCATAACGCCATCAGGGCTTAAac > 2:1794172/1‑72 (MQ=255) caA‑GG‑C‑GGGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTAAaca < 2:7252229/72‑1 (MQ=255) aa‑GG‑C‑GGGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTAAacac > 2:5583586/1‑72 (MQ=255) aa‑GG‑C‑GGGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTAAacac > 1:6003494/1‑72 (MQ=255) aa‑GG‑C‑GGGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTAAacac > 1:3635461/1‑72 (MQ=255) aa‑GG‑C‑GGGTTTTTGCTTTTCGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCACGTCTTAAacac > 2:1513028/1‑72 (MQ=255) aa‑GG‑C‑GGGGTTTTTCTTTTGGGTC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTAAacac > 2:1313752/1‑72 (MQ=255) ggggggTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTAAACACTc < 1:3339571/69‑1 (MQ=255) c‑GGGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTAAACACTccc > 1:174794/1‑72 (MQ=255) c‑GGGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGGCTTAAACACTccc > 2:4573802/1‑72 (MQ=255) gggTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTAAACACTCCCg > 1:2316818/1‑72 (MQ=255) ggTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTAAACACNCCCGa > 2:2869566/1‑72 (MQ=255) gTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTAAACACTCCCGaa > 1:7297825/1‑72 (MQ=255) tttttGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTAAACACTCCCGAat < 2:6021368/72‑1 (MQ=255) tttttGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTAAACACTCCCGAat > 1:6575609/1‑72 (MQ=255) tttttGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTAAACACTCCCGAat > 1:4881715/1‑72 (MQ=255) tttttGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCANAACGCCATCAGGTCTTAAACACTCCCGAat > 2:456524/1‑72 (MQ=255) | TAGAAGCATTCTTCGATTCCGAGTCCGGCACCACTAATTCTTAAGAACCCGCCCACAA‑GG‑C‑GGGTTTTTGCTTTTGGATC‑TTTATA‑T‑AAAATAAAC‑TTGGGAGCATAACGCCATCAGGTCTTAAACACTCCCGAAT > NC_012967/3011129‑3011264 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |