Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 910,946 G→C 8.2% L62F (TTG→TTC ycjV → predicted sugar transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE910,9460GC 8.2% 216.0 / 8.0 73L62F (TTG→TTCycjVpredicted sugar transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (32/35);  new base C (4/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAACTACG  >  minE/910877‑911016
                                                                     |                                                                      
gggCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTg                                                                        <  1:4011109/70‑1 (MQ=255)
gggCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTg                                                                        <  1:3022653/70‑1 (MQ=255)
  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)
        gTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACgg                                                                >  1:213273/1‑70 (MQ=255)
        gTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGc                                                               <  1:2877984/71‑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:819276/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)
                          gcATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGtt                                             <  1:3754410/71‑1 (MQ=255)
                            aTGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACgaatga                                                     <  1:2961061/61‑1 (MQ=255)
                            aTGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTc                                            >  1:1435543/1‑70 (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)
                                               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)
                                                       ggcggcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGgt                >  1:3100401/1‑71 (MQ=255)
                                                        gcggcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGgt                <  1:2941707/70‑1 (MQ=255)
                                                         cggcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGgtgt              >  1:2418846/1‑71 (MQ=255)
                                                           gcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTc            <  1:3927625/71‑1 (MQ=255)
                                                           gcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTc            <  1:3990331/71‑1 (MQ=255)
                                                           gcgATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTc            <  1:3243805/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:312744/50‑1 (MQ=255)
                                                              aTCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGc                              <  1:301502/50‑1 (MQ=255)
                                                              aTCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGgtgt              <  1:3423740/66‑1 (MQ=255)
                                                               tCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCa          >  1:293706/1‑69 (MQ=255)
                                                               tCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCa          >  1:4084645/1‑69 (MQ=255)
                                                               tCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAg         >  1:575362/1‑70 (MQ=255)
                                                               tCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAg         >  1:1009892/1‑70 (MQ=255)
                                                                cTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAata                          >  1:3565876/1‑52 (MQ=255)
                                                                cTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGaa       >  1:3912644/1‑71 (MQ=255)
                                                                cTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGaa       >  1:73728/1‑71 (MQ=255)
                                                                 tGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAata                          >  1:2870513/1‑51 (MQ=255)
                                                                 tGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTcc           >  1:1232187/1‑66 (MQ=255)
                                                                 tGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAAc      >  1:3246020/1‑71 (MQ=255)
                                                                 tGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAAc      >  1:1041211/1‑71 (MQ=255)
                                                                   ttGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGc                      <  1:1237288/53‑1 (MQ=255)
                                                                   ttGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGgtg               >  1:3780583/1‑60 (MQ=255)
                                                                   ttGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGa        <  1:3316595/67‑1 (MQ=255)
                                                                   ttGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAAc      <  1:1021650/69‑1 (MQ=255)
                                                                   ttGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAACTa    >  1:394562/1‑71 (MQ=255)
                                                                   ttGATCGACGGCAAACGAATGAATGACGTTCCAGCAAAAGCACGCAATATAGCGAt                   <  1:2737563/56‑1 (MQ=255)
                                                                    tGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAACTAc   >  1:373944/1‑71 (MQ=255)
                                                                     gATCGACGGCAAACGAATGAATGACGTTCCAGCCa                                      <  1:2623738/35‑1 (MQ=255)
                                                                     gATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAatat                         >  1:2104256/1‑48 (MQ=255)
                                                                     gATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAata                          >  1:2296967/1‑47 (MQ=255)
                                                                     gATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAata                          >  1:995300/1‑47 (MQ=255)
                                                                     gATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAACTAcg  >  1:2921120/1‑71 (MQ=255)
                                                                     gATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAACTAcg  >  1:701210/1‑71 (MQ=255)
                                                                     gATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAACTAcg  >  1:2475327/1‑71 (MQ=255)
                                                                     gATAGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAat                           >  1:361137/1‑46 (MQ=255)
                                                                     |                                                                      
GGGCTGCGGTAAGTCGACCACCCTGCGCATGATTGCCGGGCTTGAGGAGATCAGCGGCGGCGATCTGTTGATCGACGGCAAACGAATGAATGACGTTCCAGCCAAAGCACGCAATATAGCGATGGTGTTCCAGAACTACG  >  minE/910877‑911016

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: