Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,463 T→C 9.2% L202L (TTA→TTG) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,4630TC 9.2% 180.1 / 9.3 65L202L (TTA→TTG) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base C (4/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAT  >  minE/1015396‑1015532
                                                                   |                                                                     
gaCAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGt                                                                       <  1:1248692/68‑1 (MQ=255)
gaCAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGt                                                                       <  1:3678043/68‑1 (MQ=255)
   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: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:562725/1‑57 (MQ=255)
                         cATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTc                                                         >  1:69817/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:3507257/1‑71 (MQ=255)
                             ttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCtt                                       >  1:421302/1‑71 (MQ=255)
                             ttctCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCtt                                       >  1:3767792/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: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)
                                                                   |                                                                     
GACAGGCCGCCGTACCCAGCGTGCGCATCTTCTCAACCAGGTAATTTTTCAGTTTCCCCGGCGTCAGTAACGCTTTGGTTTCCTGATAGAGATACGTCTTGTTTGCCGAACCACCACCTTTGGCGATACAGAGGAAT  >  minE/1015396‑1015532

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: