Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 910,940 T→C 12.1% D60D (GAT→GAC ycjV → predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE910,9400TC12.1% 170.7 / 15.2 66D60D (GAT→GACycjVpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (20/38);  new base C (5/3);  total (25/41)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

CCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAG  >  minE/910871‑911009
                                                                     |                                                                     
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)
        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:819276/1‑71 (MQ=255)
                             tGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGAc                                         >  1:1196501/1‑71 (MQ=255)
                             tGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGAc                                         >  1:3636390/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)
                                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)
                                                              gcggcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGgt         <  1:2941707/70‑1 (MQ=255)
                                                               cggcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGgtgt       >  1:2418846/1‑71 (MQ=255)
                                                                 gcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTc     <  1:3243805/71‑1 (MQ=255)
                                                                 gcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTc     <  1:3927625/71‑1 (MQ=255)
                                                                 gcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTc     <  1:3990331/71‑1 (MQ=255)
                                                                   gATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCaa                     <  1:3073330/53‑1 (MQ=255)
                                                                   gATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGgtgtac     >  1:8839/1‑67 (MQ=255)
                                                                   gATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTc     >  1:3521052/1‑69 (MQ=255)
                                                                   gATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTc     >  1:1429221/1‑69 (MQ=255)
                                                                    aTCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCaaa                             >  1:1208294/1‑44 (MQ=255)
                                                                    aTCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGc                       <  1:301502/50‑1 (MQ=255)
                                                                    aTCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGc                       <  1:312744/50‑1 (MQ=255)
                                                                    aTCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGgtgt       <  1:3423740/66‑1 (MQ=255)
                                                                     tCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCa   >  1:4084645/1‑69 (MQ=255)
                                                                     tCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCa   >  1:293706/1‑69 (MQ=255)
                                                                     tCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAg  >  1:575362/1‑70 (MQ=255)
                                                                     tCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAg  >  1:1009892/1‑70 (MQ=255)
                                                                     |                                                                     
CCCGTCGGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAG  >  minE/910871‑911009

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: