Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,285,853 | 0 | G | T | 44.0% | 26.6 / 13.9 | 25 | A126S (GCC→TCC) | agaD | N‑acetylgalactosamine‑specific enzyme IID component of PTS |
Reads supporting (aligned to +/- strand): ref base G (5/9); new base T (11/0); total (16/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.01e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTT > NC_000913/3285704‑3286003 | ggTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgc < 2:222056/151‑1 (MQ=255) attaATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgccagc < 1:654770/151‑1 (MQ=255) attaATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgccag < 2:301756/150‑1 (MQ=255) ttaATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgccagc < 2:218634/150‑1 (MQ=255) aaTTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTgccagcca < 2:1007735/150‑1 (MQ=255) aaTTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGccagccag < 1:601560/151‑1 (MQ=255) aTTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAgg < 1:460738/151‑1 (MQ=255) ttCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTGCTCATCTTTTTCCAGCCCggg > 1:572658/1‑150 (MQ=255) aagaaAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAGTTTTGTCTTCCTTTgtcacaccggggaa > 2:235771/1‑135 (MQ=255) gagaAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTTTTCCTCATTTTTCACCCCGGGGAAACCGCCggggg > 2:242465/1‑146 (MQ=255) aaCCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGAATTGCTTATCATTTTCCCGCCCGGGGAAACTGCTGGGGGCGGtt > 2:122778/1‑147 (MQ=255) ccGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTAttt < 1:747818/151‑1 (MQ=255) cgACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTTCTCCTCATTTTTCCGCCCGGGGAAACTGCTGGGGGCGGTTCTttttttt > 1:538596/1‑150 (MQ=255) ggCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTTTCCTTATTTTTCCGCCCGGGGAAACTCCTGGGGGGGGTTTTTTTTTTTTCCGTTTTCCTCCtttttttttt > 1:852333/1‑150 (MQ=255) ggCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTc < 1:494864/150‑1 (MQ=255) ccAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCt < 1:926713/151‑1 (MQ=255) gCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAATTTGTCCTTCTTTTTCCGCCCGGGGAAACC‑GCTGGGGGCGGGTTTTTTTTTTTCCGTTTTTCTGCTTTTTTTTTTTCTTcgcggccggcgggccc > 2:368426/1‑139 (MQ=255) gATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGAATTTTTCCTCATTTTTCCGCCCGGGGAAACTGCTGGGGGCGGTTCTTTTTTTTTCCGTTTTCCTGCTTTTCTTTTTTTTTCGCGGCGGCTGGGccccc > 1:882472/1‑146 (MQ=255) cgatgcgatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTTCTCATCATTTTTCAACCAGGGGAACCCGCTGGGGCCGGTTCTTTTTTTTTCCGCTTTCCTTCTTTTTTTTTTTCTTCGCGGCCGCTGGACCCCCCTCgggttt > 1:229817/1‑147 (MQ=255) gatTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAAATTTGTCTTCATTTGCCACCCAGGGGAAACCGCTCGGGCGGGTTCTTTTTTTTTCCGGTTACCTTCTTATTTTTTTTCTTCCGGTCGGGTTGGCCCCCGCCgggtttt > 1:968460/1‑141 (MQ=255) ttttCTGGTTTACTTTGTTGCCGATTATGGC‑GGGGATTTTCTCATCCTTTTGCCGCCAGGGAAACCCGCTTGGGGCGGTTCTTTTTTTTTCCGGTTTCCTGCTTTTTTTTTTTCTGCGCGGCTGGTGGGCCCCCGTCGGGTTTTCAGTAgg > 2:932901/1‑151 (MQ=255) tttCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTTCTTATCATTTTCCAACCAGGGGAAACTGCTGGGGGCGGTTTTTTTTTTTTCCGTTTTCCTGCTTTTCTTTTTTCTTCGCGTCGGCTGGGCCCCCCTCGGTtttttcgt > 1:284587/1‑143 (MQ=255) gTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAg < 2:1020376/151‑1 (MQ=255) tttACTTTGTTGCCGATTATGGC‑GGGAATTTTCTCCTCATTTTCCCGCCAGGGGAAACTGCTGGGGGCGATTCTTTTTTTTTCCGTTTTCCTGCTTTTCTTTTTTCTGCGCGGCGGCTGGGGCCCCGTCGGTTTTTTTGTCGGCGTGAAgg > 1:56235/1‑151 (MQ=255) aCTTTGTTGCCGATTATGGC‑GGGAATTTGCTCCTCAATTGCCAGCCAGGGGAAACCGCTGGGGGCGCTTCTTTTTTTTTCCGTTTACCTGCTTTTCTTTTTTCTGCGCGTCGGCTGGGCCCCCGCCGGGTATTTAGTCGGCggggaagcg > 1:296225/1‑142 (MQ=255) aTGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCga < 1:687132/150‑1 (MQ=255) c‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGaa < 1:923297/151‑1 (MQ=255) tGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAg < 2:95816/150‑1 (MQ=255) tcatTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTg > 2:266419/1‑151 (MQ=255) tcatTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTg > 2:911251/1‑151 (MQ=255) atTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGcc > 2:909249/1‑151 (MQ=255) gccagccaGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGtt < 1:301355/151‑1 (MQ=255) | GGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGC‑GGGAATTTGCTCATCATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCTATTTTTCGCCGTTTACCTGCTTATCTTTTTCCTGCGCGTCGGCTGGACCCACGTCGGTTATTCAGTCGGCGTGAAGGCGATCGATAAAGTGCGAGAGAACTCGCAGATGATTGCCCGTT > NC_000913/3285704‑3286003 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |