Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 648,332 A→T 11.1% V189E (GTG→GAG)  mrdB ← cell wall shape‑determining protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606648,3320A→T11.1% 14.3 ‑1.5 ‑0.763566mrdBcell wall shape‑determining protein
Reads supporting (aligned to +/- strand):  new base (6/1):  ref base (32/24):  total (38/25)
Fisher's exact test strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.04e-01

AGAACCACAGAATCGGAATGAACGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCA  >  REL606/648297‑648367
                                   |                                   
agaaCCACAGAATCGGAATGAACGCCGCTACCAGCa                                     >  1:2715540/1‑36 (MQ=255)
agaaCCACAGAATCGGAATGAACGCCGCTACCAGCa                                     <  1:5205610/36‑1 (MQ=255)
 gaaCCACAGAATCGGAATGAACGCCGCTACCAGCac                                    <  1:5271490/36‑1 (MQ=255)
  aaCCACAGAATCGGAATGATCGCCGCTACCAGCact                                   >  1:2641008/1‑36 (MQ=255)
  aaCCACAGAATCGGAATGAACGCCGCTACCAGCact                                   >  1:1328588/1‑36 (MQ=255)
  aaCCACAGAATCGGAATGAACGCCGCTACCAGCact                                   >  1:2287517/1‑36 (MQ=255)
   aCCACAGAATCGGAATGAACGCCGCTACCAGCacta                                  >  1:4375914/1‑36 (MQ=255)
   aCCACAGAATCGGAATGAACGCCGCTACCAGCacta                                  <  1:5611517/36‑1 (MQ=255)
   aCCACAGAATCGGAATGAACGCCGCTACCAGCacta                                  >  1:5886459/1‑36 (MQ=255)
    ccACAGAATCGGAATGAACGCCGCTACCAGCactac                                 >  1:2187745/1‑36 (MQ=255)
    ccACAGAATCGGAATGAACGCCGCTACCAGCactac                                 <  1:3147675/36‑1 (MQ=255)
     cacaGAATCGGAATGAACGCCGCTACCAGCactact                                <  1:4370578/36‑1 (MQ=255)
      acaGAATCGGAATGAACGCCGCTATCAGCACTACTg                               >  1:3903793/1‑36 (MQ=255)
      acaGAATCGGAATGAACGCCGCTACCAGCACTACTg                               <  1:3295281/36‑1 (MQ=255)
      aaaGAATCGGAATGAACGCCGCTACCAGCACTACTg                               <  1:5637391/34‑1 (MQ=255)
       caGAATCGGAATGAACGCCGCTACCAGCACTACTGc                              <  1:6695908/36‑1 (MQ=255)
       caGAATCGGAATGAACGCCGCTACCAGCACTACTGc                              <  1:3706577/36‑1 (MQ=255)
        aGAATCGGAATGAACGCCGCTACCAGCTCTACTGCg                             >  1:3963204/1‑36 (MQ=255)
        aGAATCGGAATGAACGCCGCTACCAGCTCTACTGCg                             >  1:2939769/1‑36 (MQ=255)
        aGAATCGGAATGAACGCCGCTACCAGCACTTCTGCg                             >  1:3837147/1‑36 (MQ=255)
          aaTCGGAATGAACGCCGCTACCAGCTCTACTGCGAc                           <  1:3398861/36‑1 (MQ=255)
          aaTCGGAATGAACGCCGCTACCAGCACTACTGCGAc                           >  1:3892725/1‑36 (MQ=255)
          aaTCGGAATGAACGCCGCTACCAGCACTACTGCGAc                           >  1:2202204/1‑36 (MQ=255)
          aaTCGGAATGAACGCCGCTACCAGCACTACTGCGAc                           >  1:6839245/1‑36 (MQ=255)
           aTCGGAATGAACGCCGCTACCAGCTCTACTGCGACg                          >  1:5401307/1‑36 (MQ=255)
           aTCGGAATGAACGCCGCTACCAGCACTACTGCGACg                          >  1:6842023/1‑36 (MQ=255)
           aTCGGAATGAACGCCGCTACCAGCACTACTGCGACg                          >  1:217672/1‑36 (MQ=255)
           aTCGGAATGAACGCCGCTACCAGCACTACTGCGACg                          >  1:5660031/1‑36 (MQ=255)
             cGGAATGAACGCCGCTACCAGCACTACTGCGACGcc                        >  1:5316507/1‑36 (MQ=255)
              ggAATGAACGCCGCTACCAGCACTACTGCGACGCCa                       >  1:6319833/1‑36 (MQ=255)
              ggAATGAACGCCGCTACAAGCTATACTGCGACGCCa                       >  1:6638702/1‑36 (MQ=25)
               gAATGAACGCCGCTACCAGCTCTACTGCGACGCCaa                      >  1:918360/1‑36 (MQ=255)
               gAATGAACGCCGCTACCAGCACTACTGCGACGccga                      <  1:1731547/36‑3 (MQ=255)
               gAATGAACGCCGCTACCAGCACTACTGCGACGCCaa                      >  1:1640053/1‑36 (MQ=255)
               gAATGAACGCCGCTACCAGCACTACTGCGACGCCaa                      >  1:3747618/1‑36 (MQ=255)
               gAATGAACGCCGCTACCAGCACTACTGCGACGCCaa                      >  1:5209908/1‑36 (MQ=255)
                aaTGAACGCCGCTACCAGCACTACTGCGACGCCAAt                     >  1:5017219/1‑36 (MQ=255)
                aaTAAACGCCGCTACCAGCACTACTGCGACGCCAAt                     <  1:2856332/36‑1 (MQ=255)
                 aTGAACGCCGCTACCAGCACTACTGCGACGCCAATc                    <  1:4237511/36‑1 (MQ=255)
                 aTGAACGCCGCTACCAGCACTACTGCGACGCCAATc                    >  1:480605/1‑36 (MQ=255)
                 aTGAACGCCGCTACCAGCACTACTGCGACGCCAATc                    >  1:3126787/1‑36 (MQ=255)
                 aTGAACGCCGCTACCAGCACTACTGCGACGCCAATc                    <  1:2517172/36‑1 (MQ=255)
                  tGAACGCCGCTACCAGCACTACTGCGACGCCAATCa                   <  1:942969/36‑1 (MQ=255)
                  tGAACGCCGCTACCAGCACTACTGCGACGCCAATCa                   <  1:2483819/36‑1 (MQ=255)
                    aaCGCCGCTACCAGCTCTACTGCGACGCCAATCAGa                 >  1:5329720/1‑36 (MQ=255)
                    aaCGCCGCTACCAGCACTACTGCGACGCCAATCAGa                 <  1:1235368/36‑1 (MQ=255)
                    aaCGCCGCTACCAGCACTACTGCGACGCCAATCAGa                 >  1:281124/1‑36 (MQ=255)
                    aaCGCCGCTAACAACACTACTGCGACGCCAATCAGa                 <  1:722799/36‑1 (MQ=37)
                     aCGCCGCTACCAGCACTACTTCGACGCCAATCAGAc                <  1:3787633/36‑1 (MQ=255)
                     aCGCCGCTACCAGCACTACTGCGACGCCAATCAGAc                <  1:6705746/36‑1 (MQ=255)
                       gccgcTACCAGCACTACTGCGACGCCAATCAGACGc              <  1:7157647/36‑1 (MQ=255)
                         cgcTACCAGCACTACTGCGACGCCAATCAGACGCCa            <  1:7508481/36‑1 (MQ=255)
                          gcTACCAGCACTACTGCGACGCCAATCAGACGCCAg           >  1:4574216/1‑36 (MQ=255)
                          gcTACCAGCACTACTGCGACGCCAATCAGACGCCAg           >  1:3560186/1‑36 (MQ=255)
                           cTACCAGCACTACTGCGACGCCAATCAGACGCCAGc          >  1:1603841/1‑36 (MQ=255)
                           cTACCAGCACTACTGCGACGCCAATCAGACGCCAGc          <  1:724004/36‑1 (MQ=255)
                           cTACCAGCACTACTGCGACGCCAATCAGACGCCAGc          >  1:3312186/1‑36 (MQ=255)
                             aCCAGCACTACTGCGACGCCAATCAGACGCCAGCTa        >  1:7044331/1‑36 (MQ=255)
                             aCCAGCACTACTGCGACGCCAATCAGACGCCAGCTa        >  1:2203302/1‑36 (MQ=255)
                             aCCAGCACTACTGCGACGCCAATCAGACGCCAGCTa        >  1:4043887/1‑36 (MQ=255)
                             aCCAGCACTACTGCGACGCCAATCAGACGCCAGCTa        <  1:6193723/36‑1 (MQ=255)
                             aCCAGCACTACTGCGACGCCAATCAGACGCCAGCTa        <  1:1580497/36‑1 (MQ=255)
                             aCCAGCACTACTGCGACGCCAATCAGACGCCAGCTa        <  1:1104473/36‑1 (MQ=255)
                              ccAGCACTACTGCGACGCCAATCAGACGCCAGCTaa       >  1:1388809/1‑36 (MQ=255)
                              ccAGCACTACTGCGACGCCAATCAGACGCCAGCTaa       >  1:6741403/1‑36 (MQ=255)
                              ccAGCACTACTGCGACGCCAATCAGACGCCAGCTaa       <  1:2171934/36‑1 (MQ=255)
                              ccAGCACTACTGCGACGCCAATCAGACGCCAGCTaa       >  1:4263310/1‑36 (MQ=255)
                               cAGCACTACTGCGACGCCAATCAGACGCCAGCTAAg      <  1:260320/36‑1 (MQ=255)
                                aGCACTACTGCGACGCCAATCAGACGCCAGCTAAgg     >  1:2125360/1‑36 (MQ=255)
                                aGCACTACTGCGACGCCAATCAGACGCCAGCTAAgg     >  1:5498307/1‑36 (MQ=255)
                                 gCACTACTGCGACGCCAATCAGACGCCAGCTAAGGc    <  1:1486210/36‑1 (MQ=255)
                                 gCACTACTGCGACGCCAATCAGACGCCAGCTAAGGc    >  1:296326/1‑36 (MQ=255)
                                 gCACTACTGCGACGCCAATCAGACGCCAGCTAAGGc    <  1:434047/36‑1 (MQ=255)
                                 gCACTACTGCGACCCCAATCAGACGCCAGCTAAGGc    >  1:6471498/1‑36 (MQ=255)
                                  cACTACTGCGACGNCAATCAGACGCCAGCTAAGGcc   >  1:2324403/1‑36 (MQ=255)
                                   actactGCGACGCCAATCAGACGCCAGCTAAGGCCa  >  1:1771191/1‑36 (MQ=255)
                                   actactGCGACGCCAATCAGACGCCAGCTAAGGCCa  <  1:3247479/36‑1 (MQ=255)
                                   actactGCGACGCCAATCAGACGCCAGCTAAGGCCa  >  1:5005539/1‑36 (MQ=255)
                                   |                                   
AGAACCACAGAATCGGAATGAACGCCGCTACCAGCACTACTGCGACGCCAATCAGACGCCAGCTAAGGCCA  >  REL606/648297‑648367

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