Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,020,844 G→T 25.1% Q214K (CAA→AAA)  ECB_02820 ← hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,020,8440G→T25.1% 2.7 ‑3.5 0.028640ECB_02820hypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (11/9):  total (19/9)
Fisher's exact test strand distribution p-value = 2.92e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AAATATATATATTTATTATCAGCGGTAATGGCTTGCAATACATAATTATCATTGGTTAAACCTGCATCGA  >  REL606/3020810‑3020879
                                  |                                   
aaaTATATATATTTATTATCAGCGGTAATGGCTTGc                                    <  1:7891001/36‑1 (MQ=255)
 aaTATATATATTTATTATCAGCGGTAATGGCTTGCa                                   <  1:4423116/36‑1 (MQ=255)
   tatatatatTTATTATCAGCGGTAATGGCTTGCAAt                                 <  1:3153069/36‑1 (MQ=255)
     tatatatTTATTATCAGCGGTAATGGCTTTTAATaa                               >  1:2972652/1‑35 (MQ=25)
     tatatatTTATTATCAGCGGTAATGGCTTTCAATAc                               >  1:7454100/1‑36 (MQ=255)
      atatatTTATTATCAGCGGTAATGGCTTGCAATACa                              <  1:7611060/36‑1 (MQ=255)
      atatatTTATTATCAGCGGTAATGGCTTGCAATACa                              <  1:4751194/36‑1 (MQ=255)
        atatTTATTATCAGCGGTAATGGCTTTTAATACATa                            >  1:515947/1‑36 (MQ=37)
        atatTTATTATCAGCGGTAATGGCTTGCAATACATa                            <  1:4909489/36‑1 (MQ=255)
          atttattATCAGCGGTAATGGCTTTCAATACATAAt                          >  1:7413408/1‑36 (MQ=255)
          atttattATCAGCGGTAATGGCTTGCAATACATAAt                          <  1:1827706/36‑1 (MQ=255)
           tttattATCAGCGGTAATTGCTTGCAATACATAAtt                         >  1:6979168/1‑36 (MQ=255)
            ttattaTCAGCGGTAATGGCTTTTAATACATAATTa                        >  1:5573588/1‑36 (MQ=37)
            ttattaTCAGCGGTAATGGCTTTGAATATATAAttt                        >  1:3120571/1‑35 (MQ=21)
             tattatCAGCGGTAATGGCTTGCAATACATAATGAt                       >  1:7138191/1‑36 (MQ=255)
             tattatCAGCGGTAATGGCTGGCGAAACAAAATTAt                       >  1:5911316/1‑36 (MQ=16)
             tattatCAGCGGTAATGGATTTTTATAAATAATTaa                       >  1:1880561/1‑35 (MQ=12)
              attatCAGCGGGAATGGCTTTTAATACATAATTAtt                      >  1:6949518/1‑35 (MQ=25)
               ttatCAGCGGTAATGGCTTGCAATACATAATTATCa                     <  1:1167054/36‑1 (MQ=255)
                 atCAGCGGTAATGGCTTGCAATACATAATTATCAtt                   <  1:5533852/36‑1 (MQ=255)
                  tCAGCGGTAATGGCTTGCAATACATAATTATCATTg                  >  1:1086557/1‑36 (MQ=255)
                   cAGCGGTAATGGCTTGCAATACATCATTATCATTgg                 >  1:493895/1‑36 (MQ=255)
                     gCGGTAATGGCTTGCCATACATAATTATCATTGGtt               >  1:5260594/1‑36 (MQ=255)
                      cGGTAATGGCTTGCAATACATAATTATCATTGGTTa              >  1:3678250/1‑36 (MQ=255)
                         tAATGGCTTGCAACACATAATTATCATTGGTTAAAc           >  1:2966065/1‑36 (MQ=255)
                           aTGGCTTGCAATACATAATTATCATTGGTTAAACCt         >  1:1286524/1‑36 (MQ=255)
                             ggCTTGCAATACATAATTATCATTGGTTAAACCTGc       >  1:3247439/1‑36 (MQ=255)
                                 tGCAATACATAATTATCATTGGTTAAACCTGCATCg   >  1:7146243/1‑36 (MQ=255)
                                  gCAATACATAATTATCATTGGTTAAACCTGCATCGa  <  1:5423739/36‑1 (MQ=255)
                                  |                                   
AAATATATATATTTATTATCAGCGGTAATGGCTTGCAATACATAATTATCATTGGTTAAACCTGCATCGA  >  REL606/3020810‑3020879

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