Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,309,899 G→A 53.3% R31R (CGC→CGT yefM ← antitoxin of the YoeB‑YefM toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,309,8990GA53.3% 5.2 / 79.2 60R31R (CGC→CGTyefMantitoxin of the YoeB‑YefM toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (15/13);  new base A (12/20);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 2.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGATAAGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCTT  >  minE/1309829‑1309969
                                                                      |                                                                      
gTAGATAAGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTgg                                                                        >  1:2392121/1‑71 (MQ=255)
gTAGATAAGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTgg                                                                        <  1:2303759/71‑1 (MQ=255)
  aGATAAGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGc                                                                       >  1:172098/1‑70 (MQ=255)
   gATAAGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGc                                                                       >  1:1392467/1‑69 (MQ=255)
   gATAAGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGc                                                                       >  1:3142405/1‑69 (MQ=255)
        gCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGtaat                                                                <  1:1938994/71‑1 (MQ=255)
        gCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGtaa                                                                 >  1:3373805/1‑70 (MQ=255)
            tctcTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGa                                                                     <  1:616177/62‑1 (MQ=255)
                  ccAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCggg                                                      >  1:2630721/1‑71 (MQ=255)
                    aGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGc                                                    <  1:2671103/71‑1 (MQ=255)
                    aGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGc                                                    <  1:3198811/71‑1 (MQ=255)
                      cGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGAt                                                          <  1:3375388/63‑1 (MQ=255)
                       gAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTACATTCTGGCGAGTAATAAGGATCGGGGCATg                                                 <  1:3128761/71‑1 (MQ=255)
                        aGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGa                                                <  1:692245/71‑1 (MQ=255)
                         gTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGtaata                                                               <  1:3220600/55‑1 (MQ=255)
                         gTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGtaata                                                               <  1:3215993/55‑1 (MQ=255)
                         gTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGa                                                <  1:3035727/70‑1 (MQ=255)
                         gTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGAt                                               <  1:1275924/71‑1 (MQ=255)
                           tGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGAt                                               >  1:1465473/1‑69 (MQ=255)
                           tGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCt                                             >  1:2163195/1‑71 (MQ=255)
                             tATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTc                                           <  1:3616012/71‑1 (MQ=255)
                              aTTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAg                                                             >  1:1547878/1‑52 (MQ=255)
                              aTTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                           >  1:1915203/1‑70 (MQ=255)
                              aTTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                           >  1:329914/1‑70 (MQ=255)
                              aTTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTc                                           >  1:829077/1‑70 (MQ=255)
                               ttcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCaa                                         <  1:2221060/71‑1 (MQ=255)
                                tcttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAAc                                        >  1:3280305/1‑71 (MQ=255)
                                 cttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCAt                                                  >  1:761665/1‑60 (MQ=255)
                                 cttcGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCAt                                                  >  1:2552350/1‑60 (MQ=255)
                                    cGAGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGc                                     <  1:2537929/70‑1 (MQ=255)
                                     gAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCtt                                   <  1:1225046/71‑1 (MQ=255)
                                      aGTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCttt                                  <  1:3131965/71‑1 (MQ=255)
                                       gTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGAt                                                          <  1:334803/46‑1 (MQ=255)
                                       gTGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGAt                                                          <  1:1676780/46‑1 (MQ=255)
                                        tGACATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtc                                 >  1:949143/1‑70 (MQ=255)
                                          aCATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcat                               >  1:3335116/1‑70 (MQ=255)
                                           cATCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcatca                             <  1:137649/71‑1 (MQ=255)
                                            aTCAGATCACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCAt                                                  <  1:156310/49‑1 (MQ=255)
                                             tCAGAACACTAGCCTCTCCATTCTGACGTGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcatc                              <  1:3027107/68‑1 (MQ=255)
                                             tCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcatca                             <  1:239735/69‑1 (MQ=255)
                                             tCAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTcatcat                            <  1:3259485/70‑1 (MQ=255)
                                              cAGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAttg                          >  1:2770689/1‑71 (MQ=255)
                                               aGAACACAAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAttg                          <  1:672221/70‑1 (MQ=255)
                                                 aaCACAAGCCTCTCCATTCTGGCGAGTAATAAGGAt                                                          >  1:1385683/1‑36 (MQ=255)
                                                 aaCACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGc                                                    <  1:1394240/42‑1 (MQ=255)
                                                 aaCACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcatc                              <  1:1317114/64‑1 (MQ=255)
                                                  acacAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGc                                                    <  1:2528690/41‑1 (MQ=255)
                                                  acacAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGc                      >  1:1921347/1‑71 (MQ=255)
                                                   cacaAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCAt                                                  >  1:2849870/1‑42 (MQ=255)
                                                   cacaAGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGcc                     <  1:203904/71‑1 (MQ=255)
                                                     caAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGa                   >  1:226798/1‑71 (MQ=255)
                                                      aaGCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAttgttg                       <  1:1669237/66‑1 (MQ=255)
                                                        gCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTtcatca                             >  1:109406/1‑58 (MQ=255)
                                                        gCCTCTCCATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACa                 >  1:803732/1‑70 (MQ=255)
                                                           tctcCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGcc                     <  1:3465832/63‑1 (MQ=255)
                                                             tcCATTCTGGCGAGTAATAAGGATCGGGGCATGAta                                              >  1:1616365/1‑35 (MQ=255)
                                                             tcCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAAtt             <  1:2021750/69‑1 (MQ=255)
                                                               cATTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGa         >  1:2611727/1‑71 (MQ=255)
                                                                aTTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGAc        >  1:1789573/1‑71 (MQ=255)
                                                                aTTCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCAttgttg                       <  1:3687286/56‑1 (MQ=255)
                                                                 ttCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAAc                                        <  1:3786606/38‑1 (MQ=255)
                                                                  tCTGACGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGa                   >  1:858404/1‑58 (MQ=255)
                                                                      gCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCtt  <  1:693974/71‑1 (MQ=255)
                                                                      |                                                                      
GTAGATAAGCCGTCTCTTCCAGCGAGTTGTATTCTTCGAGTGACATCAGAACACAAGCCTCTCCATTCTGGCGAGTAATAAGGATCGGGGCATGATCTTCAACGGCTTTCATCATTGTTGCCGACAAATTCTGACGCGCTT  >  minE/1309829‑1309969

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