Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,540,211 T→C 15.7% E77G (GAA→GGA)  ‡ iadA ← isoaspartyl dipeptidase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,540,2110T→C15.7% 2.2 ‑3.5 ‑0.036230iadAisoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (15/15):  total (15/21)
Fisher's exact test strand distribution p-value = 3.04e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.89e-01

GCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATC  >  REL606/4540176‑4540245
                                   |                                  
gcgcCACTTCCGGCGTGCGCGTCGTGGGACCTGCtt                                    >  1:7053108/1‑36 (MQ=255)
  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:5869611/1‑36 (MQ=255)
       ttCCGGCGTGCGCGTCGTGGGACCTGCTTCgccgcc                             <  1:5048126/36‑1 (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:6581187/36‑1 (MQ=255)
            gCGTGCGCGTCGTGGGACCTGCTTCGccgccaccgc                        <  1:2978317/36‑1 (MQ=255)
            gCGGGCGCGTCGTGGGACCTGCTTCGccgccaccgc                        <  1:4431735/36‑1 (MQ=255)
              gggcgcGTCGTGGGACCTGCTTCGccgccaccgcca                      <  1:2841312/34‑1 (MQ=255)
              gTGCGCGTCGTGGGACCTGCTTCGccgccaccgcca                      >  1:1387920/1‑36 (MQ=255)
              gTGCGCGTCGTGGGACCTGCTTCGccgccaccgcca                      >  1:1556354/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)
                                   |                                  
GCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATC  >  REL606/4540176‑4540245

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