Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,037,598 +C 12.8% 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.→C12.8% 11.6 ‑3.9 0.039coding (1397/2061 nt)ECB_02838GspD, hypothetical type II secretion protein
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (20/14):  total (20/19)
Fisher's exact test strand distribution p-value = 2.02e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GCGTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.CCGCCTGCACCAGCGCCATCCAGTCGCCTTTCACC  >  REL606/3037563‑3037633
                                    |                                   
gCGTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTaa                                      <  1:947291/36‑1 (MQ=255)
 cGTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.c                                    <  1:4093135/36‑1 (MQ=255)
  gTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.cc                                   <  1:1533716/36‑1 (MQ=255)
  gTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.cc                                   >  1:2555468/1‑36 (MQ=255)
   tGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.CCg                                  <  1:7160256/36‑1 (MQ=255)
      agagCACGTTAGAGCTGGAGTCGTTTTTAA.CCGCCt                               <  1:1209241/36‑1 (MQ=255)
      agagCACGTTAGAGCTGGAGTCGTTTTTAA.CCGCCt                               <  1:1836286/36‑1 (MQ=255)
       gagCACGTTAGAGCTGGAGTCGTTTTTAA.CCGccgg                              >  1:4946423/1‑34 (MQ=255)
       gagCACGTTAGAGCTGGAGTCGTTTTTAA.CCGCCTg                              <  1:4646541/36‑1 (MQ=255)
        agCACGTTAGAGCTGGAGTCGTTTTTAA.CCGCCTGc                             >  1:5900135/1‑36 (MQ=255)
          cACGTTAGAGCTGGAGTCGTTTTTAA.CCGCCTGCAc                           >  1:5137836/1‑36 (MQ=255)
          cACGTTAGAGCGGGAGTCTTTTTTAA.CCGCCTGCAc                           <  1:7424569/36‑1 (MQ=37)
           aCGTAAGGGCTGGAGTCGTTTTAAA.CCGCCTGCAcc                          <  1:4015052/36‑1 (MQ=25)
             gTTAGAGCTGGAGTCGTTTTTAA.CCGCCTGCACCAg                        >  1:6159837/1‑36 (MQ=255)
             gTTAGAGCTGGAGTCGTTTTTAA.CCGCCTGCACCAg                        >  1:6407908/1‑36 (MQ=255)
               tAGAGCTGGAGTCGTTTTTAA.CCGCCTGCACCAgcg                      >  1:646267/1‑36 (MQ=255)
               tAGAGCTGGAGTCGTTTTTAA.CCGCCTGCACCAgcg                      >  1:4958216/1‑36 (MQ=255)
                agagCTGGAGTCGTTTTTAA.CCGCCTGCACCAgcgc                     >  1:7497480/1‑36 (MQ=255)
                agagCGGGAGTTTTTTTAAA.CCGCCTGCACCAgcgc                     <  1:2877340/36‑1 (MQ=16)
                 gagCTGGAGTCGTTTTTAA.CCGCCTGCACCAGCGcc                    >  1:3138805/1‑36 (MQ=255)
                 gagCTGGAGTCGTTTTTAA.CCGCCTGCACCAGCGcc                    >  1:2084917/1‑36 (MQ=255)
                 gagCGGGATTCGTTTTAAA.CCGCCTGCACCAGCGcc                    <  1:628144/36‑1 (MQ=25)
                   ccgggATTCTTTTTAAC.CCGCCTGCACCAGCGCCAt                  <  1:6396253/33‑1 (MQ=9)
                    cTGGAGTCGTTTTTAA.CCGCCTGCACCAGCGCCATc                 >  1:803251/1‑36 (MQ=255)
                     tGGAGTCTTTTTTAA.CCGCCTGCACCAGCGCCATcc                <  1:4946317/36‑1 (MQ=255)
                       gAGTCTTTTTAAA.CCGCCTGCACCAGCGCCATCCAg              <  1:4053734/36‑1 (MQ=37)
                        aGTCGTTTTTAA.CCGCCTGCACCAGCGCCATCCAGt             >  1:150068/1‑36 (MQ=255)
                         gTCGTTTTTAA.CCGCCTGCACCAGCGCCATCCAGTc            >  1:6651836/1‑36 (MQ=255)
                          ctctttttAACCCGCCTGCACCAGCGCCATCCAGTc            <  1:2479887/33‑1 (MQ=21)
                           tttttttAC.CCGCCTGCACCAGCGCCATCCAGTCGc          <  1:7321300/34‑1 (MQ=25)
                           cGTTTTTAA.CCGCCTGCACCAGCGCCATCCAGTCGc          >  1:7778589/1‑36 (MQ=255)
                           cGTTTTTAA.CCGCCTGCACCAGCGCCATCCAGTCGc          >  1:6113406/1‑36 (MQ=255)
                            gTTTTTAA.CCGCCTGAACCAGCGCCATCCAGTCGcc         >  1:1688620/1‑36 (MQ=255)
                             tttttAACCCGCCTGCACCAGCGCCTTCCAGTCGcc         <  1:5096557/36‑1 (MQ=21)
                             tttttAACCCGCCTGCACCAGCGCCATCCAGTCGcc         <  1:5270872/36‑1 (MQ=25)
                             tttttAACCCGCCTGCACCAGCGCCATCCAGTCGcc         <  1:7519772/36‑1 (MQ=25)
                             tttttAACCCGCCTGCACCAGCGCCATCCAGTCGcc         <  1:3685196/36‑1 (MQ=25)
                             tttttAA.CCGCCTGCACCAGCGCCATCCAGTCGCCt        >  1:6432126/1‑36 (MQ=255)
                             tttttAA.CCGCCTGCACCAGCGCCATCCAGTCGCCt        >  1:4813217/1‑36 (MQ=255)
                             ttttgAA.CCGCCTGCACCAGCGCAATCCAGTCGCCt        >  1:4309115/1‑36 (MQ=37)
                             ttttaAC.CCGCCTGCACCAGCGCCATCCAGTCGCCt        <  1:2140643/36‑1 (MQ=25)
                              ttttAA.CCGCCTGCACCAGCGCCATCCAGTCGCCtt       >  1:5366752/1‑36 (MQ=255)
                                 tAA.CCGCCTGCACCAGCGCCATCCAGTCGCCTTTca    <  1:5217033/36‑1 (MQ=255)
                                  aa.CCGCCTGCACCAGCGCCATCCAGTCGCCTTTcac   >  1:3750213/1‑36 (MQ=255)
                                  aa.CCGCCTGCACCAGCGCCATCCAGTCGCCTTTcac   >  1:1908572/1‑36 (MQ=255)
                                   a.CCGCCTGCACCAGCGCCATCCAGTCGCCTTTcacc  >  1:7414663/1‑36 (MQ=255)
                                   a.CCGCCTGCACCAGCGCCATCCAGTCGCCTTTcacc  >  1:4658482/1‑36 (MQ=255)
                                    |                                   
GCGTGGAGAGCACGTTAGAGCTGGAGTCGTTTTTAA.CCGCCTGCACCAGCGCCATCCAGTCGCCTTTCACC  >  REL606/3037563‑3037633

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG