Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,921,110 T→C 35.9% S648G (AGC→GGC)  aas ← fused 2‑acylglycerophospho‑ethanolamine acyl transferase and acyl‑acyl carrier protein synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,921,1100TC35.9% 84.3 / 89.1 92S648G (AGC→GGC) aasfused 2‑acylglycerophospho‑ethanolamine acyl transferase and acyl‑acyl carrier protein synthetase
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base C (17/16);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

GCAACTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGC  >  minE/1921041‑1921177
                                                                     |                                                                   
gCAACTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGctc                                                                    >  1:73257/1‑71 (MQ=255)
gCAACTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGctc                                                                    >  1:3025043/1‑71 (MQ=255)
gCAACTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGct                                                                     <  1:3024600/70‑1 (MQ=255)
 cAACTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATctctta                                                                   <  1:2405209/71‑3 (MQ=255)
  aaCTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCtt                                                                  <  1:2344627/71‑1 (MQ=255)
  aaCTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCtt                                                                  <  1:1086587/71‑1 (MQ=255)
   aCTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCtt                                                                  >  1:1340861/1‑70 (MQ=255)
   aCTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCtt                                                                  >  1:3677490/1‑70 (MQ=255)
   aCTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTa                                                                 <  1:2381109/71‑1 (MQ=255)
    cTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTaa                                                                >  1:1995249/1‑71 (MQ=255)
    cTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTaa                                                                >  1:3613381/1‑71 (MQ=255)
    cTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTaa                                                                <  1:505310/71‑1 (MQ=255)
       aTCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCg                                                             >  1:3430739/1‑71 (MQ=255)
       aTCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCg                                                             >  1:802348/1‑71 (MQ=255)
          gcgcgTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCa                                                           <  1:876166/70‑1 (MQ=255)
          gcgcgTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAg                                                          >  1:2813440/1‑71 (MQ=255)
            gcgTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTg                                                        >  1:711127/1‑71 (MQ=255)
                cAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGc                                                      <  1:1219735/69‑1 (MQ=255)
                cAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCAt                                                    >  1:2689556/1‑71 (MQ=255)
                 aGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCa                                                     <  1:79006/69‑1 (MQ=255)
                 aGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATg                                                   >  1:174723/1‑71 (MQ=255)
                    tCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTa                                                                 <  1:821568/54‑1 (MQ=255)
                    tCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTa                                                                 <  1:1654462/54‑1 (MQ=255)
                    tCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTa                                                                 <  1:1727696/54‑1 (MQ=255)
                      gTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAAt                                                               >  1:84397/1‑54 (MQ=255)
                      gTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATg                                                   >  1:648093/1‑66 (MQ=255)
                      gTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACt                                               <  1:3392458/70‑1 (MQ=255)
                         aTCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACt                                               >  1:3153158/1‑67 (MQ=255)
                         aTCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTAt                                           >  1:553816/1‑71 (MQ=255)
                          tCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTAt                                           <  1:2546471/70‑1 (MQ=255)
                             gtggtgAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCa                                                     <  1:559526/57‑1 (MQ=255)
                             gtggtgAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCa                                                     <  1:3231340/57‑1 (MQ=255)
                              tggtgAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTg                                        <  1:2677465/69‑1 (MQ=255)
                              tggtgAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCAt                                                    >  1:1473424/1‑57 (MQ=255)
                              tggtgAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATc                                          >  1:276060/1‑67 (MQ=255)
                              tggtgAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATc                                          >  1:116108/1‑67 (MQ=255)
                              tggtgAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTg                                        <  1:2562495/69‑1 (MQ=255)
                              tggtgAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTg                                        <  1:1698049/69‑1 (MQ=255)
                                  gAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGaaa                                  >  1:217460/1‑71 (MQ=255)
                                   aaaaGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCg                                                             <  1:3809105/43‑1 (MQ=255)
                                   aaaaGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTg                                        <  1:543168/64‑1 (MQ=255)
                                   aaaaGCACCAGTGCCTCGCCTTTCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTg                                        <  1:1173012/63‑1 (MQ=255)
                                    aaaGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAaca                                >  1:353833/1‑71 (MQ=255)
                                    aaaGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAaca                                >  1:3102209/1‑71 (MQ=255)
                                    aaaGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAaca                                >  1:1840098/1‑71 (MQ=255)
                                    aaaGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAac                                 <  1:1316407/70‑1 (MQ=255)
                                    aaaGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTCAATCGCAGTGGCATGGACTTTATCTGCCGAAac                                 <  1:1589723/70‑1 (MQ=255)
                                      aGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACAcc                              <  1:1992213/71‑1 (MQ=255)
                                      aGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACAcc                              <  1:1308268/71‑1 (MQ=255)
                                       gCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCa                             <  1:3509101/71‑1 (MQ=255)
                                       gCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCa                             <  1:374042/71‑1 (MQ=255)
                                        cACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAAt                                                               >  1:1237405/1‑36 (MQ=255)
                                        cACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAAt                                                               >  1:3711633/1‑36 (MQ=255)
                                        cACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACAcc                              >  1:3563187/1‑69 (MQ=255)
                                           cAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCaa                            <  1:2398951/68‑1 (MQ=255)
                                           cAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTg                         >  1:792197/1‑71 (MQ=255)
                                           cAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTg                         >  1:1831719/1‑71 (MQ=255)
                                            aGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGc                        <  1:256560/71‑1 (MQ=255)
                                            aGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGc                        <  1:1462917/71‑1 (MQ=255)
                                             gTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCATTGGCATGGACTTTATCTg                                        <  1:629658/54‑1 (MQ=255)
                                             gTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGc                        >  1:1080390/1‑70 (MQ=255)
                                             gTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGc                        >  1:2889844/1‑70 (MQ=255)
                                             gTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTGGCGa                                    <  1:2691287/58‑1 (MQ=255)
                                             gTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGcc                       <  1:3596863/71‑1 (MQ=255)
                                               gCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAg                     >  1:1518789/1‑71 (MQ=255)
                                                ccTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAgt                    <  1:1637560/71‑1 (MQ=255)
                                                  tCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAgttg                  >  1:629014/1‑71 (MQ=255)
                                                  tCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTg                                        >  1:602468/1‑49 (MQ=255)
                                                   cGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGTATGGACTTTATCTGGCGAAACACCAAGTGCCAgttg                  >  1:1611417/1‑70 (MQ=255)
                                                   cGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTgg                                       >  1:2034996/1‑49 (MQ=255)
                                                   cGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACAcc                              <  1:3211853/58‑1 (MQ=255)
                                                   cGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCa                             <  1:1115037/59‑1 (MQ=255)
                                                    gCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAgttgtt                <  1:103263/71‑1 (MQ=255)
                                                    gCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGAc                                                >  1:1653020/1‑39 (MQ=255)
                                                    gCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACttt                                             <  1:3491413/42‑1 (MQ=255)
                                                    gCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAgttgtt                >  1:3126994/1‑71 (MQ=255)
                                                     ccTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAgttgtt                <  1:2694214/70‑1 (MQ=255)
                                                     ccTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTc               >  1:336072/1‑71 (MQ=255)
                                                      cTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGa                                                 >  1:2367919/1‑36 (MQ=255)
                                                      cTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTc               >  1:1185614/1‑70 (MQ=255)
                                                       tttGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTa                                            <  1:3772516/40‑1 (MQ=255)
                                                       tttGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTcca             >  1:1027569/1‑71 (MQ=255)
                                                       tttGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTcca             >  1:1023091/1‑71 (MQ=255)
                                                        ttGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTccac            >  1:3491708/1‑71 (MQ=255)
                                                         tGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAaca                                >  1:1101130/1‑50 (MQ=255)
                                                           cTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCAt         <  1:2696218/71‑1 (MQ=255)
                                                            tGGCATCGCTCTTAATCGCAGTGGCATGGACTTTAt                                           <  1:2581145/36‑1 (MQ=255)
                                                              gCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAg                           <  1:3799119/50‑1 (MQ=255)
                                                              gCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTc      <  1:3817448/71‑1 (MQ=255)
                                                              gCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTc      <  1:3268943/71‑1 (MQ=255)
                                                              gCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTACCACCATTTc      >  1:3637485/1‑71 (MQ=255)
                                                               cATCGCCCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACAcc                              <  1:3651295/46‑1 (MQ=255)
                                                                aTCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTg                         >  1:1238023/1‑50 (MQ=255)
                                                                   gCTCTTAATCTCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTg                         <  1:1145470/47‑1 (MQ=255)
                                                                   gCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTg                         <  1:426780/47‑1 (MQ=255)
                                                                   gCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGc  <  1:580899/70‑1 (MQ=255)
                                                                     |                                                                   
GCAACTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGC  >  minE/1921041‑1921177

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: