Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,118,903 T→A 14.7% Y44F (TAT→TTT)  ECB_02005 ← putative phage protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,118,9030T→A14.7% 3.6 ‑3.1 ‑0.036131ECB_02005putative phage protein
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (17/13):  total (17/19)
Fisher's exact test strand distribution p-value = 4.55e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.91e-01

TCAGGCCGGAAACGAGCTGGCTTAATTTCGGGTAATAACCAA.CCACCGCCAGCCATTCCTGACCGGCG  >  REL606/2118868‑2118935
                                   |                                 
tCAGGCCGGAAACGAGCTGGCTTAATTTCGGGtaat                                   <  1:5802526/36‑1 (MQ=255)
 cAGGCCGGAAACGAGCTGGCTTAATTTCGGGtaata                                  >  1:6602519/1‑36 (MQ=255)
 cAGGCCGGAAACGAGCTGGCTTAATTTCGGGtaata                                  <  1:6007017/36‑1 (MQ=255)
  aGGCCGGAAACGAGCTGGCTTAATTTCGGGtaataa                                 >  1:5783360/1‑36 (MQ=255)
   ggCCGGAAACGAGCTGGCTTAATTTCGGGTAATaac                                <  1:3278135/36‑1 (MQ=255)
      cGGAAACGAGCTGGCTTAATTTCGGGTAATaaccaa                             <  1:8388767/36‑1 (MQ=255)
      cGGAAACGAGCTGGCTTAATTTCGGGTAATaaccaa                             >  1:6114015/1‑36 (MQ=255)
       ggAAACGAGCTGGCTTAATTTCGGGTAATaaccaa.c                           <  1:8017832/36‑1 (MQ=255)
       ggAAACGAGCTGGCTTAATTTCGGGTAATaaccaa.c                           <  1:7676338/36‑1 (MQ=255)
        gAAACGAGCTGGCTTAATTTCGGGTAATaaccaa.cc                          >  1:8348067/1‑36 (MQ=255)
        gAAACGAGCTGGCTTAATTTCGGGTAATaaccaa.cc                          >  1:617148/1‑36 (MQ=255)
        gAAACGAGCTGGCTTAATTTCGGGTAATaaccaa.cc                          <  1:4957407/36‑1 (MQ=255)
          aaCGATCTGGCTTAATTTCGGGTAATAACCAa.ccac                        <  1:204257/36‑1 (MQ=255)
          aaCGAGCTGGCTTAATTTCGGGTAATAACCAa.ccac                        >  1:4510972/1‑36 (MQ=255)
          aaCGAGCTGGCTTAATTTCGGGTAATAACCAa.ccac                        <  1:4682837/36‑1 (MQ=255)
          aaCGAGCTGGCTTAATTTCGGGTAATAACCAa.ccac                        <  1:5118792/36‑1 (MQ=255)
          aaCGAGCTGGCTTAATTTCGGGTAATAACCAa.ccac                        <  1:8431602/36‑1 (MQ=255)
           aCGAGCTGGCTTAATTTCGGGTAATAACCAa.ccacc                       <  1:941954/36‑1 (MQ=255)
            cGAGCTGGCTTAATTTCGGGTAATAACCAA.CCACCg                      >  1:3956800/1‑36 (MQ=255)
                cTGGCTTAATTTCGGGTAATAACCAA.CCACCgccag                  <  1:7474414/36‑1 (MQ=255)
                   gCTTAATTTCGGGTAATAACCAA.CCACCgccagcca               >  1:5123678/1‑36 (MQ=255)
                   gCTTAATTTCGGGTAATAACAAA.CCACCgccagcca               <  1:3355313/36‑1 (MQ=255)
                    cTTAATTTCGGGTAATAACCAA.CCACCGCCAGCCAt              >  1:3690909/1‑36 (MQ=255)
                     ttAATTTCGGTTAATAACAAA.CCCCCGCCAGCCAtt             <  1:2271648/36‑1 (MQ=25)
                     ttAATTTCGGGTAATAACCAA.CCACCGCCAGCCAtt             >  1:7705638/1‑36 (MQ=255)
                     ttAATTCCGGGTAAAAACCAA.CCCCCGCCAGCCAtt             <  1:1585783/36‑1 (MQ=25)
                      ttaTTTCCGGAAAAAAACAACCCACCGCCAGCCAtt             <  1:2959336/34‑1 (MQ=2)
                       aaTTTCGGGTAATAACCAA.CCACCGCCAGCCATTcc           >  1:91224/1‑36 (MQ=255)
                       aaTTTCGGGTAATAACCAA.CCACCGCCAGCCAGTcc           >  1:337637/1‑36 (MQ=255)
                       aaTTCCGGCAAAAAACCAA.CCACCGCCAGCCATTcc           <  1:1278371/36‑1 (MQ=16)
                         tttCGGGAAAAAACCAA.CCCCCGCCAGCCATTCCTg         <  1:4381132/36‑1 (MQ=25)
                         ttccGGGTAACACCCAG.CCCCCGCCAGCCATTCCTg         <  1:6541368/33‑1 (MQ=16)
                          ttCGGGTAATAACCAA.CCACCGCCAGCCCTTCCTGa        >  1:5406676/1‑36 (MQ=255)
                          ttCGGGTAATAACCAA.CCACCGCCAGCCATTCCTGa        >  1:1548259/1‑36 (MQ=255)
                           tCGGGTAAAAACCAA.CCACCGCCAGCCATTCCTGAc       <  1:6414347/36‑1 (MQ=255)
                             gggTAATAACCAA.CCACCGCCAGCCATTCCTGACCg     >  1:2292304/1‑36 (MQ=255)
                             gggTAATAACCAA.CCACCGCCAGCCATTCCTGACCg     >  1:2239436/1‑36 (MQ=255)
                             gggTAATAACCAA.CCACCGCCAGCCATTCCTGACCg     >  1:3451678/1‑36 (MQ=255)
                              ggTAATAACCAA.CCACCGCCAGCCATTCCTGACCgg    >  1:2705230/1‑36 (MQ=255)
                               gTAATAACCAA.CCACCGCCAGCCATTCCTGACCGGc   >  1:1599563/1‑36 (MQ=255)
                                taataaCCAA.CCACCGCCAGCCATTCCTGACCGGCg  >  1:7864905/1‑36 (MQ=255)
                                   |                                 
TCAGGCCGGAAACGAGCTGGCTTAATTTCGGGTAATAACCAA.CCACCGCCAGCCATTCCTGACCGGCG  >  REL606/2118868‑2118935

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG