Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,306,816 G→A 6.3% F260F (TTC→TTT yeeF ← predicted amino‑acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,306,8160GA 6.3% 233.1 / 4.1 79F260F (TTC→TTTyeeFpredicted amino‑acid transporter
Reads supporting (aligned to +/- strand):  ref base G (28/45);  new base A (2/3);  total (30/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CTTTACCCGCAACGTACAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATCAGCGCGGTCAGGAAAA  >  minE/1306746‑1306883
                                                                      |                                                                   
cTTTACCCGCAACGTACAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCggg                                                                     <  1:4552386/71‑1 (MQ=255)
cTTTACCCGCAACGTACAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCggg                                                                     <  1:822299/71‑1 (MQ=255)
   tACCCGCAACGTACAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAg                                                                  >  1:4307789/1‑71 (MQ=255)
    aCCCGCAACGTACAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGa                                                                    >  1:2237705/1‑68 (MQ=255)
     cccGCAACGTACAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGt                                                                 <  1:3252854/70‑1 (MQ=255)
     cccGCAACGTACAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGt                                                                 <  1:2705401/70‑1 (MQ=255)
     cccGCAACGTACAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGt                                                                 <  1:5312918/70‑1 (MQ=255)
      ccGCAACGTACAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTAc                                                               <  1:1880987/71‑1 (MQ=255)
        gCAACGTACAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAg                                                             >  1:2175576/1‑71 (MQ=255)
        ccgaCGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAg                                                             >  1:822451/4‑71 (MQ=255)
          gaCGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCt                                                           <  1:2339057/70‑1 (MQ=255)
          gaCGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACAGCt                                                           <  1:5173562/70‑1 (MQ=255)
           aCGTACAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTg                                                          <  1:5211381/71‑1 (MQ=255)
            cGTACAGCGCAATTTCAGGCAGTGCGGCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAACACCt                                                           <  1:5491982/69‑1 (MQ=255)
             gTACAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTg                                                          >  1:932665/1‑69 (MQ=255)
             gTACAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTg                                                          >  1:4251213/1‑69 (MQ=255)
              tACAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTa                                                                >  1:2815789/1‑62 (MQ=255)
                cAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAgg                                                      <  1:5340198/70‑1 (MQ=255)
                cAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGa                                                     <  1:592404/71‑1 (MQ=255)
                   cATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCggg                                                                     <  1:1238806/52‑1 (MQ=255)
                   cATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGaagt                                                  <  1:2356876/71‑1 (MQ=255)
                    aTGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGaagta                                                 <  1:1585551/71‑1 (MQ=255)
                    aTGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGaagta                                                 <  1:3349690/71‑1 (MQ=255)
                     tGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGaagtaa                                                <  1:1787441/71‑1 (MQ=255)
                          tCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTg                                            >  1:2886708/1‑70 (MQ=255)
                          tCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTg                                            >  1:3531869/1‑70 (MQ=255)
                            aGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCaa                                         <  1:837780/71‑1 (MQ=255)
                             ggCTGTGACGCATGCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCaaa                                        <  1:3197313/71‑1 (MQ=255)
                                tgtgACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGaa                                                    >  1:5084163/1‑56 (MQ=255)
                                tgtgACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGaa                                                    >  1:3547518/1‑56 (MQ=255)
                                  tgACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGAt                                     <  1:1297932/69‑1 (MQ=255)
                                  tgACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGAt                                     <  1:1530049/69‑1 (MQ=255)
                                  tgACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGAt                                     <  1:1573416/69‑1 (MQ=255)
                                  tgACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATg                                    <  1:4858311/70‑1 (MQ=255)
                                   gACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGaa                                                                   >  1:5426419/1‑38 (MQ=255)
                                      gCATCCGGATCTTTAAAGCGGCTGATATCGGGAAAGAAc                                                               >  1:5155977/1‑39 (MQ=12)
                                      gCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGaa                                                    >  1:895273/1‑50 (MQ=255)
                                         tCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCaaa                                        <  1:5295727/59‑1 (MQ=255)
                                         tCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCa                                          <  1:1245659/57‑1 (MQ=255)
                                         tCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCa                                          <  1:5386380/57‑1 (MQ=255)
                                         tCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAg                                 <  1:633966/66‑1 (MQ=255)
                                         tCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCAt                            <  1:4635857/71‑1 (MQ=255)
                                           cGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAg                                                             >  1:5350938/1‑36 (MQ=255)
                                           cGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCaaa                                        >  1:4232053/1‑57 (MQ=255)
                                           cGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGa                                   >  1:3518418/1‑62 (MQ=255)
                                           cGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATg                           <  1:1048858/70‑1 (MQ=255)
                                             gATCTTTAAAGCGAGAGATATCCGGGAAGTACAGc                                                            >  1:343508/1‑35 (MQ=255)
                                               tCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATgccgc                       >  1:2417629/1‑70 (MQ=255)
                                                cTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATgccgcc                      >  1:4550932/1‑70 (MQ=255)
                                                 tttAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGAt                                     >  1:3015784/1‑54 (MQ=255)
                                                 tttAAAGAGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGAt                                     >  1:5669920/1‑54 (MQ=255)
                                                   tAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATgccgcc                      <  1:1528210/67‑1 (MQ=255)
                                                   tAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATgccgc                       >  1:3816531/1‑66 (MQ=255)
                                                   tAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATgccgc                       >  1:3553765/1‑66 (MQ=255)
                                                   tAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGccgccg                     <  1:867860/68‑1 (MQ=255)
                                                   tAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGccgccg                     <  1:932968/68‑1 (MQ=255)
                                                   tAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGccgccg                     <  1:4123230/68‑1 (MQ=255)
                                                   tAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGccgccg                     <  1:3588404/68‑1 (MQ=255)
                                                   tAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGa                    <  1:3007535/69‑1 (MQ=255)
                                                   tAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGAt                   <  1:4230961/70‑1 (MQ=255)
                                                   tAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGAt                   <  1:5459875/70‑1 (MQ=255)
                                                   tAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGAt                   <  1:3846685/70‑1 (MQ=255)
                                                   tAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGAt                   >  1:4633378/1‑70 (MQ=255)
                                                    aaaGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATc                  >  1:2415711/1‑70 (MQ=255)
                                                     aaGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGa                                <  1:412353/55‑1 (MQ=255)
                                                     aaGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATCa                 <  1:1473492/70‑1 (MQ=255)
                                                     aaGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATCa                 <  1:4086339/70‑1 (MQ=255)
                                                     aaGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATCAg                <  1:2568976/71‑1 (MQ=255)
                                                     aaGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATCAg                <  1:2851323/71‑1 (MQ=255)
                                                            agaTATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAg                                 >  1:3462770/1‑47 (MQ=255)
                                                                atCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCa                             >  1:4056702/1‑47 (MQ=255)
                                                                atCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCa                             >  1:914147/1‑47 (MQ=255)
                                                                 tCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGAt                               <  1:3336755/44‑1 (MQ=255)
                                                                 tCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATCAGCGCGGt          <  1:5154329/65‑1 (MQ=255)
                                                                 tCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATCAGCGCGGt          <  1:4413776/65‑1 (MQ=255)
                                                                 tCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATCAGCGCGGTCAGGaa    <  1:4170785/71‑1 (MQ=255)
                                                                 tCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATCAGCGCGGTCAGGaa    <  1:4991797/71‑1 (MQ=255)
                                                                   cGGTAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGAt                   <  1:2620159/54‑1 (MQ=255)
                                                                   cGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATCAGCGCGGTCa        >  1:1930623/1‑65 (MQ=255)
                                                                   cGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATCAGCGCGGTCAGGaaaa  <  1:3332064/71‑1 (MQ=255)
                                                                   cGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATCAGCGCGGTCAGGaaaa  <  1:2906619/71‑1 (MQ=255)
                                                                   cGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATCAGCGCGGTCAGGaaa   <  1:2765055/70‑1 (MQ=255)
                                                                     ggAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATCAg                >  1:5307412/1‑55 (MQ=255)
                                                                      |                                                                   
CTTTACCCGCAACGTACAGCATGATTTCAGGCTGTGACGCATCCGGATCTTTAAAGCGAGAGATATCCGGGAAGTACAGCTGCAGGAAGTAAGTTGCAAAGATGAAGATCATGCCGCCGATCAGCGCGGTCAGGAAAA  >  minE/1306746‑1306883

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: