Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,694,258 C→T 41.7% E124K (GAA→AAA)  yrfF ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,694,2580CT41.7% 49.5 / 100.6 84E124K (GAA→AAA) yrfFpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (23/26);  new base T (23/12);  total (46/38)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

AAATAACCAGCGGCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTCGGATCGTCGGTAGA  >  minE/2694192‑2694327
                                                                  |                                                                     
aaaTAACCAGCGGCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTc                                                                       >  1:1892661/1‑67 (MQ=255)
 aaTAACCAGCGGCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGa                                                                  <  1:2244884/71‑1 (MQ=255)
 aaTAACCAGCGGCAGCGAATCAGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGa                                                                  <  1:812739/71‑1 (MQ=255)
  aTAACCAGCGGCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGa                                                                  >  1:2809193/1‑70 (MQ=255)
  aTAACCAGCGGCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGa                                                                  >  1:562708/1‑70 (MQ=255)
  aTAACCAGCGGCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGa                                                                  >  1:97513/1‑70 (MQ=255)
    aaCCAGCGGCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAg                                                               >  1:753006/1‑71 (MQ=255)
     aCCAGCGGCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGaaa                                                                >  1:528971/1‑69 (MQ=255)
     aCCAGCGGCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGaaa                                                                >  1:1645815/1‑69 (MQ=255)
     aCCAGCGGCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGaa                                                                 >  1:2450331/1‑68 (MQ=255)
     aCCAGCGGCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAgg                                                              >  1:1982441/1‑71 (MQ=255)
         gCGGCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAggggg                                                           >  1:1661255/1‑70 (MQ=255)
           ggCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGAAAGGGGGCAg                                                        <  1:3675182/71‑1 (MQ=255)
            gCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGAAAGGGGGCAgg                                                       <  1:244265/71‑1 (MQ=255)
            gCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGAAAGGGGGCAgg                                                       >  1:2736972/1‑71 (MQ=255)
            gCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGa                                                                  <  1:2038554/60‑1 (MQ=255)
            gCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAgg                                                       >  1:666811/1‑71 (MQ=255)
              aGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTg                                                     <  1:1340754/71‑1 (MQ=255)
              aGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTg                                                     <  1:2154592/71‑1 (MQ=255)
               gCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTg                                                     >  1:990529/1‑70 (MQ=255)
                cGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGa                                                   >  1:2059313/1‑71 (MQ=255)
                 gAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGAAAGGGGGCAGGTGGa                                                   >  1:1591368/1‑70 (MQ=255)
                 gAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGAc                                                  >  1:3416137/1‑71 (MQ=255)
                  aaTCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACt                                                 >  1:3597109/1‑71 (MQ=255)
                   aTCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAg                                                               >  1:1732298/1‑56 (MQ=255)
                   aTCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAg                                                               >  1:686939/1‑56 (MQ=255)
                    tCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAggggg                                                           <  1:638844/59‑1 (MQ=255)
                      ggTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTc                                                                       >  1:742373/1‑45 (MQ=255)
                      ggTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTa                                             >  1:3414894/1‑71 (MQ=255)
                        tATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAATGGGGCAGGTGGACTTCTAcc                                           <  1:3028824/71‑1 (MQ=255)
                          tGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGaa                                                                 <  1:3478996/47‑1 (MQ=255)
                          tGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCg                                          <  1:3305551/70‑1 (MQ=255)
                            aaTCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAAt                                       >  1:840066/1‑71 (MQ=255)
                               cAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCg                                     <  1:385438/70‑1 (MQ=255)
                                aGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGaaa                                                                >  1:1049634/1‑42 (MQ=255)
                                aGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGaaa                                                                >  1:1287041/1‑42 (MQ=255)
                                aGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGaaa                                                                >  1:734605/1‑42 (MQ=255)
                                aGTTCAACGGTATTCTCATCGTTGATGTACTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAg                                   >  1:2001489/1‑71 (MQ=255)
                                aGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCa                                                                    >  1:1021109/1‑38 (MQ=255)
                                 gTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCg                                          <  1:3591184/63‑1 (MQ=255)
                                 gTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAgg                                  <  1:3974459/71‑1 (MQ=255)
                                 gTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAgg                                  <  1:3834560/71‑1 (MQ=255)
                                  ttCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAgg                                                              <  1:3138413/42‑1 (MQ=255)
                                    gaaCGGTATTCTCATCGTTGATGTACTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGgtag                               >  1:3042192/2‑71 (MQ=255)
                                    cAACGGTATTCTCATCGTTGATGTACTGTTTCCAGAAAggg                                                             >  1:1322733/1‑41 (MQ=255)
                                    cAACGGTATTCTCATCGTTGATGTACTGTTTCCAGAAAGGGGGCAGGTGGAc                                                  >  1:3159500/1‑52 (MQ=255)
                                    cAACGGTATTCTCATCGTTGATGTACTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGgta                                <  1:505403/70‑1 (MQ=255)
                                    cAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAgg                                  >  1:3258748/1‑68 (MQ=255)
                                    cAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGgtag                               >  1:1337984/1‑71 (MQ=255)
                                     aaCGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCg                                          <  1:3595948/59‑1 (MQ=255)
                                       cGGTATTCTCATCGTTGATGTACTGTTTCCAGaaa                                                                >  1:2470081/1‑35 (MQ=255)
                                       cGGTATTCTCATCGTTGATGTACTGTTTCCAGaaa                                                                >  1:3785055/1‑35 (MQ=255)
                                            ttCTCATCGTTGATGTACTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGgtagta                             <  1:326741/65‑1 (MQ=255)
                                            ttCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGa                                         <  1:1539968/53‑1 (MQ=255)
                                              ctcATCGTTGATGTACTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCAtt                     >  1:3413063/1‑71 (MQ=255)
                                              ctcATCGTTGATGTACTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCAtt                     >  1:2562158/1‑71 (MQ=255)
                                              ctcATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGgt                                 >  1:3162253/1‑59 (MQ=255)
                                              ctcATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGgt                                 >  1:1000358/1‑59 (MQ=255)
                                              ctcATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCAtt                     >  1:143239/1‑71 (MQ=255)
                                               tcATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTaa                            >  1:3089855/1‑63 (MQ=255)
                                                 aTCGTTGATGTACTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGa                                         <  1:3947122/48‑1 (MQ=255)
                                                 aTCGTTGATGTACTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTa                   <  1:2052324/70‑1 (MQ=255)
                                                 aTCGTTGATGTACTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTAt                  >  1:1039223/1‑71 (MQ=255)
                                                  tCGTTGATGTACTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTAtt                 >  1:1201338/1‑71 (MQ=255)
                                                  tCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTAtt                 <  1:354524/71‑1 (MQ=255)
                                                   cGTTGATGTACTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTc                >  1:1936108/1‑71 (MQ=255)
                                                    gTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTc                <  1:3688976/70‑1 (MQ=255)
                                                    gTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGAAGTAACGCCATTTATTc                <  1:3828716/70‑1 (MQ=255)
                                                     ttGATGTACTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTCg               <  1:622119/70‑1 (MQ=255)
                                                       gATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTCGGa             <  1:3997730/70‑1 (MQ=255)
                                                         tGTACTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGa                                         <  1:3110362/40‑1 (MQ=255)
                                                         tGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGa                                    >  1:36979/1‑45 (MQ=255)
                                                         tGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTaa                            <  1:2768851/53‑1 (MQ=255)
                                                          gTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATc                                      <  1:3381061/42‑1 (MQ=255)
                                                           tACTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTCGGAtcgtc        >  1:2575376/1‑71 (MQ=255)
                                                           tACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGgta                                <  1:2838418/47‑1 (MQ=255)
                                                            aCTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTCGGAtc           >  1:1216741/1‑67 (MQ=255)
                                                            aCTGTTCCCAGAAAGGGGGCAGGTGGACTTATACCGAATCGAGGTAGTAACGCCATTTATTCGGAtcgtcg       >  1:2565569/1‑71 (MQ=255)
                                                             cTGTTTCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTCGGATCGTCgg      <  1:2793854/71‑1 (MQ=255)
                                                             cTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTCGGATCGTCgg      <  1:2289528/71‑1 (MQ=255)
                                                             cTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTCGGATCGTCgg      >  1:1537991/1‑71 (MQ=255)
                                                             cTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTCGGATCGTCgg      <  1:392486/71‑1 (MQ=255)
                                                                ttCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGgt                                 <  1:556490/41‑1 (MQ=255)
                                                                ttCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTCGGATCGTCGGTag   <  1:647874/71‑1 (MQ=255)
                                                                ttCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTCGGATCGTCGGTag   <  1:3876857/71‑1 (MQ=255)
                                                                ttCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTCGGATCGTCGGTa    <  1:926584/70‑1 (MQ=255)
                                                                  cccATAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTCGGATCGTCGGTaga  <  1:2848583/70‑1 (MQ=255)
                                                                  cccAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTCGGATCGTCGGTaga  <  1:3938370/70‑1 (MQ=255)
                                                                  |                                                                     
AAATAACCAGCGGCAGCGAATCGGTATGAATCAGTTCAACGGTATTCTCATCGTTGATGTACTGTTCCCAGAAAGGGGGCAGGTGGACTTCTACCGAATCGAGGTAGTAACGCCATTTATTCGGATCGTCGGTAGA  >  minE/2694192‑2694327

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: