Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,483,047 C→A 100% R455S (CGC→AGC)  malT → transcriptional regulator MalT

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,483,0470C→A100.0% 185.9 561363malTtranscriptional regulator MalT
Reads supporting (aligned to +/- strand):  new base (34/22):  ref base (0/0):  total (34/22)

GAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCG  >  REL606/3483014‑3483079
                                 |                                
gaagaCACCATGCACGCAGAATTTAACGCTCTgagc                                >  1:3081425/1‑36 (MQ=255)
 aagaCACCATGCACGCAGAATTTAACGCTCTgagcg                               >  1:7155609/1‑36 (MQ=255)
  agaCACCATGCACGCAGAATTTAACGCTCTgagcgc                              <  1:8432966/36‑1 (MQ=255)
  agaCACCATGCACGCAGAATTTAACGCTCTgagcgc                              >  1:1833518/1‑36 (MQ=255)
  agaCACCATGCACGCAGAATTTAACGCTCTgagcgc                              >  1:2545665/1‑36 (MQ=255)
  agaCACCATGCACGCAGAATTTAACGCTCTgagcgc                              >  1:4421630/1‑36 (MQ=255)
   gaCACCATGCACGCAGAATTTAACGCTCTGAGCGcc                             >  1:3234289/1‑36 (MQ=255)
    acacCATGCACGCAGAATTTAACGCTCTGAGCGccc                            >  1:8001967/1‑36 (MQ=255)
    acacCATGCACGCAGAATTTAACGCTCTGAGCGccc                            >  1:3801344/1‑36 (MQ=255)
    acacCATGCACGCAGAATTTAACGCTCTGAGCGccc                            <  1:8184844/36‑1 (MQ=255)
      acCATGCACGCAGAATTTAACGCTCTGAGCGCCCAg                          <  1:5758109/36‑1 (MQ=255)
      acCATGCACGCAGAATTTAACGCTCTGAGCGCCCAg                          <  1:7216205/36‑1 (MQ=255)
      acCATGCACGCAGAATTTAACGCTCTGAGCGCCCAg                          >  1:3797550/1‑36 (MQ=255)
      acCATGCACGCAGAATTGAAAGCTCTGAGCGCCCAg                          >  1:660026/1‑36 (MQ=37)
       ccATGCACGCAGAATTTAACGCTCTGAGCGCCCAgg                         <  1:7458690/36‑1 (MQ=255)
       ccATGCACGCAGAATTTAACGCTCTGAGCGCCCAgg                         >  1:4901156/1‑36 (MQ=255)
        cATGCACGCAGAATTTAACGCTCTGAGCGCCCAGGt                        >  1:569528/1‑36 (MQ=255)
        cATGCACGCAGAATTTAACGCTCTGAGCGCCCAGGt                        <  1:163089/36‑1 (MQ=255)
         aTGCACGCAGAATTTAACGCTCTGAGCGCCCCGGTg                       >  1:2782627/1‑36 (MQ=37)
         aTGCACGCAGAATTTAACGCTCTGAGCGCCCAGGTg                       >  1:675560/1‑36 (MQ=255)
         aTGCACGCAGAATTTAACGCTCTGAGCGCCCAGGTg                       >  1:759074/1‑36 (MQ=255)
          tGCACGCAGAATTTAACGCTCTGAGCGCCCAGGTgg                      <  1:6703176/36‑1 (MQ=255)
           gCACGCAGAATTTAACGCTCTGAGCGCCCAGGTGGc                     <  1:6416807/36‑1 (MQ=255)
            cACGCAGAATTTAACGCTCTGAGCGCCCAGGTGGCg                    >  1:1306837/1‑36 (MQ=255)
            cACGCAGAATTTAACGCTCTGAGCGCCCAGGTGGCg                    >  1:7553721/1‑36 (MQ=255)
             aCGCAGAATTTAACGCTCTGAGCGCCCAGGTGGCGa                   >  1:19950/1‑36 (MQ=255)
             aCGCAGAATTTAACGCTCTGAGCGCCCAGCTGGCGa                   >  1:8410853/1‑36 (MQ=37)
               gCAGAATTTAACGCTCTGAGCGCCCAGGTGGCGAtt                 <  1:7068665/36‑1 (MQ=255)
               gCAGAATTTAACGCTCTGAGCGCCCAGGTGGCGAtt                 <  1:5689770/36‑1 (MQ=255)
                 aGAATTTAACGCTCTGAGCGCCCAGGTGGCGATTaa               <  1:4618144/36‑1 (MQ=255)
                 aGAATTTAACGCTCTGAGCGCCCAGGTGGCGATTaa               >  1:6046501/1‑36 (MQ=255)
                  gAATTTAACGCTCTGAGCGCCCAGGTGGCGATTAAc              <  1:3906867/36‑1 (MQ=255)
                  gAATTTAACGCTCTGAGCGCCCAGGTGGCGATTAAc              <  1:372890/36‑1 (MQ=255)
                  gAATTTAACGCTCTGAGCGCCCAGGTGGCGATTAAc              <  1:5378840/36‑1 (MQ=255)
                   aaTTTAACGCTCTGAGCGCCCAGGTGGCGATTAACg             <  1:2590790/36‑1 (MQ=255)
                    aTTTAACGCTCTGAGGGCCCCGGTGGCGATTAACGa            >  1:7829557/1‑36 (MQ=25)
                    aTTTAACGCTCTGAGCGCCCAGGTGGCGATTAAGGa            >  1:5448687/1‑36 (MQ=38)
                    aTTTAACGCTCTGAGCGCCCAGGTGGCGATTAACGa            >  1:5274850/1‑36 (MQ=255)
                    aTTTAACGCTCTGAGCGCCCAGGTGGCGATTAACGa            >  1:7999941/1‑36 (MQ=255)
                    aTTTAACGCTCTGAGCGCCCAGGTGGCGATTAACGa            >  1:2611487/1‑36 (MQ=255)
                     tttAACGCTCTGAGCGCCCCGGGGGCGGTTTACGAt           >  1:2756956/1‑36 (MQ=16)
                     tttAACGCTCTGAGCGCCCAGGTGGCGATTAACGAt           >  1:1635733/1‑36 (MQ=255)
                     tttAACGCTCTGAGCGCCCAGGTGGCGATTAACGAt           <  1:3212536/36‑1 (MQ=255)
                     tttAACGCTCTGAGCGCCCAGGTGGCGATTAACGAt           <  1:3898922/36‑1 (MQ=255)
                      ttAACGCTCTGAGCGCCCAGGTGGCGATTAACGATg          >  1:3751389/1‑36 (MQ=255)
                      ttAACGCTCTGAGCGCCCAGGTGGCGATTAACGATg          <  1:7807808/36‑1 (MQ=255)
                       tAACGCTCTGAGCGCCCAGGTGGCGATTAACGATgg         >  1:3932167/1‑36 (MQ=255)
                         aCGCTCTGAGCGCCCAGGTGGCGATTAACGATGGTa       >  1:3203212/1‑36 (MQ=255)
                         aCGCTCTGAGCGCCCAGGTGGCGATTAACGATGGTa       >  1:2717691/1‑36 (MQ=255)
                         aCGCTCTGAGCGCCCAGGTGGCGATTAACGATGGTa       >  1:4849493/1‑36 (MQ=255)
                          cGCTCTGAGCGCCCAGGTGGCGATTAACGATGGTaa      <  1:4172861/36‑1 (MQ=255)
                           gCTCTGAGCGCCCAGGTGGCGATTAACGATGGTAAt     <  1:7901043/36‑1 (MQ=255)
                           gCTCTGAGCGCCCAGGTGGCGATTAACGATGGTAAt     >  1:2034233/1‑36 (MQ=255)
                           gCTCTGAGCGCCCAGGTGGCGATTAACGATGGTAAt     <  1:1935740/36‑1 (MQ=255)
                            ctctGAGCGCCCAGGTGGCGATTAACGATGGTAATc    >  1:8138126/1‑36 (MQ=255)
                            ctctGAGCGCCCAGGTGGCGATTAACGATGGTAATc    >  1:4734124/1‑36 (MQ=255)
                             tctGAGCGCCCAGGTGGCGATTAACGATGGTAATcc   >  1:4886504/1‑36 (MQ=255)
                             tctGAGCGCCCAGGTGGCGATTAACGATGGTAATcc   >  1:3156166/1‑36 (MQ=255)
                             actGAGCGCCCAGGTGGCGATTAACGATGGTAATcc   <  1:1276931/35‑1 (MQ=255)
                              ctGAGCGCCCAGGTGGCGATTAACGATGGTAATCCg  <  1:7379444/36‑1 (MQ=255)
                              ctGAGCGCCCAGGTGGCGATTAACGATGGTAATCCg  >  1:4197737/1‑36 (MQ=255)
                              ctGAGCGCCCAGGTGGCGATTAACGATGGTAATCCg  >  1:2571359/1‑36 (MQ=255)
                                 |                                
GAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCG  >  REL606/3483014‑3483079

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG