Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 359,423 T→A 7.2% W218R (TGG→AGG)  sfmD → predicted outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE359,4230TA 7.2% 234.4 / 8.7 83W218R (TGG→AGG) sfmDpredicted outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base T (39/38);  new base A (4/2);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

GATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGC  >  minE/359354‑359491
                                                                     |                                                                    
gATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTg                                                                     >  1:3903480/1‑71 (MQ=255)
  tAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGa                                                                   <  1:1345127/71‑1 (MQ=255)
     cTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACt                                                                <  1:2845743/71‑1 (MQ=255)
     cTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACt                                                                <  1:3310381/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:2677009/1‑60 (MQ=255)
          tttttAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCt                                                                      >  1:2594310/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:3075681/50‑1 (MQ=255)
                      tCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTgg                                                                    <  1:3045602/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:3087913/1‑71 (MQ=255)
                                     ttGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAg                                >  1:1798006/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:3366346/37‑1 (MQ=255)
                                          ccGTGGCGTTTACGCAACAATGGTTCCTGGAACTAtt                                                             <  1:1482775/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:675094/1‑70 (MQ=255)
                                           cGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaa                           >  1:2917323/1‑70 (MQ=255)
                                            gTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCAtt                                        >  1:951174/1‑56 (MQ=255)
                                               gCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGatat                                            <  1:2949179/49‑1 (MQ=255)
                                               gCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAataata                      <  1:1103300/71‑1 (MQ=255)
                                               gCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAataata                      <  1:587262/71‑1 (MQ=255)
                                               gCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAataata                      <  1:2780740/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:1499663/1‑70 (MQ=255)
                                                    tACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTgg                  >  1:1988102/1‑70 (MQ=255)
                                                    tACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTgg                  >  1:2350868/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:1388056/1‑67 (MQ=38)
                                                            aaTAGTTCCAGGAACCATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCt             >  1:3434011/1‑67 (MQ=38)
                                                            aaTAGTTCCAGGAACCATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCt             >  1:3779825/1‑67 (MQ=38)
                                                            aaTAGTTCCAGGAACCATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCt             >  1:1832313/1‑67 (MQ=38)
                                                             aTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaataa                        <  1:2849541/55‑1 (MQ=255)
                                                             aTGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTa        <  1:483774/71‑1 (MQ=255)
                                                                  tCCTGTAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATAtt                    >  1:3177618/1‑54 (MQ=255)
                                                                  tCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATAtt                    >  1:1883310/1‑54 (MQ=255)
                                                                   ccTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAgcgc  >  1:860939/1‑71 (MQ=255)
                                                                    cTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATAtt                    >  1:2233101/1‑52 (MQ=255)
                                                                     tGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCt             <  1:1441385/58‑1 (MQ=255)
                                                                     |                                                                    
GATAGCTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGC  >  minE/359354‑359491

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: