Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,728,6480C→A 6.1% 2.7 ‑6.6 0.01051313norVanaerobic nitric oxide reductase flavorubredoxin
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (0/8):  ref base (57/40):  total (57/48)
Fisher's exact test strand distribution p-value = 1.35e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

GCAGTGCAGCGTCTGCCAGTGGATTTACGATCCGGCAAAAGGCGAGCCAATGCAGGACGTTGCGCCAGGAA  >  REL606/2728613‑2728683
                                   |                                   
gcagtgcagCGTCTGCCAGTGGATTTACGATCCGGc                                     <  1:5223319/36‑1 (MQ=255)
gcagtgcagCGTCTGCCAGTGGATTTACGATCCGGc                                     <  1:6252118/36‑1 (MQ=255)
gcagtgcagCGTCTGCCAGTGGATTTACGATCCGGc                                     >  1:3414912/1‑36 (MQ=255)
gcagtgcagCGTCTGCCAGTGGA.TTACGATCCGGCa                                    <  1:1141353/36‑1 (MQ=38)
 cagtgcagCGTCTGCCAGTGGATTTACGATCCGGCa                                    <  1:2410959/36‑1 (MQ=255)
  agtgcagCGTCTGCCAGTGGATTTACGATCCGGCaa                                   >  1:4157485/1‑36 (MQ=255)
   gtgcagCGTCTGCCAGTGGATTTACGATCCGGCaaa                                  >  1:4447296/1‑36 (MQ=255)
   gtgcagCGTCTGCCAGTGGATTTACGATCCGGCaaa                                  >  1:1697297/1‑36 (MQ=255)
   gcgcagCGTCTGCCAGTGGATTTACGATCCGGCaaa                                  <  1:3932983/34‑1 (MQ=255)
    tgcagCGTCTGCCTGTGGATTTACGATCCGGCaaaa                                 <  1:7871035/36‑1 (MQ=255)
    tgcagCGTCTGCCAGTGGATTTACGATCCGGCaaaa                                 <  1:7709697/36‑1 (MQ=255)
    tgcagCGTCTGCCAGTGGATTTACGATCCGGCaaaa                                 >  1:2758312/1‑36 (MQ=255)
    tgcagCGTCTGCAAGTGGATTTACGATCCGGCaaaa                                 >  1:1297847/1‑36 (MQ=255)
    tgcagCGTCCGCCAGTGGATTTACGATCCGGCaaaa                                 <  1:4541459/36‑1 (MQ=255)
     gcagCGTCTGCCAGTGGATTTACGATCCGGCAAAAg                                <  1:3099447/36‑1 (MQ=255)
     gcagCGTCTGCCAGTGGATTTACGATCCGGCAAAAg                                >  1:1097335/1‑36 (MQ=255)
     gcagCGTCTGCCAGTGGATTTACGATCCGGCAAAAg                                >  1:8269779/1‑36 (MQ=255)
      cagCGTCTGCCAGTGGATTTACGATCCGGCAAAAgg                               >  1:6358917/1‑36 (MQ=255)
      cagCGTCTGCCAGTGGATTTACGATCCGGCAAAAgg                               >  1:6200234/1‑36 (MQ=255)
      cagCGTCTGCCAGTGGATTTACGATCCGGCAAAAgg                               <  1:4542694/36‑1 (MQ=255)
      cagCGTCTGCCAGTGGATTTACGATCCGGCAAAAgg                               <  1:4614551/36‑1 (MQ=255)
       agCGTCTGCCAGTGGATTTACGATCCGGCAAAAGGc                              <  1:1180066/36‑1 (MQ=255)
       agCGTCTGCCAGTGGATTTACGATCCGGCAAAAGGc                              <  1:3320929/36‑1 (MQ=255)
        gCGTCTGCCAGTGGATTTACGATCCGGCAAAAGGCg                             <  1:5858377/36‑1 (MQ=255)
          gTCTGCCAGTGGATTTACGATCCGGCAAAAGGCGAg                           <  1:5685698/36‑1 (MQ=255)
          gTCTGCCAGTGGATTTACGATCCGGCAAAAGGCGAg                           >  1:7270957/1‑36 (MQ=255)
          gTCTGCCAGTGGATTAACGCTCCGGCATAAAGCgcg                           >  1:7282796/1‑34 (MQ=16)
           tCTGCCAGTGGATTTACGATCCGGCAAAAGGCGAGc                          <  1:5041207/36‑1 (MQ=255)
           tCTGCCAGTGGATTTACGATCCGGCAAAAGGCGAGc                          <  1:5732133/36‑1 (MQ=255)
            cttccAGTGGATTTACGATCCGGCAAAAGGCGAGcc                         <  1:4656587/33‑1 (MQ=255)
            cTGCCAGTGGATTTACGATCCGGCAAAAGGCGAGcc                         <  1:365077/36‑1 (MQ=255)
            cTGCCAGTGGATTTACGATCCGCCAAAAGGCGAGcc                         >  1:7720202/1‑36 (MQ=255)
             tGCCAGTGGATTTACGATCCGGCAAAAGGCGAGCCa                        >  1:647595/1‑36 (MQ=255)
             tGCCAGTGGATTTACGATCCGGCAAAAGGCGAGCCa                        <  1:6954728/36‑1 (MQ=255)
             tGCCAGTGGATTTACGATCCGGCAAAAGGCGAGCCa                        <  1:1993697/36‑1 (MQ=255)
              gCCAGTGGATTTACGATCCGGCAAAAGGCGAGCCaa                       <  1:893343/36‑1 (MQ=255)
              gCCAGTGGATTTACGATCCGGCAAAAGGCGAGCCaa                       <  1:2244019/36‑1 (MQ=255)
              gCCAGTGGATTTACGATCCGGCAAAAGGCGAGCCaa                       <  1:157800/36‑1 (MQ=255)
               ccAGTGGATTTACGATCCGGCAAAAGGCGAGCCAAt                      >  1:4713202/1‑36 (MQ=255)
               ccAGTGGATTTACGATCCGGCAAAAGGCGAGCCAAt                      <  1:658604/36‑1 (MQ=255)
                cAGTGGATTTACGATCCGGCAAAAGGCGAGCCAATg                     >  1:1136859/1‑36 (MQ=255)
                cAGTGGATTTACGATCCGGCAAAAGGCGAGCCAATg                     <  1:7549469/36‑1 (MQ=255)
                cAGTGGATTTACGATCCCGCAAAAGGCGAGCCAATg                     >  1:4300755/1‑36 (MQ=255)
                cAGTGGATTTACGATCCAGCAAAAGGCGAGCCAATg                     >  1:3781302/1‑36 (MQ=255)
                 aGTGGATTTACGATCCGGCAAAAGGCGAGCCAATGc                    >  1:2446892/1‑36 (MQ=255)
                 aGTGGATTTACGATCCGGCAAAAGGCGAGCCAATGc                    >  1:261196/1‑36 (MQ=255)
                 aGTGGATTTACGATCCGGCAAAAGGCGAGCCAATGc                    >  1:3215837/1‑36 (MQ=255)
                 aGTGGATTTACGATCCGGCAAAAGGCGAGCCAATGc                    >  1:1406491/1‑36 (MQ=255)
                  gTGGATTTACGATCCGGCAAAAGGCGAGCCAATGCa                   <  1:2759707/36‑1 (MQ=255)
                  gTGGATTTACGATCCGGCAAAAGGCGAGCCAATGCa                   >  1:7173810/1‑36 (MQ=255)
                  gTGGATTTACGATCCGGCAAAAGGCGAGCCAATGCa                   >  1:271785/1‑36 (MQ=255)
                  gTGGATTTACGATCCGGCAAAAGGCGAGCCAATGCa                   >  1:209032/1‑36 (MQ=255)
                   tGGATTTACGATCCGGCAAAAGGCGAGCCAATGCAg                  >  1:4279689/1‑36 (MQ=255)
                   tGGATTTACGATCCGGCAAAAGGAGAGACAATGCAg                  >  1:3745953/1‑36 (MQ=37)
                    ggATTTACGATCCGGCAAAAGGCGAGCCAATGCAgg                 <  1:4197412/36‑1 (MQ=255)
                    ggATTTACGATCCGGCAAAAGGCGAGCCAATGCAgg                 >  1:1089276/1‑36 (MQ=255)
                    ggATTTACGATCCGGCAAAAGGCGAGCCAAGGCAgg                 >  1:4058698/1‑36 (MQ=255)
                    ggATTTACGATCCGGAAAAAGTGGATCCAATGCAgg                 <  1:2879997/36‑1 (MQ=16)
                     gATTTACGATCCGGCAAAAGGCGAGCCAATGCAGGa                <  1:2333473/36‑1 (MQ=255)
                      ttttACAATCCGGAAAAAGGGGAGCCAATGCAGGAc               <  1:3882067/35‑1 (MQ=25)
                      aTTTACGATCCGGCAAAAGGCGAGCCAATGCAgggc               >  1:2011322/1‑34 (MQ=255)
                      aTTTACGATCCGGCAAAAGGCGAGCCAATGCAGGAc               >  1:3631803/1‑36 (MQ=255)
                      aTTTACGATCCGGCAAAAGGCGAGCCAATGCAGGAc               >  1:259376/1‑36 (MQ=255)
                      aTTTACGATCCGGCAAAAGGCGAGCCAATGCAGGAc               <  1:8437147/36‑1 (MQ=255)
                      aTTTACGATCCGGCAAAAGGCGAGCCAATGCAGGAc               >  1:5626329/1‑36 (MQ=255)
                      aTTTACGATCCGGAAAAAGGCGAGCCAATGCAGGAc               <  1:5633990/36‑1 (MQ=255)
                        ttACGATCCGGCAAAAGGCGAGCCAATGCAGGACGt             >  1:5930725/1‑36 (MQ=255)
                        ttACGATCCGGCAAAAGGCGAGCCAATGCAGGACGt             >  1:2661905/1‑36 (MQ=255)
                        ttACGATCCGGAAAAAGGCGAGCCAATGCAGGACGt             <  1:3802258/36‑1 (MQ=255)
                        taaCGATCCGGAAAAAGGCGAGCCAATGCAGGACGt             <  1:2431642/34‑1 (MQ=37)
                        taaCGACCCGGAAAAAGGCGAGCCAATGCAGGACGt             <  1:3723950/34‑1 (MQ=25)
                         tACGATCCGGCAAAAGGCGCGCCAATGCAGGACGtt            >  1:7473630/1‑36 (MQ=255)
                         tACGATCCGGCAAAAGGCGAGCCAATGCAGGACGtt            <  1:1929039/36‑1 (MQ=255)
                         tACGATCCGGCAAAAGGCGAGCCAATGCAGGACGtt            <  1:6504459/36‑1 (MQ=255)
                          acaaTCCGGCAAAAGGCGAGCCAATGCAGGACGTTg           <  1:234126/33‑1 (MQ=255)
                          aCGATCCGGCAAAAGGCGAGCCAATGCAGGACGTTg           <  1:301181/36‑1 (MQ=255)
                           cGATCCGGCAAAAGGCGAGCCAATGCAGGACGTTgc          <  1:7320116/36‑1 (MQ=255)
                           cGATCCGGCAAAAGGCGAGCCAATGCAGGACGTTgc          <  1:1017481/36‑1 (MQ=255)
                            gATCCGGCAAAAGGCGAGCCAATGCAGGACGTTgcg         >  1:6327623/1‑36 (MQ=255)
                            gATCCGGCAAAAGGCGAGCCAATGCAGGACGTTgcg         >  1:4523881/1‑36 (MQ=255)
                            gATCCGGCAAAAGGCGAGCCAATGCAGGACGTTgcg         >  1:6152922/1‑36 (MQ=255)
                            gATCCGGCAAAAGCCGAGCCAATGCAGGACGTTgcg         >  1:3586750/1‑36 (MQ=255)
                            gATCCGGAAAAAGGCGAGCCAATGCAGGACGTTgcg         <  1:7215140/36‑1 (MQ=255)
                            gATCCGGAAAAAGGCGAGCAAATGCAGGACGTTgcg         <  1:2496452/36‑1 (MQ=25)
                             aTCCGGCAAAAGGCGAGCCAATGCAGGACGTTgcgc        >  1:2440394/1‑36 (MQ=255)
                             aTCCGGCAAAAGGCGAGCCAATGCAGGACGTTgcgc        >  1:2730965/1‑36 (MQ=255)
                              tCCGGCAAAAGGCGAGCCAAGGCAGGACGTTGCGcc       <  1:7849690/36‑1 (MQ=255)
                              tCCGGCAAAAGGCAAGCCAATGCAGGACGTTGCGcc       >  1:4026313/1‑36 (MQ=255)
                               ccGGCAAAAGGCGAGCCAATGCAGGACGTTGCGCCa      >  1:4986540/1‑36 (MQ=255)
                                cGGCAAAAGGCGAGCCAATGCAGGACGTTGCGCCAg     >  1:2326944/1‑36 (MQ=255)
                                cGGCAAAAGGCGAGCCAATGCAGGACGTTGCGCCAg     <  1:4733915/36‑1 (MQ=255)
                                 ggCAAAAGGCGAGCCAATGCAGGACGTTGCGCCAgg    >  1:3711123/1‑36 (MQ=255)
                                 ggCAAAAGGCGAGCCAATGCAGGACGTTGCGCCAgg    >  1:7342433/1‑36 (MQ=255)
                                 ggCAAAAGGCGAGCCAATGCAGGACGTTGCGCCAgg    >  1:2475384/1‑36 (MQ=255)
                                 ggCAAAAGGCGAGCCAATGCAGGACGTTGCGCCAgg    <  1:6449223/36‑1 (MQ=255)
                                  gCAAAAGGCGAGCCAATGCAGGACGTTGCGCCAGGa   >  1:7103494/1‑36 (MQ=255)
                                  gCAAAAGGCGAGCCAATGCAGGACGTTGCGCCAGGa   >  1:535368/1‑36 (MQ=255)
                                  gCAAAAGGCGAGCCAATGCAGGACGGTGCGCCAGGa   >  1:8439243/1‑36 (MQ=255)
                                   cAAAAGGCGAGCCAATGCAGGACGTTGCGCCAGGaa  <  1:1899576/36‑1 (MQ=255)
                                   cAAAAGGCGAGCCAATGCAGGACGTTGCGCCAGGaa  >  1:4810849/1‑36 (MQ=255)
                                   |                                   
GCAGTGCAGCGTCTGCCAGTGGATTTACGATCCGGCAAAAGGCGAGCCAATGCAGGACGTTGCGCCAGGAA  >  REL606/2728613‑2728683

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG