Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 518,720 G→T 12.9% intergenic (‑191/‑126) allD ← / → fdrA ureidoglycolate dehydrogenase/membrane protein FdrA

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606518,7200G→T12.9% 3.8 ‑3.4 ‑0.031intergenic (‑191/‑126)allD/fdrAureidoglycolate dehydrogenase/membrane protein FdrA
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (10/17):  total (14/17)
Fisher's exact test strand distribution p-value = 3.18e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.62e-01

TTAAAAAATAACACGTAATAGCATGGCGGTTATTTGAATTCTTTTGTTAATAATCCCTGTGTGATATTCAT  >  REL606/518685‑518755
                                   |                                   
ttaAAAAATAACACGTAATAGCATGGCGGTTATTTg                                     <  1:6509004/36‑1 (MQ=255)
 taAAAAATAACACGTAATAGCATGGCGGTTATTTGa                                    <  1:6494888/36‑1 (MQ=255)
    aaaaTAACACGTAATAGCATGGCGGTTATTTGAAtt                                 <  1:4279747/36‑1 (MQ=255)
    aaaaTAACACGTAATAGCATGGCGGTTATTTGAAtt                                 <  1:6129296/36‑1 (MQ=255)
      aaTAACACGTAATAGCATGGCGGTTATTTGAATTCt                               <  1:5887832/36‑1 (MQ=255)
        tAACACGTAATAGCATGGCGGTTATTTGAATTCttt                             <  1:4036641/36‑1 (MQ=255)
        tAACACGTAATAGCATGGCGGTTATTTGAATTCttt                             <  1:3155466/36‑1 (MQ=255)
          acacGTAATAGCATGGCGGGTTTTTTTATTTttttt                           >  1:3561197/1‑35 (MQ=12)
           cacGTAATAGCATGGCGGGTGTTTTTATTTtttttt                          >  1:4918897/1‑34 (MQ=12)
            acGTAATAGCATGGCGGTTATTTGAATTCTTTTGtt                         <  1:205326/36‑1 (MQ=255)
              gTAATAGCATGGCGGTTTTTTTTATTTTTTTTtttt                       >  1:1558833/1‑34 (MQ=9)
               tAATAGCATGGCGGTTATTTGAATTCTTTTGTtaat                      <  1:5564436/36‑1 (MQ=255)
                    gCATGGCGGTTATTTTTAGTATTTTTTTtataatta                 >  1:4874622/1‑34 (MQ=2)
                    gCATGGCGGTTATTTGAATTCTTTTGTTAATAATcc                 <  1:1853005/36‑1 (MQ=255)
                    gCATGGCGGTTATTTGAATTCTTTTGTTAATAATcc                 <  1:502485/36‑1 (MQ=255)
                         gCGGTTATTTGAATTCTTTTTTTAATAATCCCtgtg            >  1:1714730/1‑36 (MQ=255)
                         gCGGTTATTTGAATTCTTTTGTTAATAATCCCtgtg            <  1:6880788/36‑1 (MQ=255)
                           ggTTATTTGAATTCTTTTGTTAATAATCCCtgtgtg          <  1:6019557/36‑1 (MQ=255)
                           ggTTATTTGAATTCTTTTGTTAATAATCCCtgtgtg          >  1:1631288/1‑36 (MQ=255)
                            gTTATTTGAATTCTTTTGTTAATAATCCCTGTGTGa         >  1:3399173/1‑36 (MQ=255)
                            gTTATTTGAATTCTTTTGTTAATAATCCCTGTGTGa         <  1:110884/36‑1 (MQ=255)
                              tATTTGAATTCTTTTTTTAATAATCCCTTTTTGata       >  1:519509/1‑36 (MQ=25)
                              tATTTGAATTCTTTTGTTAATAATCCCTGTGTGata       >  1:2519651/1‑36 (MQ=255)
                               aTTTGAATTCTTTTGTTAATAAGCCCTGTGTGatat      <  1:4283188/36‑1 (MQ=255)
                                tttGAATTCTTTTGTTAATAATCCCTGTGTGATAtt     <  1:3722065/36‑1 (MQ=255)
                                tttGAATTCTTTTGTTAATAATCCCTGTGTGATAtt     >  1:2953548/1‑36 (MQ=255)
                                tttGAATTCTTTTGTTAATAATCCCTGTGTGATAtt     <  1:632457/36‑1 (MQ=255)
                                tttGAATTCTTTTGTTAATAATCCCTGTGTGATAtt     >  1:1516381/1‑36 (MQ=255)
                                 ttGAATTCTTTTGTTAATAATCCCTGTGTGATATtc    <  1:2931152/36‑1 (MQ=255)
                                 ttGAATTCTTTTGTTAATAATCCCTGTGTGATATtc    >  1:1792793/1‑36 (MQ=255)
                                  tGAATTCTTTTGTTAATAATCCCTGTGTTATAttaa   >  1:385109/1‑34 (MQ=38)
                                  tGAATTCTTTTGTTAATAATCCCTGTGTGATATtca   >  1:3505367/1‑36 (MQ=255)
                                   gAATTCTTTTGTTAATAATCCCTGTGTGATATtcat  <  1:2121610/36‑1 (MQ=255)
                                   |                                   
TTAAAAAATAACACGTAATAGCATGGCGGTTATTTGAATTCTTTTGTTAATAATCCCTGTGTGATATTCAT  >  REL606/518685‑518755

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 26 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG