Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,254,4920G→A 6.6% 0.2 ‑5.1 0.086622argGargininosuccinate synthase
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (0/8):  ref base (40/37):  total (40/46)
Fisher's exact test strand distribution p-value = 6.08e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCCACTCCAGCGTCAAAATCGTCAACCC  >  REL606/3254457‑3254526
                                   |                                  
aCATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                    >  1:4601418/1‑36 (MQ=255)
aCATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                    <  1:6593910/36‑1 (MQ=255)
aCATGCTTGGTGCAACGCATGAAGCGAAGGATCTgg                                    <  1:7204375/36‑1 (MQ=255)
 cATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGa                                   <  1:2384450/36‑1 (MQ=255)
 cATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGa                                   >  1:943458/1‑36 (MQ=255)
  aTGCTTGGTGCAACGCATGAAGCGAAGGATCTGGaa                                  <  1:2748863/36‑1 (MQ=255)
   tGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAAt                                 <  1:2344684/36‑1 (MQ=255)
   tGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAAt                                 <  1:2416808/36‑1 (MQ=255)
   tGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAAt                                 >  1:1813140/1‑36 (MQ=255)
    gCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATa                                >  1:3976739/1‑36 (MQ=255)
    gCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATa                                >  1:779642/1‑36 (MQ=255)
     cTTGGTGCAACGCATGAAGCGAAGGATCTGGAATAc                               <  1:4058048/36‑1 (MQ=255)
      ttGGTGCAACGCATGAAGCGAAGGATCTGGAATAcc                              <  1:842623/36‑1 (MQ=255)
      ttGGTGCAACGCATGAAGCGAAGGATCTGGAATAcc                              <  1:337029/36‑1 (MQ=255)
       tGGTGCAACGCATGAAGCGAAGGATCTGGAATACct                             >  1:1156612/1‑36 (MQ=255)
       tGGTGCAACGCATGAAGCGAAGGATCTGGAATACct                             <  1:2569897/36‑1 (MQ=255)
        ggTGCAACGCATGAAGCGAAGGATCTGGAATACctc                            <  1:2728395/36‑1 (MQ=255)
        ggTGCAACGCATGAAGCGAAGGATCTGGAATACctc                            <  1:3281561/36‑1 (MQ=255)
        ggTGCAACGCATGAAGCGAAGGATCTGGAATACctc                            >  1:5833342/1‑36 (MQ=255)
         gTGCAACGCATGAAGCGAAGGATCTGGAATACctcc                           >  1:6000548/1‑36 (MQ=255)
         gTGCAACGCATGAAGCGAAGGATCTGGAATACctcc                           >  1:1047201/1‑36 (MQ=255)
          tGCAACGCATGAAGCGAAGGATCTGGAATACctcca                          <  1:4323716/36‑1 (MQ=255)
          tGCAACGCATGAAGCGAAGGATCTGGAATACctcca                          <  1:179440/36‑1 (MQ=255)
          tGCAACGCATGAAGCGAAGGATCTGGAATACctcca                          >  1:1526286/1‑36 (MQ=255)
          tGCAACGCATGAAGCGAAGGATCTGGAATACctcca                          <  1:5518963/36‑1 (MQ=255)
           gCAACGCATGAAGCGAAGGATCTGGAATACctccac                         >  1:4716144/1‑36 (MQ=255)
            cAACGAATAAAGCGAAGGATCTGGAATACctccact                        <  1:2061237/36‑1 (MQ=37)
             aaCGCATGAAGCGAAGGATCTGGAATACctccactc                       >  1:2487865/1‑36 (MQ=255)
             aaCGCATGAAGCGAAGGATCTGGAATACctccactc                       <  1:3402205/36‑1 (MQ=255)
              aCGCATGAAGCGAAGGATCTGGAATACctccactcc                      >  1:6519635/1‑36 (MQ=255)
               cGCATGAAGCGAAGGATCTGGAATACctccactcca                     <  1:1340799/36‑1 (MQ=255)
                gCATGAAGCGAAGTATCTGGAATACCTCCACTCCAg                    >  1:3212488/1‑36 (MQ=255)
                 cATGAAGCGAAGGATCTGGAATACCTCCACTCCAGc                   >  1:6883400/1‑36 (MQ=255)
                 cATGAAGCGAAGGATCTGGAATACCTCCACTCCAGc                   <  1:7310308/36‑1 (MQ=255)
                 cATGAAGCGAAGGATCTGGAATACCTCCACTCCAGc                   <  1:5346470/36‑1 (MQ=255)
                 cATGAAGCGAAGGATCNGGAATACCTCCACTCCAGc                   <  1:2148436/36‑1 (MQ=255)
                  aTGAAGCGAAGGATCTGGAATACCTCCACTCCAGCg                  >  1:6443987/1‑36 (MQ=255)
                  aTGAAGCGAAGGATCTGGAATACCTCCACTCCAGCg                  >  1:7378303/1‑36 (MQ=255)
                   tGAAGCGAAGGATCTGGAATACCTCCACTCCAGCGt                 <  1:5340029/36‑1 (MQ=255)
                   tGAAGCGAAGGATCTGGAATACCTCCACTCCAGCGt                 <  1:1061314/36‑1 (MQ=255)
                    gAAGCGAAGGATCTGGAATACCTCCACTCCAGCGTc                <  1:1054987/36‑1 (MQ=255)
                    gAAGCGAAGGATCTGGAATACCTCCACTCCAGCGTc                <  1:637472/36‑1 (MQ=255)
                    gAAGCGAAGGATCTGGAATACCTCCACTCCAGCGTc                >  1:5837456/1‑36 (MQ=255)
                    gAAGCGAAGGATCTGGAATACCTCCACTCCAGCGTc                >  1:2364012/1‑36 (MQ=255)
                    gAAGCGAAGGATCTGGAATACCTCCACTCCAGCGTc                >  1:4796025/1‑36 (MQ=255)
                     aaGCGAAGGATCTGGAATACCTCCACTCCAGCGTCa               <  1:174331/36‑1 (MQ=255)
                     aaGCGAAGGATCTGGAATACCTCCACTCCAGCGTCa               <  1:6577076/36‑1 (MQ=255)
                      aGCGAAGGATCTGGAATACCTCCACTCCAGCGTCaa              <  1:5930971/36‑1 (MQ=255)
                      aGCGAAGGATCTGGAATACCTCCACTCCAGCGTCaa              <  1:942889/36‑1 (MQ=255)
                      aGCGAAGGATCTGGAATACCTCCACTCCAGCGTCaa              >  1:2440444/1‑36 (MQ=255)
                      aGCGAAGGATCTGGAATACCTCCACTCCACCGCCaa              >  1:7148640/1‑36 (MQ=37)
                      aGCGAAGGATCTGAAATACCTCCACTCCAGCGTCaa              <  1:6959411/36‑1 (MQ=255)
                      aGCGAAGGATCTGAAATACCTCCACTCCAGCGTCaa              <  1:1588941/36‑1 (MQ=255)
                        cGAAGGATCTGGAATACCTCCACTCCAGCGTCaaaa            <  1:6018449/36‑1 (MQ=255)
                        cGAAGGATCTGGAATACCTCCACTCCAGCGTCaaaa            >  1:746210/1‑36 (MQ=255)
                        cGAAGCATCTGCAATACCTCCACTCCAGCGTCaaaa            <  1:4988028/36‑1 (MQ=37)
                         gAAGGATCTGGAATACCTCCACTCCAGCGTCAAAAt           >  1:4705681/1‑36 (MQ=255)
                         gAAGGATCTGGAATACCTCCACTCCAGCGTCAAAAt           >  1:658856/1‑36 (MQ=255)
                         gAAGGATCTGGAATACCTCCACTCCAGCGTCAAAAt           >  1:1708913/1‑36 (MQ=255)
                         gAAGGATCTGGAATACCTCCACTCCAGCGTCAAAAt           >  1:2799284/1‑36 (MQ=255)
                         gAAGGATCTGGAATACCTCCACTCCAGCGTCAAAAt           >  1:5986635/1‑36 (MQ=255)
                         gAAGGATCTGGAATACCTCCACTCCAGCGTCAAAAt           <  1:7439993/36‑1 (MQ=255)
                          aaGGATCTGGAATACCTCCACGCCAGCGTCAAAATc          >  1:6788937/1‑36 (MQ=255)
                          aaGGATCTGGAAAACCTCCACTCCAGCGTCAAAATc          <  1:95271/36‑1 (MQ=255)
                          aaGGATCTGAAATACCTCCACTCCAGCGTCAAAATc          <  1:4461267/36‑1 (MQ=255)
                          aaGAATCTGAAATACCTCCACTCCAGCGTCAAAATc          <  1:3705179/36‑1 (MQ=37)
                           aGGATCTGGAATACCTCCACTCCAGCGTCACAATCg         >  1:1394198/1‑36 (MQ=255)
                           aGGATCTGGAATACCTCCACTCCAGCGTCAAAATCg         >  1:1860065/1‑36 (MQ=255)
                           aGGATCTGGAATACCTCCACTCCAGCGTCAAAATCg         >  1:6772366/1‑36 (MQ=255)
                           aGGATCTGGAATACCTCCACAAAAGGGTCAAAATCg         >  1:3066766/1‑36 (MQ=16)
                            ggATCTGGAATACCTCCACTCCAGCGTCAAAATCGt        <  1:831535/36‑1 (MQ=255)
                            ggATCTAGAATACCTCCACTCCAGCGTCAAAATCGt        <  1:141811/36‑1 (MQ=255)
                             gATCTGGAATACCTCCACTCCAGCGTCAAAATCGTc       >  1:496408/1‑36 (MQ=255)
                             gATCTGGAATACCTCCACTCCAGCGTCAAAATCGTc       <  1:152332/36‑1 (MQ=255)
                             aaTCTGAAATACCTCCACTCCAGCGTCAAAATCGTc       <  1:3980924/35‑1 (MQ=255)
                              ctCTGGAATACCTCCACTCCAGCGTCAAAATCGTCa      <  1:5057116/35‑1 (MQ=255)
                              aTCTGGAATCCCTCCACTCCAGCGTCAAAATCGTCa      <  1:4407175/36‑1 (MQ=255)
                              aTCTGGAATACCTCCACTCCAGCGTCAAAATCGTCa      >  1:2126570/1‑36 (MQ=255)
                              aTCTGGAATACCTCCACGCCGGCGTCACAATCGTCa      >  1:4210310/1‑36 (MQ=37)
                               tCTGGAATACCTCCACTCCAGCGTCAAAATCGTCaa     >  1:6504032/1‑36 (MQ=255)
                               tCTGGAATACCTCCACTCCAGCGTCAAAATCGTCaa     <  1:2676645/36‑1 (MQ=255)
                               tCTGAAATACCTCCACTCCAGCGTCAAAATCGTCaa     <  1:4560749/36‑1 (MQ=255)
                               tCTGAAATACCTCCACTCCAGCGTCAAAATCGTCaa     <  1:1737121/36‑1 (MQ=255)
                               tCTGAAATACCTCCACTCCAGCGTCAAAATCGTCaa     <  1:5628334/36‑1 (MQ=255)
                                cTGGAATACCTCCACTCCAGCGTCAAAATCGTCAAc    >  1:5387297/1‑36 (MQ=255)
                                 tGGAATACCTCCACTCCAGCGTCAAAATCGTCAAcc   >  1:4912289/1‑36 (MQ=255)
                                 tGGAATACCTCCACTCCAGCGTCAAAATCGTCAAcc   <  1:373297/36‑1 (MQ=255)
                                 tGGAATACCTCCACTCCAGCGTCAAAATCGTCAAcc   >  1:2624090/1‑36 (MQ=255)
                                 tGGAATACCTCCACTCCAGCGTCAAAATAGTAAAcc   >  1:5873128/1‑36 (MQ=37)
                                  ggAATACCTCCACTCCAGCGTCAAAATCGTCAAccc  <  1:1136826/36‑1 (MQ=255)
                                   |                                  
ACATGCTTGGTGCAACGCATGAAGCGAAGGATCTGGAATACCTCCACTCCAGCGTCAAAATCGTCAACCC  >  REL606/3254457‑3254526

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG