Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,794,591 T→C 36.7% A250A (GCA→GCG plsB ← glycerol‑3‑phosphate O‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,794,5910TC36.7% 57.1 / 61.2 98A250A (GCA→GCGplsBglycerol‑3‑phosphate O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (39/23);  new base C (29/7);  total (68/30)
Fisher's exact test for biased strand distribution p-value = 7.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

ATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCT  >  minE/2794532‑2794657
                                                           |                                                                  
aTACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATcgcg                                                           <  1:2493952/69‑1 (MQ=255)
aTACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATcgcg                                                           >  1:3202963/1‑69 (MQ=255)
aTACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATcgcg                                                           <  1:1720096/69‑1 (MQ=255)
aTACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATcgcg                                                           <  1:49226/69‑1 (MQ=255)
   cGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTc                                                        >  1:3266505/1‑69 (MQ=255)
   cGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTc                                                        >  1:1748501/1‑69 (MQ=255)
   cGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTc                                                        >  1:3182644/1‑69 (MQ=255)
   cGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTc                                                        >  1:82901/1‑69 (MQ=255)
   cGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTc                                                        >  1:2514196/1‑69 (MQ=255)
      tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTTCATCAGCGCAATCGCGTTctgc                                                     >  1:2470639/1‑69 (MQ=255)
      tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCTCGTTctgc                                                     >  1:464302/1‑69 (MQ=255)
      tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc                                                     >  1:2257310/1‑69 (MQ=255)
      tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc                                                     >  1:100190/1‑69 (MQ=255)
      tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc                                                     >  1:3386694/1‑69 (MQ=255)
      tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc                                                     >  1:438512/1‑69 (MQ=255)
      tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc                                                     >  1:2990191/1‑69 (MQ=255)
      tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc                                                     >  1:1925577/1‑69 (MQ=255)
      tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc                                                     >  1:1968414/1‑69 (MQ=255)
      tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc                                                     >  1:2166920/1‑69 (MQ=255)
      tcagtcagGCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTctgc                                                     >  1:2286494/1‑69 (MQ=255)
           cagGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc                                                >  1:1868386/1‑69 (MQ=255)
           cagGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc                                                >  1:1867776/1‑69 (MQ=255)
           cagGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc                                                >  1:3166649/1‑69 (MQ=255)
           cagGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc                                                >  1:3222560/1‑69 (MQ=255)
           cagGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc                                                >  1:733123/1‑69 (MQ=255)
           cagGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTgcgc                                                >  1:1179197/1‑69 (MQ=255)
           cagGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTACATCAGTGCAATCGCGTTCTGCTgcgc                                                >  1:2348052/1‑69 (MQ=255)
              gCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCtt                                              <  1:71310/68‑1 (MQ=255)
              gCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCtt                                              <  1:1019032/68‑1 (MQ=255)
              gCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCtt                                              <  1:2503181/68‑1 (MQ=255)
              gCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCtt                                              <  1:595794/68‑1 (MQ=255)
              gCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCtt                                              <  1:1669036/68‑1 (MQ=255)
              gCGGATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCtt                                              <  1:165557/68‑1 (MQ=255)
                gAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttttt                                           <  1:650296/69‑1 (MQ=255)
                gAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttttt                                           <  1:1057630/69‑1 (MQ=255)
                gAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttttt                                           <  1:266845/69‑1 (MQ=255)
                gAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCttttt                                           <  1:2856228/69‑1 (MQ=255)
                 gatcatcTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCtt                                              <  1:1517126/64‑1 (MQ=255)
                   tcatcTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATcgcg                                                           <  1:3316750/50‑1 (MQ=255)
                   tcatcTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCAt                                        >  1:1040728/1‑69 (MQ=255)
                               agaAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAtttttt                            <  1:2283768/69‑1 (MQ=255)
                               agaAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAtttttt                            <  1:428604/69‑1 (MQ=255)
                               agaAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAtttttt                            <  1:727560/69‑1 (MQ=255)
                               agaAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAtttttt                            <  1:3354468/69‑1 (MQ=255)
                               agaAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAtttttt                            <  1:3081675/69‑1 (MQ=255)
                               agaAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAtttttt                            <  1:2073019/69‑1 (MQ=255)
                                 aaaTTCGCCGCAATCTCTTCCATCAGTGCAATCGCGt                                                          >  1:1184999/1‑37 (MQ=255)
                                 aaaTTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTtg                          >  1:1298779/1‑69 (MQ=255)
                                  aaTTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTtgc                         <  1:2496614/69‑1 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:55422/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:1462056/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:2090637/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:1441733/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:519148/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:521480/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:544321/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:2354762/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:2822695/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:419923/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:1700829/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:1044680/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:1030726/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:885602/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:415874/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:906925/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttttttt                           >  1:3402877/1‑64 (MQ=255)
                                     tCGCCGCAATCTCTTCCATCAGCGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTtgcttg                      >  1:1658818/1‑67 (MQ=255)
                                           cAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCtt                <  1:598250/69‑1 (MQ=255)
                                           cAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCtt                <  1:2462645/69‑1 (MQ=255)
                                             aTCTCTTCCATCAGTGCAATCTCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCa              <  1:2468274/69‑1 (MQ=255)
                                             aTCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTgcgcgc                  >  1:1331692/1‑65 (MQ=255)
                                             aTCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTgcgcgc                  >  1:3405799/1‑65 (MQ=255)
                                             aTCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTgcgcgc                  >  1:91018/1‑65 (MQ=255)
                                             aTCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTgcgcgc                  >  1:2862974/1‑65 (MQ=255)
                                             aTCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTgcgcgc                  >  1:1685083/1‑65 (MQ=255)
                                             aTCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTgcgcgc                  >  1:1783491/1‑65 (MQ=255)
                                             aTCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCa              <  1:2523958/69‑1 (MQ=255)
                                             aTCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCa              <  1:142120/69‑1 (MQ=255)
                                             aTCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCa              <  1:2714145/69‑1 (MQ=255)
                                             aTCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCa              <  1:3229517/69‑1 (MQ=255)
                                             aTCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCa              <  1:1677107/69‑1 (MQ=255)
                                             aTCACTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTgcgcgc                  >  1:2642429/1‑65 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGgag                                   >  1:2942785/1‑36 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAttt                               >  1:1518993/1‑40 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCt  >  1:2513534/1‑69 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCt  >  1:2476721/1‑69 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCt  >  1:2377471/1‑69 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCt  >  1:1069959/1‑69 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCt  >  1:1207787/1‑69 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCt  >  1:1617202/1‑69 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCt  >  1:2755309/1‑69 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCt  >  1:2811310/1‑69 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCt  >  1:5504/1‑69 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCt  >  1:2342622/1‑69 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCt  >  1:2035825/1‑69 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCt  >  1:1480693/1‑69 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCt  >  1:416236/1‑69 (MQ=255)
                                                         aGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCt  >  1:1568869/1‑69 (MQ=255)
                                                           tGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGAtttttt                            <  1:2952126/41‑1 (MQ=255)
                                                           |                                                                  
ATACGGTCAGTCAGGCGAATCATCTCGTAAGAGAAATTCGCCGCAATCTCTTCCATCAGTGCAATCGCGTTCTGCTGCGCTTTTTCATGGGAGATTTTTTTGCTGCGCGCTTCATCTTCTACCGCT  >  minE/2794532‑2794657

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: