Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 993,917 | 0 | A | C | 50.0% | 51.0 / 17.9 | 50 | V43G (GTG→GGG) | ssuB | aliphatic sulfonate ABC transporter ATPase |
Reads supporting (aligned to +/- strand): ref base A (19/6); new base C (0/25); total (19/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.16e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCA > NC_000913/993772‑994060 | cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCaccacc < 2:19710/151‑1 (MQ=255) ccTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTg > 1:716137/1‑151 (MQ=255) ttCGCGTTTC‑TCCTTAATTTCCCGCCGGCGG‑TGTGGTCCCCGCAAACACATCCCGGGCGGTTGGCTTTC‑CCA‑ACCCTGCCAG‑CAGGCGCAGCAGGGTTTTTTT‑C‑CCCCCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCg < 2:828300/147‑1 (MQ=255) gtttctcCCTAATTTT‑CAGCCAGGGG‑TGGGGTGCCCGCAAACACATCGCCTGCGTTGGGCGTTC‑CCA‑GCCCTGCCAG‑CAGGCGCAGCGGGTTACTTTT‑C‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGt < 2:101775/147‑1 (MQ=255) gtgtCTTCCTGAATTT‑CAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt > 2:400890/1‑150 (MQ=255) tgtttCCCGAATTT‑TCGCCCGGGGTTGGGGTGCCCGCAAACACCCCGCCGGGGTTGGGCTTTC‑CCA‑GCCCGGCCGG‑CAGGCGCAGCAGGTTCCTTTT‑C‑CCCCCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAatat < 1:373524/144‑1 (MQ=255) cTTCCTGAATTT‑CAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATc > 1:774042/1‑151 (MQ=255) ccTGATTTC‑CAGCCGGGGG‑GGGGGGCCCCGCAACCACCCCGCCTG‑GGTTGGGGTTTCCCA‑GCCCGCCCAG‑AGGGCGCGGGGGGGT‑CTTTT‑GCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAatat < 1:499874/141‑1 (MQ=255) ccTGATTTC‑CACCCAGGGG‑TGTGGTCCCCGCTAACACACCCCCGGGGGTTGGCTTTC‑CCA‑GCCCTGCCAG‑CAGGCGCAGCAGGTTCCTTTT‑C‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAg < 2:306828/151‑1 (MQ=255) ccTGATTTC‑CACCCAGCGG‑TGGGGTCCCCGCAACCACCGCCCGGGGGTTTGGCTTTT‑CCA‑GCCCTGCCAG‑AGGGCGCAGGAGGGTTCTTTT‑C‑CCCCCACCGCTGCGCCCCACCTCCGCCACAAACTGACCTGCCGGAATATGTAAATc < 2:204026/148‑1 (MQ=255) taatttt‑CAGCGAGCGG‑GGTGGTCCCCGAAAACACATCGCGGGGGGTGGGCTTTC‑CCA‑GCCCTCCCAG‑AAGGCGCAGGGGGTTTCTTTT‑G‑CCACCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCa < 1:547616/144‑1 (MQ=255) atttt‑CACCCAGGGG‑TGTGGTCCCCGCAACCACACCGCCGGGGGTTGGCTTTC‑CCA‑GCCCTGCCAG‑CGGGCGCAG‑AGGGTTCCTTT‑GCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc < 2:458013/143‑1 (MQ=255) ttt‑CGGCCGGGGG‑TGGGTCGCCCACAAACACACCGCTTGGGGTGGGC‑TTT‑CCA‑GGCCCTCCGGCCAGGCGAGGTGGGTTCCTTTT‑G‑CCACACCGGCGGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt < 1:309657/147‑1 (MQ=255) ttt‑CAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt > 1:384035/1‑151 (MQ=255) ttccagccgggGG‑TGTGGTCCCCGAAAACACAGCGCCTGGGGTGGGCTTTC‑CCA‑GCCCTCCCGG‑AAGGCGCAGCAGGTTTCTTTT‑G‑CCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg < 1:351808/145‑1 (MQ=255) ccggggg‑TGGGGTCCCCGCAAACACACCGCCTGCGGTTGGCTTTC‑CAA‑GCCCTGCCAG‑CAGGCGCAGCAGGTTACTTTT‑C‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGa < 2:87522/145‑1 (MQ=255) caggGG‑TGTGGTGCCCGTAAACACACCGCCGGGGGTTGGCCTTT‑CCA‑CCCTGGCCGG‑CAGGCCCGGCGGGTTCCTTTC‑G‑CCCCCCCCGCCGCCGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg < 2:382676/151‑1 (MQ=255) gggtcCCCGCAAACACACCTCCTGGGTTGGGCTTTT‑CCA‑ACCCTCCCAG‑CAGGCCCAGCGGGTTTCTTTT‑C‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg < 2:703920/142‑1 (MQ=255) ggtcCCCCCTAACCCGCCGCCGGTGGTTGGCTTTC‑CCA‑GCCCGGCCAG‑CGGGCGCAGCGGGGTTTTTTT‑G‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc < 1:484043/151‑1 (MQ=255) gCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg > 1:613852/1‑151 (MQ=255) cccGCTAACACATCGCTTGCGTTGGGGTTTT‑CGG‑GCCCTGCGAG‑CTGGCGGAGCGGGTTGTTTTC‑C‑CCCCCCCCGTGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgc < 2:670113/149‑1 (MQ=255) cccGCAACCACTCCGCCTGCGTTTGCCTTTC‑CCA‑GCCCTCCCAG‑CGGGCGCGGCGGGGTATTTTT‑C‑CCACCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg < 1:175622/150‑1 (MQ=255) ccGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt > 1:727573/1‑150 (MQ=255) cGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTa > 1:8976/1‑150 (MQ=255) cTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg > 2:456380/1‑151 (MQ=255) cacaTCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg > 1:260445/1‑151 (MQ=255) cccggcccTGGGGTTTGGCTTT‑CCAGGCCCTGCCAG‑CAGGCGCAG‑AGGGTTTTTTTCG‑CCACCCCCGCCGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt < 1:313636/142‑1 (MQ=255) ggtgggtgTTC‑CAA‑GCCCTGCCAG‑CGGGCGCAGCAGGGTACTTTG‑G‑CCACCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg < 2:603078/151‑1 (MQ=255) gTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg > 2:758978/1‑150 (MQ=255) ggCTTTC‑CCA‑GCCCTGCCAG‑CAGGCGCGGGGGGTTTTTTTG‑G‑CCCCCCCCGCTGCGCCCCACCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg < 1:587880/147‑1 (MQ=255) ggCTTTC‑CCA‑GCCCTGCCAG‑CAGGCGCAGCGGGGTTCTTTC‑G‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg < 2:599983/149‑1 (MQ=255) gCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca > 1:402133/1‑150 (MQ=255) ggctttcccacccccgCCAG‑CGGGCGCAGCAGGTTCCTTTT‑C‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc < 2:326081/135‑1 (MQ=255) gTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac > 1:489921/1‑150 (MQ=255) ttt‑CCG‑GCCCGGCCAG‑CAGGCGCGGCGGGGTACTTTT‑C‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac < 2:702106/149‑1 (MQ=255) ccA‑GCCCGGCCAG‑CAGGCGCAGCGGGTTACTTTT‑C‑CCCCCACCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg < 2:644085/151‑1 (MQ=255) a‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg > 2:744994/1‑151 (MQ=255) aCCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg > 2:647227/1‑151 (MQ=255) aCCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg > 1:96581/1‑151 (MQ=255) agGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc > 1:377923/1‑150 (MQ=255) gcgcAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc > 2:763913/1‑148 (MQ=255) ggTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt > 1:845559/1‑150 (MQ=255) gtgttttttgtccacccCGGCCGCCGCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa < 2:446871/138‑1 (MQ=255) gtgttttt‑C‑CCACCCCCTCTGGCACCCACCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAg < 1:340917/144‑1 (MQ=255) ggttattttccccaccgCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa < 1:188855/138‑1 (MQ=255) ttcttttcccccccaccGTTCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGa < 1:766132/132‑1 (MQ=255) ccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATa < 1:327836/141‑1 (MQ=255) ccccaccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg < 2:796320/149‑1 (MQ=255) caccacCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg > 2:464224/1‑151 (MQ=255) cccccGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg < 2:535049/146‑1 (MQ=255) ccaccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg < 2:300310/150‑1 (MQ=255) ccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct < 2:232022/144‑1 (MQ=255) ccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa < 2:855553/147‑1 (MQ=255) | CAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCA > NC_000913/993772‑994060 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |