Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,207,889	T→C	8.5%	F513S (TTC→TCC)	PP_2812 →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,207,889	0	T	C	8.5%	129.3 / 2.5	47	F513S (TTC→TCC)	PP_2812	transporter
Reads supporting (aligned to +/- strand):  ref base T (20/23);  new base C (2/2);  total (22/25)
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.65e-01

TTCAATCGCGATGCGGTCTCGATCGCCACCAGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTTCTCCAGCAACATGCGTCAGTTCAACGAGGGGTACAACGAAGGCAATCCCAAAATGGCCGCGGTGCCCATCGACCCGATGAACTATGCCGC  >  NC_002947/3207741‑3208027                                                                                                                                                     |                                                                                                                                            ttCAATCGCGATGCGGTCTCGATCGCCACCAGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCtt                                                                                                                                             >  2:134422/1‑149 (MQ=255) ttCAATCGCGATGCGGTCTCGATCGCCACCAGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGCGTCAGTACCAGGATCGCtt                                                                                                                                             >  3:129794/1‑149 (MQ=255)   cAATCGCGATGCGGTCTCGATCGCCACCAGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCtt                                                                                                                                             <  5:130507/147‑1 (MQ=255)    aaTCGCGATGCGGTCTCGATCGCCACCAGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGt                                                                                                                                          >  3:91693/1‑149 (MQ=255)        gcgATGCGGTCTCGATCGCCACCAGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTgg                                                                                                                                      <  1:134422/149‑1 (MQ=255)         cgATGCGGTCTCGATCGCCACCAGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTgg‑g                                                                                                                                    <  3:14252/149‑1 (MQ=255)              cGGTCTCGATCGCCACCAGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATg                                                                                                                               <  5:136050/149‑1 (MQ=255)                gTCTCGATCGCCACCAGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATg                                                                                                                               >  1:25406/1‑147 (MQ=255)                     gATCGCCACCAGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTc                                                                                                                                         <  5:86908/133‑1 (MQ=255)                     gATCGCCACCAGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTc                                                                                                                                         >  6:86908/1‑133 (MQ=255)                     gATCGCCACCAGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAggggg                                                                                                                        >  2:193957/1‑149 (MQ=255)                     gATCGCCACCAGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAgggg                                                                                                                         <  4:368667/148‑1 (MQ=255)                     gATCGCCACCAGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAggg                                                                                                                          <  7:187923/147‑1 (MQ=255)                           caccaGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCg                                                                                                                  >  2:358006/1‑149 (MQ=255)                             ccaGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGgt                                                                                                                <  2:25406/149‑1 (MQ=255)                              caGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGgtg                                                                                                               <  4:235133/149‑1 (MQ=255)                               aGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAgggg                                                                                                                         <  5:305532/138‑1 (MQ=255)                               aGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAgggg                                                                                                                         >  6:305532/1‑138 (MQ=255)                                  tACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGgt                                                                                                                >  7:112542/1‑144 (MQ=255)                                  tACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGgt                                                                                                                <  8:112542/144‑1 (MQ=255)                                         ttGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTg                                                                                                    <  4:91693/149‑1 (MQ=255)                                          tGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTgg                                                                                                   >  8:37910/1‑149 (MQ=255)                                           ggCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTgg                                                                                                   <  1:358006/148‑1 (MQ=255)                                            gCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGcc                                                                                                 >  2:76707/1‑149 (MQ=255)                                             ccTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGcct                                                                                                >  2:133055/1‑149 (MQ=255)                                              cTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTCCGTCTGG‑GCCATGCAGGGGGTGcc                                                                                                                    <  6:55584/128‑1 (MQ=255)                                              cTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTCCGTCTGG‑GCCATGCAGGGGGTGcc                                                                                                                    <  6:49916/128‑1 (MQ=255)                                              cTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTCCGTCTGG‑GCCATGCAGGGGGTGcc                                                                                                                    >  5:49916/1‑128 (MQ=255)                                              cTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTCCGTCTGG‑GCCATGCAGGGGGTGcc                                                                                                                    >  5:55584/1‑128 (MQ=255)                                               tCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGcctac                                                                                              >  6:267836/1‑147 (MQ=255)                                                      gCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGTGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTtctc                                                                                         <  6:82870/147‑1 (MQ=255)                                                      gCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGGCGCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTtctc                                                                                         <  1:76707/148‑1 (MQ=255)                                                      gCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTTCTcc                                                                                        <  5:267836/148‑1 (MQ=255)                                                       cTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTTCTCCAg                                                                                      >  2:202156/1‑149 (MQ=255)                                                          aGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTTCTCCAGc                                                                                     <  1:133055/147‑1 (MQ=255)                                                                    tCGAAGCCAAGCAGGTGGCCTCTGTTTAAGTCGAGGCGGCGGGCTGGTAGGATGTCGCGGCGGGTCAGTACTAGGATCGCTTCGTCTGG‑GCCATGGTGGGGGTGGCCGGTGTGCTGTCGGTGGCCTCCTTGTCCTGGAACTTGCGtgt                                                                          <  4:129794/148‑3 (MQ=255)                                                                             aCCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCgg                                                                                                      <  4:121066/111‑1 (MQ=255)                                                                             aCCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCgg                                                                                                      >  3:121066/1‑111 (MQ=255)                                                                                             tGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTTCTCCAGCAACATGCGTCAGTTCAACGAGGGGTACAACGAAGGc                                                 >  5:22598/1‑148 (MQ=255)                                                                                                  gcgaggcgGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTTCTCCAGCAACATGCGTCAGTTCAACGAGGGGTTCAACGTAGGCAATcc                                            >  5:307323/1‑148 (MQ=255)                                                                                                          gcggcCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTTCTCCAGCAACATGCGTCAGTTCAACGAGGGGTACAACGAAGGCAATCCCAAAATGGc                                   <  7:37910/149‑1 (MQ=255)                                                                                                          gcggcCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTTCTCCAGCAACATGCGTCAGTTCAACGAGGGGTACAACGAAGGCAATCCCAAAATGGc                                   >  8:52962/1‑149 (MQ=255)                                                                                                          gcggcCTCCCAGGATGTCGCGCCGGGTCAGTACCAGGTTCGCTTCGGCTGG‑GCCATGCAGGGGGTGCCCGGTGTTCTGTCGGTGGCCTCCTGCTCCACCAACATGCGACAGGTCAACGAGGGGTACAACGACGGCACTCCCAAAATGGc                                   >  8:229742/1‑149 (MQ=255)                                                                                                               cTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTTCTCCAGCAACATGCGTCAGTTCAACGAGGGGTACAACGAAGGCAATCCCAAAATGGCCGCgg                              >  3:34683/1‑149 (MQ=255)                                                                                                                  cGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTTCTCCAGCAACATGCGTCAGTTCAACGAGGGGTACAACGAAGGCAATCCCAAAATGGCCGCGGTGc                           <  6:22598/149‑1 (MQ=255)                                                                                                                            gcggcgGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTTCTCCAGCAACATGCGTCAGTTCAACGAGGGGTACAACGAAGGCAATCCCAAAATGGCCGCGGTGCCCATCGAcc                  <  1:202156/148‑1 (MQ=255)                                                                                                                             cggcggGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTtctc                                                                                         <  1:319994/76‑1 (MQ=255)                                                                                                                             cggcggGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTtctc                                                                                         >  2:319994/1‑76 (MQ=255)                                                                                                                                  ggTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTTCTCCAGCAACATGCGTCAGTTCAACGAGGGGTACAACGAAGGCAATCCCAAAATGGCCGCGGTGCCCATCGACCCGATGaa           <  4:34683/149‑1 (MQ=255)                                                                                                                                           cAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTTCTCCAGCAACATGCGTCAGTTCAACGAGGGGTACAACGAAGGCAATCCCAAAATGGCCGCGGTGCCCATCGACCCGATGAACTATGCcgc  >  6:164726/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                            TTCAATCGCGATGCGGTCTCGATCGCCACCAGCTACGACATTGGCCTCGACTGGCTCAGCGTGGTTTTCGAAGCCAACCAGGGGGCCTCTGGTGAAGGCGAGGCGGCGGCCTGCGAGGATGTCGCGGCGGGTCAGTACCAGGATCGCTTCGTCTGG‑GCCATGCAGGGGGTGCCCGGTGTGCTGTCGGTGGCCTCCTTCTCCAGCAACATGCGTCAGTTCAACGAGGGGTACAACGAAGGCAATCCCAAAATGGCCGCGGTGCCCATCGACCCGATGAACTATGCCGC  >  NC_002947/3207741‑3208027

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