Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,206,107	T→C	7.7%	W471R (TGG→CGG)	PP_2811 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,206,107	0	T	C	7.7%	144.4 / 2.7	52	W471R (TGG→CGG)	PP_2811	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/25);  new base C (2/2);  total (25/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 = 3.54e-01

ACGAGGTCAATGGCGAACCACTCGCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTA  >  NC_002947/3205961‑3206255                                                                                                                                                    |                                                                                                                                                     aCGAGGTCAATGGCGAACCACTCGCCCCCGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTgg                                                                                                                                                    <  5:112591/149‑1 (MQ=255)   gAGGTCAATGGCGAACCACTCGCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CCGCCTGACCTATGCCGGCACGTTGCTGCCGGGCTgg                                                                                                                                                    >  6:292371/1‑147 (MQ=255)   gAGGTCAATGGCGAACCACTCGCCCCAGGGCTCAATGCCTGGCAACAGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGAACTCGTCGGGCATCAACCTCGACTT‑CAGGCTGACCTATGACGGCACGTTGCTGCCGGGCTgg                                                                                                                                                    >  7:268168/1‑147 (MQ=255)    aGGTCAATGGCGAACCACTCGCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAg                                                                                                                                                 <  8:3456/149‑1 (MQ=255)    aGGTCAATGGCGAACCACGCGCCCCAGAGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAg                                                                                                                                                 <  1:244400/149‑1 (MQ=255)      gTCAATGGCGAACCACTCGCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGt                                                                                                                                               <  6:27715/149‑1 (MQ=255)           tGGCGAACCACTCGCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGc                                                                                                                                          >  1:192489/1‑149 (MQ=255)              cGAACCACTCGCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGTCACGTTGCTGCCGGGCTGGCAGGTCACGCCCg                                                                                                                                       <  4:288890/149‑1 (MQ=255)                   cACTCGCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGtc                                                                                                                                  <  4:242744/149‑1 (MQ=255)                     cTCGCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGtctt                                                                                                                                >  6:81075/1‑149 (MQ=255)                        gCCCCAGGGCTCAATTCCTGGCAACTGTAAGCTGCGCCCGCCCCCAAGTCCCGTGGCGACAAGACCTAGTCGGGCATCAACCTTGAATTAGAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGtcttc                                                                                                                               <  1:103498/148‑1 (MQ=255)                        gCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGtcttct                                                                                                                              <  8:22023/148‑1 (MQ=255)                        gCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTa                                                                                                                             >  1:296865/1‑149 (MQ=255)                               ggCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCgcgc                                                                                                                      <  4:170244/149‑1 (MQ=255)                                 cTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGcgctc                                                                                                                    >  3:217338/1‑149 (MQ=255)                                   cAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGcgctcgc                                                                                                                  >  2:22555/1‑149 (MQ=255)                                    aaTGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGcgctcgct                                                                                                                 >  7:7887/1‑149 (MQ=255)                                        ccTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTggc                                                                                                             <  1:22555/149‑1 (MQ=255)                                          tGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTggcgg                                                                                                           >  4:174780/1‑149 (MQ=255)                                             cAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTggcgggcg                                                                                                        <  2:296865/149‑1 (MQ=255)                                                     tCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGa                                                                                                <  5:81075/149‑1 (MQ=255)                                                     tCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCAGCCCGGCAgcaac                                                                                                                                             >  6:101201/1‑100 (MQ=255)                                                     tCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCAGCCCGGCAgcaac                                                                                                                                             >  6:106362/1‑100 (MQ=255)                                                     tCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCAGCCCGGCAgcaac                                                                                                                                             <  5:106362/104‑5 (MQ=255)                                                     tCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCAGCCCGGCAgcaac                                                                                                                                             <  5:101201/104‑5 (MQ=255)                                                       cTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGaa                                                                                               <  2:128714/148‑1 (MQ=255)                                                          cGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCCCGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGTGCGAACGCCGAACCt                                                                                            >  2:63373/1‑148 (MQ=255)                                                             cAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAgc                                                                                        <  2:265413/149‑1 (MQ=255)                                                             cAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTCGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAgc                                                                                        <  5:292371/149‑1 (MQ=255)                                                                  cccAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAAc                                                                                              <  6:122315/138‑1 (MQ=255)                                                                  cccAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAAc                                                                                              >  5:122315/1‑138 (MQ=255)                                                                      aGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCtt                                                                                >  5:163887/1‑148 (MQ=255)                                                                        tCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCAc                                                                             >  5:184397/1‑149 (MQ=255)                                                                             tGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAg                                                                        >  8:252634/1‑149 (MQ=255)                                                                                     aGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGagc                                                                <  8:7887/149‑1 (MQ=255)                                                                                      gACCTCGGCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGagca                                                               <  1:105765/149‑1 (MQ=255)                                                                                        ccTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGagcagc                                                             <  4:217338/149‑1 (MQ=255)                                                                                                ggCATCAAGCTAGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGACGGGCTGGCAGGTCACGCGCGGCGTCTTCTACTCGCGCTCGCTTGGCTGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGAAGCGAGCAGCATGAACCt                                                     <  8:268168/149‑1 (MQ=255)                                                                                                ggCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCg                                                                                                 >  2:296803/1‑105 (MQ=255)                                                                                                ggCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCg                                                                                                 <  1:296803/105‑1 (MQ=255)                                                                                                ggCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCt                                                     <  3:174780/149‑1 (MQ=255)                                                                                                  cATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGt                                                   <  2:192489/149‑1 (MQ=255)                                                                                                     cAACCTGGACTT‑CAGCCTGCCCTATGACGGCACATTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTAcc                                                >  6:353861/1‑149 (MQ=255)                                                                                                         cTGGACTT‑CAGCCTGACCTAAGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACTCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGaa                                            >  4:171821/1‑149 (MQ=255)                                                                                                                  cAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTa                                     <  6:184397/148‑1 (MQ=255)                                                                                                                     ccTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCACTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGa                                  >  8:78341/1‑148 (MQ=255)                                                                                                                         aCCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAAccc                             >  2:312635/1‑149 (MQ=255)                                                                                                                         aCCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAAccc                             >  6:120671/1‑149 (MQ=255)                                                                                                                             aTGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAACACCGGc                          >  8:236674/1‑148 (MQ=255)                                                                                                                                     aCGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTgg                   >  1:322880/1‑147 (MQ=255)                                                                                                                                              ccGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCt        <  7:252634/149‑1 (MQ=255)                                                                                                                                              ccGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCt        >  6:101073/1‑149 (MQ=255)                                                                                                                                                    tGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTa  >  1:60831/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ACGAGGTCAATGGCGAACCACTCGCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTCAACTA  >  NC_002947/3205961‑3206255

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