Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 2,499,315 G→A 36.6% intergenic (‑110/‑179) maeB ← / → talA malic enzyme/transaldolase A

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,499,3150G→A36.6% 51.6 ‑3.3 ‑1.150intergenic (‑110/‑179)maeB/talAmalic enzyme/transaldolase A
Reads supporting (aligned to +/- strand):  new base (4/14):  ref base (17/14):  total (22/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 = 7.13e-02

TATTACCCTGCCGGAGTCTACAGACTTTGAGCAAGTCCAAACTCTCACCATTAATATAATGTTTTGGTA  >  REL606/2499281‑2499349
                                  |                                  
tattaCCCTGCCGGAGTCTACAGACTTTGAGCAAGt                                   >  1:1475142/1‑36 (MQ=255)
tattaCCCTGCCGGAGTCTACAGACTTTGAGCAAGt                                   >  1:5767336/1‑36 (MQ=255)
tattaCCCTGCCGGAGTCTACAGACTTTGAGCAAGt                                   >  1:2656321/1‑36 (MQ=255)
 attaCCCTGCCGGAGTCTACAGACTTTGAGCAAGTc                                  >  1:2791438/1‑36 (MQ=255)
 attaCCCTGCCGGAGTCTACAGACTTTGAGCAAATc                                  <  1:325024/36‑1 (MQ=255)
  ttaCCCTGCCGGAGTCTACAGACTTTGAGCAAGTcc                                 <  1:5831472/36‑1 (MQ=255)
  ttaCCCTGCCGGAGTCTACAGACTTTGAGCAAGTcc                                 >  1:3618106/1‑36 (MQ=255)
  ttaCCCTGCCGGAGTCTACAGACTTTGAGCAAGTcc                                 <  1:5584689/36‑1 (MQ=255)
   taCCCTGCCGGAGTCTACAGACTTTGAGCAATTCCa                                >  1:4964255/1‑36 (MQ=255)
   taCCCTGCCGGAGTCTACAGACTTTGAGCAAATCCa                                <  1:4047579/36‑1 (MQ=255)
    aCCCTGCCGGAGTCTACAGACTTTGAGCAAGTCCaa                               <  1:720915/36‑1 (MQ=255)
      ccTGCCGGAGTCTACAGACTTTGAGCAAATCCAAAc                             <  1:2328268/36‑1 (MQ=255)
       cTGCCGGAGTCTACAGACTTTGAGCAAGTCCAAAct                            >  1:5608429/1‑36 (MQ=255)
       cTGCCGGAGTCTACAGACTTTGAGCAAATCCAAAct                            <  1:2712642/36‑1 (MQ=255)
         gCCGGAGTCTACAGACTTTGAGCAAGTCCAAActct                          >  1:6157217/1‑36 (MQ=255)
          ccGGAGTCTACAGACTTTGAGCAAGTCCAAACtctc                         <  1:3116433/36‑1 (MQ=255)
          ccGGAGTCTACAGACTTTGAGCAAATCCAAACcctc                         <  1:5746625/36‑1 (MQ=255)
           cGGAGTCTACAGACTTTGAGCAAGTCCAAACTCTCa                        >  1:2345843/1‑36 (MQ=255)
           cGGAGTCTACAGACTTTGAGCAAGTCCAAACTCTCa                        >  1:883510/1‑36 (MQ=255)
            ggAGTCTACAGACTTTGAGCAAGTCCAAACTCTCAc                       <  1:6567302/36‑1 (MQ=255)
             gAGTCTACAGACTTTGAGCAAATCCAAACTCTCAcc                      <  1:3433651/36‑1 (MQ=255)
               gTCTACAGACTTTGATCAAATCCAAACTCTCACCAt                    <  1:449594/36‑1 (MQ=38)
               gTCTACAGACTTTGAGCAAATCCAAACTCTCACCAt                    <  1:6065295/36‑1 (MQ=255)
               gTCTACAGACTTTGAGCAAATCCAAACTCTCACCAt                    <  1:7408112/36‑1 (MQ=255)
               gTCTACAGACTTTGAGCAAATCCAAACTCTCACCAt                    <  1:158234/36‑1 (MQ=255)
                tCTACAGACTTTGAGCAAGTCCAAACTCTCACCAtt                   <  1:5238537/36‑1 (MQ=255)
                tCTACAGACTTTGAGCAAATCCAAACTCTCACCAtt                   <  1:306478/36‑1 (MQ=255)
                 cTACAGACTTTGAGCAAATCCAAACTCTCACCATTa                  >  1:2592108/1‑36 (MQ=255)
                 cTACAGACTTTGAGCAAATCCAAACTCTCACCATTa                  <  1:7378551/36‑1 (MQ=255)
                  tACAGACTTTGAGCAAGTCCAAACTCTCACCATTaa                 <  1:4297836/36‑1 (MQ=255)
                  tACAGACTTTGAGCAAGTCCAAACTCTCACCATTaa                 >  1:7831067/1‑36 (MQ=255)
                  tACAGACTTTCAGCAAGTCCAAACTCTCACCATTaa                 <  1:44429/36‑1 (MQ=255)
                   aCAGACTTTGAGCAAATCCAAACTCTCACCATTAat                >  1:5309788/1‑36 (MQ=255)
                   aCAGACTTTGAGCAAATCCAAACGCTCACCATTAat                <  1:2661325/36‑1 (MQ=37)
                    cAGACTTTGAGCAAATCCAAACTCTCACCATTAata               <  1:2415763/36‑1 (MQ=255)
                     aGACTTTGAGCAAGTCCAAACTCTCACCATTAatat              <  1:7177009/36‑1 (MQ=255)
                     aGACTTTGAGCAAGTCCAAACTCTCACCATTAatat              <  1:1449477/36‑1 (MQ=255)
                     aGACTTTGAGCAAATCCAAACTCTCACCATTAatat              >  1:3328364/1‑36 (MQ=255)
                        cTTTGAGCAAGTCCAAACTCTCACCATTAATATAAt           >  1:1801617/1‑36 (MQ=255)
                          ttGAGCAAGTCCAAACTCTCACCATTAATATAATGt         >  1:3949431/1‑36 (MQ=255)
                          ttGAGCAAGTCCAAACTCTCACCATTAATATAATGt         <  1:28925/36‑1 (MQ=255)
                           tGAGCAAGTCCAAACTCTCACCATTAATATAATGtt        >  1:7518597/1‑36 (MQ=255)
                            gAGCAAGTCCAAACTCTCACCATTAATATAATGttt       >  1:452634/1‑36 (MQ=255)
                             aGCAAGTCCAAACTCTCACCCTTAATATAATGtttt      >  1:3914068/1‑36 (MQ=255)
                             aGCAAATCCAAACTCTCACCATTAATATAATGtttt      >  1:6777268/1‑36 (MQ=255)
                              gCAAGTCCAAACTCTCACCATTAATATAATGTTTTg     <  1:212955/36‑1 (MQ=255)
                               cAAGTCCAAACTCTCACCATTAATATAATGTTTTgg    <  1:5712616/36‑1 (MQ=255)
                                aaGTCCAAACTCTCACCATTAATATAATGTTTTGGt   >  1:521730/1‑36 (MQ=255)
                                aaGTCCAAACTCTCACCATTAATATAATGTTTTGGt   >  1:2612790/1‑36 (MQ=255)
                                 aGTCCAAACTCTCACCATTAATATAATGTTTTGGta  <  1:7891234/36‑1 (MQ=255)
                                  |                                  
TATTACCCTGCCGGAGTCTACAGACTTTGAGCAAGTCCAAACTCTCACCATTAATATAATGTTTTGGTA  >  REL606/2499281‑2499349

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG