Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,041,581	G→A	18.7%	A374A (GCG→GCA)	ydgR →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,041,581	0	G	A	18.7%	177.0 / 31.3	80	A374A (GCG→GCA)	ydgR	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base A (10/5);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CCAAGTTTGCAATCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCATCGGGGAACT  >  minE/1041511‑1041649                                                                       |                                                                     ccaaGTTTGCAATCGGCATGGTGATGTGTTCTTGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTcgcg                                                                      <  1:705258/71‑1 (MQ=255)     gTTTGCAATCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCt                                                                   >  1:1315402/1‑70 (MQ=255)     gTTTGCAATCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCt                                                                   >  1:2727916/1‑70 (MQ=255)     gTTTGCAATCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCt                                                                   >  1:781990/1‑70 (MQ=255)     gTTTGCAATCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCt                                                                   >  1:3787953/1‑70 (MQ=255)        tGCAATCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTGTGa                                                                 >  1:2703711/1‑69 (MQ=255)        tGCAATCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCATCTGACg                                                               >  1:482132/1‑71 (MQ=255)         gCAATCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACg                                                               <  1:2098362/70‑1 (MQ=255)         gCAATCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACg                                                               <  1:1370328/70‑1 (MQ=255)           aaTCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTcgcg                                                                      <  1:1242200/61‑1 (MQ=255)           aaTCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGACGCTGGGTGCGAAATTCGCGTCTGa                                                                 >  1:687348/1‑66 (MQ=255)            aTCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTgg                                                           >  1:728458/1‑71 (MQ=255)             tCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTTGGTGCGAAATTCGCGTCTGACGCTgg                                                           >  1:36670/1‑70 (MQ=255)             tCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTgg                                                           >  1:2372247/1‑70 (MQ=255)             tCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTgg                                                           >  1:3922609/1‑70 (MQ=255)               ggCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTAt                                                        >  1:2064291/1‑71 (MQ=255)               ggCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTAt                                                        >  1:3537050/1‑71 (MQ=255)               ggCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCATCTGACGCTGGTAt                                                        >  1:3394242/1‑71 (MQ=255)                gCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCATCTGACGCTGGTa                                                         >  1:3182783/1‑69 (MQ=255)                gCATGGTGATGTGTTCTGGTGCGGTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATc                                                       <  1:752770/71‑1 (MQ=255)                 cATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCg                                                      <  1:2339188/71‑1 (MQ=255)                  aTGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCgt                                                     >  1:3874474/1‑71 (MQ=255)                  aTGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCgt                                                     <  1:712157/71‑1 (MQ=255)                  aTGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCgt                                                     >  1:844262/1‑71 (MQ=255)                    ggTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCgtgt                                                   >  1:2368161/1‑71 (MQ=255)                     gTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCgtgt                                                   <  1:3460397/70‑1 (MQ=255)                     gTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCgtgt                                                   <  1:820138/70‑1 (MQ=255)                        aTGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCt                                                                   >  1:3454233/1‑51 (MQ=255)                           tgtTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAg                                            <  1:455980/71‑1 (MQ=255)                            gtTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGa                                                                 >  1:1619853/1‑49 (MQ=255)                            gtTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGa                                                                 >  1:2688796/1‑49 (MQ=255)                            gtTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATc                                                       <  1:1500072/59‑1 (MQ=255)                             ttCTGGTGCGTTCCTGATTCTGACGCTGGGTGCGAAATTCGCATCTGACGCTGGTATCgtgt                                                   <  1:655392/62‑1 (MQ=255)                              tCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAgctg                                         >  1:1008759/1‑71 (MQ=255)                              tCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCATCTGACGCTGGTATCGTGTCTGTAAgctg                                         >  1:196918/1‑71 (MQ=255)                                 ggTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAgctggct                                      >  1:3085552/1‑71 (MQ=255)                                  gTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAgctggctg                                     <  1:627722/71‑1 (MQ=255)                                  gTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAgctggc                                       <  1:3871136/69‑1 (MQ=255)                                  gTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCATCTGACGCTGGTATCGTGTCTGTAAgctggc                                       <  1:2508539/69‑1 (MQ=255)                                   tGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGctggctgg                                    <  1:1445771/71‑1 (MQ=255)                                   tGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGctggctgg                                    <  1:605980/71‑1 (MQ=255)                                   tGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGctggctgg                                    <  1:211075/71‑1 (MQ=255)                                   tGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGctggctgg                                    <  1:610645/71‑1 (MQ=255)                                   tGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCATCTGACGCTGTTATCGTGTCTGTAAGctggctgg                                    <  1:1788554/71‑1 (MQ=255)                                   tGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCATCTGACGCTGGTATCGTGTCTGTAAGctggctgg                                    <  1:1559511/71‑1 (MQ=255)                                   tGCGCTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGAAAGctggctgg                                    <  1:2204557/71‑1 (MQ=255)                                    gCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGt                                               <  1:443501/59‑1 (MQ=255)                                     cGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAgctggc                                       >  1:1171143/1‑66 (MQ=255)                                     cGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGt                                   >  1:551262/1‑70 (MQ=255)                                      gTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCg                                 <  1:2014718/71‑1 (MQ=255)                                       ttCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGc                                <  1:3629038/71‑1 (MQ=255)                                         ccTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGc                                <  1:3839589/69‑1 (MQ=255)                                         ccTGATTCTGCCGCTGGGTGCGAAATTCGCATCTGACGCTGGTATCGTGTCTGTAAgctgg                                        >  1:3012215/1‑61 (MQ=255)                                         ccTGATTCTGCCGCTGGGTGCGAAATTCGCATCTGACGCTGGTATCGTGTCTGTAAgctgg                                        >  1:1411807/1‑61 (MQ=255)                                            gATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTa                                                         >  1:1865642/1‑41 (MQ=255)                                             aTTCTGCCGCTGGGTGCGAAATTCGCATCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTa                          >  1:1893519/1‑71 (MQ=255)                                              ttCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGt                                   <  1:2336148/61‑1 (MQ=255)                                              ttCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTAt                         >  1:3735078/1‑71 (MQ=255)                                              ttCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTAt                         >  1:592877/1‑71 (MQ=255)                                                 tGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATg                        <  1:2155892/69‑1 (MQ=255)                                                  gCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCgt                                                     >  1:3447617/1‑39 (MQ=255)                                                  gCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGctggctgg                                    <  1:1542810/56‑1 (MQ=255)                                                    cGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTg                   >  1:1112580/1‑71 (MQ=255)                                                     gCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGc                  <  1:27244/71‑1 (MQ=255)                                                      cTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCa                 <  1:1828013/71‑1 (MQ=255)                                                      cTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGCCCTGCa                 <  1:510935/71‑1 (MQ=255)                                                      cTGGGTGCGAAATTCGCATCTGACGCTTGTATAGTGTCTGTAAGCTGGCTGGTc                                  >  1:1331440/1‑54 (MQ=255)                                                      cTGGGTGCGAAATTCGCATCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCg                                 >  1:1335962/1‑55 (MQ=255)                                                      cTGGGTGCGAAATTCGCATCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCa                 <  1:2706276/71‑1 (MQ=255)                                                          gTGCGAAATTCGCATCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGc                                >  1:112335/1‑52 (MQ=255)                                                            gCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTa                          <  1:989502/56‑1 (MQ=255)                                                            gCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCAt           >  1:1954076/1‑71 (MQ=255)                                                             cGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCAt           >  1:447643/1‑70 (MQ=255)                                                             cGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCAt           >  1:816440/1‑70 (MQ=255)                                                             cGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCATc          <  1:2471817/71‑1 (MQ=255)                                                             cGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCATc          >  1:3155929/1‑71 (MQ=255)                                                              gAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCATCg         <  1:1311055/71‑1 (MQ=255)                                                                aaTTCGCGTCTGACGCTGGTATCGTGTCTGTAAGctggctgg                                    <  1:473445/42‑1 (MQ=255)                                                                aaTTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCg                                 >  1:465887/1‑45 (MQ=255)                                                                aaTTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCATCggg       >  1:3460915/1‑71 (MQ=255)                                                                  ttCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCATCGGGGa     <  1:49823/71‑1 (MQ=255)                                                                  ttCGCGTCTGACGCTGGTATCGTGTCTGT‑AGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCATCGGGGaa    <  1:1132055/71‑1 (MQ=255)                                                                    cgcgTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCATCg         >  1:56074/1‑65 (MQ=255)                                                                    cgcgTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCATCGGGGAAc   >  1:2976598/1‑71 (MQ=255)                                                                    cgcgTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCATCGGGGAAc   >  1:1243516/1‑71 (MQ=255)                                                                     gggTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCa                 <  1:3675582/54‑1 (MQ=255)                                                                     gcgTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCATCGGGGAACt  >  1:3756486/1‑71 (MQ=255)                                                                     gcgTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCATCGGGGAACt  <  1:18371/71‑1 (MQ=255)                                                                      cgTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCATCGGGGAACt  <  1:1284973/70‑1 (MQ=255)                                                                       |                                                                     CCAAGTTTGCAATCGGCATGGTGATGTGTTCTGGTGCGTTCCTGATTCTGCCGCTGGGTGCGAAATTCGCGTCTGACGCTGGTATCGTGTCTGTAAGCTGGCTGGTCGCAAGCTATGGCCTGCAGAGCATCGGGGAACT  >  minE/1041511‑1041649

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