Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,064,404 +G 5.1% coding (743/825 nt) ECB_02858 ← conserved hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,064,4041.→G 5.1% 4.0 ‑3.0 ‑0.639coding (743/825 nt)ECB_02858conserved hypothetical protein
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (7/30):  total (9/30)
Fisher's exact test strand distribution p-value = 4.86e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.61e-01

CGCTATCCAGCGCCACCGGTGCATCATACGACTGGC.GGTTTTCCATATCATCCGCCACCATGTCGTCGTGT  >  REL606/3064369‑3064439
                                    |                                   
cGCTATCCAGCGCCACCGGTGCATCATACGACTGGc                                      <  1:4230586/36‑1 (MQ=255)
cGCTATCCAGCGCCACCGGTGCATCATACGACTGGc                                      >  1:2970898/1‑36 (MQ=255)
 gCTATCCAGCGCCACCGGTGCATCATACGACTGGC.g                                    >  1:992156/1‑36 (MQ=255)
 gCTATCCAGCGCCACCGGTGCATCATACGACTGGC.g                                    >  1:5751312/1‑36 (MQ=255)
 gCTATCCAGCGCCACCGGTGCATCATACGACTGGC.g                                    >  1:1662464/1‑36 (MQ=255)
   tATCCAGCGCCACCGGTGCATCATACGACTGGC.GGt                                  <  1:4202488/36‑1 (MQ=255)
    aTCCAGCGCCACCGGTGCATCATACGACTGGC.GGtt                                 <  1:4684259/36‑1 (MQ=255)
     tCCAGCGCCACCGGTGCATCATACGACTGGC.GGttt                                <  1:3412433/36‑1 (MQ=255)
     tCCAGCGCCACCGGTGCATCATACGACTGGC.GGttt                                <  1:2771377/36‑1 (MQ=255)
      ccAGCGCCACCGGTGCATCATACGACTGGC.GGtttt                               <  1:2272945/36‑1 (MQ=255)
      ccAGCGCCACCGGTGCATCATACGACTGGC.GGcttt                               >  1:1898373/1‑36 (MQ=255)
        aGCGCCACCGGTGCATCATACGACTGGC.GGTTTTcc                             <  1:5442669/36‑1 (MQ=255)
           gcCACCGGTGCATCATACGACTGGCGGGTTTTccta                           >  1:5525677/1‑34 (MQ=37)
           gcCACCGGTGCATCATACGACTGGC.GGTTTTCCata                          <  1:3777637/36‑1 (MQ=255)
            ccACCGGTGCATCATACGACTGGCGGGTTTTccttt                          >  1:2186503/1‑33 (MQ=37)
            ccACCGGTGCATCATACGACTGGC.GGTTTTCCatat                         <  1:2730352/36‑1 (MQ=255)
             cACCGGTGCATCATACGACTGGC.GGTTTTCCATatc                        <  1:5040301/36‑1 (MQ=255)
             cACCGGTGCATCATACGACTGGC.GGTTTTCCATatc                        <  1:5667937/36‑1 (MQ=255)
              aCCGGTGCATCATACGTCTGGC.GGTTTTCCATatca                       <  1:1943566/36‑1 (MQ=255)
                cGGTGCATCATACGACTGGC.GGTTTTCCATatcatc                     <  1:6228177/36‑1 (MQ=255)
                 ggTGCATCATACGACTGGC.GGTTTTTAATCTCATcc                    >  1:6477667/1‑36 (MQ=25)
                 ggTGCATCATACGACTGGC.GGGTTTTCATTTTATcc                    >  1:5513523/1‑36 (MQ=16)
                   tGCATCATACGACTGGC.GGTTTTCTCTATCATcccc                  >  1:948286/1‑34 (MQ=37)
                   tGCATCATACGACTGGC.GGTTTTCCATATCATCCGc                  <  1:202544/36‑1 (MQ=255)
                   tGCATCATACGACTGGC.GGGTTTTCATTTCATcccc                  >  1:125367/1‑34 (MQ=21)
                     catcatACGACTTGC.GGTTTTCCATATCATCCGcca                >  1:467069/1‑36 (MQ=255)
                     catcatACGACTGGC.GGTTTTCCATATCATCCGcca                <  1:1870534/36‑1 (MQ=255)
                     catcatACGACTGGC.GGTTTTCCATATCATCCGcca                <  1:2608204/36‑1 (MQ=255)
                      atcatACGACAGGC.GGTTTTCCATATCATCCGccac               <  1:2873910/36‑1 (MQ=255)
                        catACGACTGGC.GGTTTTCCATATCATCCGccacca             <  1:4343314/36‑1 (MQ=255)
                        catACGACTGGC.GGTTTTCCATATCATCCGccacca             <  1:2680784/36‑1 (MQ=255)
                         atACGACTGGC.GGTTTTCCATATCATCCGCCACCAt            <  1:1056767/36‑1 (MQ=255)
                         atACGACTGGC.GGTTTTCCATATCATCCGCCACCAt            <  1:1277547/36‑1 (MQ=255)
                         atACGACTGGC.GGTTTTCCATATCATCCGCCACCAt            <  1:2776671/36‑1 (MQ=255)
                          tACGACTGGC.GGTTTTCCATATCATCCGCCACCATg           <  1:122693/36‑1 (MQ=255)
                             gACTGGC.GGTTTTCCATATCATCCGCCACCATgtcg        <  1:4833487/36‑1 (MQ=255)
                              aCTGGC.GGTTTTCCATATCATCCGCCACCATgtcgt       <  1:4773743/36‑1 (MQ=255)
                               cTGGC.GGTTTTCCATATCATCCGCCACCATgtcgtc      <  1:2077607/36‑1 (MQ=255)
                               cTGGC.GGTTTTCCATATCATCCGCCACCATgtcgtc      <  1:2453228/36‑1 (MQ=255)
                               cTGGC.GGTTTTCCATATCATCCGCCACCATgtcgtc      <  1:1144954/36‑1 (MQ=255)
                               cTGGC.GGTTTTCCATATCATCCGCCACCATgtcgtc      <  1:1627670/36‑1 (MQ=255)
                                tGGC.GGTTTTCCATTTCATCCCCCCCCCTGtcgtcg     >  1:5640510/1‑36 (MQ=19)
                                tGGC.GGTTTTCCATATCATCCGCCACCATGtcgtcg     <  1:1535141/36‑1 (MQ=255)
                                 ggC.GGTTTTCCATATCATCCGCNACCATGTcgtcgt    <  1:4679085/36‑1 (MQ=255)
                                   c.GGTTTTCCATATCATCCGCCACCATGTCGTCgtgt  <  1:573837/36‑1 (MQ=255)
                                   c.GGTTTTCCATATCATCCGCCACCATGTCGTCgtgt  <  1:5170793/36‑1 (MQ=255)
                                    |                                   
CGCTATCCAGCGCCACCGGTGCATCATACGACTGGC.GGTTTTCCATATCATCCGCCACCATGTCGTCGTGT  >  REL606/3064369‑3064439

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