Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,590,3600T→C 6.8% 0.2 ‑1.6 ‑0.079392dnaTprimosomal protein I
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (1/5):  ref base (36/36):  total (37/42)
Fisher's exact test strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.81e-01

TTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGCCAGTTCT  >  REL606/4590325‑4590395
                                   |                                   
tttGTTGCCACTGCACATGGTGAAAGACTTTACCtt                                     >  1:5004328/1‑36 (MQ=255)
tttGTTGCCACTGCACATGGTGAAAGACTTTACCtt                                     >  1:4163259/1‑36 (MQ=255)
tttGTTGCCACTGCACATGGTGAAAGACTTTACCtt                                     >  1:1488956/1‑36 (MQ=255)
tttGTTGCCACTGCACATGGTGAAAGACTTTACCtt                                     <  1:3306485/36‑1 (MQ=255)
tttGTTGCCACTGCACATGGTGAAAGACTTTACCtt                                     <  1:7251896/36‑1 (MQ=255)
 ttgttgCCACTGCACATGGTGAAAGACTTTACCTTc                                    <  1:5898314/36‑1 (MQ=255)
 ttgttgCCACTGCACATGGTGAAAGACTTTACCTTc                                    <  1:4690620/36‑1 (MQ=255)
 ttgttgCCACTGCACATGGTGAAAGACTTTACCTTc                                    <  1:6802352/36‑1 (MQ=255)
 ttgttgCCACTGCACATGGTGAAAGACTTTACCTTc                                    <  1:699130/36‑1 (MQ=255)
 ttgttgCCACTGCACATGGTGAAAGACTTTACCTTc                                    >  1:7915898/1‑36 (MQ=255)
 ttgttgCCACTGCACATGGTGAAAGACTTTACCTTc                                    <  1:7674720/36‑1 (MQ=255)
  tgttgCCACTGCACATGGTGAAAGCCTTTACCTTcc                                   <  1:1217348/36‑1 (MQ=255)
  tgttgCCACTGCACATGGTGAAAGACTTTACCTTcc                                   <  1:3616720/36‑1 (MQ=255)
    ttgCCACTGCACATGGTGAAATACTTTACCTTCCGc                                 <  1:2948340/36‑1 (MQ=255)
    ttgCCACTGCACATGGTGAAAGACTTTACCTTCCGc                                 <  1:5822567/36‑1 (MQ=255)
    ttgCCACTGCACATGGTGAAAGACTTTACCTTCCGc                                 <  1:613934/36‑1 (MQ=255)
    ttgCCACTGCACATGGTGAAAGACTTTACCTTCCGc                                 <  1:6393218/36‑1 (MQ=255)
    ttgCCACTGCACATGGTGAAAGACTTTACCTTCCGc                                 >  1:643158/1‑36 (MQ=255)
    ttcccACTGCACATGCTGAAAGACTTTACCTTCCGc                                 <  1:466883/33‑1 (MQ=37)
     tgCCCCTGCACATGGTGAAAGACTTTACCTACCGcc                                <  1:1862990/36‑1 (MQ=37)
     tgCCACTGCACATGGTGAAAGACTTTACCTTCCGcc                                <  1:5152932/36‑1 (MQ=255)
     tgCCACTGCACATGGTGAAAGACTTTACCTTCCGcc                                >  1:3247285/1‑36 (MQ=255)
     tgCCACTGCACATGGTGAAAGACTTTACCTTCCGcc                                <  1:683360/36‑1 (MQ=255)
      gCCACTGCACATGGTGAAAGACTTTACCTTCCGCCt                               <  1:3875349/36‑1 (MQ=255)
       ccACTGCACATGGTGAAAGACTTTACCTTCCGCCTg                              <  1:6353115/36‑1 (MQ=255)
        cACTGCACATGGTGAAAGACTTTACCTTCCGCCTGc                             >  1:1495026/1‑36 (MQ=255)
        cACTGCACATGGTGAAAGACTTTACCTTCCGCCTGc                             <  1:3267601/36‑1 (MQ=255)
        cACTGCACATGGTGAAAGACTTTACCTTCCGCCTGc                             >  1:943138/1‑36 (MQ=255)
        cACTGCACATGGTGAAAGACTTTACCTTCCGCCTGc                             <  1:858171/36‑1 (MQ=255)
        cACTGCACATGGTGAAAGACTTTACCTCCCGCCTGc                             >  1:498432/1‑36 (MQ=255)
         aCTGCACATGTTGAAAGACTTTCCCTTCCGCCTGcc                            <  1:246186/36‑1 (MQ=37)
         aCTGCACATGGTGAAAGACTTTACATTCCACCTgtc                            >  1:5200670/1‑34 (MQ=37)
         aCTGCACATGGGGAAAGACTTTACCTTCCGCCTGcc                            >  1:144897/1‑36 (MQ=255)
          cTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCa                           >  1:5439736/1‑36 (MQ=255)
            gCACATGGTGAAAGACTTTCCCTTCCGCCTGCCAGt                         <  1:2703076/36‑1 (MQ=255)
            gCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGt                         >  1:759491/1‑36 (MQ=255)
             cacaTGGTGAAAGACTTTCCCTTCCGCCTGCCAGTa                        <  1:3761417/36‑1 (MQ=255)
             cacaTGGTGAAAGACTTTACCTTCCGCCTGCCAGTa                        <  1:5095797/36‑1 (MQ=255)
             cacaTGGTGAAAGACTTTACCTTCCCCCTGCCAGTa                        <  1:7159068/36‑1 (MQ=255)
              acaTGGTGAAAGACTTTACCTTCCGCCTGCCAGTAg                       >  1:5529827/1‑36 (MQ=255)
              acaTGGTGAAAGACTTTACCTTCCGCCTGCCAGTAg                       >  1:586808/1‑36 (MQ=255)
               caTGGTGAAAGACTTTACCTTCCGCCTGCCAGTAgg                      >  1:566075/1‑36 (MQ=255)
               caTGGTGAAAGACTTTACCTTCCGCCTGCCAGTAgg                      >  1:5593278/1‑36 (MQ=255)
               caTGGTGAAAGACTTTACCTTCCGCCTGCCAGTAgg                      <  1:6189764/36‑1 (MQ=255)
               caTGGTGAAAGACTTTACCTTCCGCCTGCCAGTAgg                      <  1:1342891/36‑1 (MQ=255)
               caTGGTGAAAGACTTTACCTTCCGCCTGCCAGTAgg                      >  1:874639/1‑36 (MQ=255)
                aTGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGc                     >  1:7418277/1‑36 (MQ=255)
                aTGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGc                     >  1:7103320/1‑36 (MQ=255)
                aTGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGc                     >  1:3240264/1‑36 (MQ=255)
                aTGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGc                     >  1:7507549/1‑36 (MQ=255)
                aTGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGc                     <  1:6812789/36‑1 (MQ=255)
                aTGGTGAAAAACTTTACCTCCCGCCTGCCAGTAGGc                     <  1:988981/36‑1 (MQ=37)
                 tGGTAAACGACTTTACCTTCCCCCTGCCAGTAGGCa                    <  1:7715350/36‑1 (MQ=25)
                  ggTGAAAGACTTTCCCTTCCGCCTGCCAGTAGGCaa                   <  1:3049919/36‑1 (MQ=255)
                  ggTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCaa                   >  1:3491211/1‑36 (MQ=255)
                   gTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCaat                  >  1:1615036/1‑36 (MQ=255)
                   gTGAAAGACTTTACCTCCCCCCTGCCAGTAGGCaat                  <  1:6621089/36‑1 (MQ=38)
                   atcaaaGACTTTCCCCCCCGCCTGCCAGTAGGCaat                  <  1:1820892/33‑1 (MQ=19)
                    tGAAAGACTTTACCTTCCGCCTGCCAGTAGGCaatg                 <  1:7549161/36‑1 (MQ=255)
                    tGAAAGACTTTACCTTCCGCCTGCCAGTAGGCaatg                 >  1:6995722/1‑36 (MQ=255)
                     gAAAGACTTTACCTTCCGCCTGCCAGTAGGCaatga                >  1:3110235/1‑36 (MQ=255)
                      aaaGACTTTACCTTCCGCCTGCCAGTAGGCaatgaa               <  1:3480436/36‑1 (MQ=255)
                      aaaGACTTTACCTTCCGCCTGCCAGTAGGCaatgaa               >  1:3343278/1‑36 (MQ=255)
                       aaGACTTTACCTTCCGCCTGCCAGTAGGCaatgaat              >  1:4611953/1‑36 (MQ=255)
                       aaGACTTAACCTTCCGCCTGCCAGTAGGCaatgaat              <  1:7275044/36‑1 (MQ=255)
                        aGACTTTACCTTCCGCCTGCCAGTAGGCaatgaatg             >  1:672598/1‑36 (MQ=255)
                        aGACTTTACCTTCCGCCTGCCAGTAGGCaatgaatg             >  1:1755448/1‑36 (MQ=255)
                         gACTTTACCTTCCGCCTGCCAGTAGGCAatgaatga            <  1:3940618/36‑1 (MQ=255)
                         gACTTTACCTTCCGCCTGCCAGTAGGCAatgaatga            >  1:6399726/1‑36 (MQ=255)
                         gACTTTACCTTCCGCCTGCCAGTAGGCAatgaatga            >  1:6157781/1‑36 (MQ=255)
                          aCTTTACCTTCCGCCTGCCAGTAGACAATGAATGAg           >  1:5802279/1‑36 (MQ=255)
                            tttACCTTCCGCCTGCCAGTAGGCAATGAATGAGGc         >  1:1475208/1‑36 (MQ=255)
                            tttACCTTCCGCCTGCCAGTAGGCAATGAATGAGGc         >  1:2856638/1‑36 (MQ=255)
                             ttcccTCCCACCTGCCAGTAGGCAATGAATGAGGcc        <  1:7215610/33‑1 (MQ=25)
                             ttACCTTCCGCCTGCCAGTAGGCGATGAATGAGGcc        >  1:7783946/1‑36 (MQ=255)
                             ttACCTTCCGCCTGCCAGTAGGCAATGAATGAGGcc        <  1:7447079/36‑1 (MQ=255)
                             ttACCTTCCGCCTGCCAGTAGGCAATGAATGAGGcc        >  1:7707471/1‑36 (MQ=255)
                             ttACCTTCCGCCTGCCAGTAGGCAATGAATGAGGcc        <  1:922467/36‑1 (MQ=255)
                             ttACCTTCCGCCTGCCAGTAGGCAATGAATGAGGcc        >  1:5270041/1‑36 (MQ=255)
                              tACCTTCCGCCTGCCAGTAGGCAATGAATGAGGCCa       >  1:2593262/1‑36 (MQ=255)
                              tACCTTCCCCCTGCCAGTAGGCAATGAATGAGGCCa       <  1:1321845/36‑1 (MQ=255)
                                ccTTCCGCCTGCCAGTAGGCAATGAATGAGGCCAGt     >  1:3696650/1‑36 (MQ=255)
                                ccTCCCGCCTGCCAGTAGGCAATGAATGAGGCCAGt     <  1:2234054/36‑1 (MQ=255)
                                 cTTCCGCCTGCCAGTAGGCAATGAATGAGGCCAGtt    <  1:1507455/36‑1 (MQ=255)
                                  ttCCGCCTGCCAGTAGGCAATGAATGAGGCCAGttc   <  1:169278/36‑1 (MQ=255)
                                   tCCGCCTGCCAGTAGGCAATGAATGAGGCCAGttct  <  1:5204940/36‑1 (MQ=255)
                                   |                                   
TTTGTTGCCACTGCACATGGTGAAAGACTTTACCTTCCGCCTGCCAGTAGGCAATGAATGAGGCCAGTTCT  >  REL606/4590325‑4590395

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