Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,057,586 C→T 18.3% D137D (GAC→GAT ttdA → L‑tartrate dehydratase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,057,5860CT18.3% 158.4 / 27.0 71D137D (GAC→GATttdAL‑tartrate dehydratase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (30/28);  new base T (7/6);  total (37/34)
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.45e-01

CGGTAGAAATTTTTGACGAAGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTGCACGCTACCTGG  >  minE/2057516‑2057654
                                                                      |                                                                    
cGGTAGAAATTTTTGACGAAGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGAc                                                                      <  1:443204/71‑1 (MQ=255)
  gTAGAAATTTTTGACGAAGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGAcat                                                                    <  1:151021/71‑1 (MQ=255)
  gTAGAAATTTTTGACGAAGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGAcat                                                                    <  1:3377892/71‑1 (MQ=255)
   tAGAAATTTTTGACGAAGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGAcatc                                                                   <  1:16303/71‑1 (MQ=255)
    aGAAATTTTTGACGAAGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGAcatca                                                                  >  1:1532465/1‑71 (MQ=255)
    aGAAATTTTTGACGAAGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGAcatca                                                                  >  1:3371715/1‑71 (MQ=255)
       aaTTTTTGACGAAGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATcc                                                               >  1:2242093/1‑71 (MQ=255)
        aTTTTTGACGAAGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATccc                                                              >  1:3163398/1‑71 (MQ=255)
                 gAAGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGGtatcat                                                                 >  1:2772738/1‑59 (MQ=255)
                 gAAGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGATATCATCCCCGACAATGa                                                     >  1:2078265/1‑71 (MQ=255)
                 gAAGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGATATCATCCCCGACAATGa                                                     >  1:3877462/1‑71 (MQ=255)
                 gAAGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAAt                                                       <  1:1006417/69‑1 (MQ=255)
                  aaGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGAc                                                    <  1:1753181/71‑1 (MQ=255)
                   aGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGAc                                                    >  1:875816/1‑70 (MQ=255)
                    gTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGa                                                  >  1:2685767/1‑71 (MQ=255)
                    gTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGa                                                  >  1:2681307/1‑71 (MQ=255)
                    gTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGa                                                  >  1:1589368/1‑71 (MQ=255)
                       aaCACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGATATCATCCCCGACAATGACGATGc                                               >  1:3484009/1‑71 (MQ=255)
                       aaCACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATg                                                >  1:3770979/1‑70 (MQ=255)
                       aaCACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGc                                               <  1:1526096/71‑1 (MQ=255)
                       aaCACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGc                                               >  1:1395241/1‑71 (MQ=255)
                       aaCACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGc                                               <  1:1535222/71‑1 (MQ=255)
                         cacCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGATATCATCCCCGACAATGACGATGCgg                                             <  1:3217085/71‑1 (MQ=255)
                          acCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGATATCATCCCCGACAATGACGATGCgg                                             >  1:2754577/1‑70 (MQ=255)
                                aaaaaTACCGGTAGCGGCGTACCGTGGGTCACCTGGGATATCATCCCCGACAATGACGATGCGGAAATCGa                                      <  1:3837916/71‑1 (MQ=255)
                                aaaaaTACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATcccc                                                             <  1:1491241/48‑1 (MQ=255)
                                 aaaaTACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGaa                                     <  1:587107/71‑1 (MQ=255)
                                 aaaaTACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGaa                                     <  1:3353877/71‑1 (MQ=255)
                                  aaaTACCGGTAGCGGCGTACCGTGGGTCACCTGGGACatcatc                                                                >  1:2559559/1‑43 (MQ=255)
                                  aaaTACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCg                                       >  1:3029287/1‑68 (MQ=255)
                                   aaTACCGGTAGCGGCTTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGt                                   <  1:3443623/71‑1 (MQ=255)
                                   aaTACCGGTAGCGGCGTACCGTGGGTCACCTGGGATATCATCCCCGACAATGACGATGCGGAATTCGAAGt                                   <  1:1868175/71‑1 (MQ=255)
                                   aaTACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGt                                   >  1:2769251/1‑71 (MQ=255)
                                    aTACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGtt                                  >  1:841813/1‑71 (MQ=255)
                                         ggTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACa                              >  1:2489920/1‑70 (MQ=255)
                                          gTAGCGGCGTACCGTGGGTCACCTGGGATATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATg                            >  1:2364124/1‑71 (MQ=255)
                                          gTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACaa                                                        >  1:1666257/1‑43 (MQ=255)
                                          gTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGa                                                     >  1:3477514/1‑46 (MQ=255)
                                          gTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGa                                                     >  1:2752560/1‑46 (MQ=255)
                                          gTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACg                                                   <  1:1995158/48‑1 (MQ=255)
                                          gTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACg                                                   <  1:3485036/48‑1 (MQ=255)
                                          gTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATg                            >  1:1655358/1‑71 (MQ=255)
                                           tAGCGGCGTACCGTGGGTCACCTGGGATATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATgg                           <  1:1003747/71‑1 (MQ=255)
                                           tAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATgg                           <  1:33795/71‑1 (MQ=255)
                                             gcggcgTACCGTGGGTCACCTGGGACATCATCCCCGa                                                           >  1:2412696/1‑37 (MQ=255)
                                             gcggcgTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCa                         <  1:1685163/71‑1 (MQ=255)
                                                gcgTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTAc                               <  1:1606771/62‑1 (MQ=255)
                                                  gTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGa                                      <  1:3157787/53‑1 (MQ=255)
                                                  gTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGa                                      <  1:535677/53‑1 (MQ=255)
                                                  gTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGa                                      <  1:1581983/53‑1 (MQ=255)
                                                  gTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAggcgg                    >  1:1800432/1‑71 (MQ=255)
                                                  gTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAggcgg                    <  1:3461301/71‑1 (MQ=255)
                                                     ccGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAggcggcgg                 >  1:288371/1‑71 (MQ=255)
                                                       gTGGGTCACCTGGGATATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCt               >  1:708804/1‑71 (MQ=255)
                                                       gTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACa                              >  1:537235/1‑56 (MQ=255)
                                                        tGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTg              <  1:910924/71‑1 (MQ=255)
                                                          ggTCACCTGGGACATCATCCCCGACAATGACGATGCGGa                                            <  1:2114545/39‑1 (MQ=255)
                                                          ggTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTGc             <  1:1096458/70‑1 (MQ=255)
                                                          ggTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTGCa            <  1:1382298/71‑1 (MQ=255)
                                                            tCACCTGGGATATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTGCACg          <  1:454669/71‑1 (MQ=255)
                                                            tCACCTGGGACATCATCCCCGACAATGACGATGCGGa                                            <  1:1456607/37‑1 (MQ=255)
                                                            tCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCg                                       <  1:637029/42‑1 (MQ=255)
                                                            tCACCTGGGACATCATCACCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTGCACg          >  1:1274148/1‑71 (MQ=255)
                                                              aCCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTGCACGc         >  1:2302100/1‑70 (MQ=255)
                                                                cTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAg                                    >  1:1193979/1‑41 (MQ=255)
                                                                cTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTGCa            >  1:2286716/1‑65 (MQ=255)
                                                                cTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTGCACGCTAc      <  1:1227953/71‑1 (MQ=255)
                                                                  gggACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTa                                >  1:2569956/1‑43 (MQ=255)
                                                                  gggACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACa                              >  1:1768215/1‑45 (MQ=255)
                                                                  gggACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACa                              >  1:1728416/1‑45 (MQ=255)
                                                                  gggACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACa                              >  1:1365444/1‑45 (MQ=255)
                                                                  gggACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACAt                             <  1:1306797/46‑1 (MQ=255)
                                                                  gggACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACAt                             <  1:3549069/46‑1 (MQ=255)
                                                                  gggACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAggcggc                   <  1:1016804/56‑1 (MQ=255)
                                                                  gggACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTg              >  1:1907757/1‑61 (MQ=255)
                                                                   ggATATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTGCACGCTAcc     <  1:683603/69‑1 (MQ=255)
                                                                   ggATATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTGCACGCTACCTg   >  1:1230654/1‑71 (MQ=255)
                                                                    gACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTGCACGCTACCTgg  <  1:2165486/71‑1 (MQ=255)
                                                                      |                                                                    
CGGTAGAAATTTTTGACGAAGTAAACACCGGCAAAAATACCGGTAGCGGCGTACCGTGGGTCACCTGGGACATCATCCCCGACAATGACGATGCGGAAATCGAAGTTTACATGGCAGGCGGCGGCTGCACGCTACCTGG  >  minE/2057516‑2057654

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: