Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,585,439 C→A 29.6% L171F (TTG→TTT marC ← predicted transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,585,4390C→A29.6% 2.9 ‑4.5 0.019513marCpredicted transporter
Reads supporting (aligned to +/- strand):  new base (0/6):  ref base (9/4):  total (9/10)
Fisher's exact test strand distribution p-value = 1.08e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGG  >  REL606/1585404‑1585474
                                   |                                   
cGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCcac                                     <  1:2661305/36‑1 (MQ=255)
cGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCcac                                     <  1:4063749/36‑1 (MQ=255)
 gCATAATTGCGCCGGAGCTGCGTAGGCTTCCCcaca                                    <  1:7087222/36‑1 (MQ=255)
 gCATAATTGCGCCGGAGCTGCGTAGGCTTCCCcaca                                    >  1:3041212/1‑36 (MQ=255)
 gCATAATTGCGCCGGAGCTGCGTAGGCTTCCCcaca                                    <  1:5696769/36‑1 (MQ=255)
  cATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACaa                                   >  1:7149189/1‑36 (MQ=255)
   aTAATTGCGCCGGAGCTGCGTAGGCTTCCCCACaaa                                  <  1:638316/36‑1 (MQ=255)
        tGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAAc                             <  1:6083117/36‑1 (MQ=255)
        tGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAAc                             <  1:5008642/36‑1 (MQ=255)
         gcgcCGGAGCTGCGTAGGCTTCCCCACAAAATAAcc                            >  1:1785816/1‑36 (MQ=255)
         gcgcCGGAGCTGCGTAGGCTTCCCCACAAAATAAcc                            <  1:6403708/36‑1 (MQ=255)
            ccGGAGCTGCGTAGGCTTCCCAAAAAAATAACCGcc                         <  1:853744/36‑1 (MQ=37)
               gAGCTGCGTAGGCTTCCCCACAAAATAACCGccacc                      >  1:3121769/1‑36 (MQ=255)
                aGCTGCGTAGGCTTCCCCACAAAATAACCGccacca                     >  1:4501082/1‑36 (MQ=255)
                  cTGCGTAGGCTTACCCAAAAAATAACCGCCACCAgg                   <  1:7082960/36‑1 (MQ=37)
                   tGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGa                  >  1:4873478/1‑36 (MQ=255)
                   tGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGa                  >  1:653951/1‑36 (MQ=255)
                     cGTAGGCTTCCCAAAAAAAAAACCGCCACCAGGaaa                <  1:4563327/36‑1 (MQ=21)
                       tAGGCTTCCAAAAAAAAAAACCGCCACCAGGaaaaa              <  1:5712422/36‑1 (MQ=19)
                       gaGGCTTCCCCAAAAAATAACCGCCACCAGGaaaaa              <  1:1872565/35‑1 (MQ=255)
                         ggCTTCCCCACAAAATAACCGCCACCAGGAAAAAGa            >  1:2097025/1‑36 (MQ=255)
                           cTTCCCAAAAAAAAAACCGCCACCAGGAAAAAGATc          <  1:5852775/36‑1 (MQ=21)
                            ttCCCCACAAAATAACCGCCACCAGGAAAAAGATCa         >  1:3243232/1‑36 (MQ=255)
                            ttCCCCACAAAATAACCGCCACCAGGAAAAAGATCa         >  1:1648765/1‑36 (MQ=255)
                                   caAAATAACCGCCACCAGGAAAAAGATCAgcggcgg  >  1:5352530/1‑36 (MQ=255)
                                   |                                   
CGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGCGGCGG  >  REL606/1585404‑1585474

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 18 ≤  ATCG  < 28 ≤  ATCG  < 37 ≤  ATCG