Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 910,923 G→A 19.1% E55K (GAG→AAG) ‡ ycjV → predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE910,9230GA19.1% 113.0 / 17.5 47E55K (GAG→AAG) ‡ycjVpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (13/25);  new base A (6/3);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

TGTTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATAT  >  minE/910861‑910993
                                                              |                                                                      
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)
             gTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTgaggag                                                                      <  1:3777097/52‑1 (MQ=255)
              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)
                         tAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATc                                                                   >  1:3216752/1‑43 (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:1196501/1‑71 (MQ=255)
                                       tGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGAc                         >  1:3636390/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:2769220/70‑1 (MQ=255)
                                                            aggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa     <  1:2773716/70‑1 (MQ=255)
                                                            aggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa     <  1:2992173/70‑1 (MQ=255)
                                                             ggagATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAata   <  1:1932140/71‑1 (MQ=255)
                                                              gagaTCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAatat  <  1:3074734/71‑1 (MQ=255)
                                                              |                                                                      
TGTTTGTCGGCCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATAT  >  minE/910861‑910993

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: