Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 910,930 G→C 16.3% S57T (AGC→ACC)  ycjV → predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE910,9300GC16.3% 137.8 / 16.2 55S57T (AGC→ACC) ycjVpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (15/31);  new base C (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 = 1.00e+00

TGTTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGA  >  minE/910861‑910998
                                                                     |                                                                    
tgtTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAgc                                                                     <  1:2096284/71‑1 (MQ=255)
 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:1992317/1‑71 (MQ=255)
                  gCTGCGGTAAATCGACTTTACTGCGCATGATTGCCGGGCTTGAGACGATCACCAGCGGCGACCTGTTCATc                                                   >  1:1274434/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:3636390/1‑71 (MQ=255)
                                       tGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGAc                              >  1:1196501/1‑71 (MQ=255)
                                       tGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGAc                              >  1:819276/1‑71 (MQ=255)
                                         cgcATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGt                            <  1:746103/71‑1 (MQ=255)
                                         cgcATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGt                            >  1:1972096/1‑71 (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:1330503/70‑1 (MQ=255)
                                                   ccGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCaa                   <  1:2407201/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:1534919/49‑1 (MQ=255)
                                                                gaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGtt                           <  1:1529971/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)
                                                                     |                                                                    
TGTTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGA  >  minE/910861‑910998

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: