Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,753,067 G→A 27.9% V143V (GTG→GTA yfiQ → fused predicted acyl‑CoA synthetase NAD(P)‑binding subunit and ATP‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,753,0670GA27.9% 64.4 / 26.0 43V143V (GTG→GTAyfiQfused predicted acyl‑CoA synthetase NAD(P)‑binding subunit and ATP‑binding subunit
Reads supporting (aligned to +/- strand):  ref base G (17/14);  new base A (6/6);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGG  >  minE/1753001‑1753136
                                                                  |                                                                     
gCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTACCGa                                                                   >  1:1973452/1‑71 (MQ=255)
 cTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCg                                                                    <  1:3036660/69‑1 (MQ=255)
 cTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGa                                                                   <  1:2454159/70‑1 (MQ=255)
    ggACCAAACAGTCTGGGTTTACTGGCTCCCTGCCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGAtt                                                                 <  1:3513201/69‑1 (MQ=255)
      aCCAAACAGTCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAAc                                                             <  1:3299837/71‑1 (MQ=255)
        cAAACAGTCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAAcgc                                                           >  1:117596/1‑71 (MQ=255)
        cAAACAGTCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTACCGa                                                                   >  1:2130518/1‑63 (MQ=255)
        cAAACAGTCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTACCGATTAAAcgc                                                           >  1:2664600/1‑71 (MQ=255)
          aaCAGTCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCgg                                                         >  1:3877998/1‑71 (MQ=255)
          aaCAGTCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCgg                                                         >  1:2034158/1‑71 (MQ=255)
           aCAGTCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCgg                                                         >  1:236683/1‑70 (MQ=255)
             aGTCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCaa                                                      <  1:1761612/71‑1 (MQ=255)
              gTCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAg                                                     <  1:2234221/71‑1 (MQ=255)
              gTCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTACCGATTAAACGCGGCAAg                                                     <  1:2274366/71‑1 (MQ=255)
               tCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTACCg                                                                    <  1:189365/55‑1 (MQ=255)
                cTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCt                                                   >  1:3646830/1‑71 (MQ=255)
                cTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCt                                                   >  1:935392/1‑71 (MQ=255)
                   ggTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGc                                                <  1:1463202/71‑1 (MQ=255)
                      ttACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGtt                                             >  1:484141/1‑71 (MQ=255)
                      ttACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGtt                                             >  1:1759806/1‑71 (MQ=255)
                       tACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCtgtg                                                                       <  1:1875842/44‑1 (MQ=255)
                            gCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTACCGATTAAACGCGGCAAGCTGGCGTTTAtt                                         <  1:539261/69‑1 (MQ=255)
                            gCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTACCGATTAAACGCGGCAAGCTGGCGTTTATTTc                                       <  1:183813/71‑1 (MQ=255)
                             cTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCg                                      >  1:881026/1‑71 (MQ=255)
                              tCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTg                                                  >  1:3159929/1‑58 (MQ=255)
                              tCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCg                                      >  1:1479939/1‑70 (MQ=255)
                                 cTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAAt                                  >  1:350053/1‑71 (MQ=255)
                                 cTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTACCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAAt                                  <  1:1889190/71‑1 (MQ=255)
                                   ggCAAGGTCTGAATGCCAGCTTTTCGCCTGTACCGATTAAACGCGGCAAGCTGGCg                                               >  1:1115372/1‑56 (MQ=255)
                                    gCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTa                                                                <  1:366714/38‑1 (MQ=255)
                                       aGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCt                             >  1:56571/1‑70 (MQ=255)
                                       aGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTg                            <  1:1319935/71‑1 (MQ=255)
                                        ggTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGc                                                    <  1:3580330/46‑1 (MQ=255)
                                        ggTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCt                             >  1:701764/1‑69 (MQ=255)
                                         gTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTc                                       >  1:1389824/1‑58 (MQ=255)
                                           cTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCg                         >  1:2948122/1‑70 (MQ=255)
                                              aaTGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGtctc                     >  1:331631/1‑71 (MQ=255)
                                                 gCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCa                                    <  1:2413528/53‑1 (MQ=255)
                                                        tttCGCCTGTACCGATTAAACGCGGCAAGCTGGCGTTTa                                           >  1:1351973/1‑39 (MQ=255)
                                                          tCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCTCAATCGGCTGCCGTCTCCAACACCATcc          <  1:2816161/70‑1 (MQ=255)
                                                          tCGCCTGTACCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCt         >  1:1364960/1‑71 (MQ=255)
                                                            gCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTc        <  1:2327211/70‑1 (MQ=255)
                                                            gCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCg       <  1:2162042/71‑1 (MQ=255)
                                                            gCCTGTACCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTc        <  1:2305970/70‑1 (MQ=255)
                                                                 tgCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACt    <  1:3120322/69‑1 (MQ=255)
                                                                 tgCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTgg  <  1:3718716/71‑1 (MQ=255)
                                                                  |                                                                     
GCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTGGCAAGGTCTGAATGCCAGCTTTTCGCCTGTGCCGATTAAACGCGGCAAGCTGGCGTTTATTTCGCAATCGGCTGCCGTCTCCAACACCATCCTCGACTGG  >  minE/1753001‑1753136

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: