Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,037,598 +C 10.0% coding (1397/2061 nt) ECB_02838 ← GspD, hypothetical type II secretion protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,037,5981.→C10.0% 4.7 ‑2.7 ‑0.030coding (1397/2061 nt)ECB_02838GspD, hypothetical type II secretion protein
Reads supporting (aligned to +/- strand):  new base (0/3):  ref base (17/10):  total (17/13)
Fisher's exact test strand distribution p-value = 7.04e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.71e-01

GCGTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.CCGCCTGCACCAGCGCCATCCAGTCGCCTTTCAC  >  REL606/3037563‑3037632
                                    |                                  
gCGTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTaa                                     >  1:5603246/1‑36 (MQ=255)
gCGTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTaa                                     <  1:3173111/36‑1 (MQ=255)
gCGTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTaa                                     >  1:7936567/1‑36 (MQ=255)
 cGTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.c                                   <  1:781226/36‑1 (MQ=255)
  ttGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.cc                                  <  1:2031206/35‑1 (MQ=255)
  gTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.ct                                  >  1:3570532/1‑35 (MQ=255)
  gTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.cc                                  >  1:3060695/1‑36 (MQ=255)
  gTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.cc                                  <  1:7261109/36‑1 (MQ=255)
  gTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.cc                                  >  1:1172183/1‑36 (MQ=255)
   tGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.CCg                                 <  1:4468340/36‑1 (MQ=255)
   tGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.CCg                                 <  1:5844361/36‑1 (MQ=255)
     gagaGCACGTTAGAGCTGGAGTCGTTTTTAA.CCGcc                               <  1:6839221/36‑1 (MQ=255)
      agagCACGTTAGAGTTGGAGTCGTTTTTAA.CCGCCt                              <  1:7331458/36‑1 (MQ=255)
      agagCACGTTAGAGCTGGAGTCGTTTTTAA.CCGCCt                              <  1:1419044/36‑1 (MQ=255)
       gagCACGTTAGAGCTGGAGTCGTTTTTAA.CCGCCTg                             >  1:7528773/1‑36 (MQ=255)
           aCGTTAGGGCTGGAGTCGTTTTTAA.CCGCCTGCAcc                         >  1:7160977/1‑36 (MQ=255)
           aCGTTAGAGCTGGAGTCGTTTTTAA.CCGCCTGCAcc                         >  1:83256/1‑36 (MQ=255)
           aCGTTAGAGCTGGAGTCGTTTTTAA.CCGCCTGCAcc                         >  1:6479151/1‑36 (MQ=255)
            cGTTAGAGCTGGAGTTTTTTTAAA.CCGCCTGCACCa                        <  1:3434824/36‑1 (MQ=25)
            cGTTAGAGCGGGAGTCTTTTTAAA.CCGCCTGCACCa                        <  1:1032492/36‑1 (MQ=37)
              ttAGAGCTGGAGTCGTTTTTAA.CCGCCTGCACCAgc                      >  1:220436/1‑36 (MQ=255)
              ttAGAGCTGGAGTCGTTTTTAA.CCGCCTGCACCAgc                      >  1:5504804/1‑36 (MQ=255)
                agagCTGGAGTCGTTTTTAA.CCGCCTGCACCAgcgc                    >  1:6392106/1‑36 (MQ=255)
                  agCTGGAGTCGTTTTTAA.CCGCCTGCACCAGCGCCa                  >  1:7737602/1‑36 (MQ=255)
                  agCTGGAGTCGTTTTTAA.CCGCCTGCACCAGCGCCa                  >  1:426643/1‑36 (MQ=255)
                  agCTGGAGTCGTTTTTAA.CCGCCTGCACCAGCGCCa                  >  1:3502609/1‑36 (MQ=255)
                   gCTGGAGTCGTTTTTAA.GCCCACGCACCTGCGCCAt                 >  1:1361100/1‑36 (MQ=12)
                    cTGGAGTCGTTTTTAA.CCGCCTGCACCAGCGCCATc                >  1:8048293/1‑36 (MQ=255)
                    cTGGAGTCGTTTTTAA.CCGCCTGCACCAGCGCCATc                >  1:2166376/1‑36 (MQ=255)
                     tGGAGTCGTTTTTAA.CCGCCTGCACCAGCGCCATcc               >  1:975816/1‑36 (MQ=255)
                       gAGTCTTTTTTAC.CCGCCTGCACCAGCGCCATCCAg             <  1:1655518/36‑1 (MQ=25)
                        aGTCGTTTTTAA.CCGCCTGCACCAGCGCCATCCAGt            >  1:3497008/1‑36 (MQ=255)
                          tttttttaAA.CCGCCTGCACCAGCGCCATCCAGTCg          <  1:1309023/33‑1 (MQ=37)
                           tttttttAACCCGCCTGCACCAGCGCCATCCAGTCg          <  1:3523572/34‑1 (MQ=21)
                            gTTTTTAA.CCGCCTGCACCAGCGCCATCCAGTCGcc        >  1:2816589/1‑36 (MQ=255)
                            gTTTTTAA.CCGCCTGCACCAGCGCCATCCAGTCGcc        >  1:355297/1‑36 (MQ=255)
                             tttttAACCCGCCTGCACCAGCGCCATCCAGTCGcc        <  1:5906692/36‑1 (MQ=25)
                             tttttAACCCGCCTGCACCAGCGCCATCCAGTCGcc        <  1:4568119/36‑1 (MQ=25)
                             ttttaAC.CCGCCTGCACCAGCGCCATCCAGTCGCCt       <  1:2183870/36‑1 (MQ=25)
                                 aaa.CCGCCTGCACCAGCGCCATCCAGTCGCCTTTca   <  1:5368277/35‑1 (MQ=255)
                                  aa.CCGCCTGCACCAGCGCCATCCAGTCGCCTTTcac  >  1:3974991/1‑36 (MQ=255)
                                    |                                  
GCGTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.CCGCCTGCACCAGCGCCATCCAGTCGCCTTTCAC  >  REL606/3037563‑3037632

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG