Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,019,596 | 0 | G | A | 48.6% | 24.7 / 18.1 | 39 | V146I (GTT→ATT) | ssnA | putative chlorohydrolase/aminohydrolase |
Reads supporting (aligned to +/- strand): ref base G (5/14); new base A (17/1); total (24/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.91e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑CCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGG > NC_000913/3019452‑3019737 | cTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCC‑CAATTTGCGA‑CG‑CCCTTTTTTAAATTGGc > 2:438157/1‑150 (MQ=255) tgatgaAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTg < 1:511265/151‑1 (MQ=255) ctctATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTGGCCCTCCC‑CACTGTCCGA‑AGCCCTTTTTAAAATTTGG‑CTGCGCGCGGTGAcc > 2:472941/1‑150 (MQ=255) ctctATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGGTCCCCA‑CCCACCGCGC‑CG‑CCTTTTTTTAAGTTGGGCTGCGCGCGCTGAcc > 1:12176/1‑150 (MQ=255) ctctATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGGCCTCCC‑CCTTCCGCGC‑CG‑CACTTTTTTAAATTTGCCTGggcgcggagac > 1:329092/1‑144 (MQ=255) ctctATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATTGCGCGG‑CG‑CCTTTTTTAAAATTGGCCTGCGCGCGGTGAcc > 2:59812/1‑150 (MQ=255) ctctATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCC‑CATTTCGCGA‑CG‑CATTTTTAAAAATTGGCCTGCGCGCGGTGAcc > 2:680566/1‑150 (MQ=255) ctctATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCA‑CCATGCGCGG‑CG‑CACTTTTTAAAATTGGCAGgcgcggcgtgagc > 1:674770/1‑142 (MQ=255) ctATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGGTCTCCC‑CCCTTCGCGA‑CG‑CCTTTTTTTAAGTTTGCCTGCGCGCGGTGAc > 1:612134/1‑147 (MQ=255) ctATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCTCTCCC‑CCTTTCGCGA‑CC‑CCTTTTTTAAAATTTGCCTGCGCGCGGTggc > 2:220498/1‑145 (MQ=255) ctATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCGCCCTCCC‑CATTTCGCGA‑CG‑CCTTTTTTAAAATTTGGCTGCGCGCGGTGAc > 2:412933/1‑147 (MQ=255) aTTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCC‑CAATTCGCGA‑CG‑CATTTTTTTAAATTGGCCTGCGCGCGGTGGCCTGttt > 2:657537/1‑151 (MQ=255) aTTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGttt < 1:567044/151‑1 (MQ=255) aTTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGttt < 2:400982/151‑1 (MQ=255) aTTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGGCCCTCTCCA‑CACTTCGCGG‑CG‑CCTTTTTTAAAATTTGGCTCCGCGCGGTGACCTGttt > 2:373858/1‑151 (MQ=255) gACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAacta < 2:301093/150‑1 (MQ=255) aTTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑cc < 2:506636/151‑1 (MQ=255) gATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCC‑CACTTCCCGA‑CC‑CATTTTTTTAAATTTGCCTGCGCGCGGTGGCCTGTTTTTTAAATACTGT‑CCGTAAAAAAGGCAACaaa > 1:641965/1‑151 (MQ=255) gATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CACTGCGCGC‑CG‑CCTTTTTTAAAATTGGGCCGCGCGCGATGACCTCTTTTTAAACTACTGA‑CCCTAACAAACGCATCcaa > 2:535508/1‑151 (MQ=255) gATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCGCTCCC‑CATTTCGCGA‑CG‑CCTTTTTTTAAATTGGCCCGCGCGCGGGGGACCGTTTTTTAAATAATGA‑CCGTAACAACCGCATCaaa > 1:109260/1‑151 (MQ=255) gATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CAATTCGCGG‑CG‑CCTTTTTTAAAATTGGCCTGCGCGCGGTGACCTGTTTTTTAACTACTGA‑CCGCAAAAACGGCATCAAagagaagcgggg > 1:7771/1‑143 (MQ=255) aTGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑CCGTAACAACGGCATCAAAGAGTTGCAGGa < 1:331500/150‑1 (MQ=255) gTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGATCTCTCTCCC‑AATATCGCGG‑CG‑CACTTTTTAAAATTGGCCCGCCCGCGCTGGCCTGTTTTTTAACTCCTTT‑AAGTAACAACCGCCACAAAGAGCTGTGGGaaag > 1:82590/1‑149 (MQ=255) tACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑CCGTAACAACGGCATCAAAGAGTTGCAGGAAGgt < 2:334231/151‑1 (MQ=255) cATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑CCGTAACAACGGCATCAAAGAGTTGCAGGAAGgtgt < 2:285273/151‑1 (MQ=255) ggTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCCCCC‑CCCTTCGCGC‑CG‑CCCTTTTTTAAATTTGCCTGCGCGCGGTGACCCGTTTTTTAAACACTTA‑CCGTAAAAACAGCATCAAAAAGGTGCAGGGAGGGGTtgaa > 1:683012/1‑151 (MQ=255) ttATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑CCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTagaag < 2:665829/150‑1 (MQ=255) ttATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAATTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑CCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTagaag < 2:428225/150‑1 (MQ=255) tATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑CCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAgaagaa < 1:625424/151‑1 (MQ=255) cGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑CCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAAc < 2:108901/151‑1 (MQ=255) cACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATACCGCGA‑CG‑CATTTTTTTAAAATGGCCTGCGCGCGGGGACCTGTTTTGGAACTAATGA‑CCGTAACAAAAGCACCAAAGAGATGCAGGGAGGGGTGGAAGaaaaaaac > 1:475374/1‑146 (MQ=255) cACCATGCCTCTCCGGCGTATATCGGCGGGTCCGTCTCCACCATTCCCGGC‑CC‑CCTTTTTAAAATTTGG‑CTGCGCGCGGCGGCCTGTTTTTTAAATTCTTACCCGTAAAAAC‑GCCTCAAAAAGATGGAGGAAGGGGGAGAAGaaaaaaaa > 1:7149/1‑146 (MQ=255) ccATGCCTCTCCGGCGTATATCGGCGGGGCGCTCCCCA‑CAATTCTCGACCG‑CTTTTTTAAAAATTGG‑CTGCGCGCGGTGGCCCTTTTTTTAAATAATGA‑CACTAAAAAACGCACCAAAAAAGTGGAGGGAGGGGTAGAAGaaaaaaaacag > 1:510247/1‑144 (MQ=255) tatCGGCGGGTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑CCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTtgattgat < 2:641975/151‑1 (MQ=255) ggTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑CCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCga < 2:573034/150‑1 (MQ=255) a‑caTTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑CCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACaa > 2:590470/1‑151 (MQ=255) gcgcGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑CCGTTTCAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAg > 2:113248/1‑151 (MQ=255) cATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑CCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCt < 2:586555/151‑1 (MQ=255) tttttAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑CCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCtgg > 2:541141/1‑151 (MQ=255) | CTTGATGAAGAGTCGCTCTATTACAGCGGACTGATTTGTTCCCTGGAAGCGATTAAGAGCGGATGTACATCGGTTATCGATCACCATGCCTCTCCGGCGTATATCGGCGGGTCGCTCTCCA‑CATTGCGCGA‑CG‑CATTTTTAAAAGTTGGCCTGCGCGCGATGACCTGTTTTGAAACTACTGA‑CCGTAACAACGGCATCAAAGAGTTGCAGGAAGGTGTAGAAGAAAACATCCGTTTCGCCCGTTTGATTGATGAGGCGAAGAAAGCGACAAGCGAGCCGTATCTGG > NC_000913/3019452‑3019737 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |