Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,812,708	C→T	10.7%	A330V (GCC→GTC)	PP_3368 →	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,812,708	0	C	T	10.7%	145.5 / 8.2	56	A330V (GCC→GTC)	PP_3368	MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (25/25);  new base T (4/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

CTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACTGC  >  NC_002947/3812565‑3812849                                                                                                                                                |                                                                                                                                              cTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGa                                                                                                                                         <  1:780605/150‑1 (MQ=255) cTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGa                                                                                                                                         <  1:573861/150‑1 (MQ=255)      tACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGgcgc                                                                                                                                    >  1:89216/1‑150 (MQ=255)          gATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAg                                                                                                                                >  1:221200/1‑150 (MQ=255)           ctgcaggcggtcgCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGc                                                                                                                               <  2:777088/143‑1 (MQ=255)              tgGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTg                                                                                                                            <  1:364310/150‑1 (MQ=255)                         cGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACAGCAGCTACCCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTc                                                                                                                 >  2:20836/1‑150 (MQ=255)                            gctgctGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCaa                                                                                                              >  1:649497/1‑150 (MQ=255)                                 tgTTCTTTCCCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATg                                                                                                         <  1:411218/150‑1 (MQ=255)                                     cTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCGTgcggc                                                                                                     >  1:793734/1‑150 (MQ=255)                                       ttACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCACTGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGt                                                                                                   >  2:461986/1‑150 (MQ=255)                                           cTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCaa                                                                                               >  1:282408/1‑150 (MQ=255)                                            ttCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAAc                                                                                              >  2:659289/1‑150 (MQ=255)                                             tCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAAcc                                                                                             >  1:260327/1‑150 (MQ=255)                                                ccgccCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTgg                                                                                          >  2:748007/1‑150 (MQ=255)                                                 cgccCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGTCGTCCAACCATGCGGCGTTCAACCTGGc                                                                                         <  1:20836/150‑1 (MQ=255)                                                     cAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAAt                                                                                     <  1:540126/150‑1 (MQ=255)                                                          gcTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCt                                                                                >  2:23347/1‑150 (MQ=255)                                                            tgtgGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTgg                                                                              >  2:703067/1‑150 (MQ=255)                                                                   ccTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGGCAATGTCCATGa                                                                                                                                         <  1:524988/84‑1 (MQ=38)                                                                   ccTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGGCAATGTCCATGa                                                                                                                                         >  2:524988/1‑84 (MQ=38)                                                                        tgctgGTCATCGGCCTGGTGGGCACCAGGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTgg                                                                  <  2:322483/150‑1 (MQ=255)                                                                        tgctgGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTgg                                                                  >  1:574287/1‑150 (MQ=255)                                                                          cagggCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTg                                                                               <  1:264820/133‑1 (MQ=255)                                                                          cagggCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTg                                                                               >  2:264820/3‑135 (MQ=255)                                                                              gCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCa                                                            <  2:574287/150‑1 (MQ=255)                                                                              gCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCa                                                            <  2:260327/150‑1 (MQ=255)                                                                                aTCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATg                                                          <  1:703067/150‑1 (MQ=255)                                                                                 tCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGt                                                                                                   >  2:326972/1‑108 (MQ=255)                                                                                 tCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGt                                                                                                   <  1:326972/108‑1 (MQ=255)                                                                                 tCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTgg                                                                  >  1:228892/1‑141 (MQ=255)                                                                                 tCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTgg                                                                  <  2:228892/141‑1 (MQ=255)                                                                                 tCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATgg                                                         >  1:529259/1‑150 (MQ=255)                                                                                 tCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTCCCACCGGCGCAATGCCCTGGGGCCGTGGCTGGGAGGCATgg                                                         <  1:461986/150‑1 (MQ=255)                                                                                            ggCACCATGATCGCCTTGGCGGCACCGCAACAAATGGGCCGCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGcc                                              >  1:751844/1‑150 (MQ=255)                                                                                             gCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCg                                             >  1:278387/1‑150 (MQ=255)                                                                                               aCCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGc                                                                     >  2:260029/1‑124 (MQ=255)                                                                                               aCCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGc                                                                     <  1:260029/124‑1 (MQ=255)                                                                                               aCCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTAGGAGGAAGGGCGATTACTGGCgtc                                           >  2:767787/1‑148 (MQ=255)                                                                                                   tGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATgg                                       <  1:741385/150‑1 (MQ=255)                                                                                                    gATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATggg                                      >  2:718425/1‑150 (MQ=255)                                                                                                     aTCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCATTGCCGGCATGGGc                                     <  2:793734/150‑1 (MQ=255)                                                                                                     aTCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGc                                     <  1:701712/150‑1 (MQ=255)                                                                                                       cGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTg                                   <  2:649497/150‑1 (MQ=255)                                                                                                                  cACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAATCATGCGGCGTTCAGCCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCg                        <  2:89216/150‑1 (MQ=255)                                                                                                                    ccGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGa                      <  1:23347/150‑1 (MQ=255)                                                                                                                      gCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATa                    <  1:748007/150‑1 (MQ=255)                                                                                                                           aGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCg               <  2:221200/150‑1 (MQ=255)                                                                                                                             aTCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGgc             <  2:282408/150‑1 (MQ=255)                                                                                                                                 gCCTCATGGGCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGc                                                        >  2:484110/1‑103 (MQ=39)                                                                                                                                 gCCTCATGGGCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGc                                                        <  1:484110/103‑1 (MQ=39)                                                                                                                                  ccTCATGGGCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGcgccgc        >  2:349814/1‑150 (MQ=255)                                                                                                                                  ccTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGcgccgc        >  1:131566/1‑150 (MQ=255)                                                                                                                                      aTGGGCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCAACGGATACATCGGCGCCGCGACt    >  1:411680/1‑150 (MQ=255)                                                                                                                                      aTGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGTGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACt    >  1:645164/1‑150 (MQ=255)                                                                                                                                        ggACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACTgc  >  2:610466/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                              CTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACTGC  >  NC_002947/3812565‑3812849

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑