Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 862,403 C→A 26.3% L128L (CTC→CTA yliB → predicted peptide transporter subunit: periplasmic‑binding component of ABC superfamily

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606862,4030C→A26.3% 6.1 ‑3.3 0.038384yliBpredicted peptide transporter subunit: periplasmic‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  new base (0/11):  ref base (10/17):  total (10/28)
Fisher's exact test strand distribution p-value = 3.72e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCT  >  REL606/862368‑862435
                                   |                                
gAAGCGATCGATCCGACAACGGTAAAGATTACCCTc                                  >  1:7478580/1‑36 (MQ=255)
gAAGCGATCGATCCGACAACGGTAAAGATTACCCTc                                  >  1:4219859/1‑36 (MQ=255)
 aaGCGATCGATCCGACAACGGTAAAGATTACCCTCa                                 <  1:1948076/36‑1 (MQ=255)
 aaGCGATCGATCCGACAACGGTAAAGATTACCCTCa                                 >  1:350284/1‑36 (MQ=255)
  aGCGATCGATCCGACAACGGTAAAGATTACCCTCaa                                <  1:3676516/36‑1 (MQ=255)
   gCGATCGATCCGACAACGGTAAAGATTACCCTCaaa                               <  1:2163573/36‑1 (MQ=255)
    cgatcgatCCTACAACGGTAAAGATTACCCTCAAAc                              <  1:4392200/36‑1 (MQ=255)
    cgatcgatCCGACAACGGTAAAGATTACCCTCAAAc                              >  1:5940271/1‑36 (MQ=255)
       tcgatcCGACAACGGTAAAGATTACCCTCAAACAGc                           <  1:357732/36‑1 (MQ=255)
       tcgatcCGACAACGGTAAAGATTACCCTCAAACAGc                           <  1:3554718/36‑1 (MQ=255)
        cgatcCGACAACTGTAAATATTACCCTCAAACAGcc                          <  1:7010116/36‑1 (MQ=37)
        cgatcCGACAACGGTAAAGATTACCCTCAAACAGcc                          <  1:7075634/36‑1 (MQ=255)
              gACAACGGTAAAGATTACCCTCAAACAGCCGTtctc                    <  1:7399986/36‑1 (MQ=255)
              gACAACGGTAAAGATTACCCTCAAACAGCAGTtctc                    >  1:6155577/1‑36 (MQ=255)
               aaaaCGGTAAAGATTACCCTCAAACAGCCGTTCTCa                   <  1:2033530/34‑1 (MQ=255)
               aaaaCGGAAAAGATTACCCTAAAACAGCCGTTCTCa                   <  1:3083759/34‑1 (MQ=37)
               aCAACGGTAAAGATTACCCTCAAACAGCCGTTCTCa                   <  1:3663252/36‑1 (MQ=255)
               aCAACGGAAAAGATTACCCTAAAACAGCCGTTCTCa                   <  1:5059799/36‑1 (MQ=38)
                cAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAg                  <  1:1971476/36‑1 (MQ=255)
                 aaCGGTAAAGATTACCCTCAAACAGCCGTTCTCAGc                 >  1:1902899/1‑36 (MQ=255)
                 aaCGGTAAAGATTACCCTAAACCAGCCGTTCTCAGc                 <  1:147509/36‑1 (MQ=38)
                 aaCGGTAAAGATTACCCTAAAACAGCCGTTCTCAGc                 <  1:2382608/36‑1 (MQ=255)
                   cGGAAAAGATCACCCTCAAACAGCCGTTCTCAGCGt               <  1:6990065/36‑1 (MQ=37)
                    ggaaaaGATTCCCCTAAACCAGCCGTTCTCAGCGtt              <  1:1183738/33‑1 (MQ=19)
                    ggaaaaGAATCCCCGCAAACAGCCGTTCTCAGCGtt              <  1:617786/33‑1 (MQ=19)
                     gTAAAGATTACCCTCAAACAGCCGTTCTCAGCGttt             >  1:6147214/1‑36 (MQ=255)
                     gTAAAGATTACCCTAAAACAGCCGTTCTCAGCGttt             <  1:5325851/36‑1 (MQ=255)
                      tAAAGATTACCCTCAAACAGCCGTTCTCAGCGTtta            <  1:1291880/36‑1 (MQ=255)
                      aagaGATGACCCTCAACCAGCCGTTCTCAGCGTtta            <  1:1603481/33‑1 (MQ=25)
                       aaaGATTACCCTCAAACAGCCGTTCTCAGCGctttt           >  1:6445337/1‑34 (MQ=38)
                       aaaGATTACCCTCAAACAGCCGTTCTCAGCGTttat           >  1:5181269/1‑36 (MQ=255)
                       aaaGATTACCCTCAAACAGCCGTTCTCAGCGTttat           <  1:701305/36‑1 (MQ=255)
                         aGATTACTCTAAAACAGCCGTTCTCAGCGTttatta         <  1:2764919/36‑1 (MQ=38)
                         aGATTACCCTCAAACAGCCGTTCTCAGCGTttatta         >  1:4291911/1‑36 (MQ=255)
                          gATTACCCTAAAACAGCCGTTCTCAGCGTTTATTaa        <  1:175859/36‑1 (MQ=255)
                            ttACCCTAAAACAGCCGTTCTCAGCGTTTATTAata      <  1:2816685/36‑1 (MQ=255)
                              aCCCTAAAACAGCCGTTCTCAGCGTTTATTAATAtt    <  1:3435449/36‑1 (MQ=255)
                                ccTCAAACAGCCGTTCTCAGCGTTTATTAATATTCt  >  1:4345929/1‑36 (MQ=255)
                                ccTCAAACAGCCGTTCTCAGCGTTTATTAATATTCt  >  1:5777440/1‑36 (MQ=255)
                                ccTAAAACAGCCGTTCTCAGCGTTTATTAATATTCt  <  1:619553/36‑1 (MQ=255)
                                   |                                
GAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCT  >  REL606/862368‑862435

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 17 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG