Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,598,634	A→G	5.3%	R257R (CGT→CGC)	yjiZ ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,598,634	0	A	G	5.3%	201.2 / 2.6	75	R257R (CGT→CGC)	yjiZ	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (40/31);  new base G (2/2);  total (42/33)
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.76e-01

AAGCCAGGCCAGATACAGCCACGCAGTGTAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCATTGACGCTACCTGC  >  W3110S.gb/4598565‑4598702                                                                      |                                                                     aaGCCAGGCCAGATACAGCCACGCAGTGTAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTAc                                                                     >  1:441364/1‑71 (MQ=255) aaGCCAGGCCAGATACAGCCACGCAGTGTAGTTGATGCCAATGAATCCGAGCATCATTCCCCACATTGTAc                                                                     <  1:79471/71‑1 (MQ=255)  aGCCAGGCCAGATACAGCCACGCAGTGTAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACg                                                                    <  1:1695586/71‑1 (MQ=255)          cagATACAGCCACGCAGTGTAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACgg                                                             >  1:2178106/1‑70 (MQ=255)          cagATACAGCCACGCAGTGTAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACgg                                                             >  1:1630132/1‑70 (MQ=255)          cagATACAGCCACGCAGTGTAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACgg                                                             >  1:1601745/1‑70 (MQ=255)              tACAGCCACGCAGTGTAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTa                                                                >  1:3203865/1‑63 (MQ=255)                 aGCCACGCAGTGTAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCt                                                     >  1:1303123/1‑71 (MQ=255)                  gCCACGCAGTGTAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTa                                                                <  1:2528646/59‑1 (MQ=255)                       gCAGTGTAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTGCGGTTACGGAACAAGCTGCGCCa                                               >  1:3724199/1‑71 (MQ=255)                       gCAGTGTAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTGCGGTTACGGAACAAGCTGCGCCa                                               >  1:3891866/1‑71 (MQ=255)                         aGTGTAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCAtt                                             <  1:157666/71‑1 (MQ=255)                          gtgtAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTc                                            <  1:3501570/71‑1 (MQ=255)                             tAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACg                                                              <  1:3682172/50‑1 (MQ=255)                             tAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCgg                                          >  1:1523395/1‑70 (MQ=255)                             tAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCgg                                          >  1:1518283/1‑70 (MQ=255)                             tAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCgg                                          >  1:3486729/1‑70 (MQ=255)                             tAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCgg                                          >  1:3386108/1‑70 (MQ=255)                               gTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGt                                                                  >  1:3536742/1‑44 (MQ=255)                               gTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCaa                                       <  1:2956473/71‑1 (MQ=255)                               gTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCaa                                       <  1:1918853/71‑1 (MQ=255)                               gTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCa                                        >  1:74087/1‑70 (MQ=255)                               gTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCa                                        >  1:766727/1‑70 (MQ=255)                               gTTGATGCCACTGAATCCGAGAATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCa                                        >  1:3895119/1‑70 (MQ=255)                                ttGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCaaa                                      <  1:3402019/71‑1 (MQ=255)                                ttGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCaa                                       <  1:1346211/70‑1 (MQ=255)                                  gATGCCACTGAATCCGAGCATCATTCCCCACATTGTGCGGTTACGGAACAAGCTGCGCCATTCGGCaaaa                                     <  1:1914966/70‑1 (MQ=255)                                  gATGCCACTGAATCCGAGCATCATTCCCCACATTGTGCGGTTACGGAACAAGCTGCGCCATTCGGCaaaa                                     <  1:1541866/70‑1 (MQ=255)                                  gATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTAc                                                               >  1:755330/1‑44 (MQ=255)                                    tGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACt                                   >  1:494361/1‑70 (MQ=255)                                    tGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTg                                  <  1:997344/71‑1 (MQ=255)                                    tGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTg                                  <  1:303756/71‑1 (MQ=255)                                          tGAATCCGAGCATCATTCCCCACATTGTACGGTTa                                                                <  1:3629195/35‑1 (MQ=255)                                          tGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGaa                                                           >  1:30270/1‑40 (MQ=255)                                          tGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGaa                                                           >  1:1215091/1‑40 (MQ=255)                                           gAATCCGAGCATCATTCCCCACATTGTACGGTTACTGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGAt                           >  1:2133647/1‑71 (MQ=255)                                            aaTCCGAGCATCATTCCCCACATTGTACGGTTGCGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGa                            <  1:4113704/69‑1 (MQ=255)                                            aaTCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGAtc                          <  1:2742272/71‑1 (MQ=255)                                            aaTCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGAtc                          >  1:443977/1‑71 (MQ=255)                                              tCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGAtctc                        <  1:654422/71‑1 (MQ=255)                                              tCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGAtctc                        <  1:4136620/71‑1 (MQ=255)                                              tCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGAtctc                        >  1:1233196/1‑71 (MQ=255)                                               ccGAGCATCATTCCCCACATTGTACGGTTACGGAAc                                                          >  1:651011/1‑36 (MQ=255)                                                cGATCATCATTCCCCACATTGTACGGTTACGGAACAGGCt                                                     >  1:2001827/1‑40 (MQ=255)                                                cGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTcgg                      >  1:2134560/1‑71 (MQ=255)                                                cGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTcg                       <  1:707909/70‑1 (MQ=255)                                                     atcatTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCggg                  <  1:390866/70‑1 (MQ=255)                                                       catTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCaa                                       <  1:344397/47‑1 (MQ=255)                                                        atTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGc                                         >  1:163986/1‑44 (MQ=255)                                                        atTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCAtt              <  1:2193337/71‑1 (MQ=255)                                                        atTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCAtt              <  1:1218917/71‑1 (MQ=255)                                                        atTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCAtt              <  1:1898880/71‑1 (MQ=255)                                                        atTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCAtt              >  1:269715/1‑71 (MQ=255)                                                        atTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGAAtt              >  1:783113/1‑71 (MQ=255)                                                         ttCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCAtt              <  1:1377800/70‑1 (MQ=255)                                                         ttCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCATTg             <  1:1739876/71‑1 (MQ=255)                                                         ttCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCATTg             <  1:1096878/71‑1 (MQ=255)                                                            cccACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAAc                                    >  1:2688278/1‑45 (MQ=255)                                                               acaTTGTACGGTTACGGAACAGGCTGCGCCATTCGGCaaaa                                     >  1:4156031/1‑41 (MQ=255)                                                               acaTTGTACGGTTACGGAACAGGCTGCGCCATTCGGCaaaa                                     >  1:420281/1‑41 (MQ=255)                                                               acaTTGTACGGTTACGGAACAGGCTGCGCCATTCGGCaa                                       <  1:3663876/39‑1 (MQ=255)                                                               acaTTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCATTGACGCTa       >  1:279047/1‑71 (MQ=255)                                                                caTTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCATTGACGCTAc      >  1:2141618/1‑71 (MQ=255)                                                                  ttGTACGGTTACGGAACAGGCTGCGCCATTCTGCAAAACTGAGCGGATCTCGGCGGGCATTGACGCTACCt    >  1:4363818/1‑71 (MQ=255)                                                                  ttGTACGGTTACGGAACAGGCTGCGCCATTCGGCaaaa                                     >  1:3255710/1‑38 (MQ=255)                                                                  ttGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGa                                 <  1:4367142/42‑1 (MQ=255)                                                                  ttGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCATTGACGCTAcc     >  1:2837750/1‑70 (MQ=255)                                                                  ttGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCATTGACGCTACCt    >  1:1138612/1‑71 (MQ=255)                                                                   tGTACGGTTACGGAACAGGCTGCGCCATTCGGCaaaa                                     >  1:1666322/1‑37 (MQ=255)                                                                   tGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCAtt              <  1:1944064/60‑1 (MQ=255)                                                                    gTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGc                               >  1:3166273/1‑42 (MQ=255)                                                                    gTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTcggc                     >  1:2411011/1‑52 (MQ=255)                                                                    gTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTcggc                     >  1:840919/1‑52 (MQ=255)                                                                    gTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCATTGACGCTACCTg   <  1:2786564/70‑1 (MQ=255)                                                                     tACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCATTGACGCTACCTGc  <  1:1539528/70‑1 (MQ=255)                                                                      aCGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGa                                 >  1:3651355/1‑38 (MQ=255)                                                                      |                                                                     AAGCCAGGCCAGATACAGCCACGCAGTGTAGTTGATGCCACTGAATCCGAGCATCATTCCCCACATTGTACGGTTACGGAACAGGCTGCGCCATTCGGCAAAACTGAGCGGATCTCGGCGGGCATTGACGCTACCTGC  >  W3110S.gb/4598565‑4598702

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