Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 910,932 G→A 16.4% G58S (GGC→AGC)  ycjV → predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE910,9320GA16.4% 144.9 / 17.8 55G58S (GGC→AGC) ycjVpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (15/31);  new base A (6/3);  total (21/34)
Fisher's exact test for biased strand distribution p-value = 7.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GTTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGT  >  minE/910862‑911002
                                                                      |                                                                      
gtTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAgcg                                                                        >  1:1511818/1‑71 (MQ=255)
 tttGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATTATTGCCGGGCTTGAGGAGATCAgcg                                                                        <  1:27326/70‑1 (MQ=255)
 tttGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAgcg                                                                        <  1:3433251/70‑1 (MQ=255)
   tGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAgcggcg                                                                     >  1:1242701/1‑71 (MQ=255)
       ggCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGa                                                                 >  1:3324589/1‑71 (MQ=255)
        gCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGa                                                                 <  1:119747/70‑1 (MQ=255)
         accGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATc                                                               <  1:1267345/70‑1 (MQ=255)
          ccGTCTGGCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCt                                                              <  1:1181329/71‑1 (MQ=2)
             tctggCTGCGGTAAATCGCCTTTAATGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGCCCTGtt                                                           >  1:3183041/4‑71 (MQ=2)
               gggCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTg                                                          <  1:3022653/70‑1 (MQ=255)
               gggCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTg                                                          <  1:4011109/70‑1 (MQ=255)
                 gCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGac                                                                >  1:3526274/1‑61 (MQ=16)
                 gCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATc                                                       >  1:1274434/1‑71 (MQ=255)
                 gCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATc                                                       >  1:1992317/1‑71 (MQ=255)
                   tGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGa                                                     >  1:3724427/1‑71 (MQ=255)
                    gCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCggt                                                    >  1:3202758/1‑69 (MQ=255)
                       gTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTg                                                             <  1:2326693/59‑1 (MQ=255)
                       gTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACgg                                                  >  1:213273/1‑70 (MQ=255)
                       gTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGc                                                 <  1:2877984/71‑1 (MQ=255)
                          aGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATc                                                               >  1:1641979/1‑54 (MQ=255)
                          aGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTg                                                             <  1:2261862/56‑1 (MQ=255)
                          aGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTg                                                          <  1:3968957/59‑1 (MQ=255)
                            tCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACg                                            <  1:1871673/71‑1 (MQ=255)
                                ccaccCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatg                                        <  1:3637183/71‑1 (MQ=255)
                                   tacTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCGGTGAGAAACggatga                                       <  1:3949279/67‑1 (MQ=255)
                                   cccTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatgaat                                     <  1:2784703/71‑1 (MQ=255)
                                    ccTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTg                                                          <  1:3178361/49‑1 (MQ=255)
                                    ccTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatgaat                                     <  1:1328184/70‑1 (MQ=255)
                                     cTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCGGTGAGAAACGGatgaatga                                   <  1:3337133/71‑1 (MQ=255)
                                      tGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGAc                                  >  1:819276/1‑71 (MQ=255)
                                      tGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGAc                                  >  1:3636390/1‑71 (MQ=255)
                                      tGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGAc                                  >  1:1196501/1‑71 (MQ=255)
                                        cgcATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGt                                >  1:1972096/1‑71 (MQ=255)
                                        cgcATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGt                                <  1:746103/71‑1 (MQ=255)
                                        cgcATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATCGGTGAGAAACGGATGAATGAcac                                >  1:232109/1‑69 (MQ=255)
                                         gcATGATTGCCGGGCTTGAGGAGATCAgcggcggcg                                                                  <  1:3699612/36‑1 (MQ=255)
                                         gcATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGtt                               <  1:3754410/71‑1 (MQ=255)
                                           aTGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatga                                       <  1:2961061/61‑1 (MQ=255)
                                           aTGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTc                              >  1:1435543/1‑70 (MQ=255)
                                               ttGCCGGGCTTGAGGAGATCAGCGTCGGCGATCTg                                                             >  1:3106202/1‑35 (MQ=255)
                                                 gCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCa                        <  1:1422315/70‑1 (MQ=255)
                                                  ccGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCaa                       <  1:2407201/70‑1 (MQ=255)
                                                  ccGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCaa                       <  1:1330503/70‑1 (MQ=255)
                                                     ggCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACg                                            <  1:3144211/46‑1 (MQ=255)
                                                       cTTGATGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAg                     <  1:1031705/67‑1 (MQ=255)
                                                          gaggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGc                <  1:982836/69‑1 (MQ=255)
                                                          gaggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCa               <  1:246438/70‑1 (MQ=255)
                                                           aggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa              <  1:2992173/70‑1 (MQ=255)
                                                           aggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa              <  1:2773716/70‑1 (MQ=255)
                                                           aggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa              <  1:2769220/70‑1 (MQ=255)
                                                            ggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAata            <  1:1932140/71‑1 (MQ=255)
                                                             gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAatat           <  1:3074734/71‑1 (MQ=255)
                                                              agaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAata            >  1:3769610/1‑69 (MQ=255)
                                                               gaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGtt                               <  1:1529971/49‑1 (MQ=255)
                                                               gaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGtt                               <  1:1534919/49‑1 (MQ=255)
                                                                aTCAGCGGCGGCGATCTGTTGATCGACGGCAAACg                                            <  1:503137/35‑1 (MQ=255)
                                                                aTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGc        >  1:4091666/1‑71 (MQ=255)
                                                                  cAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGa      >  1:2414120/1‑71 (MQ=255)
                                                                  cAGCGGAGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGc        >  1:160553/1‑69 (MQ=255)
                                                                   aGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCAc                  >  1:3388101/1‑58 (MQ=255)
                                                                      ggcggcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGgt  >  1:3100401/1‑71 (MQ=255)
                                                                      |                                                                      
GTTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGT  >  minE/910862‑911002

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