Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,122,7000T→C 8.3% 0.9 ‑5.8 0.090811mdoCglucans biosynthesis protein
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (0/9):  ref base (42/39):  total (42/48)
Fisher's exact test strand distribution p-value = 3.01e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGTA  >  REL606/1122665‑1122735
                                   |                                   
gAGGACCATGGTGATCACCGACTCGGTTTCGTACAt                                     <  1:6764546/36‑1 (MQ=255)
gAGGACCATGGTGATCACCGACTCGGTTTCGTACAt                                     >  1:1283955/1‑36 (MQ=255)
gAGGACCATGGTGATCACCGACTCGGTTTCGTACAt                                     >  1:303547/1‑36 (MQ=255)
gAGGACCATGGTGATCACCGACTCGGTTTCGTACAt                                     >  1:7476567/1‑36 (MQ=255)
gAGGACCATGGTGATCACCGACTCGGTTTCGTACAt                                     >  1:2871568/1‑36 (MQ=255)
 aGGACCATGGTGATCACCGACTCGGTTTCGTACATc                                    <  1:220502/36‑1 (MQ=255)
 aGGACCATGGTGATCACCGACTCGGTTTCGTACATc                                    >  1:1404674/1‑36 (MQ=255)
 aGGACCATGGTGATCACCGACTCGGTTTCGTACATc                                    >  1:3674015/1‑36 (MQ=255)
 aGGACCATGGTGATCACCGACTCGGTTTCGTACATc                                    <  1:1634222/36‑1 (MQ=255)
 aGGACAACGGTGTTCGCCGACTCGGTTTCGTACATc                                    <  1:2297601/36‑1 (MQ=16)
  ggACCATGGTGATCACCGACTCGGTTTCGTACATcc                                   <  1:6634503/36‑1 (MQ=255)
  ggACCATGGTGATCACCGACTCGGTTTCGTACATcc                                   >  1:2378873/1‑36 (MQ=255)
  ggACCATGGTGATCACCGACTCGGTTTCGTACATcc                                   <  1:2020317/36‑1 (MQ=255)
   gACCATGGTGATCACCGACTCGGTTTCGTACATCCa                                  <  1:6589348/36‑1 (MQ=255)
    aCCATGGTGATCCCCGACTCGGTTTCGTACATCCAg                                 <  1:1238524/36‑1 (MQ=255)
    aCCATGGTGATCACCGACTCGGTTTCGTACATCCAg                                 >  1:5891284/1‑36 (MQ=255)
    aCCATGGTGATCACCGACTCGGTTTCGTACATCCAg                                 <  1:1151885/36‑1 (MQ=255)
     ccATGGTGATCGCCGACTCGGTTTCGTACATCCAgg                                <  1:2834738/36‑1 (MQ=255)
     ccATGGTGATCACCGACTCGGTTTCGTACATCCAgg                                <  1:6196117/36‑1 (MQ=255)
      cATGGTGATCACCGACTCGGTTTCGTACATCCAGGc                               >  1:4334724/1‑36 (MQ=255)
       aTGGTGATCACCGACTCGGTTTCGTACATCCAGGCa                              >  1:7325442/1‑36 (MQ=255)
       aTGGTGATCACCGACTCGGTTTCGTACATCCAGGCa                              >  1:354143/1‑36 (MQ=255)
       aTGGTGATCACCGACTCGGTTTCGTACATCCAGGCa                              >  1:5816958/1‑36 (MQ=255)
        tGGTGATCACCGACTCGGTTTCGTACATCCAGGCAt                             <  1:5819573/36‑1 (MQ=255)
          gTGATCACCGACTCGGTTTCGTACATCCAGGCATCg                           <  1:241511/36‑1 (MQ=255)
          gTGATCACCGACTCGGTTTCGTACATCCAGGCATCg                           >  1:5346573/1‑36 (MQ=255)
           tGATCACCGACTCGGTTTCGTACATCCAGGCATCGc                          >  1:4891935/1‑36 (MQ=255)
           tGATCACCGACTCGGTTTCGTACATCCAGGCATCGc                          >  1:320077/1‑36 (MQ=255)
            gATCGCCGACTCGGTTTCGTGCATCCAGGCATCGcc                         <  1:6715502/36‑1 (MQ=37)
            gATCACCGACTCGGTTTCGTACATCCAGGCATCGcc                         >  1:7389853/1‑36 (MQ=255)
              tCACCGACTCGGTTTCGTACATCCAGGCATCGCCAc                       <  1:6931039/36‑1 (MQ=255)
              tCACCGACTCGGTTTCGTACATCCAGGCATCGCCAc                       <  1:7177083/36‑1 (MQ=255)
              tCACCGACTCGGTTTCGTACATCCAGGCATCGCCAc                       >  1:1248702/1‑36 (MQ=255)
              tCACCGACTCGGTTTCGTACATCCAGGCATCGCCAc                       >  1:4667323/1‑36 (MQ=255)
               ccccGACTCGGTTTCGTACACCCAGGCATCGCCACt                      <  1:4040831/34‑1 (MQ=38)
               cACCGACTCGGTTTCGTACATCCAGGCATCGCCACt                      >  1:5636904/1‑36 (MQ=255)
               cACCGACTCGGTTTCGTACATCCAGGCATCGCCACt                      >  1:4176020/1‑36 (MQ=255)
                aCCGACTCGGTTTCGTACATCCATGCATCCTCACTg                     >  1:3825495/1‑36 (MQ=25)
                aCCGACTCGGTTTCGTACATCCAGGCATCGCCACTg                     >  1:344460/1‑36 (MQ=255)
                aCCGACTCGGTTTCGTACATCCAGGCATCGCCACTg                     >  1:7461146/1‑36 (MQ=255)
                aCCGACTCGGTTTCGTACATCCAGGCATCGCCACTg                     >  1:4438816/1‑36 (MQ=255)
                 ccGACTCGGTTTCGTACATCCAGGCATCGCCACTGc                    >  1:3691865/1‑36 (MQ=255)
                 ccGACTCGGTTTCGTACATCCAGGCATCGCCACTGc                    <  1:1132818/36‑1 (MQ=255)
                  cGACTCGGTTTCGTACATCCAGGCATCGCCACTGcc                   <  1:6314972/36‑1 (MQ=255)
                  cGACTCGGTTTCGTACATCCAGGCATCGCCACTGcc                   >  1:4024992/1‑36 (MQ=255)
                  cGACTCGGTTTCGTACATCCAGGCATCGCCACTGcc                   <  1:4227756/36‑1 (MQ=255)
                  cGACTCGGTTTCGTACATCCAGGCATCGCCACTGcc                   <  1:919162/36‑1 (MQ=255)
                  cGACTCGGTTTCGTACATCCAGGCATCGCCACTGcc                   <  1:1391227/36‑1 (MQ=255)
                   gACTCGGTTTCGTACATCCAGGCATCGCCACTGCCa                  >  1:2722797/1‑36 (MQ=255)
                   gACTCGGTTTCGTACATCCAGGCATCGCCACTGCCa                  <  1:6495459/36‑1 (MQ=255)
                   gACTCGGTTTCGTACATCCAGGCATCGCCACTGCCa                  >  1:1846663/1‑36 (MQ=255)
                    aCTCGGTTTCGTACATCCAGGCATCGCCACTGCCAt                 <  1:93021/36‑1 (MQ=255)
                     cccGGTTCCGAAAACCCAGTCACCGCCACTGCCATa                <  1:5612976/34‑1 (MQ=2)
                      tCGGTTTCGTACATCCAGGCATCGCCACTGCCATAg               >  1:6422462/1‑36 (MQ=255)
                      tCGGTTTCGTACATCCAGGCATCGCCACTGCCATAg               >  1:3678683/1‑36 (MQ=255)
                      tCGGTTTCGTACATCAAGGCATCGCCACTGCCATAg               <  1:5950480/36‑1 (MQ=255)
                       cGGTTTCGTACATACAGCAATCCCCACTGCCATAgc              <  1:6337349/36‑1 (MQ=19)
                        ggTTTCGTACATCCAGGCATCGCCACTGCCATAgcg             <  1:685385/36‑1 (MQ=255)
                        ggTTTCGTACACCCAGGCATCGCCACTGCCATAgcg             <  1:2477775/36‑1 (MQ=255)
                         gTTTCGTACATCCAGGCATCGCCACTGCCATAgcgc            >  1:441074/1‑36 (MQ=255)
                         gTTTCGTACATCCAGGCATCGCCACTGCCATAgcgc            >  1:6322699/1‑36 (MQ=255)
                          tttCGTACATCCAGGCATCGCCACTGCCATAGCGCt           <  1:3499832/36‑1 (MQ=255)
                          tttCGTACACCCAGGAATCGCCACTGCCATAGCGCt           <  1:5407420/36‑1 (MQ=37)
                           ttCGTACATCCAGGCATCGCCACTGCCATAGCGCTg          >  1:1648375/1‑36 (MQ=255)
                           ttCGTACATCAAGGCATCGCCACTGCCATAGCGCTg          <  1:7055521/36‑1 (MQ=255)
                           ttCGTACACCCAGGCATCGCCACTGCCATAGCGCTg          <  1:4426743/36‑1 (MQ=255)
                           ttCGAACACCCAGGCATCTCCACTGCCATAGCGCTg          <  1:5160956/36‑1 (MQ=25)
                            tCGTACATCCAGGCATCGCCACTGCCATAGCGCTgg         >  1:6882321/1‑36 (MQ=255)
                            tCGTACATCCAGGCATCGCCACTGCCATAGCGCTgg         >  1:5940183/1‑36 (MQ=255)
                            tCGTACATCCAGGCATCGCCACTGCCATAGCGCTgg         >  1:2024485/1‑36 (MQ=255)
                            tCGTACATCCAGGCATCGCCACTGCCATAGCGCTgg         <  1:7503496/36‑1 (MQ=255)
                            tCGTACATCCAGGCATCGCCACTGCCATAGCGCTgg         >  1:3926450/1‑36 (MQ=255)
                             cGTACATCCAGGCATCGCCACTGCCCTAGCGCTGGt        >  1:4587930/1‑36 (MQ=255)
                             cGTACATCCAGGCATCGCCACTGCCATAGCGCTGGt        >  1:2616151/1‑36 (MQ=255)
                             cGTACATCCAGGCATCGCCACTGCCATAGCGCTGGt        <  1:1392442/36‑1 (MQ=255)
                             cGTACATCCAGGCATCGACACTGCCATAGAGCTGGt        >  1:1462644/1‑36 (MQ=38)
                             cGTAAATCCAGGCATCGCCACTGCCATAGCGCTGGt        <  1:4374221/36‑1 (MQ=255)
                              gTACATCCAGTCATCGCCACTGCCATAGCGCTGGtt       <  1:6376093/36‑1 (MQ=255)
                              gTACATCCAGGCATCGCCACTGCCATAGCGCTGGtt       >  1:7241382/1‑36 (MQ=255)
                              gTACATCCAGGCATCGCCACTGCCATAGCGCTGGtt       >  1:2698362/1‑36 (MQ=255)
                              gTACATCCAGGCATCGCCACTGCCATAGCGCTGGtt       <  1:218300/36‑1 (MQ=255)
                              gTACATCCAGGCATCGCCACTGCCATAGCGCTGGtt       >  1:3071627/1‑36 (MQ=255)
                              gTACATCCAGGAATCGCCACTGCCATAGCGCTGGtt       <  1:911226/36‑1 (MQ=255)
                              gTACACCCAGGCATCGCCACTGCCATAGCGCTGGtt       <  1:5602064/36‑1 (MQ=255)
                               tACATCCAGGCATCGCCACTGCCATAGCGCTGGTTg      <  1:6739698/36‑1 (MQ=255)
                               tACACCCAGGCATCGCCACTGCCATAGCGCTGGTTg      <  1:1231362/36‑1 (MQ=255)
                                aCATCCAGTAATTGCCTCTGCCATCGCGCTGGTTGa     <  1:4254422/36‑1 (MQ=16)
                                aCATCCAGGCATCGCCACTGCCATAGCGCTGGTTGa     >  1:145819/1‑36 (MQ=255)
                                aCACCCAGGCATCGCCACTGCCATAGCGCTGGTTGa     <  1:1591040/36‑1 (MQ=255)
                                  aTCCGGGCATCGCAACTGCCATAGCGCTGGTTGAGt   <  1:7379434/36‑1 (MQ=37)
                                  aTCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGt   >  1:1768515/1‑36 (MQ=255)
                                  aTCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGt   <  1:1829907/36‑1 (MQ=255)
                                  aTCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGt   <  1:4843911/36‑1 (MQ=255)
                                  aTCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGt   >  1:2495698/1‑36 (MQ=255)
                                  aTCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGt   <  1:3648535/36‑1 (MQ=255)
                                   tCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGTa  >  1:6822823/1‑36 (MQ=255)
                                   tCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGTa  >  1:3388707/1‑36 (MQ=255)
                                   |                                   
GAGGACCATGGTGATCACCGACTCGGTTTCGTACATCCAGGCATCGCCACTGCCATAGCGCTGGTTGAGTA  >  REL606/1122665‑1122735

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG