Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 138,668 G→A 7.0% S500L (TCA→TTA)  htrE ← predicted outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported138,6680GA 7.0% 336.9 / 10.9 100S500L (TCA→TTA) htrEpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base G (50/43);  new base A (4/3);  total (54/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

CCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTT  >  Exported/138602‑138733
                                                                  |                                                                 
ccGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATTCGTGAGTa                                                               <  1:3209585/71‑1 (MQ=255)
     aTCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGa                                                                  >  1:737002/1‑63 (MQ=255)
       cTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACg                                                         >  1:630053/1‑70 (MQ=255)
        ttttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCa                                                       >  1:3620740/1‑71 (MQ=255)
          ttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACg                                                         <  1:3159958/67‑1 (MQ=255)
          ttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACg                                                         <  1:479278/67‑1 (MQ=255)
          ttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATg                                                     >  1:1435694/1‑71 (MQ=255)
           tttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACg                                                         <  1:3165100/66‑1 (MQ=255)
           tttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATg                                                     >  1:4126721/1‑70 (MQ=255)
           tttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATg                                                     >  1:2420384/1‑70 (MQ=255)
           tttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATg                                                     >  1:1757488/1‑70 (MQ=255)
            ttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATgg                                                    <  1:3049796/70‑1 (MQ=255)
            ttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATgg                                                    <  1:1177864/70‑1 (MQ=255)
                 aaaTTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                              >  1:3105421/1‑71 (MQ=255)
                 aaaTTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                              >  1:3094467/1‑71 (MQ=255)
                 aaaTTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                              >  1:1004124/1‑71 (MQ=255)
                 aaaTTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                              >  1:3719719/1‑71 (MQ=255)
                    tttCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGtt                                           >  1:3597348/1‑71 (MQ=255)
                    tttCAACGGTTGGTTAATACTGACCGTAACCTGATATTTCATGCGTGAGTAATAACGCATgg                                                    >  1:2151492/1‑62 (MQ=255)
                       cAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCa                                                       >  1:4112014/1‑56 (MQ=255)
                       cAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCa                                         <  1:1403787/70‑1 (MQ=255)
                        aaCGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCaa                                        <  1:2131423/70‑1 (MQ=255)
                          cGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                              >  1:2979259/1‑62 (MQ=255)
                          cGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                              >  1:565893/1‑62 (MQ=255)
                           ggttggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                    <  1:1607005/71‑1 (MQ=255)
                            gttggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATgg                                                    <  1:2166734/54‑1 (MQ=255)
                              tggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTg                                 >  1:2954923/1‑71 (MQ=255)
                              tggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTg                                 >  1:1921383/1‑71 (MQ=255)
                               ggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTg                                 <  1:1868192/70‑1 (MQ=255)
                               ggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTg                                 <  1:110519/70‑1 (MQ=255)
                                gttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGtt                               <  1:526344/71‑1 (MQ=255)
                                   aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAAt                                                             >  1:3466460/1‑38 (MQ=255)
                                   aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                    >  1:279796/1‑63 (MQ=255)
                                   aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                    >  1:155244/1‑63 (MQ=255)
                                   aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                    >  1:383213/1‑63 (MQ=255)
                                   aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                            >  1:2334553/1‑71 (MQ=255)
                                   aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                            >  1:1623906/1‑71 (MQ=255)
                                   aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                            >  1:2924489/1‑71 (MQ=255)
                                   aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCg                             <  1:713027/70‑1 (MQ=255)
                                    aTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGAt                                     >  1:3366235/1‑61 (MQ=255)
                                    aTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGAt                                     >  1:1416778/1‑61 (MQ=255)
                                    aTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCggg                           <  1:1579855/71‑1 (MQ=255)
                                    aTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                            <  1:334817/70‑1 (MQ=255)
                                    aTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                            <  1:2489645/70‑1 (MQ=255)
                                    aTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                            <  1:1549132/70‑1 (MQ=255)
                                      aCTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                    <  1:4298124/60‑1 (MQ=255)
                                       cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCaa                                        >  1:228666/1‑55 (MQ=255)
                                       cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCaa                                        >  1:4057896/1‑55 (MQ=255)
                                       cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                    >  1:548630/1‑59 (MQ=255)
                                       cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                    >  1:984845/1‑59 (MQ=255)
                                       cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                    >  1:2110137/1‑59 (MQ=255)
                                       cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGAt                         <  1:2195673/70‑1 (MQ=255)
                                       cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATg                        >  1:2544496/1‑71 (MQ=255)
                                       cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATg                        >  1:478779/1‑71 (MQ=255)
                                         gACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGtt                      >  1:1036854/1‑71 (MQ=255)
                                          aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGttt                     <  1:3919660/71‑1 (MQ=255)
                                          aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGttt                     <  1:758373/71‑1 (MQ=255)
                                          aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGttt                     <  1:3120584/71‑1 (MQ=255)
                                          aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGttt                     <  1:1856634/71‑1 (MQ=255)
                                          aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGttt                     <  1:321674/71‑1 (MQ=255)
                                          aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGttt                     <  1:3365012/71‑1 (MQ=255)
                                             gTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTc                                          <  1:146049/47‑1 (MQ=255)
                                              tAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACt                 >  1:3627984/1‑71 (MQ=255)
                                              tAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACt                 >  1:2555248/1‑71 (MQ=255)
                                                aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGAt                         >  1:4023860/1‑61 (MQ=255)
                                                aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGAt                         >  1:1919919/1‑61 (MQ=255)
                                                aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGAt                         >  1:1147256/1‑61 (MQ=255)
                                                aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTc                    >  1:79171/1‑66 (MQ=255)
                                                aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                >  1:747287/1‑70 (MQ=255)
                                                aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                >  1:146286/1‑70 (MQ=255)
                                                    gATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtcatc            >  1:4102601/1‑70 (MQ=255)
                                                     aTTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtcatc            <  1:4118998/69‑1 (MQ=255)
                                                      tttttCATGCGTGAGTAATTACGCATGGATTTCGGTTCaa                                        <  1:1033615/40‑1 (MQ=255)
                                                      tttttCATGCGTGAGTAATTACGCATGGATTTCGGTTCaa                                        <  1:1570073/40‑1 (MQ=255)
                                                      tttttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtcatc            <  1:244486/68‑1 (MQ=255)
                                                        tttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCtgg                           <  1:3618505/51‑1 (MQ=255)
                                                        tttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCggg                           <  1:4000234/51‑1 (MQ=255)
                                                        tttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAAtt        <  1:1475779/70‑1 (MQ=255)
                                                        atccATGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAAt         >  1:616378/4‑69 (MQ=25)
                                                         ttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTa       <  1:2923398/70‑1 (MQ=255)
                                                         ttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTa       <  1:2648552/70‑1 (MQ=255)
                                                         ttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTa       <  1:1417112/70‑1 (MQ=255)
                                                          tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCtt                                  >  1:1078763/1‑42 (MQ=255)
                                                          tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCtt                                  >  1:863713/1‑42 (MQ=255)
                                                          tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCtt                                  >  1:2362988/1‑42 (MQ=255)
                                                          tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCtt                                  >  1:3215845/1‑42 (MQ=255)
                                                          tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCtt                                  >  1:1660399/1‑42 (MQ=255)
                                                          ccATGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                <  1:4251833/59‑1 (MQ=25)
                                                          ccATGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                <  1:2730351/59‑1 (MQ=25)
                                                          ccATGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                <  1:406446/59‑1 (MQ=25)
                                                            aTGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGt   <  1:1069039/71‑1 (MQ=255)
                                                            aTGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtca              >  1:3247666/1‑60 (MQ=37)
                                                            aTGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtca              >  1:1580655/1‑60 (MQ=37)
                                                            aTGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtca              >  1:2794561/1‑60 (MQ=37)
                                                             tGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGt   >  1:65728/1‑70 (MQ=255)
                                                             tGCCTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGtt  <  1:1331536/71‑1 (MQ=255)
                                                                 tGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTc                    <  1:105678/49‑1 (MQ=255)
                                                                 tGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTc                    <  1:1743654/49‑1 (MQ=255)
                                                                 tGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTc                    <  1:783074/49‑1 (MQ=255)
                                                                 tGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGt   <  1:933718/66‑1 (MQ=255)
                                                                  |                                                                 
CCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTT  >  Exported/138602‑138733

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: