Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 887,247 T→A 15.7% R452R (CGT→CGA ECB_00825 → putative replication protein for prophage

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606887,2470T→A15.7% 4.4 ‑4.3 0.0401356ECB_00825putative replication protein for prophage
Reads supporting (aligned to +/- strand):  new base (0/7):  ref base (17/16):  total (17/23)
Fisher's exact test strand distribution p-value = 1.42e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GTACTGTGCACTGGCATCTCATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGC  >  REL606/887213‑887281
                                  |                                  
gTACTGTGCACTTGCATCTCATGTGTTTCATGCGTa                                   <  1:7298357/36‑1 (MQ=255)
gTACTGTGCACTGGCATCTCATGTGTTTCATGCGTa                                   <  1:7972315/36‑1 (MQ=255)
 tACTGTGCACTGGCATCTCATGTGTTTCATGCGTaa                                  <  1:8455820/36‑1 (MQ=255)
  aCTGTGCACTGGCATCTCATGTGTTTCATGCGTaaa                                 <  1:1468121/36‑1 (MQ=255)
     gtgCACTGGCATCTCATGTGTTTCATGCGTaaaaaa                              <  1:1540761/36‑1 (MQ=255)
     gtgCACTGGCATCTCATGTGTTTCATGCGTaaaaaa                              <  1:5826070/36‑1 (MQ=255)
       gCACTGGCATCTCATGTGTTTCATGCGTAAAAAAGa                            <  1:8208275/36‑1 (MQ=255)
       gCACTGGCATCTCATGTGTTTCATGCGTAAAAAAGa                            <  1:2741457/36‑1 (MQ=255)
         aCTGGCATCTCATGTGTTTCATGCGTAAAAAAGAcc                          >  1:6563766/1‑36 (MQ=255)
          cTGGCATCTCATGTGTTTCATGCGTAAAAAAGAccg                         >  1:6109237/1‑36 (MQ=255)
          cTGGCATCTCATGTGTTTCATGCGTAAAAAAGAccg                         <  1:1738438/36‑1 (MQ=255)
            ggCATCTCATGTGTTTCATGCGTAAAAAAGAccgcc                       <  1:6839608/36‑1 (MQ=255)
            ggCATCTCATGTGTTTCATGCGTAAAAAAGAccgcc                       <  1:4139686/36‑1 (MQ=255)
            ggCATCTCATGTGTTTCATGCGTAAAAAAGAccgcc                       >  1:7681310/1‑36 (MQ=255)
             gCATCTCATGTGTTTCATGCGTAAAAAAGAccgccg                      >  1:5901273/1‑36 (MQ=255)
                   cATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCAt                >  1:1993059/1‑36 (MQ=255)
                   cATGTGTTTCATGCAAAAAAAAGCCCGCCGCGCCAt                <  1:7247085/36‑1 (MQ=25)
                   cATGTGTTTAATGCGAAAAAAAGACCGCCGCGCCAt                <  1:3460517/36‑1 (MQ=37)
                    aTGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATc               >  1:3042447/1‑36 (MQ=255)
                    aTGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATc               >  1:3983177/1‑36 (MQ=255)
                    aTGTGTTTCATCCATAAAAAAGGCCGCCGCGCCATc               <  1:1194819/36‑1 (MQ=25)
                    aTGTGTTTCAGGCGTAAAAAAGACCGCCGCGCCATc               <  1:4627074/36‑1 (MQ=255)
                     tgtgTTTCGTGAAAAAAAAACACCGCCGCGCCATCa              <  1:5238731/36‑1 (MQ=16)
                     tgtgTTTCATGCGTAAAAAAGACCGCCGCGCCATCa              >  1:4531047/1‑36 (MQ=255)
                     tgtgTTTCATGCGAAAAAAAGCCCGCCGCGCCATCa              <  1:2325084/36‑1 (MQ=37)
                     agtgtTTCATGAATAAAAAACCCCGCCGCGCCATCa              <  1:2404539/35‑1 (MQ=16)
                        gtTTCATGCGTAAAAAAGACCGCCGCGCCATCACTg           >  1:3585033/1‑36 (MQ=255)
                         tttCATGCGTAAAAAGGACCGCCGCGCCATCACTGc          <  1:6517860/36‑1 (MQ=255)
                         tttAATGCGAAAAAAAGACCGCCGCGCCATCACTGc          <  1:4145637/36‑1 (MQ=37)
                           tCATGCAAAAAAAAGACCGCCGCGCCATCACTGCAt        <  1:1932984/36‑1 (MQ=37)
                             aTGCGTAAAAAAGACCGCCGCGCCATCACTGCATTa      >  1:4822834/1‑36 (MQ=255)
                             aTGCGTAAAAAAGACCGCCGCGCCATCACTGCATTa      >  1:3372260/1‑36 (MQ=255)
                             aTGCGTAAAAAAGACCGCCGCGCCATCACTGCATTa      >  1:2869401/1‑36 (MQ=255)
                             aTGCGTAAAAAAGACCGCCGCGCCATCACTGCATTa      >  1:6739951/1‑36 (MQ=255)
                              tGCGTAAAAAAGACCGCCGCGCCATCACTGCATTAc     <  1:7934319/36‑1 (MQ=255)
                              tGCGTAAAAAAGACCGCCGCGCCATCACTGCATTAc     >  1:4699447/1‑36 (MQ=255)
                               gCGTAAAAAAGACCGCCGCGCCATCACTGCATTACt    >  1:2729215/1‑36 (MQ=255)
                               gCGAAAAAAAGCCCGCCGCGCCATCACTGCATTACt    <  1:990580/36‑1 (MQ=37)
                                cGTAAAAAAGACCGCCGCGCCATCACTGCATTACTg   >  1:643388/1‑36 (MQ=255)
                                 gTAAAAAAGACCGCCGCGCCATCACTGCATTACTGc  >  1:7179167/1‑36 (MQ=255)
                                  |                                  
GTACTGTGCACTGGCATCTCATGTGTTTCATGCGTAAAAAAGACCGCCGCGCCATCACTGCATTACTGC  >  REL606/887213‑887281

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 20 ≤  ATCG  < 31 ≤  ATCG  < 37 ≤  ATCG