Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,812,702	A→G	10.5%	D328G (GAC→GGC)	PP_3368 →	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,812,702	0	A	G	10.5%	152.5 / 8.0	57	D328G (GAC→GGC)	PP_3368	MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base G (4/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

CGCGCCGTGGGCTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACTGC  >  NC_002947/3812554‑3812849                                                                                                                                                     |                                                                                                                                                    cgcgCCGTGGGCTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGAc                                                                                                                                                    >  1:565689/1‑150 (MQ=255)            cTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGa                                                                                                                                         <  1:573861/150‑1 (MQ=255)            cTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGa                                                                                                                                         <  1:780605/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                                                                               >  2:264820/3‑135 (MQ=255)                                                                                     cagggCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTg                                                                               <  1:264820/133‑1 (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                                                        <  1:484110/103‑1 (MQ=39)                                                                                                                                            gCCTCATGGGCAATGTCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGc                                                        >  2:484110/1‑103 (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)                                                                                                                                                     |                                                                                                                                                    CGCGCCGTGGGCTTGGTACTGATGTGGTCGGTCGCCGTGCTGCTGTTCTTTACCTTCGCCGCCCAGGCGCTGTGGAGCCTGCTGCTGGGCATCGGCCTGGTGGGCACCATGATCGCCTTGGCGGCACCGCTACAGATCCGCCTCATGGACATTGCCCATGAGGCGCCTAGCCTGGCAGCGGCGTCCAACCATGCGGCGTTCAACCTGGCCAATGCCCTGGGGCCGTGGCTGGGAGGCATGGCGATCACTGCCGGCATGGGCTGGACCAGCACCGGATACATCGGCGCCGCGACTGC  >  NC_002947/3812554‑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: ‑