Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,197,456	C→A	13.8%	K376N (AAG→AAT)	PP_4578 ←	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,197,456	0	C	A	13.8%	153.5 / 15.3	58	K376N (AAG→AAT)	PP_4578	MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (23/27);  new base A (6/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

GCCGCGGGTTTCAGGCAAGCTCAAGGCCGCCACGATCACAATGCCGTAGGACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCAGCCACAGCATCACGCCATCGCTG  >  NC_002947/5197318‑5197600                                                                                                                                           |                                                                                                                                                 gCCGCGGGTTTCAGTCATGCTCAAGGCCGCCACTATCACAATGCCGTAGTACACCTCATAAAACACGCCGATACCCATGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCTGCAATGACCTTTCCGAtg                                                                                                                                         <  4:268611/148‑1 (MQ=255)   cgcgGGTTTCAGGCAAGCTCAAGGCCGCCACGATCACAATGCCGTAGGACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGAtgttg                                                                                                                                      <  1:296236/149‑1 (MQ=255)    gcgGGTTTCAGGCAAGCTCAAGGCCGCCACGATCACAATGCCGTAGGACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGAtgttgt                                                                                                                                     <  7:3471/149‑1 (MQ=255)        gTTTCAGGCAAGCTCAAGGCCGCCACGATCACAATGCCGTAGGACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACa                                                                                                                                 <  7:262047/149‑1 (MQ=255)            cAGGCAAGCTCAAGGCCGCCACGATCACAATGCCGTAGGACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGAtgttgt                                                                                                                                     <  7:42920/141‑1 (MQ=255)            cAGGCAAGCTCAAGGCCGCCACGATCACAATGCCGTAGGACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGAtgttgt                                                                                                                                     >  8:42920/1‑141 (MQ=255)             aGGCAAGCTCAAGGCCGCCACGATCACAATGCCGTAGGACACCGCAGAAAACACGCCGATACCCACGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGc                                                                                                                            <  7:282642/149‑1 (MQ=255)                     tCAAGGCCGCCACGATCACAATGCCGTAGGACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGcccc                                                                                                                    >  4:95046/1‑149 (MQ=255)                        aGGCCGCCACGATCACAATGCCGTAGGACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGccccgaa                                                                                                                 <  7:74322/149‑1 (MQ=255)                            cgccACGATCACAATGCAGTAGGACACCGCAGAAAACACGCAGATACCCAGGCCCAGCGGAACATTCTGGCAGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCTAAGTTGTAACAGAAGCAATGccacgaacccc                                                                                                             >  8:65966/1‑149 (MQ=255)                                     cacaATGCCGTAGGACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGccccgaa                                                                                                                 <  4:2113/136‑1 (MQ=255)                                     cacaATGCCGTAGGACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGccccgaa                                                                                                                 >  3:2113/1‑136 (MQ=255)                                        aaTGCCGTAGGACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGaa                                                                                                                              >  1:303278/1‑120 (MQ=255)                                        aaTGCCGTAGGACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGaa                                                                                                                              <  2:303278/120‑1 (MQ=255)                                         aTGCCGTAGGACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTggg                                                                                                <  3:62568/149‑1 (MQ=255)                                             cGTAGGACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCCAGCATGCCGATCATCAACGGAAACAGCGCGTCAATGACCTTCCCGATGTTGTAACAGAAGCCCTGACCCGAACCCCGAATGCGCGTGGGGaa                                                                                             >  2:275951/1‑148 (MQ=255)                                                 ggACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGtt                                                                                        <  2:5973/149‑1 (MQ=255)                                                     acCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGccc                                                                                                                     >  3:46556/1‑116 (MQ=255)                                                     acCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGccc                                                                                                                     <  4:46556/116‑1 (MQ=255)                                                     acCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGt                                                                                    <  5:159687/149‑1 (MQ=255)                                                     acCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGt                                                                                    <  1:198374/149‑1 (MQ=255)                                                                      gATACCCATGCCCAGCGGAACATTCTGGCCGAGCATGCAGATCCTCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGc                                                                   <  2:308056/149‑1 (MQ=255)                                                                      gATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGc                                                                   <  3:95046/149‑1 (MQ=255)                                                                           ccAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCTCGGCAATGACCTTGCAGATGTTGTAACAGAAGACCTGCCCCGAACCCCGAATGCGCGTGGGGATCAGTTCGCTCAGAAACGACCCCATGCCGcg                                                               >  8:281175/1‑147 (MQ=255)                                                                           ccAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCt                                                               >  7:259609/1‑148 (MQ=255)                                                                               gCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGt                                                                                                   >  6:152507/1‑108 (MQ=255)                                                                               gCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGt                                                                                                   <  5:152507/108‑1 (MQ=255)                                                                                cccAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGaa                                                            >  2:194106/1‑146 (MQ=255)                                                                                cccAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGaa                                                            <  1:194106/146‑1 (MQ=255)                                                                                 ccAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGTGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATa                                                        <  5:133101/149‑1 (MQ=255)                                                                                  cAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATAc                                                       >  6:268935/1‑149 (MQ=255)                                                                                       gAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAg                                                  >  1:270454/1‑149 (MQ=255)                                                                                            ttCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGc                                                                   <  2:192134/127‑1 (MQ=255)                                                                                            ttCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGc                                                                   >  1:192134/1‑127 (MQ=255)                                                                                            ttCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCg                                               >  1:334481/1‑147 (MQ=255)                                                                                               tGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCgaag                                            >  8:7031/1‑147 (MQ=255)                                                                                                ggCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGaagaag                                         >  5:85683/1‑149 (MQ=255)                                                                                                          gCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGt                                                                                    >  8:118995/1‑96 (MQ=255)                                                                                                            cGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGGTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAgg                             <  8:225241/149‑1 (MQ=255)                                                                                                                 atcaACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCa                         <  2:334481/148‑1 (MQ=255)                                                                                                                  tcaACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCAGc                       >  3:129255/1‑149 (MQ=255)                                                                                                                      cGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCAGCca                     >  4:177814/1‑147 (MQ=255)                                                                                                                       ggAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCAGCCACAg                  >  3:335184/1‑149 (MQ=255)                                                                                                                        gAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGATCAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCAGCCACAg                  >  4:167221/1‑148 (MQ=255)                                                                                                                        gAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCAGCCACAGc                 <  5:268935/149‑1 (MQ=255)                                                                                                                        gAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCAGCCACAGc                 <  8:259609/149‑1 (MQ=255)                                                                                                                            aacatcgGCAAGGTCATTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGaa                                                            <  1:110037/97‑1 (MQ=255)                                                                                                                            aacatcgGCAAGGTCATTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGaa                                                            >  2:110037/6‑102 (MQ=255)                                                                                                                            aacatcgGCAAGGTCATTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGa                                                             >  8:111174/6‑101 (MQ=255)                                                                                                                            aacatcgGCAAGGTCATTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGa                                                             <  7:111174/96‑1 (MQ=255)                                                                                                                             aGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCAGCCACAGCATCAc            <  2:270454/149‑1 (MQ=255)                                                                                                                               atcgGCAAGGTCATTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGGAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCAGCCACAGCATCACGc          >  3:116967/3‑149 (MQ=255)                                                                                                                               atcgGCAAGGTCATTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCAGCCCCAGCATCACGc          >  8:333373/3‑149 (MQ=255)                                                                                                                                tcgGCAAGGTCATTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCAGCCACAGCATCACGc          >  6:152124/2‑148 (MQ=255)                                                                                                                                 cgGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCAGCCACAGCATCACGCCa        >  6:47304/1‑149 (MQ=255)                                                                                                                                  ggCAAGGTCATTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCAGCCACAGCATCACGCCAt       >  4:77390/1‑149 (MQ=255)                                                                                                                                    cAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCAGCCACAGCATCACGCCATCg     <  4:335184/149‑1 (MQ=255)                                                                                                                                       tGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCTGTCAGAAACGCCCCCATGCCGCTGAAGATACCCTAGGCGAAGAAGCCCAACGGGAGGCCCAGCCACAGCATCACGCCATCGCTg  <  1:275951/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                 GCCGCGGGTTTCAGGCAAGCTCAAGGCCGCCACGATCACAATGCCGTAGGACACCGCAGAAAACACGCCGATACCCAGGCCCAGCGGAACATTCTGGCCGAGCATGCCGATCATCAACGGAAACAGCGCGGCAATGACCTTGCCGATGTTGTAACAGAAGCCCTGCCCCGAACCCCGAATGCGCGTGGGGAACAGTTCGGTCAGAAACGCCCCCATGCCGCTGAAGATACCCGAGGCGAAGAAGCCCAACGGGAGGCCCAGCCACAGCATCACGCCATCGCTG  >  NC_002947/5197318‑5197600

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