Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,466 C→G 8.5% A201A (GCG→GCC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,4660CG 8.5% 203.3 / 6.8 70A201A (GCG→GCCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base G (4/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTG  >  minE/1015400‑1015535
                                                                  |                                                                     
ggCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCtt                                                                   <  1:1332529/71‑1 (MQ=255)
ggCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCtt                                                                   >  1:1826154/1‑71 (MQ=255)
 gccgccGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCttt                                                                  >  1:3175734/1‑71 (MQ=255)
 gccgccGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCttt                                                                  >  1:468997/1‑71 (MQ=255)
  ccgccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTg                                                                 <  1:2762064/71‑1 (MQ=255)
     ccgTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGtt                                                              <  1:4108094/71‑1 (MQ=255)
       gTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTc                                                            <  1:2342474/71‑1 (MQ=255)
        tACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTc                                                            <  1:253518/70‑1 (MQ=255)
           ccAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTg                                                         <  1:876715/70‑1 (MQ=255)
             aGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGa                                                        >  1:392222/1‑69 (MQ=255)
               cGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATaga                                                    <  1:1523652/71‑1 (MQ=255)
               cGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATaga                                                    <  1:140574/71‑1 (MQ=255)
                 tacgcATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTg                                                         >  1:617274/3‑64 (MQ=255)
                 tGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAgaga                                                  <  1:14176/71‑1 (MQ=255)
                 tGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAgaga                                                  <  1:1466124/71‑1 (MQ=255)
                   cgcATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATa                                                >  1:666274/1‑71 (MQ=255)
                    gcATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATAc                                               >  1:3279619/1‑71 (MQ=255)
                     cATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTc                                                            >  1:69817/1‑57 (MQ=255)
                     cATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTc                                                            >  1:562725/1‑57 (MQ=255)
                       tcttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCt                                                          >  1:306130/1‑57 (MQ=255)
                       tcttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATagag                                                   <  1:2890554/64‑1 (MQ=255)
                       tcttctCAACCAGGTAATTTTTAAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACg                                              >  1:919753/1‑69 (MQ=255)
                         ttctCGACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCtt                                          >  1:1537136/1‑71 (MQ=255)
                         ttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCtt                                          >  1:3767792/1‑71 (MQ=255)
                         ttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCtt                                          >  1:421302/1‑71 (MQ=255)
                         ttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCtt                                          >  1:3507257/1‑71 (MQ=255)
                         ttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCtt                                          >  1:2054421/1‑71 (MQ=255)
                          tctcAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTg                                         >  1:3527936/1‑71 (MQ=255)
                           ctcAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTg                                         <  1:234154/70‑1 (MQ=255)
                               aCCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTg                                                                 >  1:928060/1‑42 (MQ=255)
                               aCCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGc                                    >  1:2576816/1‑71 (MQ=255)
                                ccAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGcc                                   >  1:3564336/1‑71 (MQ=255)
                                 cAGGTAATTTTTCAGTTTCCCCGGTGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGcc                                   >  1:201857/1‑70 (MQ=255)
                                 cAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTc                                            <  1:3335276/61‑1 (MQ=255)
                                  aGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGa                                 <  1:2795270/71‑1 (MQ=255)
                                    gTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGa                                                        <  1:1526825/46‑1 (MQ=255)
                                      aaTTTTTCAGTTTCCCCGGCGTCAGTAGCGCTTTGGt                                                               >  1:3311029/1‑37 (MQ=255)
                                      aaTTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAacca                             >  1:2331983/1‑71 (MQ=255)
                                       gtttttCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTgg                                                        <  1:180025/42‑2 (MQ=16)
                                       aTTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGtt                                                              <  1:3381044/37‑1 (MQ=255)
                                       aTTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccac                            >  1:830267/1‑71 (MQ=255)
                                        tttttCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGttgg                                     >  1:800671/1‑59 (MQ=255)
                                        tttttCAGTTTCCCCGGCGTCAGTAACGCTTTGGtt                                                              <  1:3273351/36‑1 (MQ=255)
                                        tttttCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGt                                        <  1:1571339/58‑1 (MQ=255)
                                         ttttCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTc                                                            >  1:617957/1‑37 (MQ=255)
                                          tttCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATa                                                      >  1:998649/1‑42 (MQ=255)
                                          tttCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccaccac                         >  1:2515483/1‑71 (MQ=255)
                                          tttCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAaccacca                          >  1:2589072/1‑70 (MQ=255)
                                          tttCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTACCTGATAGAGATAAGTCt                                           >  1:2293951/1‑53 (MQ=39)
                                              aGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCttt                     >  1:1511918/1‑71 (MQ=255)
                                               gTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTg                    >  1:3868438/1‑71 (MQ=255)
                                               gTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTg                    <  1:371532/71‑1 (MQ=255)
                                                 ttCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGc                  >  1:3917992/1‑71 (MQ=255)
                                                  tCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGc                  <  1:2411132/70‑1 (MQ=255)
                                                           tCAGTAACGCTTTGGTTTCCTGATAGAGATCCGTCTTGTTTGCCg                                  <  1:3976805/45‑1 (MQ=255)
                                                           tCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGa                                 <  1:1757944/46‑1 (MQ=255)
                                                           tCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGc                  <  1:1308070/61‑1 (MQ=255)
                                                           tCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAgg        <  1:3807248/71‑1 (MQ=255)
                                                           tCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAgg        <  1:4005873/71‑1 (MQ=255)
                                                           tCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGACGa                <  1:1420871/63‑1 (MQ=255)
                                                             aGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGa                                 <  1:2269200/44‑1 (MQ=255)
                                                              gTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCt                       <  1:3892123/53‑1 (MQ=255)
                                                              gTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATa              <  1:495208/62‑1 (MQ=255)
                                                              gTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAt     <  1:1271286/71‑1 (MQ=255)
                                                               tAACGCTTTGGTTTCCTGATAGAGATACGTCTTGttt                                      <  1:460999/37‑1 (MQ=255)
                                                                aaCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAttt   >  1:1381802/1‑71 (MQ=255)
                                                                 aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGaa                                >  1:986943/1‑41 (MQ=255)
                                                                 aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg  <  1:728345/71‑1 (MQ=255)
                                                                 aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg  >  1:2228545/1‑71 (MQ=255)
                                                                 aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg  >  1:66845/1‑71 (MQ=255)
                                                                  |                                                                     
GGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTG  >  minE/1015400‑1015535

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: