Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,754,163 G→A 9.8% A509T (GCG→ACG)  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,754,1630GA 9.8% 188.3 / 8.8 61A509T (GCG→ACG) yfiQfused predicted acyl‑CoA synthetase NAD(P)‑binding subunit and ATP‑binding subunit
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base A (3/3);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AGTTCAGCCCATCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGAAATTGCGTTCGCCGGATATTCCACATA  >  minE/1754093‑1754232
                                                                      |                                                                     
aGTTCAGCCCATCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAg                                                                       <  1:164919/71‑1 (MQ=255)
  ttCAGCCCATCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCg                                                                     >  1:1908340/1‑71 (MQ=255)
    cAGCCCATCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGt                                                                   >  1:157251/1‑71 (MQ=255)
    cAGCCCATCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGt                                                                   >  1:153258/1‑71 (MQ=255)
     aGCCCATCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTg                                                                  <  1:2838482/71‑1 (MQ=255)
       cccATCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCa                                                                >  1:2966824/1‑71 (MQ=255)
       cccATCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCa                                                                <  1:3356992/71‑1 (MQ=255)
       cccATCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCa                                                                <  1:3360698/71‑1 (MQ=255)
        ccATCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCat                                                               <  1:2270538/71‑1 (MQ=255)
        ccATCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCat                                                               <  1:809576/71‑1 (MQ=255)
          aTCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCatat                                                             <  1:53240/71‑1 (MQ=255)
          aTCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCatat                                                             <  1:2491899/71‑1 (MQ=255)
           tCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCat                                                               <  1:2518848/68‑1 (MQ=255)
           tCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCat                                                               <  1:1745154/68‑1 (MQ=255)
            ccTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTg                                                           >  1:783213/1‑71 (MQ=255)
             cTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAg                                                                       >  1:2749733/1‑58 (MQ=255)
             cTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGc                                                                      <  1:1926885/59‑1 (MQ=255)
             cTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGc                                                          >  1:1345670/1‑71 (MQ=255)
             cTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCA‑‑GATAGCACCGAAg                                                                       >  1:3050930/1‑56 (MQ=255)
              tGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCa                                                                <  1:2922249/64‑1 (MQ=255)
                cAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCata                                                              >  1:1811307/1‑64 (MQ=255)
                cAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCata                                                              >  1:505370/1‑64 (MQ=255)
                cAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCata                                                              >  1:1610653/1‑64 (MQ=255)
                cAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCata                                                              >  1:1366152/1‑64 (MQ=255)
                cAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCg                                                        >  1:2301934/1‑70 (MQ=255)
                   gCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGc                                                                      <  1:539818/53‑1 (MQ=255)
                   gCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAACGGTGCATATTGCCGAACa                                                    <  1:1922287/71‑1 (MQ=255)
                     gTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGaa                                                      >  1:2705860/1‑67 (MQ=255)
                     gTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGaa                                                      >  1:2080784/1‑67 (MQ=255)
                        tGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTg                                               >  1:1878706/1‑71 (MQ=255)
                        tGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTg                                               <  1:1940289/71‑1 (MQ=255)
                        tGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTg                                               >  1:2034114/1‑71 (MQ=255)
                               aaCACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAg                                                                       >  1:1176196/1‑40 (MQ=255)
                               aaCACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATcc                                        >  1:1396585/1‑71 (MQ=255)
                               aaCACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATcc                                        >  1:1740375/1‑71 (MQ=255)
                               aaCACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATc                                         >  1:1568170/1‑70 (MQ=255)
                               aaCACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATc                                         >  1:1680304/1‑70 (MQ=255)
                                 cacGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTg                                               <  1:45013/62‑1 (MQ=255)
                                   cGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCg                                                                     >  1:963825/1‑38 (MQ=255)
                                   cGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACa                                                    >  1:1594469/1‑55 (MQ=255)
                                   cGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACa                                                    >  1:2248264/1‑55 (MQ=255)
                                   cGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTg                                    >  1:1948173/1‑71 (MQ=255)
                                   cGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTg                                    >  1:2912246/1‑71 (MQ=255)
                                    gCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTgg                                   >  1:1622891/1‑71 (MQ=255)
                                           aCCTGGATTGCCAGCGATAGCACCGAAGCGGTGCa                                                                >  1:2536123/1‑35 (MQ=255)
                                            ccTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGa                             <  1:651013/69‑1 (MQ=255)
                                            ccTGGATTGCCAGCGATAGCACCGAAACGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGaaa                           <  1:1602296/71‑1 (MQ=255)
                                                gATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACa                                                    <  1:737623/42‑1 (MQ=255)
                                                 aTTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGAAATTGCg                      <  1:1985515/71‑1 (MQ=255)
                                                  ttGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGt                                     <  1:1448542/55‑1 (MQ=255)
                                                     ccAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGaaa                           >  1:2668958/1‑62 (MQ=255)
                                                        gCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGAAATTGCGTTCGCCg               <  1:2787627/71‑1 (MQ=255)
                                                        gCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGAAATTGCGTTCGCCg               <  1:1094552/71‑1 (MQ=255)
                                                        gCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGAAATTGCGTTCGCCg               <  1:2592398/71‑1 (MQ=255)
                                                        gCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGAAATTGCGTTCGCCg               <  1:3270021/71‑1 (MQ=255)
                                                        gCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGAAATTGCGTTCGCCg               <  1:2793237/71‑1 (MQ=255)
                                                         cGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGAAATTGCGTTCGCCg               <  1:2411067/70‑1 (MQ=255)
                                                         cGATAGCACCGAAACGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGAAATTGCGTTCGCCgg              <  1:1331937/71‑1 (MQ=255)
                                                           aTAGCACCGAAGCGGTGCATATTGCCGAACAGATTg                                               <  1:2998778/36‑1 (MQ=255)
                                                             aGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGAAATTGCGTTCGCCGGatat          >  1:1800414/1‑71 (MQ=255)
                                                                aCCGAAACGGTGCATATTGCCGAACAGATTGGTTATCCGGTgg                                   >  1:2647253/1‑43 (MQ=255)
                                                                aCCGAAACGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGAAATTGCGTTCGCCGGATATTc        >  1:2665051/1‑70 (MQ=255)
                                                                aCCGAAACGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGAAATTGCGTTCGCCGGATATTc        >  1:748769/1‑70 (MQ=255)
                                                                 ccGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGAAAtt                         <  1:1204734/52‑1 (MQ=255)
                                                                     aGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGAAATTGCGTTCGCCGGATATTCCACATa  <  1:1492794/71‑1 (MQ=255)
                                                                      |                                                                     
AGTTCAGCCCATCCTGCAAGCGTATGGCATGAACACGCTCCCTACCTGGATTGCCAGCGATAGCACCGAAGCGGTGCATATTGCCGAACAGATTGGTTATCCGGTGGCGCTGAAATTGCGTTCGCCGGATATTCCACATA  >  minE/1754093‑1754232

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: