Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,333,604 G→T 33.6% pseudogene (73/129 nt) yciX → predicted protein; b4522(b4574)

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,333,6040G→T33.6% 4.9 ‑4.0 0.022pseudogene (73/129 nt)yciXpredicted protein; b4522(b4574)
Reads supporting (aligned to +/- strand):  new base (8/0):  ref base (6/8):  total (14/8)
Fisher's exact test strand distribution p-value = 1.77e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TTAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGCACTGCTG  >  REL606/1333571‑1333639
                                 |                                   
ttAAAACGTACAGCCTGGCCGGCGGTTTTTCTTGGt                                   <  1:985513/36‑1 (MQ=37)
 tAAAACGTAAAGCCTGGCTGGCGGTTTTTTTTTGtt                                  >  1:590996/1‑36 (MQ=37)
 tAAAACGTAAAGCCTGGCTGGCGGGTTTTTTTTGtt                                  >  1:802377/1‑36 (MQ=21)
    aaCGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTg                               <  1:3535100/36‑1 (MQ=255)
     aCGTAAAGCCTTGCTGGCGGTTTTTCTTGGTTCTGc                              <  1:651358/36‑1 (MQ=255)
         aaaGCCTGGCTGGCGGTTTTTTTTTTTTCTGTACtt                          >  1:1762075/1‑36 (MQ=16)
          aaGCCTGGCTGGCGGTTTTTTTTTGTTTTTCACttt                         >  1:4605790/1‑36 (MQ=16)
            gCCTGGCTGGCGGTTTTTCTTGGTTCTGCACttttt                       >  1:1181244/1‑35 (MQ=255)
             ccTGGCTGGCGGTTTTTTTTTGTTCTTCACtttttt                      >  1:4027652/1‑34 (MQ=21)
              ctggctggCGGTTTTTTTTTTTTCTTCACTTTTCtt                     >  1:5364358/1‑35 (MQ=16)
              ctggctggCGGTTTTTCTTGGTTCTGCACTTTTCTg                     <  1:4090500/36‑1 (MQ=255)
              ctggctggCGGGTTTTTTTTGTTGTGCCCttttttt                     >  1:1612488/1‑33 (MQ=9)
               tggctggcGGTTTTTCTTGGTTCTGCACTTTTCTgg                    <  1:889600/36‑1 (MQ=255)
                  ctggcGGTTTTTTTTTGTTTTTCACTTTTCTTggtt                 >  1:5845521/1‑35 (MQ=9)
                     gcGGTTTTTCTTGGTTCTGCACTTTTCTGGGTGGtt              <  1:4535365/36‑1 (MQ=255)
                       ggTTTTTCTTGGTTCTGCACTTTTCTTGGTGGTTGt            >  1:892875/1‑36 (MQ=255)
                       ggTTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGt            >  1:3026124/1‑36 (MQ=255)
                        gTTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTc           >  1:5374164/1‑36 (MQ=255)
                          ttttCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGc         >  1:400940/1‑36 (MQ=255)
                          ttttCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGc         <  1:178216/36‑1 (MQ=255)
                             tCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGCAct      <  1:2503137/36‑1 (MQ=255)
                               ttGGTTCTGCACTTTTCTGGGTGTTTGTCGCActgc    >  1:2083825/1‑36 (MQ=255)
                                 ggTTCTGCACTTTTCTGGGTGGTTGTCGCActgctg  >  1:601300/1‑36 (MQ=255)
                                 |                                   
TTAAAACGTAAAGCCTGGCTGGCGGTTTTTCTTGGTTCTGCACTTTTCTGGGTGGTTGTCGCACTGCTG  >  REL606/1333571‑1333639

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 10 ≤  ATCG  < 26 ≤  ATCG  < 37 ≤  ATCG