Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 638,168 T→A 6.0% D141E (GAT→GAA ycbB → predicted carboxypeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported638,1680TA 6.0% 242.3 / 6.1 83D141E (GAT→GAAycbBpredicted carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base T (43/35);  new base A (2/3);  total (45/38)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

ACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATGGGCTATCTCCATTTCATTG  >  Exported/638099‑638238
                                                                     |                                                                      
aCAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATc                                                                       >  1:3191807/1‑71 (MQ=255)
 cAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATcc                                                                      <  1:1026154/71‑1 (MQ=255)
    cTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTgg                                                                   >  1:75916/1‑71 (MQ=255)
       gCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGgtgt                                                                >  1:26921/1‑71 (MQ=255)
        cAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGtt                                                               <  1:3621181/71‑1 (MQ=255)
         agagGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGtt                                                               >  1:1118174/1‑70 (MQ=255)
         agagGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGtt                                                               >  1:687747/1‑70 (MQ=255)
         agagGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTa                                                              <  1:2110291/71‑1 (MQ=255)
         agagGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTa                                                              <  1:2114550/71‑1 (MQ=255)
           agGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAAc                                                            >  1:731058/1‑71 (MQ=255)
           agGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGAACCTGGTGTTAAc                                                            <  1:2926793/71‑1 (MQ=255)
            ggTGGCGATTGCCGGTTTCCAGCCTCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACg                                                           <  1:3178462/71‑1 (MQ=255)
              tGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACggg                                                         <  1:1043267/71‑1 (MQ=255)
              tGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACggg                                                         <  1:1047503/71‑1 (MQ=255)
                  gATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGGTAACGGGATgg                                                     >  1:2981477/1‑71 (MQ=255)
                   aTTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATgg                                                     >  1:2733007/1‑70 (MQ=255)
                    ttGCCGGTTTCCAGCTGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAAc                                                            <  1:3696303/62‑1 (MQ=255)
                     tGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGa                                                        >  1:3168847/1‑65 (MQ=255)
                      gCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTgg                                                                   >  1:19091/1‑53 (MQ=255)
                       ccGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCAcgc                                                <  1:1921196/71‑1 (MQ=255)
                       ccGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCAcg                                                 <  1:2858711/70‑1 (MQ=255)
                       ccGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCAcg                                                 <  1:1730340/70‑1 (MQ=255)
                       ccGGTTTCCAGCCGCAGTTTAATAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCAcgc                                                <  1:414990/70‑1 (MQ=255)
                         ggTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACg                                                           >  1:2082962/1‑58 (MQ=255)
                         ggTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGa                                              >  1:2069163/1‑71 (MQ=255)
                         ggTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGa                                              <  1:298502/71‑1 (MQ=255)
                          gTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGAc                                             >  1:1653049/1‑71 (MQ=255)
                          gTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGAACCTGGTGTTAACGGGATGGCACGCGAc                                             >  1:2156107/1‑71 (MQ=255)
                             tCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACgt                                           <  1:3273763/70‑1 (MQ=255)
                               cAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACgtggt                                        >  1:787926/1‑71 (MQ=255)
                               cAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACgtggt                                        >  1:1132252/1‑71 (MQ=255)
                               cAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACgtgg                                         <  1:1703634/70‑1 (MQ=255)
                                  ccGCAGTTTAATAAATGGGTAGAGTTACTGACCGAACCTGGTGTTAACGGGATGGCACGCGACGTGGTGct                                     >  1:2146218/1‑71 (MQ=255)
                                   cGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACgg                                                          <  1:996665/49‑1 (MQ=255)
                                   cGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACgg                                                          <  1:3560891/49‑1 (MQ=255)
                                   cGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACgg                                                          <  1:1314110/49‑1 (MQ=255)
                                    gCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAAc                                                            >  1:56064/1‑46 (MQ=255)
                                    gCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAAc                                                            >  1:3341748/1‑46 (MQ=255)
                                    gCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAAc                                                            >  1:2197688/1‑46 (MQ=255)
                                      aGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGctctct                                 >  1:3133943/1‑71 (MQ=255)
                                      aGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGctctct                                 <  1:2429812/71‑1 (MQ=255)
                                      aGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGctctc                                  >  1:2845895/1‑70 (MQ=255)
                                      aGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGAGTTAACGGGATGGCACGCGACGTGGTGctctct                                 <  1:3365232/71‑1 (MQ=255)
                                       gTTTAATAAATGGGTAGAGTTACTGACCGAACCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTg                                <  1:450324/71‑1 (MQ=255)
                                        tttAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATgg                                                     <  1:1814229/49‑1 (MQ=255)
                                        tttAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATgg                                                     <  1:901181/49‑1 (MQ=255)
                                        tttAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACgt                                           >  1:1236204/1‑59 (MQ=255)
                                        tttAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACgt                                           >  1:2863420/1‑59 (MQ=255)
                                        tttAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGctct                                   <  1:122123/67‑1 (MQ=255)
                                          taataaATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGAc                                             >  1:2021440/1‑55 (MQ=255)
                                           aataaATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATg                                                      <  1:1220849/45‑1 (MQ=255)
                                           aataaATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATg                                                      <  1:840060/45‑1 (MQ=255)
                                            ataaATGGGTAGAGTTACTGACCGAACCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGc                            <  1:3176296/70‑1 (MQ=255)
                                                  gggTAGAGTTACTGACCGATCCTGGTGTTAACGTGATGGCACGCGACGTGGTGCTCTCTGATGCGatgatg                     >  1:1540872/1‑71 (MQ=255)
                                                  gggTAGAGTTACTGACCGATCCTGGTGTTAACGTGATGGCACGCGACGTGGTGCTCTCTGATGCGatgatg                     >  1:3162586/1‑71 (MQ=255)
                                                  gggTAGAGTTACTGACCGATCCTGGTGTTAACGTGATGGCACGCGACGTGGTGCTCTCTGATGCGatgatg                     >  1:962557/1‑71 (MQ=255)
                                                   ggTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGc                                                    >  1:417711/1‑39 (MQ=255)
                                                   ggTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGatgatg                     >  1:845383/1‑70 (MQ=255)
                                                   ggTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATgg                    >  1:1349900/1‑71 (MQ=255)
                                                    gTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACgtggtg                                       >  1:652724/1‑51 (MQ=255)
                                                    gTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGa                               >  1:2098340/1‑59 (MQ=255)
                                                    gTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATggg                   <  1:2852282/71‑1 (MQ=255)
                                                    gTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATggg                   <  1:2848039/71‑1 (MQ=255)
                                                     tAGAGTTACTGACCGATCCTGGTGTTAACGGGATg                                                      <  1:3622866/35‑1 (MQ=255)
                                                     tAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGa                               <  1:557324/58‑1 (MQ=255)
                                                     tAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGa                               <  1:1945687/58‑1 (MQ=255)
                                                     tAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATGGGc                  >  1:3647644/1‑71 (MQ=255)
                                                     tAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATGGGc                  >  1:3436790/1‑71 (MQ=255)
                                                     tAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATGGGc                  >  1:1330276/1‑71 (MQ=255)
                                                       gagTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATGGGCt                 <  1:1545633/70‑1 (MQ=255)
                                                         gTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATGGGCTAt               <  1:1615417/70‑1 (MQ=255)
                                                           tACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATGGGCTAtctc            >  1:521488/1‑71 (MQ=255)
                                                           tACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATGGGCTAtctc            >  1:1529245/1‑71 (MQ=255)
                                                            aCTGACCGATCCTGGTGTTAACGGGATGGCACGCGAc                                             >  1:3735539/1‑37 (MQ=255)
                                                            aCTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATGGGCTAtctc            >  1:2857700/1‑70 (MQ=255)
                                                            aCTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATGGGCTAtct             >  1:1468093/1‑69 (MQ=255)
                                                                  cGATCCTGGTGTTAACGGGATGGCAGGCGACGTGGTGCTCTCTGATGCGATGATGGGCTATCTCCATTTCa     >  1:171557/1‑71 (MQ=255)
                                                                  cGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATGGGCTATCTCCATTTCa     >  1:688889/1‑71 (MQ=255)
                                                                  cGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATGGGCTATCTCCATTTCa     >  1:2204862/1‑71 (MQ=255)
                                                                    aTCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCgat                         <  1:1086082/49‑1 (MQ=255)
                                                                    aTCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCgat                         <  1:3275434/49‑1 (MQ=255)
                                                                    aTCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCgat                         <  1:3047324/49‑1 (MQ=255)
                                                                    aTCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATGGGCTATCTCCATTTCAtt   >  1:646071/1‑71 (MQ=255)
                                                                     tCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATGGGCTATCTCCATTTCATTg  >  1:565969/1‑71 (MQ=255)
                                                                     |                                                                      
ACAGCTGGCAGAGGTGGCGATTGCCGGTTTCCAGCCGCAGTTTAATAAATGGGTAGAGTTACTGACCGATCCTGGTGTTAACGGGATGGCACGCGACGTGGTGCTCTCTGATGCGATGATGGGCTATCTCCATTTCATTG  >  Exported/638099‑638238

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: