Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 359,417 G→A 7.1% G216S (GGT→AGT)  sfmD → predicted outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE359,4170GA 7.1% 260.0 / 6.9 84G216S (GGT→AGT) sfmDpredicted outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base G (39/39);  new base A (4/2);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

TAACGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTA  >  minE/359350‑359485
                                                                   |                                                                    
tAACGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGt                                                                    >  1:1237440/1‑70 (MQ=255)
  aCGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTcc                                                                 >  1:2038952/1‑71 (MQ=255)
  aCGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTcc                                                                 <  1:1862958/71‑1 (MQ=255)
    gATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTg                                                               >  1:3903480/1‑71 (MQ=255)
    gATAGCTATTTTTTTAGTGAGCTCAGGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATgg                                                                     >  1:1976055/1‑64 (MQ=255)
      tAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGa                                                             <  1:1345127/71‑1 (MQ=255)
         cTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACt                                                          <  1:3310381/71‑1 (MQ=255)
         cTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACt                                                          <  1:2845743/71‑1 (MQ=255)
            tttttttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTAtt                                                       >  1:943216/1‑71 (MQ=255)
            tttttttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTAtt                                                       >  1:2800852/1‑71 (MQ=255)
             ttttttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTAttt                                                      >  1:3097954/1‑71 (MQ=255)
              tttttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCt                                                                >  1:2594310/1‑60 (MQ=255)
              tttttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCt                                                                >  1:2677009/1‑60 (MQ=255)
               ttttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTc                                                    >  1:3697154/1‑71 (MQ=255)
                tttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTc                                                    <  1:1952200/70‑1 (MQ=255)
                tttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTc                                                    <  1:2370498/70‑1 (MQ=255)
                tttAGTGAGCTAAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTc                                                    <  1:1201137/70‑1 (MQ=255)
                   aGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTg                                                               >  1:583824/1‑56 (MQ=255)
                   aGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTa                                                         <  1:2406506/62‑1 (MQ=255)
                   aGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTcgcg                                                 <  1:3149208/70‑1 (MQ=255)
                    gTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAAc                                                           <  1:1753766/59‑1 (MQ=255)
                    gTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGa                                               <  1:699160/71‑1 (MQ=255)
                     tGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGaa                                              <  1:306744/71‑1 (MQ=255)
                       aGCTCAGCGTGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAAt                                            <  1:2675801/71‑1 (MQ=255)
                       aGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTc                                                    <  1:3205524/63‑1 (MQ=255)
                       aGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGaaa                                             <  1:3555589/70‑1 (MQ=255)
                       aGCTAAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTg                                                               <  1:2537574/52‑1 (MQ=255)
                         cTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATgg                                          <  1:2468665/71‑1 (MQ=255)
                          tCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTgg                                                              <  1:3045602/50‑1 (MQ=255)
                          tCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTgg                                                              <  1:3075681/50‑1 (MQ=255)
                          tCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGa                                         >  1:2853392/1‑71 (MQ=255)
                           cAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGa                                         >  1:2623568/1‑70 (MQ=255)
                            aGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTa                                                         <  1:1799914/53‑1 (MQ=255)
                            aGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGata                                       >  1:303934/1‑71 (MQ=255)
                            aGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGat                                        <  1:2332296/70‑1 (MQ=255)
                              cGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGatat                                      >  1:518109/1‑70 (MQ=255)
                                     aaTATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACt                                                          <  1:1604878/43‑1 (MQ=255)
                                     aaTATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAg                               >  1:960523/1‑70 (MQ=255)
                                      atatTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGaa                             >  1:3390950/1‑71 (MQ=255)
                                         ttGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAg                          >  1:1798006/1‑71 (MQ=255)
                                         ttGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAg                          >  1:3087913/1‑71 (MQ=255)
                                           ggCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTcgc                                                  >  1:2818704/1‑45 (MQ=255)
                                            gCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTcg                                                   <  1:2869156/43‑1 (MQ=255)
                                            gCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTg                        >  1:2443327/1‑70 (MQ=255)
                                             cccGTGGCGTTTACGCAACAATGGTTCCTGGAACCATTTTCGCGGAAATGGATATCATTCAGAACAGTgg                       >  1:3140493/1‑70 (MQ=255)
                                             cccGTGGCGTTTACGCAAAAATAGTTCCAGGAACCAttg                                                      <  1:2424985/39‑2 (MQ=11)
                                             cccGTGGCGTTTACGCAAAAATAGTTCCAGGAACCAttg                                                      <  1:697921/39‑2 (MQ=11)
                                              ccGTGGCGTTTACGCAACAATGGTTCCTGGAACTAtt                                                       <  1:1482775/37‑1 (MQ=255)
                                              ccGTGGCGTTTACGCAACAATGGTTCCTGGAACTAtt                                                       <  1:3366346/37‑1 (MQ=255)
                                              ccGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTcgcg                                                 >  1:2234246/1‑43 (MQ=255)
                                              ccGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTcgcg                                                 >  1:2174738/1‑43 (MQ=255)
                                               cGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATgg                                          >  1:3887220/1‑49 (MQ=255)
                                               cGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaa                     >  1:2917323/1‑70 (MQ=255)
                                               cGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaa                     >  1:675094/1‑70 (MQ=255)
                                                gTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCAtt                                  >  1:951174/1‑56 (MQ=255)
                                                   gCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGatat                                      <  1:2949179/49‑1 (MQ=255)
                                                   gCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAataata                <  1:587262/71‑1 (MQ=255)
                                                   gCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAataata                <  1:2780740/71‑1 (MQ=255)
                                                   gCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAataata                <  1:1103300/71‑1 (MQ=255)
                                                   gCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAataata                >  1:1568784/1‑71 (MQ=255)
                                                     gTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATTGATATCATTCAGAACAGTGGAATAATAtt              <  1:2253011/71‑1 (MQ=255)
                                                       ttACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTgg            >  1:268936/1‑71 (MQ=255)
                                                       ttACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTgg            >  1:80757/1‑71 (MQ=255)
                                                       ttACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTg             >  1:3409074/1‑70 (MQ=255)
                                                        tACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATAtt              >  1:228967/1‑68 (MQ=255)
                                                        tACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTgg            >  1:2350868/1‑70 (MQ=255)
                                                        tACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTgg            >  1:1499663/1‑70 (MQ=255)
                                                        tACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTgg            >  1:1988102/1‑70 (MQ=255)
                                                        tACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTg             >  1:1049441/1‑69 (MQ=255)
                                                          cGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCa          >  1:727291/1‑70 (MQ=255)
                                                          cGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCAc         <  1:3557522/71‑1 (MQ=255)
                                                           gCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAatat               <  1:3729971/64‑1 (MQ=255)
                                                            caacaaTGGTTCCTGGAACTATTTTCGCGGAAATTGATATCATTCAGAACAGTGGAATAatat               <  1:282161/63‑1 (MQ=255)
                                                            caacaaTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCt       <  1:2059998/71‑1 (MQ=255)
                                                             aacaaTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaataat                 <  1:2185194/60‑1 (MQ=255)
                                                              acaaTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaat                    <  1:1506794/56‑1 (MQ=255)
                                                               caaTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaataat                 >  1:110996/1‑58 (MQ=255)
                                                               caaTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTggg    <  1:2045159/71‑1 (MQ=255)
                                                               caaTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTgg     >  1:2165305/1‑70 (MQ=255)
                                                                aaTAGTTCCAGGAACCATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCt       >  1:1832313/1‑67 (MQ=38)
                                                                aaTAGTTCCAGGAACCATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCt       >  1:3779825/1‑67 (MQ=38)
                                                                aaTAGTTCCAGGAACCATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCt       >  1:3434011/1‑67 (MQ=38)
                                                                aaTAGTTCCAGGAACCATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCt       >  1:1388056/1‑67 (MQ=38)
                                                                 aTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaataa                  <  1:2849541/55‑1 (MQ=255)
                                                                 aTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTa  <  1:483774/71‑1 (MQ=255)
                                                                   |                                                                    
TAACGATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTA  >  minE/359350‑359485

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: