Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,000,768	C→G	7.7%	T76T (ACG→ACC)	PP_2624 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,000,768	0	C	G	7.7%	163.3 / 4.2	52	T76T (ACG→ACC)	PP_2624	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/24);  new base G (2/2);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

GGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCA  >  NC_002947/3000620‑3000915                                                                                                                                                     |                                                                                                                                                       gggCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTAcc                                                                                                                                                        >  6:146974/1‑149 (MQ=255)       ccAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGAGCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACg                                                                                                                                                  >  8:162100/1‑149 (MQ=255)         aGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTc                                                                                                                                                <  6:298582/149‑1 (MQ=255)           tGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTc                                                                                                                                                >  4:187528/1‑147 (MQ=255)                   gcAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTcaca                                                                                                                                      >  5:220489/1‑149 (MQ=255)                   gcAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTcaca                                                                                                                                      <  1:5264/149‑1 (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                                                                                                                                             <  8:155499/134‑1 (MQ=255)                           agcagcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAg                                                                                                                                             >  7:155499/1‑134 (MQ=255)                              agcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTg                                                                                                                                 >  3:262856/1‑141 (MQ=255)                              agcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTg                                                                                                                                 <  4:262856/141‑1 (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                                                                                                     >  4:254313/1‑138 (MQ=255)                                                            aTCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAat‑gta                                                                                                     <  3:254313/138‑1 (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                                                            >  1:12686/1‑149 (MQ=255)                                                                                          cTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACg                                                            >  8:205284/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                  >  6:168938/1‑149 (MQ=255)                                                                                                                                    ccGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTtgtg                  <  5:168938/149‑1 (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)                                                                                                                                                     |                                                                                                                                                       GGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAAT‑GTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCA  >  NC_002947/3000620‑3000915

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