Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,540,213 G→C 13.8% G76G (GGC→GGG)  ‡ iadA ← isoaspartyl dipeptidase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,540,2130G→C13.8% 2.2 ‑3.1 ‑0.034228iadAisoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (17/12):  total (17/17)
Fisher's exact test strand distribution p-value = 4.45e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.90e-01

GCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATC  >  REL606/4540178‑4540245
                                   |                                
gcCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTcg                                  <  1:4802982/36‑1 (MQ=255)
gcCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTcg                                  >  1:4279395/1‑36 (MQ=255)
   aCTTCCGGCGTGCGCGTCGTGGGACCTGCTTcgccg                               <  1:4965451/36‑1 (MQ=255)
   aCTTCCGGCGTGCGCGTCGTGGGACCTGCTTcgccg                               <  1:4390972/36‑1 (MQ=255)
     ttCCGGCGTGCGCGTCGTGGGACCTGCTTCgccgcc                             <  1:5048126/36‑1 (MQ=255)
     ttCCGGCGTGCGCGTCGTGGGACCTGCTTCgccgcc                             >  1:5869611/1‑36 (MQ=255)
       ccGGCGGGCGCGTCGTGGGACCTGCTTCGccgccac                           <  1:6494865/36‑1 (MQ=255)
        cGGCGTGCGCGTCGTGGGACCTGCTTCGccgccacg                          <  1:3399552/36‑2 (MQ=255)
        cGGCGTGCGCGTCGTGGGACCTGCTTCGccgccacc                          <  1:5932760/36‑1 (MQ=255)
         ggCGTGCGCGTCGTGGGACCTGCTTCGccgccaccg                         >  1:3068045/1‑36 (MQ=255)
          gCGTGCGCGTCGTGGGACCTGCTTCGccgccaccgc                        <  1:5845138/36‑1 (MQ=255)
          gCGTGCGCGTCGTGGGACCTGCTTCGccgccaccgc                        <  1:402886/36‑1 (MQ=255)
          gCGTGCGCGTCGTGGGACCTGCTTCGccgccaccgc                        <  1:2978317/36‑1 (MQ=255)
          gCGTGCGCGTCGTGGGACCTGCTTCGccgccaccgc                        <  1:6581187/36‑1 (MQ=255)
          gCGGGCGCGTCGTGGGACCTGCTTCGccgccaccgc                        <  1:4431735/36‑1 (MQ=255)
            gggcgcGTCGTGGGACCTGCTTCGccgccaccgcca                      <  1:2841312/34‑1 (MQ=255)
            gTGCGCGTCGTGGGACCTGCTTCGccgccaccgcca                      >  1:1556354/1‑36 (MQ=255)
            gTGCGCGTCGTGGGACCTGCTTCGccgccaccgcca                      >  1:1387920/1‑36 (MQ=255)
             tGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCaa                     <  1:1345779/36‑1 (MQ=255)
              gcgcGTCGTGGGACCTGCTTCGCCGCCACCGCCAAt                    >  1:1699223/1‑36 (MQ=255)
               cgcgTCGTGGGCCTGGCTCCGCCGCCACCGCCAATc                   <  1:1864350/36‑1 (MQ=19)
               cgcgGCGGGGGCCCTGCGCCGCCGCCACCGCCAATc                   <  1:562821/36‑1 (MQ=12)
                 cgtcgtGGGACCTGCTTCGCCGCCACCGCCAATCaa                 >  1:3152689/1‑36 (MQ=255)
                 cgtcggGGGCCTTGCTTCCCCCCCACCGCCAATCaa                 <  1:1650867/36‑1 (MQ=16)
                  cgcgtGGACCCTGCTCCCCCCCCACCGCCAATCaaa                <  1:4718265/34‑1 (MQ=12)
                   tggtGGCACCTGCTCCCCCCCCACCGCCAATCAAAt               <  1:4836235/34‑1 (MQ=16)
                    cgtGGGACCTGCTTCGCCGCCACCGCCAATCAAATg              >  1:6591049/1‑36 (MQ=255)
                    cgggggCCCTGCTCCCCCCCCACCGCCAATCAAATg              <  1:7470660/33‑1 (MQ=16)
                      tGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGa            >  1:4871541/1‑36 (MQ=255)
                       gggACCTGCTTCGCCGCCACCGCCAATCAAATGGAc           >  1:3690869/1‑36 (MQ=255)
                       ggcaCCTGCCCCCCCGCCACCGCCAATCAAATGGAc           <  1:333242/33‑1 (MQ=16)
                          aCCTGCTTCGCCGCCACCGCCAATCAAATGGACgtg        >  1:2230059/1‑36 (MQ=255)
                          aCCTGCTTCGCCGCCACCGCCAATCAAATGGACgtg        >  1:105417/1‑36 (MQ=255)
                            cTGCTTCGCCGCCACCGCCATTCAAATACACGTGtt      >  1:432175/1‑36 (MQ=21)
                             tGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTg     >  1:5696863/1‑36 (MQ=255)
                             tGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTg     >  1:2592800/1‑36 (MQ=255)
                                ttCGCCGCCACCGCCAATCACATGGACGTGTCGATc  >  1:6267594/1‑36 (MQ=37)
                                ttCGCCGCCACCGCCAATCAAATGGACGTGTTGATc  >  1:631916/1‑36 (MQ=255)
                                ttCGCCGCCACCGCCAATCAAATGGACGTGTTGATc  >  1:3801717/1‑36 (MQ=255)
                                ttCGCCGCCACCGCCAATCAAATGGACGTGTTGATc  >  1:7769215/1‑36 (MQ=255)
                                   |                                
GCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATC  >  REL606/4540178‑4540245

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 17 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG