Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,417,054 A→G 14.1% Y416C (TAT→TGT)  yejM → predicted hydrolase, inner membrane

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,417,0540AG14.1% 165.5 / 18.0 64Y416C (TAT→TGT) yejMpredicted hydrolase, inner membrane
Reads supporting (aligned to +/- strand):  ref base A (37/18);  new base G (4/5);  total (41/23)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GTTCTCGTGGGTTTCTTTCAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGTGATTCTGGCAAA  >  Exported/1416986‑1417124
                                                                    |                                                                      
gTTCTCGTGGGTTTCTTTCAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAata                                                                        >  1:1195119/1‑69 (MQ=255)
   ctcGTGGGTTTCTTTCAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGcc                                                                   >  1:2548237/1‑71 (MQ=255)
     cGTGGGTTTCTTTCAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCgg                                                                 >  1:678184/1‑71 (MQ=255)
      gTGGGTTTCTTTCAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCggg                                                                >  1:3658127/1‑71 (MQ=255)
       tGGGTTTCTTTCAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAatat                                                                       <  1:2277805/63‑1 (MQ=255)
       tGGGTTTCTTTCAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAatat                                                                       <  1:2441275/63‑1 (MQ=255)
       tGGGTTTCTTTCAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAatat                                                                       <  1:1620695/63‑1 (MQ=255)
        gggTTTCTTTCAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGcc                                                                   <  1:2394460/66‑1 (MQ=255)
        gggTTTCTTTCAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGcc                                                                   <  1:1569953/66‑1 (MQ=255)
         ggTTTCTTTCAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGggcg                                                              <  1:2780878/70‑1 (MQ=255)
          gTTTCTTTCAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGggcggc                                                            <  1:3126101/71‑1 (MQ=255)
          gTTTCTTTCAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGggcg                                                              >  1:398475/1‑69 (MQ=255)
           tttctttcAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATGTAGCCGGGCggca                                                           <  1:2223785/71‑1 (MQ=255)
             tctttcAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATGTAGCCGGGCggcagg                                                         >  1:256091/1‑71 (MQ=255)
             tctttcAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCggcagg                                                         >  1:3475184/1‑71 (MQ=255)
                  cAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAAt                                                     <  1:3541877/70‑1 (MQ=255)
                  cAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAAt                                                     >  1:1640135/1‑70 (MQ=255)
                   aaTGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGt                                                   >  1:3446790/1‑71 (MQ=255)
                   aaTGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGt                                                   >  1:1003008/1‑71 (MQ=255)
                    aTGTTAGTACCATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTc                                                  <  1:3419839/71‑1 (MQ=255)
                    aTGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTc                                                  <  1:3552544/71‑1 (MQ=255)
                     tGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATGTAGCCGGGCGGCAGGCAATGTCg                                                 <  1:1627305/71‑1 (MQ=255)
                     tGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTc                                                  >  1:1129485/1‑70 (MQ=255)
                       gTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGAt                                               <  1:3148883/71‑1 (MQ=255)
                        tACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATg                                              <  1:1981291/71‑1 (MQ=255)
                             aCATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGAc                                            >  1:2856818/1‑68 (MQ=255)
                               aTTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGAc                                            >  1:341784/1‑66 (MQ=255)
                               aTTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGAc                                            >  1:318113/1‑66 (MQ=255)
                                 tGACGACAGCAATCAGCCGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATc                                      >  1:2447754/1‑70 (MQ=255)
                                 tGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATc                                      >  1:1604536/1‑70 (MQ=255)
                                  gacgacAGCAATCAGCAGGCATTTGCACGGAAATGTAGc                                                                    >  1:456257/1‑39 (MQ=255)
                                  gacgacAGCAATCAGCAGGCATTTGCACGGAAATGTAGc                                                                    >  1:3120405/1‑39 (MQ=255)
                                  gacgacAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCggcaggc                                                        >  1:2387287/1‑51 (MQ=255)
                                  gacgacAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCggcaggc                                                        >  1:2580858/1‑51 (MQ=255)
                                  gacgacAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCggcaggc                                                        >  1:3227209/1‑51 (MQ=255)
                                  gacgacAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCggcaggc                                                        >  1:1931701/1‑51 (MQ=255)
                                  gacgacAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCggcaggc                                                        >  1:2949321/1‑51 (MQ=255)
                                    cgacAGCAATCAGCAGGCATTTGCACGGAAATGTAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAAcc                                  <  1:3123868/71‑1 (MQ=255)
                                     gacAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATg                                              <  1:2895456/58‑1 (MQ=255)
                                     gacAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATg                                              <  1:2911668/58‑1 (MQ=255)
                                     gacAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACcg                                 <  1:776546/71‑1 (MQ=255)
                                       cAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGAc                                            >  1:640518/1‑58 (MQ=255)
                                       cAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACcgcg                               <  1:3639550/71‑1 (MQ=255)
                                       cAGCAATCAGCAGGCAGTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACcgcg                               >  1:302120/1‑71 (MQ=255)
                                          cAATCAGCAGGCATTTGCACGGAAATGTAGCCGGGCGGCAGGCAATGTCGATGACCAGATc                                      >  1:222869/1‑61 (MQ=255)
                                          cAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGc                            >  1:753547/1‑71 (MQ=255)
                                          cAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGc                            >  1:2561748/1‑71 (MQ=255)
                                             tCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGa                                                >  1:1461176/1‑48 (MQ=255)
                                             tCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGt                              <  1:2323854/66‑1 (MQ=255)
                                             tCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCt                           >  1:401017/1‑69 (MQ=255)
                                                gcagGCATTTGCACGGAAATGTAGCCGGGCGGCAGGCAATGTCGATGACCAg                                         <  1:2219347/52‑1 (MQ=255)
                                                     cATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTg                 <  1:1968236/71‑1 (MQ=255)
                                                      aTTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGc                >  1:1461387/1‑71 (MQ=255)
                                                        ttGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGt              >  1:3105810/1‑71 (MQ=255)
                                                        ttGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGt              >  1:3687935/1‑71 (MQ=255)
                                                         tGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGt              <  1:1558312/70‑1 (MQ=255)
                                                           cACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCa                                          >  1:3523959/1‑40 (MQ=255)
                                                           cACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCa                                          >  1:394529/1‑40 (MQ=255)
                                                            aCGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGTGAtt          >  1:1543679/1‑71 (MQ=255)
                                                            aCGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGTGAtt          >  1:30939/1‑71 (MQ=255)
                                                            aCGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGTGAtt          >  1:194177/1‑71 (MQ=255)
                                                            aCGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGTGAtt          >  1:3567832/1‑71 (MQ=255)
                                                            aCGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGTGAtt          >  1:1605780/1‑71 (MQ=255)
                                                            aCGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGTGAtt          >  1:3336303/1‑71 (MQ=255)
                                                             cGGAAATGTAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGTGATTc         <  1:3526236/71‑1 (MQ=255)
                                                             cGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGTGATTc         <  1:403774/71‑1 (MQ=255)
                                                               gAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGc                            >  1:1766080/1‑50 (MQ=255)
                                                               gAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGTGATTCTg       <  1:2801892/71‑1 (MQ=255)
                                                                  atatAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTc                          <  1:2283453/49‑1 (MQ=255)
                                                                  atatAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTc                          <  1:1859439/49‑1 (MQ=255)
                                                                   tataGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGTGATTCTGGCaa   <  1:2790223/71‑1 (MQ=255)
                                                                   tataGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGTGATTCTGGCa    >  1:500098/1‑70 (MQ=255)
                                                                    ataGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGTGATTCTGGCaaa  >  1:2911039/1‑71 (MQ=255)
                                                                    |                                                                      
GTTCTCGTGGGTTTCTTTCAATGGTACTAACATTGACGACAGCAATCAGCAGGCATTTGCACGGAAATATAGCCGGGCGGCAGGCAATGTCGATGACCAGATCAACCGCGTGCTCAATGCACTGCGTGATTCTGGCAAA  >  Exported/1416986‑1417124

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