Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,141,700	T→C	7.5%	V12A (GTG→GCG)	PP_4526 →	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,141,700	0	T	C	7.5%	149.8 / 2.3	53	V12A (GTG→GCG)	PP_4526	MFS transporter
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base C (2/2);  total (26/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 = 8.60e-01

TGAATTCCACGGCATTCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTATTTCGTGGGCGCGGTGGACGCC  >  NC_002947/5141552‑5141848                                                                                                                                                      |                                                                                                                                                      tGAATTCCACGGCATTCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGt                                                                                                                                                       >  1:87166/1‑149 (MQ=255) tGAATTCCACGGCATTCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGt                                                                                                                                                       >  2:63788/1‑149 (MQ=255)     ttCCACGGCATTCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGgcg                                                                                                                                                   >  1:229210/1‑149 (MQ=255)        cACGGCATTCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCtg                                                                                                                                                <  2:217574/149‑1 (MQ=255)        cACGGCATTCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCt                                                                                                                                                 >  6:101972/1‑148 (MQ=255)         aCGGCATTCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGTTCGCCAGCGCCGTGGCGCtga                                                                                                                                               >  3:64355/1‑149 (MQ=255)         aCGGCATTCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCtga                                                                                                                                               <  2:87166/149‑1 (MQ=255)            tcATTCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCtgatga                                                                                                                                            <  6:236976/148‑1 (MQ=255)            gCATTCTCCAAAACAGAACAACTTGCTCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCtgatga                                                                                                                                            <  6:29073/149‑1 (MQ=255)            gCATTCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCTTGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCACGGCGCtggcga                                                                                                                                            <  1:217743/149‑5 (MQ=255)              aTTCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGatgatg                                                                                                                                          <  8:122281/149‑1 (MQ=255)               ttCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATgg                                                                                                                                         <  3:496/149‑1 (MQ=255)                    cAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACAGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGAACTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGAGCCGTGGCTCTGATGATGG‑ACATg                                                                                                                                   >  6:249812/1‑149 (MQ=255)                            aacaacTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCgg                                                                                                                           >  5:56731/1‑149 (MQ=255)                                acTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAACTTGCTCGCCAGCGCCGTGGCGCTGATGATGGCCCATGGCCATcatca                                                                                                                         <  1:63788/149‑5 (MQ=255)                                     atcatcATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGc                                                                                                                  <  5:101972/149‑1 (MQ=255)                                        atcatTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGccc                                                                                                                >  5:149707/1‑148 (MQ=255)                                        atcatTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCATACCTGCGCAGTAGACCTCTCTTGTAGCCACTCATACAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGCGCATCGGACGTTTCGCCCt                                                                                                               >  7:128343/1‑149 (MQ=255)                                                tCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCccg                                                                                                        >  3:225316/1‑148 (MQ=255)                                                  aaTTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCGCCccgcaa                                                                                                     >  2:76836/1‑149 (MQ=255)                                                        aTTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCccgcaactgcc                                                                                                <  2:229210/148‑1 (MQ=255)                                                                  aaaCCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCtg                                                                                                                                                <  3:106597/91‑1 (MQ=255)                                                                  aaaCCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCtg                                                                                                                                                >  4:106597/1‑91 (MQ=255)                                                                     ccGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCTTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCg                                                                                  <  8:184753/149‑1 (MQ=255)                                                                      cGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCTACTGCCGCAACTGATTGCCGa                                                                                 <  6:144337/149‑1 (MQ=255)                                                                        cccTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATTG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAgg                                                                               <  5:249812/149‑1 (MQ=255)                                                                        cccTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAgg                                                                               >  7:79496/1‑149 (MQ=255)                                                                                tGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAgg                                                                               >  6:173984/1‑141 (MQ=255)                                                                                tGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAgg                                                                               <  5:173984/141‑1 (MQ=255)                                                                                       tATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGa                                                                                 >  1:225889/1‑132 (MQ=255)                                                                                       tATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGa                                                                                 <  2:225889/132‑1 (MQ=255)                                                                                       tATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGa                                                                >  1:19592/1‑149 (MQ=255)                                                                                       tATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGa                                                                >  3:499/1‑149 (MQ=255)                                                                                        aTCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGa                                                                >  3:257509/1‑148 (MQ=255)                                                                                         tCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGa                                                                <  8:226761/147‑1 (MQ=255)                                                                                            cTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCccgcaactgccgc                                                                                              <  6:47216/114‑1 (MQ=255)                                                                                            cTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCccgcaactgccgc                                                                                              >  5:47216/1‑114 (MQ=255)                                                                                                  tGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGa                                                                     >  2:97170/1‑133 (MQ=255)                                                                                                  tGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCgctggg                                                     >  5:60750/1‑149 (MQ=255)                                                                                                  tGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCCCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGa                                                                     <  1:97170/133‑1 (MQ=255)                                                                                                    cgcAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCgctgggct                                                   <  4:259069/149‑1 (MQ=255)                                                                                                        ggAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTc                                               >  1:36407/1‑149 (MQ=255)                                                                                                         gagaCCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTcg                                              <  2:19592/149‑1 (MQ=255)                                                                                                                   aTGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCaa                                    <  4:225316/149‑1 (MQ=255)                                                                                                                   aTGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCaa                                    >  2:257776/1‑149 (MQ=255)                                                                                                                    tGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAAt                                   <  5:290903/149‑1 (MQ=255)                                                                                                                             tCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTg                             <  5:86368/146‑1 (MQ=255)                                                                                                                             tCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTg                             >  6:86368/1‑146 (MQ=255)                                                                                                                                ttCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTAt                       >  4:120211/1‑149 (MQ=255)                                                                                                                                  cAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCAGCGGCCAATTACCTGGGCTAtt                      >  4:15693/1‑148 (MQ=255)                                                                                                                                       ccatcatcAGCGCCACGGCGCTGGCGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTATTTCGTg                 >  3:95266/8‑149 (MQ=255)                                                                                                                                       ccatcatcAGCGCCACGGCGCTGGCGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTATTTCGTg                 <  4:95266/142‑1 (MQ=255)                                                                                                                                         cTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCAGCGGCCAATTACCTGGGCTATTTCGTg                 <  3:15693/147‑1 (MQ=255)                                                                                                                                            atcAGCGCCACGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTATTTCGTGGgcgc            >  4:173129/3‑149 (MQ=255)                                                                                                                                                  gcCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTATTTCGTGGGCGCGGTGGa      <  4:257509/149‑1 (MQ=255)                                                                                                                                                      tGGCGATGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCGCCCCGCAACTGCCTCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTATTTCGTGGGCGCGGTGGACGcc  <  1:76836/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      TGAATTCCACGGCATTCTCCAAAACAGAACAACTTGATCATCATTATTCAATTTTATTGAAAAATAAACCGCCCTAACCTGCCATCTATCGCTTACCTGCGCAGGAGACCTGTCATGTCGCCACTCATTCAAC‑TGCTCGCCAGCGCCGTGGCGCTGATGATGG‑CCATGGGCATCGGCCGTTTCGCCCTCACCCCGCAACTGCCGCAACTGATTGCCGAGGGCCAGTTCGACCTGACAGTCGCTGGGCTGGTCGCGGCGGCCAATTACCTGGGCTATTTCGTGGGCGCGGTGGACGCC  >  NC_002947/5141552‑5141848

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