Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,000,765	T→G	7.4%	V77V (GTA→GTC)	PP_2624 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,000,765	0	T	G	7.4%	167.3 / 3.4	54	V77V (GTA→GTC)	PP_2624	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base G (2/2);  total (27/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.97e-01

CAGGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACC  >  NC_002947/3000618‑3000912                                                                                                                                                    |                                                                                                                                                      cAGGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGt                                                                                                                                                       <  3:215004/148‑1 (MQ=255)  aGGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTa                                                                                                                                                      <  4:232112/148‑1 (MQ=255)   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                                                                                                                                  <  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                                                                                                                             >  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‑‑TTGGCGTCGTTCAACACGCTGTTGAatg                                                                                                    >  5:141643/1‑148 (MQ=255)                                                    cTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAatgtat                                                                                                 <  3:187528/149‑1 (MQ=255)                                                    cTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAatgtat                                                                                                 >  3:280174/1‑149 (MQ=255)                                                       gcgcgcGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAatgtatgt                                                                                               <  5:262612/148‑1 (MQ=255)                                                         gcgcgATCCGGTCCGGCCCGAAATCGTTGATGGATCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTa                                                                                            <  7:162100/149‑1 (MQ=255)                                                              aTCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAatgta                                                                                                  >  4:254313/1‑138 (MQ=255)                                                              aTCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAatgta                                                                                                  <  3:254313/138‑1 (MQ=255)                                                                 cggtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATgc                                                                                     <  1:290748/148‑1 (MQ=255)                                                                   gtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCa                                                                                  >  8:130654/1‑149 (MQ=255)                                                                                tCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTc                                                                     <  6:220489/149‑1 (MQ=255)                                                                                  gTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTcc                                                                    <  5:126208/148‑1 (MQ=255)                                                                                      aTGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAg                                                               >  3:144870/1‑149 (MQ=255)                                                                                          aCCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCa                                                           <  4:280174/149‑1 (MQ=255)                                                                                            cTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACg                                                         >  8:205284/1‑149 (MQ=255)                                                                                            cTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACg                                                         >  1:12686/1‑149 (MQ=255)                                                                                             tGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACgg                                                        >  1:48133/1‑149 (MQ=255)                                                                                               aaCTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACgg                                                        <  2:12686/147‑1 (MQ=255)                                                                                                aCTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTg                                                     <  2:48133/149‑1 (MQ=255)                                                                                                       gTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCTCATGACCCGACGGTCCTCAAGCGCACGGCTGTCATCCa                                              <  6:141643/149‑1 (MQ=255)                                                                                                             ggCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACAAGGTGGCGTCGTTCAACACGCTGTTGAATTTATTTTTATCGATGCGCATTACCTGACGGTCCTAAAGCGCACGGCTGTAATCCAGGc                                           >  1:93596/1‑148 (MQ=255)                                                                                                              gCCAGCTCTGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGTCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCTCATGACCTGACGGTCCTCAAGCGCACTGCTGTGATCCAGGCCGaa                                       <  6:311611/149‑1 (MQ=255)                                                                                                               ccAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCAAA‑‑TTGGCGTCGTTCAACACGCTGTTGATTGTATGTTTTTCGATGCGCATGACCTGCCGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAg                                      >  8:282127/1‑149 (MQ=255)                                                                                                               ccAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACCGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAg                                      >  5:247648/1‑149 (MQ=255)                                                                                                                        tATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCccagc                             <  7:205284/149‑1 (MQ=255)                                                                                                                               ggCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCccagcaccagca                      >  6:286677/1‑149 (MQ=255)                                                                                                                                gCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCAt                     <  5:286677/149‑1 (MQ=255)                                                                                                                                   gTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCtt                  <  7:22685/149‑1 (MQ=255)                                                                                                                                      ccGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTtgtg               >  6:168938/1‑149 (MQ=255)                                                                                                                                      ccGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTtgtg               <  5:168938/149‑1 (MQ=255)                                                                                                                                      ccGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTtgtg               >  2:259362/1‑149 (MQ=255)                                                                                                                                       tgAGCCTGGACGTGACGGTACCGTCCAGGCTCGCA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTtgttc              >  7:254451/2‑147 (MQ=21)                                                                                                                                       cGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGc              <  7:130654/149‑1 (MQ=255)                                                                                                                                       cGAGCCTGGACCGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGc              <  6:247648/149‑1 (MQ=255)                                                                                                                                         aGCCTGGACGTGACGGTACCGTCCAGGCTCGCA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTTCgg            <  8:254451/149‑1 (MQ=21)                                                                                                                                             tGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGc        >  6:80369/1‑149 (MQ=255)                                                                                                                                                   gTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCAcc  >  6:221129/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CAGGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACA‑‑TTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACC  >  NC_002947/3000618‑3000912

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