Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	296,977	G→A	E191K (GAA→AAA)	ykgG →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	296,977	0	G	A	98.6%	273.7 / ‑4.2	70	E191K (GAA→AAA)	ykgG	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (37/32);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTC  >  REL606/296904‑297048                                                                          |                                                                        ccTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTaaa                                                                       >  2:519374/1‑76 (MQ=255)   tGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACg                                                                     <  1:891486/76‑1 (MQ=255)   tGCGTAAAAGCACTACCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACg                                                                     <  1:1395818/76‑1 (MQ=255)    gCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGa                                                                    <  1:12924/76‑1 (MQ=255)     cGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGaa                                                                   <  1:1022631/76‑1 (MQ=255)     cGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGaa                                                                   <  1:1009031/76‑1 (MQ=255)     cGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGaa                                                                   >  1:875442/1‑76 (MQ=255)      gTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAAt                                                                  >  2:1146140/1‑76 (MQ=255)       tAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATg                                                                 >  2:533665/1‑76 (MQ=255)       tAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATg                                                                 >  2:420644/1‑76 (MQ=255)       tAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATg                                                                 >  1:975044/1‑76 (MQ=255)        aaaaGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAGTGc                                                                >  2:106330/1‑76 (MQ=255)        aaaaGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGc                                                                >  1:1824668/1‑76 (MQ=255)           aGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCtt                                                             <  2:1391110/76‑1 (MQ=255)           aGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCtt                                                             <  2:130035/76‑1 (MQ=255)             cACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTcc                                                           >  1:75194/1‑76 (MQ=255)              aCTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCt                                                          >  1:1990515/1‑76 (MQ=255)               cTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTg                                                         >  2:1978686/1‑76 (MQ=255)                 aTCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCa                                                       >  1:1752462/1‑76 (MQ=255)                 aTCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCa                                                       >  1:1094166/1‑76 (MQ=255)                 aTCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCa                                                       >  1:470603/1‑76 (MQ=255)                   ccTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCAtt                                                     >  1:625164/1‑76 (MQ=255)                      gCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAAc                                                  >  1:570896/1‑76 (MQ=255)                      gCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAAc                                                  >  1:120729/1‑76 (MQ=255)                       ccGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAAca                                                 <  2:1935771/76‑1 (MQ=255)                          cgTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAAcatca                                              <  2:1386061/76‑1 (MQ=255)                            tgtAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCAtt                                            >  2:1605731/1‑76 (MQ=255)                             gtAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTa                                           <  1:1776399/76‑1 (MQ=255)                             gtAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTa                                           >  2:1682771/1‑76 (MQ=255)                             gtAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTa                                           <  1:175264/76‑1 (MQ=255)                             gtAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTa                                           <  2:2015582/76‑1 (MQ=255)                              tAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTACACGCATGCCTTCCTGCATTACCATCATTAg                                          >  2:180714/1‑76 (MQ=255)                              tAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAg                                          <  2:584097/76‑1 (MQ=255)                              tAGCGCAACTCGCAGAAAAATTGCAACAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAg                                          <  1:861573/76‑1 (MQ=255)                               aGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGc                                         <  2:1015536/76‑1 (MQ=255)                               aGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGc                                         <  2:1605615/76‑1 (MQ=255)                               aGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGCAAACGAATGCCTTCCTGCATTAACATCATTAGc                                         >  2:788376/1‑76 (MQ=255)                                gcgcAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCg                                        <  1:928508/76‑1 (MQ=255)                                 cgcAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCgg                                       >  2:1710918/1‑76 (MQ=255)                                   cAACTCGCAGAAAAATGGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGcc                                     <  2:414639/76‑1 (MQ=255)                                      cTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCa                                  <  2:1413479/76‑1 (MQ=255)                                        cGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCGTTCCTGCATTAACATCATTAGCGGCCCCCGt                                >  1:364643/1‑76 (MQ=255)                                         gCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGtt                               <  2:1066365/76‑1 (MQ=255)                                          gaGAAAAATTGCATCAGAAAGCGCCGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTc                              <  1:478005/75‑1 (MQ=255)                                          cAGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTc                              <  1:1210273/76‑1 (MQ=255)                                           aGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCa                             <  2:1518448/76‑1 (MQ=255)                                           aGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCa                             >  2:493628/1‑76 (MQ=255)                                           aGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCa                             <  1:1492199/76‑1 (MQ=255)                                           aGAAAAATTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCa                             >  1:1261026/1‑76 (MQ=255)                                              aaaaTTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAAcn                          <  1:248405/76‑2 (MQ=255)                                              aaaaTTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAAcg                          >  2:683459/1‑76 (MQ=255)                                               aaaTTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAAcgg                         <  2:1834697/76‑1 (MQ=255)                                                aaTTGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAAcggc                        >  1:1590025/1‑76 (MQ=255)                                                  ttGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAAcggcgg                      <  2:939055/76‑1 (MQ=255)                                                  tcncATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAAcggcgg                      >  2:1422196/4‑76 (MQ=255)                                                   tGCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGa                     >  2:1999616/1‑76 (MQ=255)                                                    gCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGat                    <  2:1665787/76‑1 (MQ=255)                                                    gCATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGat                    >  2:234108/1‑76 (MQ=255)                                                     cATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGata                   >  1:744221/1‑76 (MQ=255)                                                     cATCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGata                   >  1:2009425/1‑76 (MQ=255)                                                      aTCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGatat                  >  2:1939158/1‑76 (MQ=255)                                                       tCAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATAtt                 <  2:524414/76‑1 (MQ=255)                                                        cAGAAAGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTg                >  2:361271/1‑76 (MQ=255)                                                          gAAAGCGCAGGCCGGTAAACGGATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAg              <  2:117005/76‑1 (MQ=255)                                                           aaaGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGc             >  2:443274/1‑76 (MQ=255)                                                            aaGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGATCt            >  2:874078/1‑76 (MQ=255)                                                            aaGCGCAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCt            >  1:1198662/1‑76 (MQ=255)                                                                 cAGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCa       <  2:1703243/76‑1 (MQ=255)                                                                  aGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCaa      >  1:1174620/1‑76 (MQ=255)                                                                  aGGCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCaa      >  1:999229/1‑76 (MQ=255)                                                                   ggCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCaaa     <  1:572992/76‑1 (MQ=255)                                                                    gCCGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAg    <  1:882329/76‑1 (MQ=255)                                                                      cGGTAAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAgtc  >  1:1389971/1‑76 (MQ=255)                                                                          |                                                                        CCTGCGTAAAAGCACTATCCTGCCGCGTGTAGCGCAACTCGCAGAAAAATTGCATCAGAAAGCGCAGGCCGGTGAACGAATGCCTTCCTGCATTAACATCATTAGCGGCCCCAGTTCAACGGCGGATATTGAGCTTATCAAAGTC  >  REL606/296904‑297048

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑