Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,042,503	A→C	10.4%	V18G (GTG→GGG)	PP_1816 ←	alcohol dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,042,503	0	A	C	10.4%	130.0 / 6.3	47	V18G (GTG→GGG)	PP_1816	alcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (22/20);  new base C (2/3);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

GGCGGGATGTAGGACTTGCCTTCGTTCATCCAGCCGCGCATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAGGCGTTCAATCTACCTTGTTGATGAAGGCGGTTCCTTAAAATCACTGACAAGCATGGTAGC  >  NC_002947/2042365‑2042647                                                                                                                                             |                                                                                                                                                 ggcgggATGTAGGACTTGCCTTCGTTCATCCAGCCGCGCATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCCGTGCGCCgtg                                                                                                                                        <  4:365690/149‑3 (MQ=255)       aTGTAGGACTTGCCTTCGTTCATCCAGCCGCGCATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCg                                                                                                                                  >  4:144343/1‑149 (MQ=255)        tGTAGGACTTGCCTTCGTTCATCCAGCCGCGCATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGt                                                                                                                                 >  1:144292/1‑149 (MQ=255)         gTAGGACTTGCCTTCGTTCATCCAGCCGCGCATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCCGTGCGCCgtga                                                                                                                                       >  1:285823/1‑139 (MQ=255)         gTAGGACTTGCCTTCGTTCATCCAGCCGCGCATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCCGTGCGCCgtga                                                                                                                                       <  2:285823/142‑4 (MQ=255)          tAGGACTTGCCTTCGTTCATCCAGCCGCGCATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGttt                                                                                                                               <  2:307339/149‑1 (MQ=255)              aCTTGCCTTCGTTCATCCAGCCGCGCATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGttt                                                                                                                               >  8:251926/1‑145 (MQ=255)              aCTTGCCTTCGTTCATCCAGCCGCGCATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGttt                                                                                                                               <  7:251926/145‑1 (MQ=255)                     ttCGTTCATCCAGCCGCGCATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTCCCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCCCGGCGCCGACCGGGCGTTTGGCAAGCAgg                                                                                                                     >  3:218784/1‑148 (MQ=255)                      tCGTTCATCCAGCCGCACATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAgg                                                                                                                     >  6:342055/1‑147 (MQ=255)                        gTTCATCCAGCCGCGCATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCCGTGCGCCgtgac                                                                                                                                      <  5:61367/128‑5 (MQ=255)                                aGCCGCGCATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTg                                                                                                         <  1:314464/149‑1 (MQ=255)                                      ggaTGGCCGGGTCCAGCGACAGGTACAGGTTGCGG‑ACCAGTACCTGGCCTTCGCCGGGGTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTgg                                                                                                   <  7:97730/147‑1 (MQ=255)                                      gcATGGCCTTGTCCAGCTTCAGTTACAGGTTGCGC‑ACCAGTACCTGGGCTTCGGCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCTT‑CACGGCGGACGGCGCCGACCGGGCGTTTGGCAAGCAGTAAGCGGCGGTTTGTCTgg                                                                                                   <  8:316582/149‑1 (MQ=255)                                      gcATGGCCGGTTCCAGCGACAGTTACAGGTTGCGC‑ACCATTACCTGGCCTTCGCCGGGTTTTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTgg                                                                                                   <  8:232896/149‑1 (MQ=255)                                      gcATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGTCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTgg                                                                                                   <  2:167099/149‑1 (MQ=255)                                      gcATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTgg                                                                                                   <  2:144292/149‑1 (MQ=255)                                            ccGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATg                                                                                             <  7:152959/149‑1 (MQ=255)                                             cGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑cacggcgcac                                                                                                                                                 >  7:188652/1‑96 (MQ=255)                                             cGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑cacggcgcac                                                                                                                                                 <  8:188652/96‑1 (MQ=255)                                                gTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCgg                                                                                         <  4:307772/149‑1 (MQ=255)                                                      ggccaggtactgg‑tGCGCAACCTGTACCTCGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCtt                                                                                   >  3:87319/4‑149 (MQ=255)                                                       gACAGGTACAGGTTGAGC‑ACCAGTACCTGGCCTTCGCCGGGGTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGGCATGGCGGGTCCTtg                                                                                  <  4:218784/149‑1 (MQ=255)                                                        acaggtacaggtTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGCTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGCCAAGCAGGAAGCGGCGGTTGGTCAGGGTCATGGCGGGTCCTtgt                                                                                 >  2:356611/1‑149 (MQ=255)                                                         caggtacaggtTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTtgtg                                                                                >  8:134248/1‑149 (MQ=255)                                                            gtacaggtTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTtgt                                                                                 <  5:167406/145‑1 (MQ=255)                                                            gtacaggtTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTtgt                                                                                 >  6:167406/1‑145 (MQ=255)                                                             tacaggtTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCg                                                                            <  5:342055/149‑1 (MQ=255)                                                              acaggtTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTg                                                                                                         >  4:273304/1‑119 (MQ=255)                                                                       gc‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCg                                                                  >  7:226844/1‑149 (MQ=255)                                                                            cAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGttt                                                                                                                               <  2:143621/84‑1 (MQ=255)                                                                            cAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGttt                                                                                                                               >  1:143621/1‑84 (MQ=255)                                                                            cAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAg                                                              <  3:144343/149‑1 (MQ=255)                                                                              gTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCAAGGCGCCGACCGGGCGTTTCGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAGGc                                                            >  6:54859/1‑149 (MQ=255)                                                                                  cTGGCTTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAGGCGtt                                                         >  7:261770/1‑148 (MQ=255)                                                                                  cTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTgg                                                                                                                             <  5:85085/80‑1 (MQ=255)                                                                                  cTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTgg                                                                                                                             >  6:85085/1‑80 (MQ=255)                                                                                  cTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAGGCGtt                                                         >  2:123390/1‑148 (MQ=255)                                                                                   tGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAGGCGTTc                                                        >  2:337512/1‑148 (MQ=255)                                                                                           gCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAGGCGTTCAATCTAcc                                                <  7:134248/148‑1 (MQ=255)                                                                                                   gttCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAGGCGTTCAATCTACCTTGTtgat                                        >  6:318210/1‑148 (MQ=255)                                                                                                   gttCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAGGCGTTCAATCTACCTTGTtgat                                        >  7:223912/1‑148 (MQ=255)                                                                                                    ttCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGGTGGTCTGGGTCATGGCGGGTCCTGGTGGGCGGTGGAGGGCGTAAGGc                                                            >  2:269110/1‑127 (MQ=255)                                                                                                       ggcggGTACGGTCTCGAAGCTGAAGTCGTGCACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAGGCGTTCAATCTACCTTGTTGATGAAg                                    >  1:27035/1‑149 (MQ=255)                                                                                                               cGGTCTCGAAGCTGCAGTCGG‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGGTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAGGCGTTCAATCTACCTTGTTGATGAAGGCGGTTCCt                           <  1:356611/149‑1 (MQ=255)                                                                                                                      gAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAGGCGTTCAATCTACCTTGTTGATGAAGGCGGTTCCTTaaaa                      <  2:27035/147‑1 (MQ=255)                                                                                                                      gAAGCTGAAGACGT‑CACGGCGCACGGCGCCGCCCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGCTCCTTGTGGGCGGTGGACGGCGTAACGCGTTCAATCTACCTTGTAGATCAAGGCTGTTCCTTAAAAt                     >  6:347047/1‑148 (MQ=255)                                                                                                                           tGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAGGCGTTCAATCTACCTTGTTGATGAAGGCGGTTCCTTAAAATCACTg                >  3:142620/1‑148 (MQ=255)                                                                                                                             aaGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAGGCGTTCAATCTACCTTGTTGATGAAGGCGGTTCCTTAAAATCACTGAc              <  8:223912/148‑1 (MQ=255)                                                                                                                                     cacggcgcacggcgCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAGGCGTTCAATCTACCTTGTTGATGAAGGCGGTTCCTTAAAATCACTGACAAGCATgg      >  8:111108/1‑149 (MQ=255)                                                                                                                                      ccggcgcacggcgcCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCGGGGTCATGGCGGGTCCTTGTGGGGGGTGGAGGGCGTAAGGGGTTCAATCTAGCTTGTTGATGAAGGCGGTTCCTTAAAATCACTGACAAGCATGGt     <  2:352420/148‑1 (MQ=255)                                                                                                                                      acggcgcacggcgcGGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTTGTCTGGGTCATGGCGGGTCCTTGTGGGCGTTGTAGGGCGTAATGCGTTCAATCTACCTTGTTGATGAAGGCGGTTCCTTAAAATCACTGACAAGCATgg      <  5:54859/148‑1 (MQ=255)                                                                                                                                         gcgcacggcgcCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGTGGGCGTAAGGCGTTCAATCTACGTTGTTGATGAAGGCGGTTCCTTAAAATCACTGACAAGCATGGTAGc  <  4:87319/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                 GGCGGGATGTAGGACTTGCCTTCGTTCATCCAGCCGCGCATGGCCGGGTCCAGCGACAGGTACAGGTTGCGC‑ACCAGTACCTGGCCTTCGCCGGGTTGTTCGGCGGGTACGGTCTCGAAGCTGAAGTCGT‑CACGGCGCACGGCGCCGACCGGGCGTTTGGCAAGCAGGAAGCGGCGGTTGGTCTGGGTCATGGCGGGTCCTTGTGGGCGGTGGAGGGCGTAAGGCGTTCAATCTACCTTGTTGATGAAGGCGGTTCCTTAAAATCACTGACAAGCATGGTAGC  >  NC_002947/2042365‑2042647

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑