Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,363,338 C→A 11.5% G263G (GGC→GGA yjeM → predicted transporter

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,363,3380C→A11.5% 27.2 ‑0.5 ‑0.0154G263G (GGC→GGAyjeMpredicted transporter
Reads supporting (aligned to +/- strand):  new base (9/9):  ref base (58/78):  total (67/87)
Fisher's exact test strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.58e-01

Only 100 of 160 total aligned reads displayed.

CAATCGGTTATTCGCTGGCAATATTTTTATGGGGCGTCAGCACAAACTGGCAGCAGGTATTAAGTAATG  >  REL606/4363304‑4363372
                                  |                                  
cAATCGGTTATTCGCTGGCAATATTTTTTTGGGGc                                    >  1:19012924/1‑35 (MQ=255)
cAATCGGTTATTCGCTGGCAATATTTTTATGGGGc                                    <  1:19969693/35‑1 (MQ=255)
 aaTCGGTTATTCGCTGGCAATATTTTTTTGGGGCg                                   >  1:15103297/1‑35 (MQ=255)
  atgggTTATTCGCTGGCAATATTTTTATGGGGCGt                                  <  1:14810841/32‑1 (MQ=255)
  aTCGGTTATTCGCTGGCAATATTTTTATGGGGCGt                                  <  1:9387237/35‑1 (MQ=255)
  aTCGGTTATTCGCTGGCAATATTTTTATGGGGCGt                                  <  1:20204328/35‑1 (MQ=255)
  aTCGGTTATTCGCTGGCAATATTTTTATGGGGCGt                                  <  1:19694301/35‑1 (MQ=255)
  aTCGGTTATTCGCTGGCAATATTTTTATGGGGCGt                                  <  1:5817491/35‑1 (MQ=255)
     ggTTATTCGCTGGCAATATTTTTATGGGGCGTCAg                               <  1:14653458/35‑1 (MQ=255)
     ggTTATTCGCTGGCAATATTTTTATGGGGCGTCAg                               <  1:13837572/35‑1 (MQ=255)
     ggTTATTCGCTGGCAATATTTTTATGGGGCGTCAg                               >  1:1610354/1‑35 (MQ=255)
     ggTTATTCGCTGGCAATATTTTTATGGGGCGTCAg                               <  1:320520/35‑1 (MQ=255)
      gTTATTCGCTGGCAATATTTTTATGGGGCGTCAGc                              >  1:6600245/1‑35 (MQ=255)
      gTTATTCGCTGGCAATATTTTTATGGGGCGTCAGc                              >  1:6870170/1‑35 (MQ=255)
      gTTATTCGCTGGCAATATTTTTATGGGGCGTCAGc                              >  1:15161553/1‑35 (MQ=255)
      gTTATTCGCTGGCAATATTTTTATGGGGCGTCAGc                              <  1:3608177/35‑1 (MQ=255)
      gTTATTCGCTGGCAATATTTTTATGGGGAGTCAGc                              >  1:7054221/1‑35 (MQ=255)
       ttATTCGCTGGCAATATTTTTATGGGGCGTCAGca                             >  1:5588952/1‑35 (MQ=255)
       ttATTCGCTGGCAATATTTTTATGGGGCGTCAGca                             <  1:16314029/35‑1 (MQ=255)
       ttATTCGCTGGCAATATTTTTATGGGGCGTCAGca                             <  1:4699152/35‑1 (MQ=255)
       ttATTCGCTGGCAATATTTTTATGGGGCGTCAGca                             <  1:2052341/35‑1 (MQ=255)
       ttATTCGCTGGCAATATTTTTATGGGGCGTCAGca                             >  1:14516825/1‑35 (MQ=255)
        tATTCGCTGGCAATATTTTTATGGGGCGTCAGcac                            >  1:17527473/1‑35 (MQ=255)
        tATTCGCTGGCAATATTTTTATGGGGCGTCAGcac                            >  1:224783/1‑35 (MQ=255)
        tATTCGCTGGCAATATTTTTATGGGGCGTCAGcac                            >  1:14007338/1‑35 (MQ=255)
        tATTCGCTGGCAATATTTTTATGGGGAGTCAGcac                            >  1:9513626/1‑35 (MQ=255)
         aTTCGCTGGCAATATTTTTATGGGGCGTCAGcaca                           <  1:17110336/35‑1 (MQ=255)
         aTTCGCTGGCAATATTTTTATGGGGCGTCAGcaca                           >  1:6534436/1‑35 (MQ=255)
         aTTCGCTGGCAATATTTTTATGGGGCGTCAGcaca                           >  1:903505/1‑35 (MQ=255)
         aTTCGCTGGCAATATTTTTATGGGGCGTCAGcaca                           >  1:1803348/1‑35 (MQ=255)
         aTTCGCTGGCAATATTTTTATGGGGAGTCAGcaca                           >  1:11903269/1‑35 (MQ=255)
          ttCGCTGGCAATATTTTTATGGGGCGTCAGCACaa                          >  1:10345085/1‑35 (MQ=255)
          ttCGCTGGCAATATTTTTATGGGGCGTCAGCACaa                          <  1:3976977/35‑1 (MQ=255)
          ttCGCTGGCAATATTTTTATGGGGCGTCAGCACaa                          <  1:1130194/35‑1 (MQ=255)
           tCGCTGGCAATATTTTTATGGGGCGTCAGCACaaa                         <  1:18443146/35‑1 (MQ=255)
           tCGCTGGCAATATTTTTATGGGGCGTCAGCACaaa                         <  1:4185408/35‑1 (MQ=255)
           tCGCTGGCAATATTTTTATGGGGCGTCAGCACaaa                         <  1:8731119/35‑1 (MQ=255)
            cGCTGGCAATATTTTTATGGGGCGTCAGCACAAAc                        <  1:15796889/35‑1 (MQ=255)
            cGCTGGCAATATTTTTATGGGGCGTCAGCACAAAc                        <  1:19162781/35‑1 (MQ=255)
            cGCTGGCAATATTTTTATGGGGCGTCAGCACAAAc                        <  1:16510174/35‑1 (MQ=255)
             gCTGGCAATATTTTTATGGGGCGTCAGCACAAACt                       >  1:16408615/1‑35 (MQ=255)
             gCTGGCAATATTTTTATGGGGAGTCAGCACAAACt                       <  1:6310903/35‑1 (MQ=255)
              cTGGCAATATTTTTATGGGGCGTCAGCACAAACTg                      <  1:3763709/35‑1 (MQ=255)
              cTGGCAATATTTTTATGGGGCGTCAGCACAAACTg                      >  1:14985250/1‑35 (MQ=255)
              cTGGCAATATTTTTATGGGGCGTCAGCACAAACTg                      >  1:601485/1‑35 (MQ=255)
               tGGCAATATTTTTATGGGGCGTCAGCACAAACTgg                     <  1:15512232/35‑1 (MQ=255)
                ggCAATATTTTTATGGGGCGTCAGCACAAACTGgc                    >  1:18203394/1‑35 (MQ=255)
                ggCAATATTTTTATGGGGAGTCAGCACAAACTGgc                    >  1:10231799/1‑35 (MQ=255)
                 gCAATATTTTTATGGGGCGTCAGCACAAACTGgca                   >  1:12202863/1‑35 (MQ=255)
                 gCAATATTTTTATGGGGAGTCAGCACAAACTGgca                   >  1:13664215/1‑35 (MQ=255)
                  cAATATTTTTATGGGGAGTCAGCACAAACTGgcag                  <  1:7921532/35‑1 (MQ=255)
                   aaTATTTTTATGGGGCGTCAGCACAAACTGgcagc                 >  1:4582657/1‑35 (MQ=255)
                   aaTATTTTTATGGGGCGTCAGCACAAACTGgcagc                 >  1:11081077/1‑35 (MQ=255)
                    atatTTTTATGGGGCGTCAGCACAAACTGgcagca                <  1:19996537/35‑1 (MQ=255)
                    atatTTTTATGGGGCGTCAGCACAAACTGgcagca                >  1:13364469/1‑35 (MQ=255)
                    atatTTTTATGGGGAGTCAGCACAAACTGgcagca                <  1:4097541/35‑1 (MQ=255)
                    atatTTTTATGGGGAGTCAGCACAAACTGgcagca                <  1:10414160/35‑1 (MQ=255)
                     tatTTTTATGGGGCGTCAGCACAAACTGGcagcag               >  1:17067423/1‑35 (MQ=255)
                     tatTTTTATGGGGAGTCAGCACAAACTGGcagcag               <  1:7470616/35‑1 (MQ=255)
                      atTTTTATGGGGCGTCAGCACAAACTGGCAGCAgg              >  1:11277145/1‑35 (MQ=255)
                      atTTTTATGGGGCGTCAGCACAAACTGGCAGCAgg              >  1:9036647/1‑35 (MQ=255)
                      atTTTTATGGGGCGTCAGCACAAACTGGCAGCAgg              <  1:1276015/35‑1 (MQ=255)
                      atTTTTATGGGGCGTCAGCACAAACTGGCAGCAgg              <  1:8100326/35‑1 (MQ=255)
                       tttttATGGGGCGTCAGCACAAACTGGCAGCAggg             >  1:10072186/1‑34 (MQ=255)
                       tttttATGGGGCGTCAGCACAAACTGGCAGCAGGt             <  1:9821566/35‑1 (MQ=255)
                       tttttATGGGGCGTCAGCACAAACTGGCAGCAGGt             <  1:9155700/35‑1 (MQ=255)
                       tttttATGGGGCGTCAGCACAAACTGGCAGCAGGt             >  1:20218677/1‑35 (MQ=255)
                       tttttATGGGGCGTCAGCACAAACTGGCAGCAGGt             <  1:19679213/35‑1 (MQ=255)
                       tttttATGGGGCGTCAGCACAAACTGGCAGCAGGt             <  1:6064400/35‑1 (MQ=255)
                       tttttATGGGGCGTCAGCACAAACTGGAAGCAGGt             <  1:9471143/35‑1 (MQ=255)
                        ttttATGGGGCGTCAGCACAAACTGGCAGCAGGTa            <  1:10159370/35‑1 (MQ=255)
                        ttttATGGGGCGTCAGCACAAACTGGCAGCAGGTa            <  1:16587156/35‑1 (MQ=255)
                        ttttATGGGGCGTCAGCACAAACTGGCAGCAGGTa            >  1:17931514/1‑35 (MQ=255)
                         tttATGGGGCGTCAGCACAAACTGGCAGCAGGGAt           >  1:2561296/1‑35 (MQ=255)
                          ttATGGGGCGTCAGCACAAACTGGCAGCAGGTAtt          <  1:12366246/35‑1 (MQ=255)
                          ttATGGGGCGTCAGCACAAACTGGCAGCAGGTAtt          <  1:13232446/35‑1 (MQ=255)
                           tATGGGGCGTCAGCACAAACTGGCAGGAGGGATTa         >  1:16691506/1‑35 (MQ=37)
                           tATGGGGCGTCAGCACAAACTGGCAGCAGGTATTa         >  1:4384097/1‑35 (MQ=255)
                           tATGGGGCGTCAGCACAAACTGGCAGCAGGTATTa         >  1:4972372/1‑35 (MQ=255)
                           tATGGGGCGTCAGCACAAACTGGCAGCAGGTATTa         <  1:9414123/35‑1 (MQ=255)
                           tATGGGGAGTCAGCACAAACTGGCAGCAGGGAttt         >  1:20292696/1‑34 (MQ=25)
                            aTGGGGCGTCAGCACAAACTGGCAGCAGGTATTaa        <  1:18010610/35‑1 (MQ=255)
                            aTGGGGCGTCAGCACAAACTGGCAGCAGGTATTaa        <  1:19672229/35‑1 (MQ=255)
                            aTGGGGCGTCAGCACAAACTGGCAGCAGGTATTaa        >  1:10817427/1‑35 (MQ=255)
                            aTGGGGAGTCAGCACAAACTGGCAGCAGGTATTaa        <  1:8314285/35‑1 (MQ=255)
                            aTGGGGAGTCAGCACAAACTGGCAGCAGGGATTaa        >  1:15737313/1‑35 (MQ=25)
                             tGGGGCGTCAGCACAAACTGGCAGCAGGTATTAAg       <  1:20539144/35‑1 (MQ=255)
                             tGGGGCGTCAGCACAAACTGGCAGCAGGGATTAAg       >  1:12581834/1‑35 (MQ=255)
                              ggggCGTCAGCACAAACTGGCAGCAGGTATTAAGt      <  1:1245035/35‑1 (MQ=255)
                              ggggCGTCAGCACAAACTGGCAGCAGGTATTAAGt      <  1:13734546/35‑1 (MQ=255)
                              ggggCGTCAGCACAAACTGGCAGCAGGTATTAAGt      >  1:16624762/1‑35 (MQ=255)
                               gggCGTCAGCACAAACTGGCAGCAGGTATTAAGTa     >  1:1614201/1‑35 (MQ=255)
                               gggCGTCAGCACAAACTGGCAGCAGGTATTAAGTa     >  1:17090510/1‑35 (MQ=255)
                               gggCGTCAGCACAAACTGGCAGCAGGTATTAAGTa     <  1:2065955/35‑1 (MQ=255)
                               gggCGTCAGCACAAACTGGCAGCAGGGATTAAGTa     >  1:14159096/1‑35 (MQ=255)
                                ggCGTCAGCACAAACTGGCAGCATGTATTAAGTaa    >  1:8586556/1‑35 (MQ=255)
                                ggCGTCAGCACAAACTGGCAGCAGGTATTAAGTaa    <  1:18331901/35‑1 (MQ=255)
                                ggCGTCAGCACAAACTGGCAGCAGGTATTAAGTaa    >  1:3248842/1‑35 (MQ=255)
                                 gCGTCAGCACAAACTGGCAGCAGGTATTAAGTAAt   <  1:4946584/35‑1 (MQ=255)
                                  cGTCAGCACAAACTGGCAGCAGGTATTAAGTAATg  <  1:15849370/35‑1 (MQ=255)
                                  |                                  
CAATCGGTTATTCGCTGGCAATATTTTTATGGGGCGTCAGCACAAACTGGCAGCAGGTATTAAGTAATG  >  REL606/4363304‑4363372

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.