Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,159,0850C→A 9.6% 0.3 ‑11.5 0.068596rplA50S ribosomal protein L1
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (0/8):  ref base (49/11):  total (49/19)
Fisher's exact test strand distribution p-value = 1.02e-05
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGA  >  REL606/4159051‑4159120
                                  |                                   
aaCCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCa                                    >  1:5576679/1‑36 (MQ=255)
aaCCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCa                                    >  1:1407003/1‑36 (MQ=255)
 aCCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCaa                                   >  1:5007411/1‑36 (MQ=255)
  ccTTGAAGCTCTGCTGGTTGCGCTGAAAAAAGCaaa                                  <  1:6037926/36‑1 (MQ=255)
   cTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCaaaa                                 <  1:6576102/36‑1 (MQ=255)
   cTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCaaaa                                 <  1:2614739/36‑1 (MQ=255)
    tGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAAc                                <  1:1058494/36‑1 (MQ=255)
    tGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAAc                                >  1:4910552/1‑36 (MQ=255)
    tGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAAc                                >  1:384029/1‑36 (MQ=255)
    tGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAAc                                >  1:3833465/1‑36 (MQ=255)
    tGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAAc                                <  1:3814595/36‑1 (MQ=255)
     ggAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAAcc                               >  1:4388001/1‑36 (MQ=255)
     ggAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAAcc                               >  1:1860406/1‑36 (MQ=255)
     ggAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAAcc                               >  1:155865/1‑36 (MQ=255)
       aaGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGa                             >  1:2498110/1‑36 (MQ=255)
       aaGCTCTGCTGGCTGCGCTGAAAAAAGCAAAACCGa                             >  1:7016736/1‑36 (MQ=255)
        aGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGAc                            <  1:5992426/36‑1 (MQ=255)
         gCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACt                           >  1:6709210/1‑36 (MQ=255)
          ctctGCTGGTTGCGCTGAAAAAAGCAAAACCGACTc                          >  1:5948365/1‑36 (MQ=255)
           tctGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCa                         >  1:1812002/1‑36 (MQ=255)
            ctgctgGTTGCGCTGAAAAAAGCAAAACCGACTCAg                        <  1:1748304/36‑1 (MQ=255)
            ctgctgGTTGCGCTGAAAAAAGCAAAACCGACTCAg                        >  1:2900312/1‑36 (MQ=255)
            ctgctgGTTGCGCTGAAAAAAGCAAAACCGACTCAg                        >  1:123169/1‑36 (MQ=255)
               ctgGTTGCGCTGACAAAAGCAAAACCGACTCAGGCg                     >  1:3239139/1‑36 (MQ=255)
                tgGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGa                    >  1:7342907/1‑36 (MQ=255)
                tgGTTGCGCTGAAAAAAGAAAAACCGACTCAGGCGa                    <  1:2468724/36‑1 (MQ=255)
                 ggTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGaa                   >  1:6243152/1‑36 (MQ=255)
                  gTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGaaa                  >  1:920656/1‑36 (MQ=255)
                   ttGCGCTGAAAAAAGCAAAAACGACTCAGGCGAAAg                 >  1:810635/1‑36 (MQ=255)
                   ttGCGCTGAAAAAAGAAAAACCGACTCAGGCGAAAg                 <  1:6673116/36‑1 (MQ=255)
                    tGCGCTGAAAAAAGAAAAACCGACTCAGGCGAAAgg                <  1:1235162/36‑1 (MQ=255)
                    tGCGATGAAAAAAGAAAAACCGACTCAGGCGAAAgg                <  1:1163998/36‑1 (MQ=37)
                     gcgcTGAAAAAAGCAAAGCCGACGCAGGCGAAAGGc               <  1:431582/36‑1 (MQ=37)
                     gcgcTGAAAAAAGCAAAACCGACTCAGGCGAAAGGc               >  1:218417/1‑36 (MQ=255)
                     gcgcTGAAAAAAGAAAAACCGACTCAGGCGAAAGGc               <  1:3212076/36‑1 (MQ=255)
                     gcgcGGAAAAAAGAAAAACCGACTCAGGCGAAAGGc               <  1:1485477/36‑1 (MQ=37)
                      cgcTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCg              >  1:6202230/1‑36 (MQ=255)
                      cgcTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCg              >  1:3609439/1‑36 (MQ=255)
                      cgcTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCg              >  1:2402328/1‑36 (MQ=255)
                       gcTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCgt             >  1:5115089/1‑36 (MQ=255)
                       gcTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCgt             >  1:2793890/1‑36 (MQ=255)
                       gcTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCgt             >  1:3815202/1‑36 (MQ=255)
                       gcTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCgt             >  1:2485867/1‑36 (MQ=255)
                        cTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCgtg            >  1:5307840/1‑36 (MQ=255)
                        cTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCgtg            >  1:4180999/1‑36 (MQ=255)
                         tGAAAAAAGCAAAACCGACTCAGGCGAAAGGCgtgt           >  1:3505316/1‑36 (MQ=255)
                         tGAAAAAAGCAAAACCGACTCAGGCGAAAGGCgtgt           >  1:7112081/1‑36 (MQ=255)
                         tGAAAAAAGCAAAACCGACTCAGGCGAAAGGCgtgt           >  1:7188984/1‑36 (MQ=255)
                         tGAAAAAAGCAAAACCGACTCAGGCGAAAGGCgtgt           >  1:2479067/1‑36 (MQ=255)
                          gAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTa          >  1:6224576/1‑36 (MQ=255)
                          aaaaaaaCAAAAACCGACTCAGGCGAAAGGCGTGTa          <  1:2538847/35‑1 (MQ=37)
                            aaaaaGCAAAACCGACTCAGGCGAAAGGCGTGTACa        <  1:4583302/36‑1 (MQ=255)
                            aaaaaGCAAAACCGACTCAGGCGAAAGGCGTGTACa        >  1:6520946/1‑36 (MQ=255)
                            aaaaaGCAAAACCGACTCAGGCGAAAGGCGTGTACa        >  1:5653360/1‑36 (MQ=255)
                            aaaaaGCAAAACCGACTCAGGCGAAAGGCGTGTACa        <  1:4903338/36‑1 (MQ=255)
                             aaaaGCAAAACCGTCTCAGGCGAAAGGAGTGTACAt       >  1:4285715/1‑36 (MQ=37)
                             aaaaGCAAAACCGACTCAGGCGAAAGGCGTGTACAt       >  1:4844943/1‑36 (MQ=255)
                             aaaaGCAAAACCGACTCAGGCGAAAGGCGTGTACAt       <  1:2050972/36‑1 (MQ=255)
                             aaaaGCAAAACCGACTCAGGCGAAAGGCGTGCACAt       >  1:6378878/1‑36 (MQ=255)
                             aaaaGAAAAACCGACTCAGGCGAAAGGCGTGTACAt       <  1:976367/36‑1 (MQ=255)
                              aaaGCAAAACCGACTCAGGCGAAAGGCGTGTACATc      >  1:4372587/1‑36 (MQ=255)
                              aaaGCAAAACCGACTCAGGCGAAAGGCGTGTACATc      >  1:579390/1‑36 (MQ=255)
                              aaaGCAAAACCGACTCAGGCGAAAGGCGTGTACATc      >  1:1544838/1‑36 (MQ=255)
                               aaGCAAAACCGACTCAGGCGAAAGGCGTGTACATCa     >  1:978909/1‑36 (MQ=255)
                                aGCAAAACCGACTCAGGCGAAAGGCGTGTACATCaa    >  1:6694875/1‑36 (MQ=255)
                                aGCAAAACCGACTCAGGCGAAAGGCGTGTACATCaa    >  1:6447761/1‑36 (MQ=255)
                                aGCAAAACCGACTCAGGCGAAAGGCGTGTACATCaa    >  1:6054150/1‑36 (MQ=255)
                                 gCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAg   >  1:6153914/1‑36 (MQ=255)
                                  cAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGa  >  1:5858241/1‑36 (MQ=255)
                                  |                                   
AACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGA  >  REL606/4159051‑4159120

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG