Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 359,429 T→C 9.2% Y220H (TAT→CAT)  sfmD → predicted outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE359,4290TC 9.2% 210.8 / 10.5 76Y220H (TAT→CAT) sfmDpredicted outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base T (36/33);  new base C (5/2);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

CTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTAT  >  minE/359359‑359499
                                                                      |                                                                      
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)
      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)
          aGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTa                                                                       <  1:2406506/62‑1 (MQ=255)
          aGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTcgcg                                                               <  1:3149208/70‑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)
                cTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATgg                                                        <  1:2468665/71‑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:2174738/1‑43 (MQ=255)
                                     ccGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTcgcg                                                               >  1:2234246/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: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:3434011/1‑67 (MQ=38)
                                                       aaTAGTTCCAGGAACCATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCt                     >  1:1388056/1‑67 (MQ=38)
                                                       aaTAGTTCCAGGAACCATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCt                     >  1:3779825/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)
                                                                  gAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAgcgcgc        <  1:625772/69‑1 (MQ=255)
                                                                  gAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGcc       <  1:579548/70‑1 (MQ=255)
                                                                      tATTTTCGCGGAAATGGATATCATTCAGAACAGTGGaataa                                <  1:3810824/41‑1 (MQ=255)
                                                                      tATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCattat  >  1:2079629/1‑71 (MQ=255)
                                                                      tATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCattat  >  1:1190528/1‑71 (MQ=255)
                                                                      |                                                                      
CTATTTTTTTAGTGAGCTCAGCGGGATTAATATTGGCCCGTGGCGTTTACGCAACAATGGTTCCTGGAACTATTTTCGCGGAAATGGATATCATTCAGAACAGTGGAATAATATTGGCACCTGGGTACAGCGCGCCATTAT  >  minE/359359‑359499

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: