Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,285,853 | 0 | G | T | 42.9% | 62.8 / 24.1 | 42 | A126S (GCC→TCC) | agaD | N‑acetylgalactosamine‑specific enzyme IID component of PTS |
Reads supporting (aligned to +/- strand): ref base G (2/22); new base T (18/0); total (20/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.37e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01 |
GGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGT > NC_000913/3285703‑3286002 | gggTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTg < 2:703824/151‑1 (MQ=255) taATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgccag‑cc < 1:626746/150‑1 (MQ=255) aaTTTCGATGGGAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgccag‑cca < 2:263700/150‑1 (MQ=255) aaTTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgccag‑cca < 2:721008/150‑1 (MQ=255) ggAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACct < 2:942698/150‑1 (MQ=255) aaaaCCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTc < 1:926179/150‑1 (MQ=255) aCCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTTTCCTCCTTTTTCCG‑CCAGGGGAAACTGCTGGGGGCGGTTCtttt > 2:865970/1‑151 (MQ=255) gcgACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTAttttt < 2:116135/151‑1 (MQ=255) cacCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAAATTTTTTACCATTTTCCCG‑CCCGGGGAACCCGCTTGGGCCGGTTCTAtttttttcc > 2:892969/1‑147 (MQ=255) acCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGTCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCg < 1:339522/151‑1 (MQ=255) cATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGGATTTG‑TCCTCCTTTTTCAGCCCAGGAAAACCTCTTGGGGCGGTTCCTtttttttccgtt > 1:229758/1‑145 (MQ=255) cATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTTGTCCTCCTTTTTCAG‑CCACGGGAAACTGCTGGGGGCGGTTCTtttttt > 1:236736/1‑144 (MQ=255) aTTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTGCTCCTCATTTTCCCG‑CCCGGGGAACCTGCTGGGGGCGGTTTTTTTTTTTGCCGttt > 1:369461/1‑151 (MQ=255) aTTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGttt < 2:716319/151‑1 (MQ=255) tAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCTGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTa < 2:385384/150‑1 (MQ=255) aaaGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTTTTCATCATTTTTCCC‑CCCGGGGAACCCGCTTGGGGCGGTTTTTTTTTTTTCCGTTTTcc > 2:659540/1‑151 (MQ=255) aaGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTCTCATCATTTTTCAG‑CCAGGGGAAACCGCTGGGGGCGGTTTTATTTTTTTCCGttt > 2:14858/1‑147 (MQ=255) aCTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTg < 1:647323/151‑1 (MQ=255) cTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGAATTGGTTCTCCTTTTCCCG‑CCAGGGGAAACCGCTGGGGGCGGTTCTTTTTTTTTCCGCTTTCCTTCttttttttttt > 1:47450/1‑146 (MQ=255) cTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCt < 1:129706/149‑1 (MQ=255) ttGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTTCCCG‑CCAGGGGAAACCGCTGGGGGCGGTTCTATTTTTTTCCGGTTTCCTGCTTTTCttttt > 2:713886/1‑142 (MQ=255) cAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTg < 1:46017/151‑1 (MQ=255) tCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAAATTTTTTCTCCTTTTCCCG‑CCCGGGAAAACCGCTCGGGGCGGTTCTATTTTTTCCCGTTTTCCTGCCTTTTTTTTTTCTTCGCCccggccgg > 2:570641/1‑149 (MQ=255) cGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTCTCTTCATTTTTCAG‑CCAGGGGAACCTGCTGGGGGCGGTTTTATTTTTTTCCGCTTTCCTGCTTATTTTTTTTCTTCGCGTCGGCTGGGcc > 2:758338/1‑151 (MQ=255) gATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTCTCATCCTTTTTCAG‑CCCGGGAAAACTCCTGGGGGCGGTTCTTTTTTTTTCCCTTTTCCTCCTTATTTTTTTTCTTCGCGCCGGCCGGGCCCCCgcc > 2:254693/1‑149 (MQ=255) ttGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGCATTTTGTCCTCATTTTTCAG‑CCAGGGGAAACTTCTGGGGGCGGGTCTTTTTTTTTCCGTTTTCCTGCTTTTTTTTTTTCTTCGCGTCGGCCgggacccc > 1:734663/1‑140 (MQ=255) tGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCgg < 2:244462/150‑1 (MQ=255) ggCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGtt < 1:484428/151‑1 (MQ=255) atgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTc < 1:82474/151‑1 (MQ=255) gcgatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCCTTATTTTTCAC‑CCCGGGGAACCCCCTGGGGGCGGTTTTCTTTTTTTCCGTTTTCCTGCTTTTTTTTTTTCTTCGCGTCGGCTGGGCCCCCGCCGGTTAtt > 1:331237/1‑148 (MQ=255) cgatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTTTCATCATTTTCCAG‑CCAGGGGAAACTGCTGGGGCCGGTTCTATTTTTTTCCGCTTTCCTGCTTTTCTTTTTTCTTCGCGTCGGCTGGGCCCCCCCCGGGTATTTAGt > 2:58560/1‑151 (MQ=255) cgatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGt < 2:510425/151‑1 (MQ=255) atTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGAATTTTTCCTCATTTTCCAG‑CCCGGGGAAACCGCTGGGGGCGATTCTTTTTTTTGCCGTTTTTCTGCTTTTTTTTTTTCTTCGCGGCGGGTGGGCCCCCCGCGGGTATTTAGt > 2:475406/1‑149 (MQ=255) atTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCg < 2:333267/151‑1 (MQ=255) tCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCg < 2:684845/150‑1 (MQ=255) tCTGGTTTACTTTGGTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCAGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCg < 1:681319/150‑1 (MQ=255) cTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTCTCCTCATTTTCCCG‑CCCGGGGAAACTTCTGGGGGCGATTCTTTTTTTTTCCGTTTTCCTGCTTATCTTTTTTCTTCGCGTCGGCCGGGCCCCCCTCGGTTTTTTAGTCGGCGTg > 2:191663/1‑151 (MQ=255) ggTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGt < 1:567686/148‑1 (MQ=255) aCTTTGTTGCCGATTATGGC‑GGGAATTTTCTCATTATTTTCCAC‑CCAGGGGAAACTGCTGGGGGCGGGTTTTTTTTTTTCCGCTTTTCTGCTTTTTTTTTTTCTGCGCGTCGGCTGGGCCCACGCCGGGTTTTTTGTCGGCGTGAAGGcgg > 1:150338/1‑150 (MQ=255) tgttgCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGcgatcg < 2:421165/150‑1 (MQ=255) gCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATaa < 1:87840/150‑1 (MQ=255) ccGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAg < 1:808554/151‑1 (MQ=255) gATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTg < 2:479704/151‑1 (MQ=255) ggC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCgaga < 1:892940/150‑1 (MQ=255) ttGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGc > 1:583960/1‑149 (MQ=255) tGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGc > 1:561169/1‑148 (MQ=255) gccag‑ccaGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGt < 2:962582/150‑1 (MQ=255) | GGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAG‑CCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGT > NC_000913/3285703‑3286002 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |