Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,531,994 T→C 50.7% L314L (TTA→TTG kdtA ← 3‑deoxy‑D‑manno‑octulosonic‑acid transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,531,9940TC50.7% 8.9 / 221.5 154L314L (TTA→TTGkdtA3‑deoxy‑D‑manno‑octulosonic‑acid transferase
Reads supporting (aligned to +/- strand):  ref base T (37/39);  new base C (38/40);  total (75/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 157 total aligned reads displayed.

CAGCGGATTATGCCCACCACGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCTGGTGGAGGGGACTTCCCCTGAAGAGC  >  minE/2531925‑2532064
                                                                     |                                                                      
cAGCGGATTATGCCCACCACGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGt                                                                        >  1:711515/1‑70 (MQ=255)
  gCGGATTATGCCCACCACGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAAc                                                                     <  1:2027605/71‑1 (MQ=255)
   cGGATTATGCCCACCACGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACa                                                                    >  1:3776368/1‑71 (MQ=255)
   cGGATTATGCCCACCACGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACa                                                                    >  1:1881103/1‑71 (MQ=255)
   cGGATTATGCCCACCACGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACa                                                                    >  1:1137642/1‑71 (MQ=255)
     gATTATGCCCACCACGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATc                                                                  <  1:182715/71‑1 (MQ=255)
     gATTATGCCCACCACGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATc                                                                  <  1:3158350/71‑1 (MQ=255)
       ttATGCCCACCACGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACa                                                                    >  1:3997874/1‑67 (MQ=255)
        tATGCCCACCACGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCaa                                                                >  1:3891478/1‑70 (MQ=255)
                ccacGTTCAACCAGTGAACGGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCAt                                                       >  1:3641849/1‑71 (MQ=255)
                ccacGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGc                                                           <  1:2600964/67‑1 (MQ=255)
                ccacGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCAt                                                       >  1:2399120/1‑71 (MQ=255)
                ccacGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCAt                                                       >  1:1123782/1‑71 (MQ=255)
                 cacGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATc                                                      >  1:3421135/1‑71 (MQ=255)
                 cacGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATc                                                      >  1:3121734/1‑71 (MQ=255)
                 cacGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATc                                                      >  1:2319502/1‑71 (MQ=255)
                  acGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATc                                                      >  1:1569376/1‑70 (MQ=255)
                   cGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACa                                                                    >  1:506054/1‑55 (MQ=255)
                   cGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTc                                                             <  1:1899778/62‑1 (MQ=255)
                    gTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTa                                                   <  1:3046682/71‑1 (MQ=255)
                       cAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCg                                                >  1:3688538/1‑71 (MQ=255)
                       cAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCAATACAGTAACATCAACTCGCCCATCGTATCg                                                >  1:4102533/1‑71 (MQ=255)
                        aaCCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCGTACAGTAACATCa                                                                 <  1:171921/53‑1 (MQ=255)
                          ccAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGt                                                    <  1:2645015/64‑1 (MQ=255)
                           cAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCaa                                                                >  1:1839474/1‑51 (MQ=255)
                           cAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCaa                                                                >  1:111486/1‑51 (MQ=255)
                           cAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCaa                                                                >  1:4543511/1‑51 (MQ=255)
                           cAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCaa                                                                >  1:1609911/1‑51 (MQ=255)
                            aGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCg                                                     <  1:15806/61‑1 (MQ=255)
                             gTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCaa                                                                      >  1:4449031/1‑43 (MQ=255)
                               gAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACca                                         >  1:1020958/1‑70 (MQ=255)
                               gAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACca                                         >  1:4321675/1‑70 (MQ=255)
                                aaCCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCAt                                                       >  1:4241972/1‑55 (MQ=255)
                                aaCCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACcaca                                       <  1:3529213/71‑1 (MQ=255)
                                aaCCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACcaca                                       <  1:1738471/71‑1 (MQ=255)
                                aaCCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACcaca                                       <  1:2972861/71‑1 (MQ=255)
                                aaCCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACcac                                        >  1:4114520/1‑70 (MQ=255)
                                aaCCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACcac                                        <  1:3901735/70‑1 (MQ=255)
                                 aCCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACcaca                                       <  1:1783174/70‑1 (MQ=255)
                                 aCCGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACcaca                                       <  1:722036/70‑1 (MQ=255)
                                   cGCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGt                                                    <  1:1791278/55‑1 (MQ=255)
                                    gCCAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCt                                   <  1:1024036/71‑1 (MQ=255)
                                     ccAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCg                                                <  1:584882/57‑1 (MQ=255)
                                     ccAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTg                                  >  1:796306/1‑71 (MQ=255)
                                     ccAACAAAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCa                                             >  1:3698335/1‑60 (MQ=255)
                                      caacaaAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGc                                 <  1:4030820/71‑1 (MQ=255)
                                      caacaaAGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGc                                 <  1:3066702/71‑1 (MQ=255)
                                       aacaaAGGCGAGATCGGCAATGCCATACAGCAACAt                                                                   >  1:944335/1‑36 (MQ=255)
                                        acaaAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGt                               >  1:2426182/1‑71 (MQ=255)
                                         caaAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTg                              >  1:3954883/1‑71 (MQ=255)
                                         caaAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTg                              >  1:86960/1‑71 (MQ=255)
                                         caaAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTg                              >  1:2164829/1‑71 (MQ=255)
                                         caaAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTg                              <  1:2791198/71‑1 (MQ=255)
                                          aaaGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTAt                                                  >  1:15204/1‑50 (MQ=255)
                                          aaaGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTAt                                                  >  1:2428921/1‑50 (MQ=255)
                                           aaGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCt                            >  1:1178104/1‑71 (MQ=255)
                                           aaGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCt                            <  1:4021700/71‑1 (MQ=255)
                                            aGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCg                                                     <  1:1917290/45‑1 (MQ=255)
                                            aGGCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCt                                   <  1:2060332/63‑1 (MQ=255)
                                             ggCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCtgg                          <  1:1208418/71‑1 (MQ=255)
                                              gCGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCtggt                         <  1:2768775/71‑1 (MQ=255)
                                               cGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGt                                                    <  1:3867727/43‑1 (MQ=255)
                                               cGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCtggtg                        <  1:3683089/71‑1 (MQ=255)
                                               cGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCtggtg                        <  1:4434590/71‑1 (MQ=255)
                                               cGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCtggtg                        <  1:4295044/71‑1 (MQ=255)
                                               cGAGATCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCtggt                         >  1:3194331/1‑70 (MQ=255)
                                                gagaTCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCtggtgg                       <  1:1815889/71‑1 (MQ=255)
                                                gagaTCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCtggtgg                       <  1:3540938/71‑1 (MQ=255)
                                                   aTCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCtggt                         >  1:4285401/1‑66 (MQ=255)
                                                   aTCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCt                            <  1:1512847/63‑1 (MQ=255)
                                                   aTCGGCAATGCCATACAGTAACATAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCt                            <  1:4153627/62‑1 (MQ=255)
                                                   aTCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCtggtgg                       <  1:4251291/68‑1 (MQ=255)
                                                   aTCGGCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCtggtgg                       <  1:3183166/68‑1 (MQ=255)
                                                    tCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCTGGTGGAggg                   <  1:3592157/71‑1 (MQ=255)
                                                      ggCAATGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCTGGTGGAgggg                  <  1:1926811/70‑1 (MQ=255)
                                                          aTGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCTGGTGGAGGGGACTTc             <  1:4580199/71‑1 (MQ=255)
                                                          aTGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCTGGTGGAGGGGACTTc             >  1:3190999/1‑71 (MQ=255)
                                                           tGCCATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCTGGTGGAgggg                  <  1:3443497/65‑1 (MQ=255)
                                                            gCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACaa                                      >  1:855855/1‑44 (MQ=255)
                                                            gCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAAc                                     >  1:802301/1‑45 (MQ=255)
                                                             ccATACAGCAACATCAACTCGCCCATCGTATCGCCAAcc                                          >  1:2668380/1‑39 (MQ=255)
                                                             ccATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTg                              >  1:1365504/1‑51 (MQ=255)
                                                             ccATACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTg                              >  1:4400523/1‑51 (MQ=255)
                                                              cATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTg                                  <  1:212137/46‑1 (MQ=255)
                                                              cATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTTCGTGCt                            >  1:3302455/1‑52 (MQ=255)
                                                              cATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCt                            >  1:3769776/1‑52 (MQ=255)
                                                              cATACAGTAACATCAACTCCCCCATCGTATCGCCAACCACAACCTg                                  <  1:2212210/46‑1 (MQ=255)
                                                               aTACAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCt                            <  1:2081645/51‑1 (MQ=255)
                                                                tACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGc                             <  1:394085/49‑1 (MQ=255)
                                                                tACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCTGGTGGAGGGGACTTCCCCTGa       <  1:2739631/71‑1 (MQ=255)
                                                                tACAGCAACATCAACTCGCCCATCGTATCGCCAACcaca                                       <  1:1344269/39‑1 (MQ=255)
                                                                 aCAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCtggtg                        <  1:1089470/53‑1 (MQ=255)
                                                                 aCAGCAACATCAACTCGCCCATCGTATCGCCAACcac                                        <  1:1305791/37‑1 (MQ=255)
                                                                  cAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCTGGTGGAGGGGACTTCCCCTg        >  1:4361236/1‑68 (MQ=255)
                                                                  cAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCTGGTGGAGGGGACTTCCCCTGAag     <  1:216347/71‑1 (MQ=255)
                                                                  cAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCTGGTGGAg                     >  1:489492/1‑55 (MQ=255)
                                                                  cAGCAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCTGGTGGAGGGGACTTCCCCTg        >  1:2492078/1‑68 (MQ=255)
                                                                    gTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCTGGTGGAGGGGACTTCCCCTGAagag   <  1:2129634/71‑1 (MQ=255)
                                                                     tAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCtggtg                        <  1:2824271/49‑1 (MQ=255)
                                                                     tAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCTGGTGGAGGGGACTTCCCCTGAAGAGc  >  1:4083113/1‑71 (MQ=255)
                                                                     |                                                                      
CAGCGGATTATGCCCACCACGTTCAACCAGTGAACCGCCAACAAAGGCGAGATCGGCAATGCCATACAGTAACATCAACTCGCCCATCGTATCGCCAACCACAACCTGCGTGCTGGTGGAGGGGACTTCCCCTGAAGAGC  >  minE/2531925‑2532064

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: