Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,566,735 C→A 12.9% L123F (TTG→TTT yneL ← predicted transcriptional regulator

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,566,7350C→A12.9% 2.4 ‑3.0 ‑0.138369yneLpredicted transcriptional regulator
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (18/15):  total (18/20)
Fisher's exact test strand distribution p-value = 4.80e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.80e-01

GCCCCATCAAGTATTAATACCTCACCCGTACGGTCCAAAATTTGAATACCGACATTTTGGGCAGAGCCTGC  >  REL606/1566700‑1566770
                                   |                                   
gCCCCATCAAGTATTAATACCTCACCCGTACGGTcc                                     >  1:2248333/1‑36 (MQ=255)
gCCCCATCAAGTATTAATACCTCACCCGTACGGTcc                                     >  1:625479/1‑36 (MQ=255)
 ccccATCAAGTATTAATACCTCACCCGTACGGTCCa                                    <  1:488495/36‑1 (MQ=255)
  cccATCAAGTATTAATACCTCACCCGTACGGTCCaa                                   <  1:4371116/36‑1 (MQ=255)
  cccATCAAGTATTAATACCTCACCCGTACGGTCCaa                                   <  1:4445782/36‑1 (MQ=255)
    cATCAAGTATTAATACCTCACCCGTACGGTCCaaaa                                 <  1:6664541/36‑1 (MQ=255)
       cAAGTATTAATACCTCACCCGTACGGTCCAAAAttt                              <  1:1053015/36‑1 (MQ=255)
        aaGTATTAATACCTCACCCGTACGGTCCAAAATTTg                             >  1:6944170/1‑36 (MQ=255)
        aaGTATTAATACCTCACCCGTACGGGCCAAAATTTg                             >  1:2390079/1‑36 (MQ=255)
          gTATTAATACCTCACCCGCACGCTCCAAAATTTGaa                           >  1:6737962/1‑36 (MQ=37)
            aTTAATACCTCACCCGTACGGTCCAAAATTTGAATa                         <  1:5933772/36‑1 (MQ=255)
             ttAATACCTCACCCGTACGGTCCAAAATTTGAATAc                        <  1:2323769/36‑1 (MQ=255)
             ttAATACCTCACCCGTACGGTCCAAAATTTGAATAc                        >  1:5610562/1‑36 (MQ=255)
              tAATACCTCACCCGTACGGTCCAAAATTTGAATAcc                       >  1:3151202/1‑36 (MQ=255)
              tAATACCTCACCCGTACGGTCCAAAATTTGAATAcc                       <  1:7123749/36‑1 (MQ=255)
               aaTACCTCACCCGTACGGTCCAAAATTTGAATACCg                      <  1:4716247/36‑1 (MQ=255)
               aaTACCTCACCCGTACGGTCCAAAATTTGAATACCg                      <  1:957180/36‑1 (MQ=255)
                   ccTCACCCGTACGGTCCAAAATTTGAATACCGACAt                  >  1:421930/1‑36 (MQ=255)
                   ccTCACCCGTACGGTCCAAAATTTGAATACCGACAt                  >  1:4377997/1‑36 (MQ=255)
                   ccTCACCCGTACGGTCCAAAATTTGAATACCGACAt                  <  1:906318/36‑1 (MQ=255)
                    ctaaCCCGTACGGTAAAAAATTTGAATACCGACAtt                 <  1:4830800/33‑1 (MQ=25)
                     tCACCCGTACGGTCCAAAATTTGAATACCGACAttt                <  1:373321/36‑1 (MQ=255)
                     tCACCCGTACGGTCCAAAATTTGAATACCGACAttt                >  1:7762837/1‑36 (MQ=255)
                     tCACCCGTACGGTCCAAAATTTGAATACCGACAttt                <  1:2860708/36‑1 (MQ=255)
                     tCACCCGTACGGTCAAAAATTTGAATACCGACAttt                <  1:5401980/36‑1 (MQ=255)
                      cACCCGTACGGTCCAAAATTTGAATACCGACAtttt               >  1:2044558/1‑36 (MQ=255)
                      cACCCGTACGGTCCAAAATTTGAATACCGACAtttt               <  1:5184013/36‑1 (MQ=255)
                       aCCCGTACGGTCCAAAATTTGAATACCGACATTTTg              >  1:244522/1‑36 (MQ=255)
                        cccGTACGGTCCAAAATTTGAATACCGACATTTTgg             <  1:7760011/36‑1 (MQ=255)
                        cccGTACGGTCCAAAATTTGAATACCGACATTTTgg             >  1:3543513/1‑36 (MQ=255)
                        cccGTACGGTCAAAAATTTGAATACCGACATTTTgg             <  1:6515003/36‑1 (MQ=255)
                         ccGTACGGTCCAAAATTTGAATACCGACATTTTggg            >  1:3486287/1‑36 (MQ=255)
                         ccGTACGGTAAAAAATTTGAATACCGACATTTTggg            <  1:2980827/36‑1 (MQ=37)
                          cGTACGGTCCAAAATTTGAATACCGACATTTTGGGc           >  1:1019497/1‑36 (MQ=255)
                          cGTACGGTCCAAAATTTGAATACCGACATTTTGGGc           >  1:2038213/1‑36 (MQ=255)
                           gTACGGTCCAAAATTTGAATACCGACATTTTGGGCa          >  1:5516671/1‑36 (MQ=255)
                             aCGGTAAAAAATTTGAATACCGACATTTTGGGCaga        <  1:2510074/36‑1 (MQ=37)
                               ggTCCAAAATTTGAATACCGACATTTTGGGCAGAGc      >  1:2255321/1‑36 (MQ=255)
                                 tCCAAAATTTGAATACCGACATTTTGGGCAGAGCCt    >  1:2342145/1‑36 (MQ=255)
                                 tCCAAAATTTGAATACCGACATTTTGGGCAGAGCCt    >  1:563141/1‑36 (MQ=255)
                                   cAAAATTTGAATACCGACATTTTGGGCAGAGCCCGc  >  1:7766415/1‑36 (MQ=255)
                                   |                                   
GCCCCATCAAGTATTAATACCTCACCCGTACGGTCCAAAATTTGAATACCGACATTTTGGGCAGAGCCTGC  >  REL606/1566700‑1566770

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG