Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,206,103	G→A	7.7%	P469P (CCG→CCA)	PP_2811 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,206,103	0	G	A	7.7%	146.4 / 2.4	52	P469P (CCG→CCA)	PP_2811	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/25);  new base A (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 = 6.30e-01

ACGAGGTCAATGGCGAACCACTCGCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT  >  NC_002947/3205961‑3206249                                                                                                                                                |                                                                                                                                                   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:106362/1‑100 (MQ=255)                                                     tCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCAGCCCGGCAgcaac                                                                                                                                       >  6:101201/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                                                                                        >  5:122315/1‑138 (MQ=255)                                                                  cccAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAAc                                                                                        <  6:122315/138‑1 (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                                                                                           <  1:296803/105‑1 (MQ=255)                                                                                                ggCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCg                                                                                           >  2:296803/1‑105 (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)                                                                                                                                                |                                                                                                                                                   ACGAGGTCAATGGCGAACCACTCGCCCCAGGGCTCAATGCCTGGCAACTGGATCCTGCGACAGCGCCCAAGTCCCGTGGCGACAAGACCTCGTCGGGCATCAACCTGGACTT‑CAGCCTGACCTATGACGGCACGTTGCTGCCGGGCTGGCAGGTCACGCCCGGCGTCTTCTACTCGCGCTCGCTTGGCGGGCGAACGCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT  >  NC_002947/3205961‑3206249

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