Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,227,339	C→T	T369I (ACC→ATC)	ygjI →	putative transporter YgjI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,227,339	0	C	T	85.7%	84.6 / ‑1.3	35	T369I (ACC→ATC)	ygjI	putative transporter YgjI
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (22/8);  total (24/11)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGTTATCTTCCTTGGCTTTG  >  NC_000913/3227098‑3227581                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   ggTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCCACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGAAGGATCTGATGAATACTATTATTAATATGATAGCCGCAgg                                                                                                                                                                                                                                           >  2:368974/1‑250 (MQ=255)          tCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCt                                                                                                                                                                                                                                  >  2:166916/1‑251 (MQ=255)              ttACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTcc                                                                                                                                                                                                                               <  1:368974/250‑1 (MQ=255)                    cgATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGtt                                                                                                                                                                                                                         >  1:235868/1‑250 (MQ=255)                             gTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAAGCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATg                                                                                                                                                                                                               <  1:568337/251‑1 (MQ=255)                               tCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATCCCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTCTTCATCATGCt                                                                                                                                                                                                             >  2:26572/1‑251 (MQ=255)                                                cGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTAc                                                                                                                                                                                            >  1:313629/1‑251 (MQ=255)                                                         ggTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGTTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTGACCTGAATTTACGCGCCAAAt                                                                                                                                                                                   <  2:181437/251‑1 (MQ=255)                                                                    tcttctCCGAAATCCCGTAAGGCATCTTTGGTAAGAAACCCGTCGAACTTAACGAAAACGGCGTTCCGGCTCGCGCAGCGTTGATCCATTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGGTATTCATCATGCTGGCTTACCTGAATTTACGCGCCTAATTAGATCACCt                                                                                                                                                                         <  1:555151/250‑1 (MQ=255)                                                                              aaTCCCGGAAGGCATCTTTGGTAAGAAACCCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCGGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGGGCAGGATCTGATGAATCCTATTATTAATATGATCGCCGCAGCGTCTATGCTGCCGCCGTTATTCATCATTCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGAtt                                                                                                                                                              <  2:235868/251‑1 (MQ=255)                                                                                      aaGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCTTCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATAGGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCGGGCTTACCTGAATTTACGCGACAAATTAGATCACCTGCCACGCGATTTCCGTATg                                                                                                                                                      <  2:251353/251‑1 (MQ=255)                                                                                            tCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTcc                                                                                                                                                >  2:559094/1‑251 (MQ=255)                                                                                                            cGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGtatt                                                                                                                                 >  1:85996/1‑250 (MQ=255)                                                                                                                         gAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAAt                                                                                                                   >  1:409559/1‑251 (MQ=255)                                                                                                                           aaaCGCCGTTCAGGCGCACCCAGCGAGGCTCCCGTTCCTGATCGTCATCACGCTGATGATTATCGCGATGCACGGCTCCAATACCGTGCAGGGTCAGATGGATACTATTATTAATCTGATCGCCGCAGCGTCCATGCTTCCGCAGTTATTCATCATGCTGGCTTCCCTGAATTTACGCGCCAAATTAGATCACCTGCCACGAGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGc                                                                                                                 >  2:31387/1‑251 (MQ=255)                                                                                                                                    tCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGAACCGGTATTATCGTTGTTTCAATGCTGATTg                                                                                                           >  1:301912/1‑248 (MQ=255)                                                                                                                                     ccGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGat                                                                                                       >  1:280889/1‑251 (MQ=255)                                                                                                                                           cgcgcAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGc                                                                                                 >  2:272334/1‑251 (MQ=255)                                                                                                                                                  cggggATCCAGTTCCTGATCGCCACCCCGCTGATGATTATCCCGATGCCCGGTTCCAATCCCGTGAGGGCTCTGATAAATACTATTATTAATATAACCGCCGCAGCGTCCATGCTTCCGCCTTTATTCATCAGGCTGTTTTCCCTAAATTTACGCGCCAAATTAGACCACTTGCCGCGCGATTTCCGAATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATAGTTGCCGTAggg                                                                                          <  1:557441/248‑1 (MQ=255)                                                                                                                                                  cGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGatat                                                                                                     >  1:219861/1‑240 (MQ=255)                                                                                                                                                   gTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGt                                                                                         >  1:366911/1‑251 (MQ=255)                                                                                                                                                              tCCTGATCGTCATCCCTCTGTTGATTATCCCGATGCTCGGTTCCAATCCCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCAGGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGACCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTc                                                                               <  1:202462/250‑1 (MQ=255)                                                                                                                                                                     cGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTcc                                                                       >  1:84701/1‑251 (MQ=255)                                                                                                                                                                      gTCACCCCGCTGATGATTACCCCGAGGCCCGGTTCAAATACCGTGCAGGATCGGATAAATATTATTATTAATAGGTCCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTccc                                                                      <  1:559094/251‑1 (MQ=255)                                                                                                                                                                        cATCCCGCTGATGCTTATCCCGCTGCCCGGCTCCAATCCCGTGCACGATATGATGCATCCTCTTATTATTCTCATCGCCGCAGCGTCCATGCCTCCGCCGTTATTCCTCTTGCTGCCTTACCTGAATTTACGCGCCAACTTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGCTATTATCGTTGTTTCACTGCTGATTGCGCTATTTGCCGTAGGGTTTGTCGCTTCGCCCTTCCCgc                                                                    >  2:286027/1‑250 (MQ=255)                                                                                                                                                                                         tCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTg                                                   >  2:118273/1‑251 (MQ=255)                                                                                                                                                                                           ccGATGCTCGTTTCAAATCCCGTGCAGGATCTGATGAATACTATTATAAATTTGATCGCGGCAGCGTCCGTGCTTCCGCGGTTTTTCATCATGCTGGTTTACCTGAATTTCCGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATGTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGAc                                                 <  1:249234/251‑1 (MQ=255)                                                                                                                                                                                                   cGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCAttt                                         >  1:207803/1‑251 (MQ=255)                                                                                                                                                                                                       cccAATCCCGTGCAGGATCTGATGAATCATATTATTAATATGATCGCCGCAGCGTCCAGGCTTCCCCCGTTATTCATCTTGTTGGCTTACCTGAATTTCCGCGCCAAATTAGACCCCCTGCCACGCGATTTCCGTAGGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCAtttttt                                      <  1:499729/249‑1 (MQ=255)                                                                                                                                                                                                                      gATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGCATATTCTGACCATCATTTTTTATAACGTCGGCGGTAt                      >  2:540765/1‑251 (MQ=255)                                                                                                                                                                                                                             tgaATACTATTATTAATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAAt                                                                                                                                                                                   >  1:547969/1‑87 (MQ=255)                                                                                                                                                                                                                                 tACTATTATTAATATGATCGCCGCAGAGTCAATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTCAATTTACGCGCCACATTAGATCACCTGCCACGCGATTTCCGTCTGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTCATTGCGATATTTGCCGTAGGGTTTGGCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCCTTTTATATaaagtaggatgtgttgttttcttcc           >  2:531847/1‑228 (MQ=255)                                                                                                                                                                                                                                          taATATGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGTTATCTTCCTTGGCTTTg  >  2:239099/1‑251 (MQ=255)                                                                                                                                                                                                                                          taATATGACCGCCGCAGAGTCCATACTTCCGACGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTCTCGTTGGTTCAATGCTGATTGCGATATTAGCCGTAGGGGTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGGTATCTTCCTTGGGtttt  >  2:285858/1‑250 (MQ=255)                                                                                                                                                                                                                                            atatGATCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTa                       >  1:106703/1‑228 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   GGTCTGGTGTCCTTTACCGCGATGTTCGGTTCCCTGCTGATGTGGACCGCAACGCCGGTGAAAATTTTCTTCTCCGAAATCCCGGAAGGCATCTTTGGTAAGAAAACCGTCGAACTGAACGAAAACGGCGTTCCGGCGCGCGCAGCGTGGATCCAGTTCCTGATCGTCATCCCGCTGATGATTATCCCGATGCTCGGTTCCAATACCGTGCAGGATCTGATGAATACTATTATTAATATGACCGCCGCAGCGTCCATGCTTCCGCCGTTATTCATCATGCTGGCTTACCTGAATTTACGCGCCAAATTAGATCACCTGCCACGCGATTTCCGTATGGGCTCCCGCCGCACCGGTATTATCGTTGTTTCAATGCTGATTGCGATATTTGCCGTAGGGTTTGTCGCTTCGACATTCCCGACTGGCGCGAATATTCTGACCATCATTTTTTATAACGTCGGCGGTATTGTTATCTTCCTTGGCTTTG  >  NC_000913/3227098‑3227581

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