Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 146,093 Δ1 bp 4.8% coding (512/927 nt) yadG → predicted transporter subunit: ATP‑binding component of ABC superfamily

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606146,0930T→. 4.8% 5.1 ‑1.4 ‑0.463coding (512/927 nt)yadGpredicted transporter subunit: ATP‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (30/30):  total (33/30)
Fisher's exact test strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.01e-01

TCTCGACGAACCGACCGCAGGCGTGGATATTGAACTTCGCCGCTCAATGTGGGGCTTTTTGAAGGATTTAA  >  REL606/146058‑146128
                                   |                                   
tctcGACGAACCGACCGCAGGCGTGGATATTGAACt                                     <  1:8423234/36‑1 (MQ=255)
 ctcGACGAACCGACCGCAGGCGTGGATATTGAACtt                                    >  1:3531297/1‑36 (MQ=255)
  tcGACGAACCGACCGCAGGCGTGGATATTGAACTTc                                   <  1:5880427/36‑1 (MQ=255)
  tcGACGAACCGACCGCAGGCGTGGATATTGAACTTc                                   <  1:5900309/36‑1 (MQ=255)
  tcGACGAACCGACCGCAGGCGTGGATATTGAACTTc                                   >  1:3393771/1‑36 (MQ=255)
  tcGACGAACCGACCGCAGGCGTGGATATTGAACTTc                                   <  1:2887609/36‑1 (MQ=255)
   cgacgaACCGACCGCAGGCGTGGATATTGAACTTcg                                  >  1:486334/1‑36 (MQ=255)
    gacgaACCGACCGCAGGCGTGGATATTGAACTTcgc                                 <  1:4956126/36‑1 (MQ=255)
    gacgaACCGACCGCAGGCGTGGATATTGAACTTcgc                                 >  1:4742182/1‑36 (MQ=255)
    gacgaACCGACCGCAGGCGTGGATATTGAACTTcgc                                 <  1:8349051/36‑1 (MQ=255)
      cgaACCGACCGCAGGCGTGGATATTGAACTTcgccg                               <  1:7221453/36‑1 (MQ=255)
      cgaACCGACCGCAGGCGTGGATATTGAACTTcgccg                               <  1:141562/36‑1 (MQ=255)
       gaACCGACCGCAGGCGTGGATATTGAACTTcgccgc                              >  1:528949/1‑36 (MQ=255)
       gaACCGACCGCAGGCGTGGATATTGAACTTcgccgc                              >  1:2469460/1‑36 (MQ=255)
       gaACCGACCGCAGGCGTGGATATTGAACTTcgccgc                              >  1:7599280/1‑36 (MQ=255)
        aaCCGACCGCAGGCGTGGATATTGAACTTCGCCGCt                             <  1:5017173/36‑1 (MQ=255)
        aaCCGACCGCAGGCGTGGATATTGAACTTCGCCGCt                             <  1:8409430/36‑1 (MQ=255)
         accgaccgCAGGCGTGGATATTGAACTTCGCCGCTc                            <  1:7191828/36‑1 (MQ=255)
         accgaccgCAGGCGTGGATATTGAACTTCGCCGCTc                            >  1:3376391/1‑36 (MQ=255)
         accgaccgCAGGCGTGGATATTGAACTTCGCCGCTc                            <  1:1705230/36‑1 (MQ=255)
          ccgaccgCAGGCGTGGATATTGAACTTCGCCGCTCa                           <  1:2388729/36‑1 (MQ=255)
            gaccgCAGGCGTGGATATTGAACTTCGCCGCTCAAt                         <  1:3983154/36‑1 (MQ=255)
            gaccgCAGGCGTGGATATTGAACTTCGCCGCTCAAt                         >  1:5021688/1‑36 (MQ=255)
             accgCAGGCGTGGATATTGAACTTCGCCGCTCAAtg                        >  1:6207745/1‑36 (MQ=255)
             accgCAGGCGTGGATATTGAACTTCGCCGCTCAAtg                        <  1:7719809/36‑1 (MQ=255)
             accgCAGGCGTGGATATTGAACTTCGCCGCGCAAtg                        >  1:7743296/1‑36 (MQ=255)
              ccgCAGGCGTGGATATTGAACTTCGCCGCTCAAtgt                       <  1:7987372/36‑1 (MQ=255)
              ccgCAGGCGTGGATATTGAACTTCGCCGCTCAAtgg                       >  1:1906822/1‑35 (MQ=255)
               cgCAGGCGTGGATATTGAACTTCGCCGCTCAAtgtg                      <  1:5082452/36‑1 (MQ=255)
                gCAGGCGTGGATCTTGAACTTCGCTGCTCAATGTgg                     >  1:7764408/1‑36 (MQ=38)
                gCAGGCGTGGATATTGAACTTCGCCGCTCCATGTgg                     >  1:4827822/1‑36 (MQ=255)
                gCAGGCGTGGATATTGAACTTCGCCGCTCAATGTgg                     >  1:4483034/1‑36 (MQ=255)
                gCAGGCGTGGATATTGAACTTCGCCGCTCAATGTgg                     <  1:5418932/36‑1 (MQ=255)
                gCAGGCGTGGATATTGAACTTCGCCGCTCAATGCgg                     >  1:6926809/1‑36 (MQ=255)
                 cAGGCGTGGATATTGACCTTCGCCGCTCAATGTggg                    <  1:1056386/36‑1 (MQ=255)
                  aGGCGTGGATATTGAACTTCGCCGCTCAATGTgggg                   >  1:7744031/1‑36 (MQ=255)
                     cGTGGATATTGAACTTCGCCGCTCAATGTGGGGCtt                >  1:4129274/1‑36 (MQ=255)
                      gTGGATATTGAACTTCGCCGCTCAATGTGGGGCttt               >  1:7842594/1‑36 (MQ=255)
                       tGGATATTGAACTTCGCCGCTCCATGTGGGGCtttt              >  1:8201538/1‑36 (MQ=255)
                       tGGATATTGAACTTCGCCGCTCAATGTGGGGCtttt              <  1:8406934/36‑1 (MQ=255)
                       tGGATATTGAACTTCGCCGCTCAATGTGGGGCtttt              <  1:2860202/36‑1 (MQ=255)
                       tGGATATTGAACTTCGCCGCTCAATGTGGGGCtttt              >  1:6516722/1‑36 (MQ=255)
                        ggATATTGAACTTCGCCGCTCAATGTGGGGCttttt             <  1:7729688/36‑1 (MQ=255)
                        ggATATTGAACTTCGCCGCTCAATGTGGGGCttttt             <  1:3544713/36‑1 (MQ=255)
                         gATATTGAACTTCGCCGCTCAATGTGGGGCTTTTTg            >  1:598257/1‑36 (MQ=255)
                          atatTGAACTTCGCCGCTCAATGTGGGGCTTTTTGa           >  1:1160025/1‑36 (MQ=255)
                           tatTGAACTTCGCCGCTCAATGTGGGGCTTTTTGaa          >  1:7786064/1‑36 (MQ=255)
                           tatTGAACTTCGCCGCTCAATGTGGGGCTTTTTGaa          <  1:3256931/36‑1 (MQ=255)
                           tatTGAACTTCGCCGATCAATGTGGGGCTTTTTGaa          >  1:2412878/1‑36 (MQ=255)
                            cttGAAC.TCGCCGCTCAATGTGGGGCTTTTTGAAgg        >  1:6704612/2‑36 (MQ=38)
                            cttGAAC.TCGCCGCTCAATGTGGGGCTTTTTGAAgg        >  1:5167537/2‑36 (MQ=38)
                            cttGAAC.TCGCCGCTCAATGTGGGGCTTTTTGAAgg        >  1:5249888/2‑36 (MQ=38)
                            atTGAACTTCGCCGCTCAATGTGGGGCTTTTTGAAg         >  1:3203856/1‑36 (MQ=255)
                            atTGAACTTCGCCGCTCAATGTGGGGCTTTTTGAAg         <  1:572567/36‑1 (MQ=255)
                             ttGAACTTCGCCGCTCAATGTGGGGCTTTTTGAAgg        >  1:7869928/1‑36 (MQ=255)
                              tGAACTTCGCCGCTCAATGTGGGGCTTTTTGAAGGa       >  1:7883569/1‑36 (MQ=255)
                              tGAACTTCGCCGCTCAATGTGGGGCTTTTTGAAGGa       >  1:2882305/1‑36 (MQ=255)
                              tGAACTTCGCCGCTCAATGTGGGGCTTTTTGAAGGa       <  1:373484/36‑1 (MQ=255)
                              tGAACTTCGCCGCTCAATGTGGGGCTTTTTGAAGGa       >  1:5664748/1‑36 (MQ=255)
                                aaCTTCGCCGCTCAATGTGGGGCTTTTTGAAGGAtt     <  1:4586351/36‑1 (MQ=255)
                                 aCTTCGCCGCTCAATGTGGGGCTTTTTGAAGGAttt    <  1:8065228/36‑1 (MQ=255)
                                 aCTTCGCCGCTCAATGTGGGGCTTTTTGAAGGAttt    <  1:5176297/36‑1 (MQ=255)
                                  cTTCGCCGCTCAATGTGGGGCTTTTTGAAGGATTTa   <  1:2666190/36‑1 (MQ=255)
                                  cTTCGCCGCTCAATGTGGGGCTTTTTGAAGGATTTa   >  1:3181659/1‑36 (MQ=255)
                                  cTTCGCCGCTCAATGTGGGGCTTTTTGAAGGATTTa   <  1:140739/36‑1 (MQ=255)
                                   ttCGCCGCTCAATGTGGGGCTTTTTGAAGGATTTaa  >  1:4548700/1‑36 (MQ=255)
                                   |                                   
TCTCGACGAACCGACCGCAGGCGTGGATATTGAACTTCGCCGCTCAATGTGGGGCTTTTTGAAGGATTTAA  >  REL606/146058‑146128

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG