Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,465 A→G 10.7% L202L (TTA→CTA) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,4650AG10.7% 174.0 / 11.1 65L202L (TTA→CTA) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (28/30);  new base G (5/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTG  >  minE/1015399‑1015535
                                                                  |                                                                      
aGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTaa                                                                        <  1:805717/67‑1 (MQ=255)
 ggCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCtt                                                                   <  1:1332529/71‑1 (MQ=255)
 ggCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCtt                                                                   >  1:1826154/1‑71 (MQ=255)
  gccgccGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCttt                                                                  >  1:468997/1‑71 (MQ=255)
  gccgccGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCttt                                                                  >  1:3175734/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:140574/71‑1 (MQ=255)
                cGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATaga                                                    <  1:1523652/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:421302/1‑71 (MQ=255)
                          ttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCtt                                          >  1:3507257/1‑71 (MQ=255)
                          ttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCtt                                          >  1:2054421/1‑71 (MQ=255)
                          ttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCtt                                          >  1:3767792/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:4005873/71‑1 (MQ=255)
                                                            tCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAgg        <  1:3807248/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:2228545/1‑71 (MQ=255)
                                                                  aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg  <  1:728345/71‑1 (MQ=255)
                                                                  aCGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTg  >  1:66845/1‑71 (MQ=255)
                                                                  |                                                                      
AGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAATTTG  >  minE/1015399‑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: