Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,9170AC50.0% 51.0 / 17.9 50V43G (GTG→GGG) ssuBaliphatic 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: ‑