Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,000,772	A→C	7.3%	V75G (GTG→GGG)	PP_2624 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,000,772	0	A	C	7.3%	169.0 / 4.4	55	V75G (GTG→GGG)	PP_2624	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base C (2/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

CCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCC  >  NC_002947/3000626‑3000918                                                                                                                                                   |                                                                                                                                                      ccAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGAGCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACg                                                                                                                                                     >  8:162100/1‑149 (MQ=255)   aGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTc                                                                                                                                                   <  6:298582/149‑1 (MQ=255)     tGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTc                                                                                                                                                   >  4:187528/1‑147 (MQ=255)             gcAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTcaca                                                                                                                                         <  1:5264/149‑1 (MQ=255)             gcAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTcaca                                                                                                                                         >  5:220489/1‑149 (MQ=255)                 cgcCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTgg                                                                                                                                   >  2:290748/1‑149 (MQ=255)                  gcCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGc                                                                                                                                  >  8:22685/1‑149 (MQ=255)                    cagcagCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGCGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGcgt                                                                                                                                >  7:81035/1‑149 (MQ=255)                     agcagcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAg                                                                                                                                                >  7:155499/1‑134 (MQ=255)                     agcagcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAg                                                                                                                                                <  8:155499/134‑1 (MQ=255)                        agcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTg                                                                                                                                    <  4:262856/141‑1 (MQ=255)                        agcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTg                                                                                                                                    >  3:262856/1‑141 (MQ=255)                                cAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGTGACGGTACCGTCCAg                                                                                                                                                >  6:140277/1‑123 (MQ=255)                                cAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGTGACGGTACCGTCCAg                                                                                                                                                <  5:140277/123‑1 (MQ=255)                                        gCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAat                                                                                                            <  5:146974/149‑1 (MQ=255)                                        gCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGGGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAat                                                                                                            <  4:283221/149‑1 (MQ=255)                                          aCCTGGCGAGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAat‑g                                                                                                          >  5:141643/1‑148 (MQ=255)                                            cTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAat‑gtat                                                                                                       <  3:187528/149‑1 (MQ=255)                                            cTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAat‑gtat                                                                                                       >  3:280174/1‑149 (MQ=255)                                               gcgcgcGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAat‑gtatgt                                                                                                     <  5:262612/148‑1 (MQ=255)                                                 gcgcgATCCGGTCCGGCCCGAAATCGTTGATGGATCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTa                                                                                                  <  7:162100/149‑1 (MQ=255)                                                      aTCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAat‑gta                                                                                                        <  3:254313/138‑1 (MQ=255)                                                      aTCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAat‑gta                                                                                                        >  4:254313/1‑138 (MQ=255)                                                         cggtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATgc                                                                                           <  1:290748/148‑1 (MQ=255)                                                           gtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCa                                                                                        >  8:130654/1‑149 (MQ=255)                                                                        tCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTc                                                                           <  6:220489/149‑1 (MQ=255)                                                                          gTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTcc                                                                          <  5:126208/148‑1 (MQ=255)                                                                              aTGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAg                                                                     >  3:144870/1‑149 (MQ=255)                                                                                  aCCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCa                                                                 <  4:280174/149‑1 (MQ=255)                                                                                    cTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACg                                                               >  8:205284/1‑149 (MQ=255)                                                                                    cTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACg                                                               >  1:12686/1‑149 (MQ=255)                                                                                     tGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACgg                                                              >  1:48133/1‑149 (MQ=255)                                                                                       aaCTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACgg                                                              <  2:12686/147‑1 (MQ=255)                                                                                        aCTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTg                                                           <  2:48133/149‑1 (MQ=255)                                                                                               gTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCTCATGACCCGACGGTCCTCAAGCGCACGGCTGTCATCCa                                                    <  6:141643/149‑1 (MQ=255)                                                                                                     ggCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACAAGGTGGCGTCGTTCAACACGCTGTTGAAT‑TTATTTTTATCGATGCGCATTACCTGACGGTCCTAAAGCGCACGGCTGTAATCCAGGc                                                 >  1:93596/1‑148 (MQ=255)                                                                                                      gCCAGCTCTGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGTCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCTCATGACCTGACGGTCCTCAAGCGCACTGCTGTGATCCAGGCCGaa                                             <  6:311611/149‑1 (MQ=255)                                                                                                       ccAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCAAA‑‑TTGGCGTCGTTCAACACGCTGTTGATT‑GTATGTTTTTCGATGCGCATGACCTGCCGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAg                                            >  8:282127/1‑149 (MQ=255)                                                                                                       ccAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACCGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAg                                            >  5:247648/1‑149 (MQ=255)                                                                                                                tATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCccagc                                   <  7:205284/149‑1 (MQ=255)                                                                                                                       ggCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCccagcaccagca                            >  6:286677/1‑149 (MQ=255)                                                                                                                        gCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCAt                           <  5:286677/149‑1 (MQ=255)                                                                                                                           gTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCtt                        <  7:22685/149‑1 (MQ=255)                                                                                                                              ccGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTtgtg                     >  2:259362/1‑149 (MQ=255)                                                                                                                              ccGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTtgtg                     <  5:168938/149‑1 (MQ=255)                                                                                                                              ccGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTtgtg                     >  6:168938/1‑149 (MQ=255)                                                                                                                               tgAGCCTGGACGTGACGGTACCGTCCAGGCTCGCA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTtgttc                    >  7:254451/2‑147 (MQ=21)                                                                                                                               cGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGc                    <  7:130654/149‑1 (MQ=255)                                                                                                                               cGAGCCTGGACCGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGc                    <  6:247648/149‑1 (MQ=255)                                                                                                                                 aGCCTGGACGTGACGGTACCGTCCAGGCTCGCA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTTCgg                  <  8:254451/149‑1 (MQ=21)                                                                                                                                     tGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGc              >  6:80369/1‑149 (MQ=255)                                                                                                                                           gTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCAcc        >  6:221129/1‑149 (MQ=255)                                                                                                                                               cGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGGTATGTTTATCGATGCGCATGACCTGACGGTCATCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACAAGAATATTGTGCGGCAGCTCCACCTCa     >  1:71226/1‑149 (MQ=18)                                                                                                                                                gTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGc   <  1:259362/149‑1 (MQ=255)                                                                                                                                                 tCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGcc  >  3:273467/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                      CCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGCC  >  NC_002947/3000626‑3000918

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