Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,366,935 G→T 24.2% A197A (GCC→GCA yjeP ← predicted mechanosensitive channel

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,366,9350G→T24.2% 4.4 ‑6.0 0.040591yjePpredicted mechanosensitive channel
Reads supporting (aligned to +/- strand):  new base (9/2):  ref base (7/22):  total (16/24)
Fisher's exact test strand distribution p-value = 2.53e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACG  >  REL606/4366900‑4366968
                                   |                                 
tctGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTg                                   <  1:2838539/36‑1 (MQ=255)
 gtGAGCGGAAGCGCGCTAATTCCTGGCGGTTATTgg                                  <  1:7778366/35‑1 (MQ=255)
  tGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGc                                 <  1:623835/36‑1 (MQ=255)
        tAAGCGCGCTAATTCCTGGCGGTTTTTTGCAGAccg                           >  1:2833135/1‑34 (MQ=37)
          aGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCt                         <  1:5017578/36‑1 (MQ=255)
          aGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCt                         >  1:1769183/1‑36 (MQ=255)
          aGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCt                         <  1:6710227/36‑1 (MQ=255)
           gcgcgcTAATTCCTGGCGGTTATTGGCAGACAGCTg                        <  1:5130665/36‑1 (MQ=255)
             gcgcTAATTCCTGGCGGTTATTTGCAGACAGCTgcg                      >  1:2975340/1‑36 (MQ=255)
             gcgcTAATTCCTGGCGGTTATTTGCAGACAGATgcg                      >  1:3373202/1‑36 (MQ=37)
               gcTAATTCCTGGCGGTTATTGGCAGACAGCTGCGcc                    <  1:4773011/36‑1 (MQ=255)
                cTAATTCCTGGCGGTTATTTGCAGACAGCTGCGCCa                   >  1:4209264/1‑36 (MQ=255)
                cTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCa                   <  1:4553093/36‑1 (MQ=255)
                 tAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAg                  <  1:6692140/36‑1 (MQ=255)
                 tAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAg                  <  1:1170525/36‑1 (MQ=255)
                 tAATTCCTGGCGGGTATTTGGAGAGAGGTGTGccct                  >  1:1173503/1‑34 (MQ=2)
                  aaTTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGc                 <  1:5103097/36‑1 (MQ=255)
                   aTTCCTGGCGGTTATTTGCAGACAGCTGCGCCAGCt                >  1:7647005/1‑36 (MQ=255)
                   aTTCCTGGCGGTTATTGGCAGACAGCTGCGCCCGCt                >  1:85647/1‑36 (MQ=255)
                   aTTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCt                <  1:1981390/36‑1 (MQ=255)
                    ttCCTGGCGGTTATTTGGAGACAGCTGTGCCCGCTc               >  1:1973913/1‑36 (MQ=19)
                    ttCCTGGCGGTTATTTGCAGAGAGCTGCGCCAGATc               >  1:481303/1‑36 (MQ=25)
                     tCCTGGCGGTTATTTGCAGAAAGCTGCGCAAGCCcc              >  1:2433008/1‑36 (MQ=19)
                      ccTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCa             <  1:937116/36‑1 (MQ=255)
                       cTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAg            <  1:426550/36‑1 (MQ=255)
                          gCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTtc         <  1:1891704/36‑1 (MQ=255)
                            gtttAGTGGGAGACAGCTGCGCCAGCTCCAGTtcat       <  1:890473/34‑1 (MQ=25)
                            ggTTATTTACAGACAGCTGCGGCCGCTCCAGTtcat       <  1:5558762/36‑1 (MQ=19)
                            ggTTATTGGCAGATAGCTGCGCCAGCTCCAGTtcat       >  1:7210693/1‑36 (MQ=255)
                            ggTTATTGGCAGACAGCTGCGCCAGCTCCAGTtcat       <  1:75280/36‑1 (MQ=255)
                            ggTTATTGGCAGACAGCTGCGCCAGCTCCAGTtcat       >  1:25612/1‑36 (MQ=255)
                            gattATTGGCAGACAGCTTCGCCAGCTCCAGTtcat       <  1:7560816/34‑1 (MQ=37)
                             gTTATTGGCAGACAGCTGCGCCAGCTCCAGTtcatc      <  1:5577172/36‑1 (MQ=255)
                             gTTATTGGCAGACAGCTGCGCCAGCTCCAGTtcatc      <  1:7668834/36‑1 (MQ=255)
                              ttATTGGCAGACAGCTGCGCCAGCTCCAGTtcatca     <  1:3158364/36‑1 (MQ=255)
                              ttATTGGCAGACAGCTGCGCCAGCTCCAGTtcatca     >  1:5528304/1‑36 (MQ=255)
                               tATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCaa    >  1:4545001/1‑36 (MQ=255)
                               tATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCaa    >  1:4295481/1‑36 (MQ=255)
                               tATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCaa    <  1:4866499/36‑1 (MQ=255)
                                aTTTGCATACAGCTGCGCCAGCTCCAGTTCATCAAc   <  1:467187/36‑1 (MQ=25)
                                 ttGGCAGACATCTGCGCCAGCTCCAGTTCATCAACg  <  1:5794869/36‑1 (MQ=255)
                                 ttGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACg  <  1:6594543/36‑1 (MQ=255)
                                   |                                 
TCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACG  >  REL606/4366900‑4366968

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 14 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG