Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,913	0	G	C	47.1%	67.1 / 20.1	52	G44G (GGC→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (20/7);  new base C (0/24);  total (20/32)
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

ACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGG  >  NC_000913/993765‑994051                                                                                                                                                            |                                                                                                                                           aCGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGcc                                                                                                                                          >  2:805518/1‑151 (MQ=255)        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)                                                                                                                                                            |                                                                                                                                           ACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAG‑CAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGG  >  NC_000913/993765‑994051

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: ‑