Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,640,981 | 0 | C | T | 44.4% | 33.4 / 17.2 | 36 | F40F (TTC→TTT) | dtpB | dipeptide and tripeptide permease B |
Reads supporting (aligned to +/- strand): ref base C (3/17); new base T (16/0); total (19/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.35e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTT > NC_000913/3640836‑3641119 | aGGATGGGATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGTttatttct > 1:238633/1‑148 (MQ=255) gTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCTTTtaa > 2:431345/1‑149 (MQ=255) tGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCCGTTaacccct > 1:391832/1‑145 (MQ=255) ggTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTGCGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGAta < 1:127998/151‑2 (MQ=255) ggTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCAGTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGAtt < 1:164102/151‑1 (MQ=255) aataaataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAAAAACTTTGATtct > 2:373666/1‑149 (MQ=255) aaataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCa < 2:699114/150‑1 (MQ=255) aataTGACTACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTaaaaaaactgggttt > 1:296220/1‑132 (MQ=255) tatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGGGTTTTTTGTTAAACactctgggttct > 1:200692/1‑132 (MQ=255) cATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACtttt < 2:604009/151‑1 (MQ=255) gggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTg < 2:260408/151‑1 (MQ=255) gATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGGTAAACAACTTGGGTTCTCCCAAGAGGAGGCTTTTTTCACtttttg > 1:175440/1‑145 (MQ=255) aTGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCtt < 1:286354/151‑1 (MQ=255) ctgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCTTTTAAAAGCTTTGGATCTCGCAAGAGCAGGCTTTTTTCACTTTTTGTGGttttt > 2:289101/1‑150 (MQ=255) tgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTTAAAAGCCTGGATTCTTGCCAGAGGCGGCTTTTTTCACTTTTTGGGGtttt > 1:288142/1‑149 (MQ=255) gcagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTtgc < 1:447222/150‑1 (MQ=255) gcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTTAACAACCTTGGTTCTCGCGAGAGCAGGCTTTTTTTACTTTTTGTGCttttttt > 1:47358/1‑145 (MQ=255) caACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTgcgc < 2:472457/151‑1 (MQ=255) aaCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTTTTTTAAAAGCTTTGGTTTTTGCCAAAGCAGGGTTTTTTTTCTTTTTGGGGttttttt > 1:240606/1‑143 (MQ=255) cccATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTAt < 2:290799/151‑1 (MQ=255) tttCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAAAAGCTTGGGTTTTCGCCAGAGCAGGGTTTTTTCACTTTTTGTGGTTTTTTTTCGCTGGGCTATgggccc > 2:638141/1‑147 (MQ=255) ttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAACAACTTTGATTTTCGCAAGAGCAGGCTTTTTTCACTTTTTGTGCTTTTTCTGCGCTGGGCTATGGCCTCCttttcctt > 2:695997/1‑143 (MQ=255) ttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTCGTTAAACCGCTTGGGTTTTCGCAAGAGCAGGCTTTTTTTCCTTTTGGTGGTTTTTTTGTGCCGGTCTTTGGCCTCCttttccct > 1:378141/1‑143 (MQ=255) gAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTAAAAAGCCTGGGTTTTCCCAAGAGCAGGCCTTTTTCTCTTTTTGGGGTTTTGCTGCGCTGGGTTATGGCCTCATTTCCAtttgggg > 2:411833/1‑141 (MQ=255) aGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTAAAAAGCTTGGATTTTCGCAAGAGGAGGGTTTTTTCACTTTTGGTGCTTTTTCTGCGCTGGGCTATGGCCCCATTTTCATTTGGGGGTATGt > 2:258245/1‑151 (MQ=255) aGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTCTTAAAAAGCTTGGATTATCGCCAGAGCAGGCTTTTTTTCCCTTTTGTGCTTTTTCTGCGGTTGGCTTTGGCCTCCTTTTCCTTtgcgg > 1:420656/1‑145 (MQ=255) aGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGt < 1:143992/151‑1 (MQ=255) tATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTCTTTAACCCCATGGATTCTACCCAAAGGCGGCTTTTTGGCCTTTTTGGGCTTTTTTTGCGGCGGCTTATGGGCCCCTTTCCCTTggcgg > 2:579616/1‑142 (MQ=255) gggAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCg < 1:431887/150‑1 (MQ=255) actacGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTgggg < 1:713379/149‑1 (MQ=255) tacGGCGTGCAGGGCGTACTGGC‑GGTTTTTTTTGTTTAACAGCCTTGATTTTCGCAAGAGCAGGCTTTTTTTACTTTTTGTGGTTTTTCTGCGCTGGTCTTTTGCCTCATTTTCCTTTGGGGGTATGTTGGCGgccccccccttggggac > 2:535607/1‑139 (MQ=255) acGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCaa < 2:593432/151‑1 (MQ=255) acGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCa < 2:378360/150‑1 (MQ=255) ggCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTTACTTTTTGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTTGCGGCTATGTCGGCGACCCCCCGCTGGGGACCaaa > 2:79106/1‑150 (MQ=255) gCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCAcc < 2:438998/151‑1 (MQ=255) gCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGt < 2:93425/151‑1 (MQ=255) aCTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCATGAGCAGGCTTGTTTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCtt < 1:535825/151‑1 (MQ=255) | AGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTT > NC_000913/3640836‑3641119 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |