Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,331,656	T→A	R385R (CGT→CGA)	proP →	osmolyte:H(+) symporter ProP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,331,656	0	T	A	84.4%	52.8 / 10.2	33	R385R (CGT→CGA)	proP	osmolyte:H(+) symporter ProP
Reads supporting (aligned to +/- strand):  ref base T (0/5);  new base A (10/17);  total (10/23)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-03

CATTATGATCGGTATGCTGTTTGTCCAGCCGGTGATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTC  >  NC_000913/4331410‑4331893                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               cATTATTATCGGTATTCTGTTTGTCCAGCCGGTGATGGGCTGGCTGAGTGACCGTTTGGGCCTTCGTCCGTTTGTGCTACTTGTTAGTGTTGCCCTGTTTGTGTTTGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTATCTGCTTTACGGGCGGTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCAGATCCGATACa                                                                                                                                                                                                                                            <  2:379164/251‑1 (MQ=255)    tATGATCGGTATGCTGTTTGTCCAGCCGGTGATGGGCTTGCTGAGTTTCCGTATTGGCCGTCGTCCGTTTTTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCGTGCCAGCGATGTTCCCGACGCATATCCGATACAgc                                                                                                                                                                                                                                          <  1:383517/250‑1 (MQ=255)          cGGTATGCTGTTTGTCCAGCCGGTGATGGGCTTGCTGAGGGACCGTTTTGGCCGTCGGCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGGTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTGGCCGGGTTACTGATGCTGGCGGTGATCCTTAACGGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCTATGTTCCCGACGCATATCCGTTACAGCGCGCTg                                                                                                                                                                                                                                    <  2:398251/250‑1 (MQ=255)              aTGCTGTTTGTCCAGCCGGTTATTGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGGCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTAGGGCGTCTACCTTGCCAGGGATGTTCCCGACGCATATCCGATACAGCGCGCTggcggc                                                                                                                                                                                                                               <  2:200929/251‑1 (MQ=255)                     ttggccAGCCGTTGATGGGCTTGCTGAGGGGCCGTTTTGTCCGTCGGCCGTTTGTGCTACTTGGTGGTGTGGCTCTGTTTGTGTGGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGAAGCTGGCGGTTCTCCTTATCTGCTTTACGGGCGTTATGGCTTCTCCCTTGCCAGCGCTGTGCCCGACGCATATCCGATTCAGCGCGCTGGCGGCGGCAtt                                                                                                                                                                                                                         <  2:200652/246‑1 (MQ=255)                     ttGTCCAGCCGGTGATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCAttt                                                                                                                                                                                                                        >  1:319473/1‑251 (MQ=255)                       gTCCAGCCGGGGAGGGGCTTGCTGAGTGACCGTTTTGGCCTTCGTCCGTGTGTGCTGCGTGGTAGTGTTGCTCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGAAGCATATCCGATACAGCGCGTTGGCGGCGGCATTTa                                                                                                                                                                                                                       <  2:303893/250‑1 (MQ=255)                             ccGGTGATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTTCCCGCGATGTTCCCGACGCATATCCGATACAGCGCGCCTGCGGCGGCCTTTAATAtttt                                                                                                                                                                                                                >  2:312142/1‑250 (MQ=255)                             ccGGTGATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTc                                                                                                                                                                                                                <  2:6509/251‑1 (MQ=255)                              gggggATGGGCTGGCTGGGTGACCGTTTTGGCCGGCGTCCGTTTGTGCTCCTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCGTCTACCTTGGCAGCGATGTTCCCGACGCATAGCCGATACAGCGCGTTGGCGGCGGCATGTAATATTTCg                                                                                                                                                                                                               <  2:351127/247‑1 (MQ=255)                                gTGATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTGGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTGGCCGGGTTGCTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCGTGCCATCGATGTTCCCGACTCAGATCCGAGACAGCGCGCGGGCGGCGGCATTTAATAGTTCGGt                                                                                                                                                                                                             <  2:183357/251‑1 (MQ=255)                                   aTGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCt                                                                                                                                                                                                          >  1:348280/1‑251 (MQ=255)                                         ttGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTTCCGACGCATATCCGATACAGCGCGCTGGCGGCGGGATTTAATATTTCGGGGCTGGTTGc                                                                                                                                                                                                    >  1:290812/1‑251 (MQ=255)                                          tGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCCACGCATATCCGATCCAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGcc                                                                                                                                                                                                   >  1:57249/1‑251 (MQ=255)                                                              cgtcggcagttggggttACTTGGTAGTGGTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGTCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTggc                                                                                                                                                                              <  2:360875/235‑1 (MQ=255)                                                                        ttGTGCTCCTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCctggct                                                                                                                                                                     <  2:319473/250‑1 (MQ=255)                                                                                 ttGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATTGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCCTTTAATATTTCGGTGCTGGGTGCCGGGTTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGc                                                                                                                                                           >  1:392847/1‑251 (MQ=255)                                                                                     tAGTGTTACCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTTTGGCTTCTACATTTCCCGGGATGTTCCCGACGCATACTCGATAAAACGCGCTGGCGGCGGGCTTTTATATTTCGGTGCTGGGTGCCGGTCTGCCGCCCAAC‑CTggcgggcagggcgggcggaagcccgc                                                                                                                                                       >  2:350320/1‑230 (MQ=255)                                                                                                                     cGTTTATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCAGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATg                                                                                                                       <  1:24921/251‑1 (MQ=255)                                                                                                                       tttATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCCGCGATGTTCCCGCCGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCgggg                                                                                                                                                                                                            >  1:200888/1‑163 (MQ=255)                                                                                                                       tttATTCTGATTAACAGTAACGTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTg                                                                                                                                                                                                            <  2:200888/165‑1 (MQ=255)                                                                                                                               gATTAACATGAACGTCATCGGTCTGATTTTGGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTggagg                                                                                                             <  2:104737/251‑4 (MQ=255)                                                                                                                                            gTCATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATTCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTTGGTGCTGGTTGCCGGTCTGACGCC‑CACGCTGGCGGCCTGGCTGGGCGGAAGCTCGCAGAATCTGATGATGCCTGCCTATTTCCTGATGGTAGTGGCgggtggtggggtttt                                                                                                >  1:17265/1‑238 (MQ=255)                                                                                                                                              cATCGGTCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGGTGCCTGCCTATTACCTGATGGTAGTGGCgggggggggtgttttt                                                                                               >  1:417536/1‑236 (MQ=255)                                                                                                                                                       gttttttGCCGGGTTCCTGATGTTGGCGGTGTTCCTTAACTGCTTGTCGTGCGTTAGGGCTTTTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTGCAGCGCGCTGGTGGCGGGATTTAATATGTCGGTGCTGGTTGCCGGTCTGTCGCC‑AACGCTGGCGGATGGGCTGGTCGAAAGGTCGCAGATTTTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTAGTTTAATCACCGGCGTaa                                                                                     <  2:404643/249‑1 (MQ=255)                                                                                                                                                                                 cGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCg                                                           <  1:402898/251‑1 (MQ=255)                                                                                                                                                                                                     aCGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACcggc                                       <  1:350320/251‑1 (MQ=255)                                                                                                                                                                                                         gCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAAGCACCGGCGTAACCATGAAAGAGACGGTAAAGCGTCCGTTGAAAGGTGCGACACCGgcggcg                                   <  2:57249/251‑1 (MQ=255)                                                                                                                                                                                                            tttgggctgtttccTTGCCAGAGAGTTTCCCGGCGCTTATCCGGTACTGTGTTCTGGCGGCGGCTTTTAATATTTCGGTGCTGGTTTCCGGTTTGTAGCC‑AACGCTGGCGGCTTGGCTGTGCGAAAGCTCGCAGATTCTGATGATGCCTGACTATTACCTGATGGTATTGGCGGTGGGTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACCCCGgtggcggc                                 <  2:318973/238‑3 (MQ=255)                                                                                                                                                                                                                           tGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGc                  <  1:369418/250‑1 (MQ=255)                                                                                                                                                                                                                                  gATGTTCCCGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGaaa          <  1:12644/251‑1 (MQ=255)                                                                                                                                                                                                                                          cGACGCATATCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCTGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGGGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGGCAGATATACAGGAATCGAAGGAAAGGCtcgtc  <  2:272492/251‑1 (MQ=255)                                                                                                                                                                                                                                                atatCCGATACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGa                                                                               <  2:123762/168‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               CATTATGATCGGTATGCTGTTTGTCCAGCCGGTGATGGGCTTGCTGAGTGACCGTTTTGGCCGTCGTCCGTTTGTGCTACTTGGTAGTGTTGCCCTGTTTGTGTTGGCGATCCCGGCGTTTATTCTGATTAACAGTAACGTCATCGGCCTGATTTTTGCCGGGTTACTGATGCTGGCGGTGATCCTTAACTGCTTTACGGGCGTTATGGCTTCTACCTTGCCAGCGATGTTCCCGACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCC‑AACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGTC  >  NC_000913/4331410‑4331893

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