Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,087,666 C→T 50.0% T54I (ACC→ATC)  yhaH → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,087,6660CT50.0% ‑4.1 / 81.9 58T54I (ACC→ATC) yhaHpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (13/16);  new base T (15/14);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTG  >  minE/2087599‑2087730
                                                                   |                                                                
atcTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCt                                                               <  1:1160262/71‑1 (MQ=255)
atcTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCt                                                               <  1:1718860/71‑1 (MQ=255)
  cTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTga                                                             >  1:2107860/1‑71 (MQ=255)
      aCGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTgacga                                                          >  1:1870671/1‑70 (MQ=255)
       cGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACa                                                        <  1:3506990/71‑1 (MQ=255)
        gTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGTGAAGGCACCCTGACGACaa                                                       <  1:2933982/71‑1 (MQ=255)
          tCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATc                                                     >  1:1884559/1‑71 (MQ=255)
           cGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATc                                                     <  1:2078660/70‑1 (MQ=255)
             tGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTAc                                                  <  1:3219448/71‑1 (MQ=255)
              gCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACg                                                 >  1:866703/1‑71 (MQ=255)
                  gggCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCAtcc                                                                >  1:1140223/1‑52 (MQ=255)
                  gggCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTAt                                             >  1:1847813/1‑71 (MQ=255)
                      tgctgGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATc                                                     <  1:2522821/59‑1 (MQ=255)
                      tgctgGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTg                                         >  1:1375785/1‑71 (MQ=255)
                      tgctgGATAAAATGTTAGGCTGGCAACGCGCTGGGGGCGAAGGCACCCTGACGACAATCTACGGTATTc                                           <  1:1294701/69‑1 (MQ=255)
                       gctgGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTa                                                   <  1:47548/60‑1 (MQ=255)
                         tgGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGgtg                                      >  1:1471453/1‑71 (MQ=255)
                         tgGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTga                                                             >  1:2448509/1‑48 (MQ=255)
                          ggATAAGATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGgtgt                                     >  1:3818/1‑71 (MQ=255)
                           gATAAAATGTTAGGCTTGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGgtg                                      <  1:1098831/69‑1 (MQ=255)
                           gATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGtt                                    >  1:508157/1‑71 (MQ=255)
                            aTAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGttt                                   >  1:3499719/1‑71 (MQ=255)
                               aaaTGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTga                                                             <  1:991775/42‑1 (MQ=255)
                               aaaTGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTa                                >  1:3533259/1‑71 (MQ=255)
                                aaTGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGttttt                                 <  1:1114170/69‑1 (MQ=255)
                                aaTGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGttttt                                 <  1:2589095/69‑1 (MQ=255)
                                 aTGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTAcc                              <  1:3018038/71‑1 (MQ=255)
                                    ttAGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCgtg                           >  1:2952382/1‑71 (MQ=255)
                                     tAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACa                                                        <  1:3783319/41‑1 (MQ=255)
                                      aGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCgtgg                          <  1:1441087/70‑1 (MQ=255)
                                      aGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTa                                                   <  1:2625455/45‑1 (MQ=255)
                                      aGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTa                                                   <  1:3777413/45‑1 (MQ=255)
                                            gCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTg                                         <  1:3342894/49‑1 (MQ=255)
                                            gCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCgtg                           >  1:172173/1‑63 (MQ=255)
                                             cAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGt                  >  1:2527548/1‑71 (MQ=255)
                                             cAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGt                  <  1:2697472/71‑1 (MQ=255)
                                              aaCGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGtt                 >  1:1441953/1‑71 (MQ=255)
                                               aCGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTAc                               >  1:433353/1‑56 (MQ=255)
                                               aCGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCgg                   >  1:3827355/1‑68 (MQ=255)
                                               aCGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCgg                   >  1:1370214/1‑68 (MQ=255)
                                               aCGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCgg                   >  1:702810/1‑68 (MQ=255)
                                                  cgcTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTa                                <  1:177989/52‑1 (MQ=255)
                                                    cTGGTGGCGAAGGCACCCTGACGACAATCTACTGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGtt            <  1:3399400/70‑1 (MQ=255)
                                                     tggtggCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGtttt                                  >  1:3429698/1‑47 (MQ=255)
                                                     tggtggCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGtttt                                  >  1:2639621/1‑47 (MQ=255)
                                                       gtggCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTccgc        >  1:1249732/1‑71 (MQ=255)
                                                       gtggCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTccg         <  1:3309672/70‑1 (MQ=255)
                                                       gtggCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTccg         <  1:8150/70‑1 (MQ=255)
                                                         ggCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCgtg                           >  1:473384/1‑50 (MQ=255)
                                                          gCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTAcc                              <  1:2887599/46‑1 (MQ=255)
                                                          gCGAAGGCACCCTGACGACAATCTACGGGATTCTGGTGTTTTTAcc                              <  1:2892368/46‑1 (MQ=255)
                                                           cGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCgccgcc    <  1:718166/71‑1 (MQ=255)
                                                           cGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCgccgcc    <  1:939856/71‑1 (MQ=255)
                                                           cGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCgccgcc    >  1:2947938/1‑71 (MQ=255)
                                                           agAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCgg                   >  1:1352860/2‑56 (MQ=255)
                                                            gAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTc                <  1:1877075/58‑1 (MQ=255)
                                                            gAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCa               <  1:1452325/59‑1 (MQ=255)
                                                              aGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTg  >  1:3430486/1‑70 (MQ=255)
                                                              aGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGtt            >  1:2787890/1‑60 (MQ=255)
                                                                   |                                                                
ATCTTTACGTTCGTGCTGGGGCTGCTGGATAAAATGTTAGGCTGGCAACGCGCTGGTGGCGAAGGCACCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTG  >  minE/2087599‑2087730

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: