Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,724,534 C→T G121S (GGT→AGT)  secB ← export protein SecB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,724,5340CT93.4% 271.4 / 7.5 91G121S (GGT→AGT) secBexport protein SecB
Reads supporting (aligned to +/- strand):  ref base C (6/0);  new base T (39/46);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 1.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

ACAGCGCATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACCGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGAACAGAA  >  REL606/3724486‑3724581
                                                |                                               
aCAGCGCATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACtgc                                               >  1:620879/1‑51 (MQ=255)
aCAGCGCATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACtgc                                               <  1:6004004/51‑1 (MQ=255)
 cAGCGCATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACtgcg                                              >  1:3553232/1‑51 (MQ=255)
    cgcATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGa                                           >  1:3479397/1‑51 (MQ=255)
    cgcATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGa                                           <  1:137330/51‑1 (MQ=255)
      cATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATa                                         <  1:5932053/51‑1 (MQ=255)
      cATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATa                                         >  1:244786/1‑51 (MQ=255)
      cATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATa                                         >  1:1989470/1‑51 (MQ=255)
      cATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATa                                         <  1:5144793/51‑1 (MQ=255)
       aTCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATAc                                        <  1:5198420/51‑1 (MQ=255)
       aTCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATAc                                        >  1:3466983/1‑51 (MQ=255)
       aTCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATAc                                        >  1:1500980/1‑51 (MQ=255)
       aTCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATAc                                        <  1:7249220/51‑1 (MQ=255)
        tCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATAcc                                       >  1:2163004/1‑51 (MQ=255)
        tCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATAcc                                       <  1:110621/51‑1 (MQ=255)
         cGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCa                                      <  1:1614333/51‑1 (MQ=255)
         cGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCa                                      <  1:5284411/51‑1 (MQ=255)
         cGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCa                                      >  1:6624644/1‑51 (MQ=255)
          gAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCAt                                     <  1:1427537/51‑1 (MQ=255)
          gAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCAt                                     <  1:4390007/51‑1 (MQ=255)
          gAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCAt                                     <  1:2271787/51‑1 (MQ=255)
          gAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCAt                                     <  1:5098390/51‑1 (MQ=255)
          gAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACCGCGGGATACCAt                                     >  1:3359958/1‑51 (MQ=255)
           aaGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATg                                    <  1:1034430/51‑1 (MQ=255)
           aaGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATg                                    <  1:1287218/51‑1 (MQ=255)
           aaGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATg                                    <  1:6601053/51‑1 (MQ=255)
           aaGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATg                                    >  1:5506121/1‑51 (MQ=255)
            aGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGc                                   >  1:3026867/1‑51 (MQ=255)
            aGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGc                                   >  1:7110419/1‑51 (MQ=255)
              ttAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTg                                 >  1:4617462/1‑51 (MQ=255)
               tAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTgg                                <  1:3185900/51‑1 (MQ=255)
               tAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTgg                                <  1:6698852/51‑1 (MQ=255)
               tAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTgg                                >  1:585069/1‑51 (MQ=255)
               tAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTgg                                >  1:3147443/1‑51 (MQ=255)
               tAACCGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTgg                                <  1:4991908/51‑1 (MQ=255)
                aaCCGGCGCAAGGTTCAGTTGCGGGAATGTACCGCGGGATACCATGCTGGt                               >  1:4117331/1‑51 (MQ=255)
                 aCCGGCGCAAGGTTCAGTTGCGGGAATGTACCGCGGGATACCATGCTGGTg                              >  1:1956430/1‑51 (MQ=255)
                  ccGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGa                             >  1:3795945/1‑51 (MQ=255)
                  ccGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGa                             >  1:4832860/1‑51 (MQ=255)
                   cGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGAt                            <  1:5117704/51‑1 (MQ=255)
                   cGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGAt                            >  1:6600480/1‑51 (MQ=255)
                   cGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGAt                            <  1:6016272/51‑1 (MQ=255)
                   cGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGAt                            <  1:3437427/51‑1 (MQ=255)
                   cGGCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGAt                            >  1:6962455/1‑51 (MQ=255)
                    ggCGCAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATg                           <  1:2993805/51‑1 (MQ=255)
                    ggCGCAAGGTTCAGTTGCGGGAATGTACCGCGGGATACCATGCTGGTGATg                           >  1:6132270/1‑51 (MQ=255)
                     gcgcAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGc                          >  1:6830774/1‑51 (MQ=255)
                     gcgcAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGc                          >  1:4326473/1‑51 (MQ=255)
                     gcgcAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGc                          <  1:6321389/51‑1 (MQ=255)
                       gcAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCAc                        >  1:3501804/1‑51 (MQ=255)
                       gcAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCAc                        >  1:1035055/1‑51 (MQ=255)
                       gcAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCAc                        >  1:3267797/1‑51 (MQ=255)
                        cAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACt                       >  1:6541466/1‑51 (MQ=255)
                        cAAGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACt                       <  1:6031757/51‑1 (MQ=255)
                         aaGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTc                      >  1:4786345/1‑51 (MQ=255)
                         aaGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTc                      <  1:3171561/51‑1 (MQ=255)
                          aGGTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCa                     <  1:2504627/51‑1 (MQ=255)
                           ggTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCAc                    >  1:1067871/1‑51 (MQ=255)
                            gTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACg                   <  1:2819875/51‑1 (MQ=255)
                            gTTCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACg                   <  1:5538166/51‑1 (MQ=255)
                             ttCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGa                  <  1:4817066/51‑1 (MQ=255)
                             ttCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGa                  <  1:1915151/51‑1 (MQ=255)
                             ttCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGa                  <  1:4526804/51‑1 (MQ=255)
                              tCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCa                     >  1:113560/1‑47 (MQ=255)
                              tCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAg                 <  1:3782379/51‑1 (MQ=255)
                              tCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAg                 <  1:3832358/51‑1 (MQ=255)
                              tCAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAg                 <  1:3028207/51‑1 (MQ=255)
                               cAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGc                <  1:1387611/51‑1 (MQ=255)
                               cAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGc                >  1:5063398/1‑51 (MQ=255)
                               cAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGc                <  1:3253341/51‑1 (MQ=255)
                               cAGTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGc                >  1:3755716/1‑51 (MQ=255)
                                 gTTGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCAt              <  1:5512507/51‑1 (MQ=255)
                                  ttGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATa             >  1:1578216/1‑51 (MQ=255)
                                   tGCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATAc            <  1:6867980/51‑1 (MQ=255)
                                    gCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACg           <  1:4990586/51‑1 (MQ=255)
                                    gCGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACg           >  1:2505363/1‑51 (MQ=255)
                                     cGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACgg          <  1:526675/51‑1 (MQ=255)
                                     cGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACgg          >  1:6242792/1‑51 (MQ=255)
                                     cGGGAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACgg          >  1:6188168/1‑51 (MQ=255)
                                      gggAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACggg         >  1:6430238/1‑51 (MQ=255)
                                      gggAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACggg         <  1:5274790/51‑1 (MQ=255)
                                      gggAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACggg         >  1:5316233/1‑51 (MQ=255)
                                      gggAATGTACCGCGGGATACCATGCTGGTGATGCACTCACGAGCATACggg         >  1:3641746/1‑51 (MQ=255)
                                       ggAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGa        >  1:6259834/1‑51 (MQ=255)
                                       ggAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGa        >  1:4818710/1‑51 (MQ=255)
                                       ggAATGTACCGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGa        >  1:4816546/1‑51 (MQ=255)
                                        gAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGaa       <  1:1056563/51‑1 (MQ=255)
                                        gAATGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGaa       <  1:6901649/51‑1 (MQ=255)
                                          aTGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGAACa     <  1:1863623/51‑1 (MQ=255)
                                          aTGTACTGCGGGATACCATACTGGTGATGCACTCACGAGCATACGGGAACa     <  1:4546027/51‑1 (MQ=255)
                                           tGTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGAACAg    >  1:693714/1‑51 (MQ=255)
                                            gTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGAACAGa   >  1:141000/1‑51 (MQ=255)
                                            gTACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGAACAGa   <  1:2018159/51‑1 (MQ=255)
                                             tACTGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGAACAGaa  >  1:6378489/1‑51 (MQ=255)
                                                |                                               
ACAGCGCATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACCGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGAACAGAA  >  REL606/3724486‑3724581

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: