Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,273,680 Δ1 bp 11.4% coding (17/552 nt) yrbD ← predicted ABC‑type organic solvent transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,273,6800T→.11.4% 13.5 ‑0.4 ‑0.444coding (17/552 nt)yrbDpredicted ABC‑type organic solvent transporter
Reads supporting (aligned to +/- strand):  new base (1/4):  ref base (14/25):  total (15/29)
Fisher's exact test strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.20e-01

TGCTGCTAATAAAAAGATACCCACCCAAATTTCATTTTTTTTCGTTTGCATGAACTCAATTCCCAAACA  >  REL606/3273646‑3273714
                                  |                                  
tgctgcTAATAAAAAGATACCCACCCAAATTTCAtt                                   <  1:7817936/36‑1 (MQ=255)
tgctgcTAATAAAAAGATACCCACCCAAATTTCAtt                                   >  1:6679242/1‑36 (MQ=255)
tgctgcTAATAAAAAGATACCCACCCAAATTTCAtt                                   >  1:2834918/1‑36 (MQ=255)
 gctgctAATAAAAATATACCCACCCAAATTTCAttt                                  >  1:1867114/1‑36 (MQ=255)
  ctgctAATAAAAAGATACCCACCCAAATTTCAtttt                                 >  1:5485177/1‑36 (MQ=255)
  ctgctAATAAAAAGATACCCACCCAAATTTCAtttt                                 >  1:4034919/1‑36 (MQ=255)
    gctAATAAAAAGATACCCACCCAAATTTCAtttttt                               >  1:6761995/1‑36 (MQ=255)
     ctAATAAAAAGATCCCCACCCAAATTTCAttttttt                              >  1:3190077/1‑36 (MQ=255)
      taataaAAAGATACCCACCCAAATTTCAtttttttt                             <  1:5952253/36‑1 (MQ=255)
       aataaAAAGATACCCACCCAAATTTCATTTTTTTTc                            <  1:1678291/36‑1 (MQ=255)
         taaAAAGATACCCACCCAAATTTCATTTTTTTTCGt                          >  1:6872342/1‑36 (MQ=255)
         taaAAAGATACCCACCCAAATTTCATTTTTTTTCGt                          >  1:4690053/1‑36 (MQ=255)
          aaaaaGATACCCACCCAAATTTCATTTTTTTTCGtt                         <  1:3369137/36‑1 (MQ=255)
          aaaaaGATACCCACCCAAATTTCATTTTTTTTCGtt                         >  1:4760402/1‑36 (MQ=255)
           aaaaGATACCCACCCAAATTTCATTTTTTTTCGttt                        >  1:2956596/1‑36 (MQ=255)
           aaaaGATACCCACCCAAATTTCATTTTTTTTCGttt                        >  1:5703096/1‑36 (MQ=255)
            aaaGATACCCACCCAAATTTCATTTTTTTTCGTTTg                       >  1:3545851/1‑36 (MQ=255)
            aaaGATACCCACCCAAATTTCATTTTTTTTCGTTTg                       <  1:389812/36‑1 (MQ=255)
            aaaGATACCCACCCAAATTTCA.TTTTTTTCGTTTGc                      >  1:5535918/1‑36 (MQ=37)
               gATACCCACCCAAATTTCATTTTTTTTCGTTTGCAt                    <  1:2268272/36‑1 (MQ=255)
               gATACCCACCCAAATTTCATTTTTTTTCGTTTGCAt                    <  1:3465224/36‑1 (MQ=255)
               gATACCCACCCAAATTTCATTTTTTTTCGTTTGCAt                    <  1:923356/36‑1 (MQ=255)
                attcccacccaAATTTCATTTTTTTTCGTTTGCATg                   >  1:4762090/4‑36 (MQ=255)
                aTACCCACCCAAATATCATTTTTTTTCGTTTGCAtt                   <  1:5840151/36‑2 (MQ=255)
                  gcccacccaAATCTCATTTTTTTTCTTTTGCATGaa                 >  1:8049493/2‑36 (MQ=37)
                   cccacccaAATTTCATTTTTTTTCGTTTGCATGAAc                <  1:5554772/36‑1 (MQ=255)
                    ccacccaAATTTCATTTTTTTTCGTTTGCATGAACt               >  1:216656/1‑36 (MQ=255)
                    ccacccaAATTTCATTTTTTTTCGTTTGCATGAACt               <  1:6587057/36‑1 (MQ=255)
                    ccacccaAATTTCA.TTTTTTTCGTTTGCATGAACTc              <  1:702419/36‑1 (MQ=37)
                    ccacccaAATTTCA.TTTTTTTCGTTTGCATGAACTc              <  1:697749/36‑1 (MQ=37)
                     cacccaAATTTCATTTTTTTTCGTTTGCATGAACTc              <  1:3082569/36‑1 (MQ=255)
                     cacccaAATTTCATTTTTTTTCGTTTGCATGAACTc              <  1:857176/36‑1 (MQ=255)
                     cacccaAATTTCATTTTTTTTCGTTTGCATGAACTc              <  1:7263281/36‑1 (MQ=255)
                      acccaAATTTCATTTTTTGTCGTGTGCTTGAAATCa             >  1:6813319/1‑36 (MQ=16)
                       cccaAATTTCATTTTTTTTCGTTTGCATGAACTCaa            >  1:2155203/1‑36 (MQ=255)
                         caAATTTCATTTTTTTTCGTTTGCATGAACTCAAtt          <  1:659109/36‑1 (MQ=255)
                         caAATTTCATTTTTTTTCGTTTGCATGAACTCAAtt          >  1:4637762/1‑36 (MQ=255)
                          aaaTTTCATTTTTTTTCGTTTTCATGAACTCAATTc         >  1:7570526/1‑36 (MQ=255)
                          aaaTTTCATTTTTTTTCGTTTGCATGAACTCAATTc         <  1:5452424/36‑1 (MQ=255)
                          aaaTTTCA.TTTTTTTCGTTTGCATGAACTCAATTcc        <  1:7512251/36‑1 (MQ=38)
                          aaaTTTCA.TTTTTTTCGTTTGCATGAACTCAATTcc        <  1:8254286/36‑1 (MQ=38)
                           aaTTTCATTTTTTTTCGTTTGCATGAACTCAATTcc        <  1:3790371/36‑1 (MQ=255)
                            aTTTCATTTTTTTTCTTTTGCATGAACTCAATTccc       >  1:7350727/1‑36 (MQ=255)
                            aTTTCATTTTTTTTCGTTTGCATGAACTCAATTccc       <  1:8091714/36‑1 (MQ=255)
                             tttCATTTTTTTTCGTTTGCATGAACTCAATTCCCa      <  1:1705184/36‑1 (MQ=255)
                             tttCATTTTTTTTCGTTTGCATGAACTCAATTCCCa      <  1:7884694/36‑1 (MQ=255)
                             tttCATTTTTTTTCGTTTGCATAAACTCAATTCCCa      <  1:4777128/36‑1 (MQ=255)
                              ttCATTTTTTTTCGTTTGCATGAACTCAATTCCCaa     <  1:2738135/36‑1 (MQ=255)
                                cATTTTTTTTCTTTTGCATGAACTCAATTCCCAAAc   <  1:2970113/36‑1 (MQ=255)
                                cATTTTTTTTCGTTTGCATGAACTCAATTCCCAAAc   <  1:5459542/36‑1 (MQ=255)
                                cATTTTTTTTCCTTTCCATGAGCTCAATTCCCAAAc   <  1:3624726/36‑1 (MQ=21)
                                 aTTTTTTTTCGTTTGCATGAACTCAATTCCCAAACa  <  1:7963755/36‑1 (MQ=255)
                                 aTTTTTTTTCGTTTGCATGAACTCAATTCCCAAACa  <  1:7386085/36‑1 (MQ=255)
                                  |                                  
TGCTGCTAATAAAAAGATACCCACCCAAATTTCATTTTTTTTCGTTTGCATGAACTCAATTCCCAAACA  >  REL606/3273646‑3273714

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG