Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 509,811 T→A 13.5% intergenic (+36/‑147) ybhI → / → ybhJ predicted transporter/predicted hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE509,8110TA13.5% 127.7 / 12.1 52intergenic (+36/‑147)ybhI/ybhJpredicted transporter/predicted hydratase
Reads supporting (aligned to +/- strand):  ref base T (22/23);  new base A (4/3);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

GTGGAATATGCTGATCGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCT  >  minE/509742‑509879
                                                                     |                                                                    
gtggAATATGCTGATCGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGctgcgtc                                                                     >  1:2255932/1‑71 (MQ=255)
 tggAATATGCTGATCGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGctgcgtct                                                                    >  1:2917766/1‑71 (MQ=255)
 tggAATATGCTGATCGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGctgcgtct                                                                    <  1:5213687/71‑1 (MQ=255)
 tggAATATGCTGATCGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGctgcggct                                                                    <  1:3499539/71‑1 (MQ=255)
   gAATATGCTGATCGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGctgcgtctgc                                                                  >  1:5389750/1‑71 (MQ=255)
   gAATATGCTGATCGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGctgcgtctg                                                                   <  1:4563685/70‑1 (MQ=255)
     atatGCTGATCGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGctgcgtctgcgt                                                                <  1:3158806/71‑1 (MQ=255)
      tatGCTGATCGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGctgcgtctgcgt                                                                <  1:5195992/70‑1 (MQ=255)
         gCTGATCGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGGGTTCGc                                                            >  1:1816370/1‑71 (MQ=255)
         gCTGATCGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGc                                                            >  1:4312324/1‑71 (MQ=255)
            gATCGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGctgcgtct                                                                    <  1:4504083/60‑1 (MQ=255)
             aTCGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAc                                                        >  1:4114385/1‑71 (MQ=255)
              tCGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAcc                                                       <  1:5402810/71‑1 (MQ=255)
               cGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTc                                                              >  1:4279611/1‑63 (MQ=255)
               cGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTc                                                              >  1:1504264/1‑63 (MQ=255)
               cGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAccc                                                      >  1:5059687/1‑71 (MQ=255)
               cGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAccc                                                      >  1:3245123/1‑71 (MQ=255)
               cGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAccc                                                      >  1:1641784/1‑71 (MQ=255)
               cGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAccc                                                      >  1:2338810/1‑71 (MQ=255)
                ggCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAccc                                                      <  1:3040718/70‑1 (MQ=255)
                     gAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGt                                                  <  1:5541444/69‑1 (MQ=255)
                     gAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGt                                                  <  1:3622603/69‑1 (MQ=255)
                     gAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGt                                                  <  1:390914/69‑1 (MQ=255)
                     gAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGt                                                  <  1:5508383/69‑1 (MQ=255)
                     gAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCa                                                <  1:5594135/71‑1 (MQ=255)
                     gAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCa                                                >  1:5298674/1‑71 (MQ=255)
                       aCATGCTGTAAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCact                                              <  1:5525119/71‑1 (MQ=255)
                       aCATGCTGTAAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCact                                              <  1:3786743/71‑1 (MQ=255)
                       aCATGCTGTAAATATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCact                                              <  1:2539823/71‑1 (MQ=255)
                         aTGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCactta                                            <  1:4185844/71‑1 (MQ=255)
                         aTGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCactta                                            <  1:3844019/71‑1 (MQ=255)
                         aTGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCactt                                             >  1:3590426/1‑70 (MQ=255)
                         aTGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAACCCAGTCactta                                            <  1:5693374/71‑1 (MQ=255)
                          tGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAcccc                                                     >  1:3472906/1‑61 (MQ=255)
                            cTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCacttact                                          <  1:5056631/70‑1 (MQ=255)
                               tAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAccc                                                      <  1:1002421/55‑1 (MQ=255)
                               tAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCActtactta                                        <  1:1834590/69‑1 (MQ=255)
                                 aaTATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCg                                                             >  1:5362582/1‑46 (MQ=255)
                                 aaTATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGc                                                            >  1:861187/1‑47 (MQ=255)
                                 aaTATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGc                                                            >  1:257486/1‑47 (MQ=255)
                                 aaTATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAccc                                                      >  1:3156029/1‑53 (MQ=40)
                                 aaTATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAccc                                                      >  1:2259062/1‑53 (MQ=40)
                                 aaTATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGt                                     >  1:3830634/1‑70 (MQ=39)
                                 aaTATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGt                                     <  1:3301971/70‑1 (MQ=39)
                                 aaTATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGt                                     <  1:4915030/70‑1 (MQ=39)
                                    atTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGc                                 >  1:3927941/1‑71 (MQ=34)
                                     ttCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATg                                      <  1:3969571/65‑1 (MQ=33)
                                      tctcGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGt                                                  <  1:4111234/52‑1 (MQ=31)
                                      tctcGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGt                                                  <  1:2971770/52‑1 (MQ=31)
                                                      gTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGc               <  1:3917907/71‑1 (MQ=32)
                                                        tGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCtt             >  1:1800773/1‑71 (MQ=34)
                                                        tGCATGTACTGCGACTGCGTTCGTTCACCCCAGTCACTTACTTATGTAAGCTCCTGGCGATTCACTCGCtt             >  1:2246424/1‑71 (MQ=11)
                                                         gCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCtt             <  1:4398392/70‑1 (MQ=35)
                                                          cATGTACTGCGACTGCGTTCGTTCACCCCAGTCACTTACTTATGTAAGCTCCTGGCGATTCACTCGCTTg            <  1:1322603/70‑1 (MQ=11)
                                                          cATGTACTGCGACTGCGTTCGTTCACCCCAGTCACTTACTTATGTAAGCTCCTGGCGATTCACTCGCTTg            <  1:4728494/70‑1 (MQ=11)
                                                          cATGTACTGCGACTGCGTTCGTTCACCCCAGTCACTTACTTATGTAAGCTCCTGGCGATTCACTCGCTTg            <  1:4247990/70‑1 (MQ=11)
                                                           aTGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTc          >  1:299285/1‑71 (MQ=34)
                                                            tgtaCTGCGACTGCGTTCGTTCACCCCAGTCACTTACTTATGTAAGCTCCTGGCGATTCACTCGCTTGTc          >  1:4104332/1‑70 (MQ=11)
                                                            tgtaCTGCGACTGCGTTCGTTCACCCCAGTCACTTACTTATGTAAGCTCCTGGCGATTCACTCGCTTGTCg         >  1:3145766/1‑71 (MQ=11)
                                                            tgtaCTGCGACTGCGTTCGTTCACCCCAGTCACTTACTTATGTAAGCTCCTGGCGATTCACTCGCTTGTCg         >  1:2072234/1‑71 (MQ=11)
                                                              tgctgcGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCg         >  1:1710652/1‑69 (MQ=35)
                                                               actgcgactgcgtTCGTTCACCCCAGTCACTTACTTATGTAAGCTCCTGGCGATTCACTCGCTTGTCGCCt      >  1:2545224/2‑71 (MQ=11)
                                                                ctgcgtctgcgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCt      <  1:5105963/70‑1 (MQ=35)
                                                                ctgcgtctgcgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCt      <  1:2113211/70‑1 (MQ=35)
                                                                 tgcgtctgcgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGc               >  1:2382228/1‑60 (MQ=33)
                                                                   cgtctgcgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTcc   <  1:2278782/70‑1 (MQ=35)
                                                                   cgtctgcgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCt  >  1:3089674/1‑71 (MQ=34)
                                                                   cgtctgcgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCt  >  1:2018292/1‑71 (MQ=34)
                                                                   cgtctgcgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCt  >  1:4230004/1‑71 (MQ=34)
                                                                     |                                                                    
GTGGAATATGCTGATCGGCTGGAACATGCTGTAAATATTCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCGCTTGTCGCCTTCCT  >  minE/509742‑509879

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: