Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,728,1861.→G 1.2% 0.3 ‑4.9 ‑0.0169coding (851/1440 nt)norVanaerobic nitric oxide reductase flavorubredoxin
Rejected: E-value exceeds prediction threshold.
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (13/154):  total (15/154)
Fisher's exact test strand distribution p-value = 7.40e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.28e-01

Only 100 of 220 total aligned reads displayed.

GCCCAGGGAATTGCGGAAACCGACCCACGCGTGGCGGTGAAAATTTTCAACGTCGCCCGAAGCGATAAA  >  REL606/2728152‑2728220
                                                                     
gCCCAGGGAATTGCGGAAACCGACCCACGCGTGGc                                    <  1:3302589/35‑1 (MQ=255)
gCCCAGGGAATTGCGGAAACCGACCCACGCGTGGc                                    <  1:8622854/35‑1 (MQ=255)
gCCCAGGGAATTGCGGAAACCGACCCACGCGTGGc                                    >  1:19950873/1‑35 (MQ=255)
gCCCAGGGAATTGCGGAAACCGACCCACGCGTGGc                                    <  1:19297910/35‑1 (MQ=255)
gCCCAGGGAATTGCGGAAACCGACCCACGCGTGGc                                    <  1:12436494/35‑1 (MQ=255)
gCCCAGGGAATTGCGGAAACCGACCCACGCGTGGc                                    <  1:11484795/35‑1 (MQ=255)
 cccGGGGAATTGCGGAAACCGACCCACGCGTGGCg                                   <  1:18369278/35‑1 (MQ=255)
 cccAGGGAATTGCGGAAACCGACCCACGCGTGGCg                                   <  1:18850236/35‑1 (MQ=255)
 cccAGGGAATTGCGGAAACCGACCCACGCGTGGCg                                   <  1:8127357/35‑1 (MQ=255)
  ccAGGGAATTGCTGAAACCGACCCACGCGTGGCgg                                  >  1:9044983/1‑35 (MQ=255)
  ccAGGGAATTGCGGAAACCGACCCACGCGTGGCgg                                  >  1:16680743/1‑35 (MQ=255)
  ccAGGGAATTGCGGAAACCGACCCACGCGTGGCga                                  >  1:5243747/1‑34 (MQ=255)
  ccAGGGAATTGCGGAAACCCACCCACGCGTGGCgg                                  >  1:4864648/1‑35 (MQ=255)
   cAGGGAATTGCGGAAACCGACCCACGCGTGGCggg                                 >  1:11837833/1‑34 (MQ=255)
   cAGGGAATTGCGGAAACCGACCCACGCGTGGCggg                                 >  1:2189330/1‑34 (MQ=255)
   cAGGGAATTGCGGAAACCGACCCACGCGTGGCGGt                                 <  1:6176006/35‑1 (MQ=255)
   cAGGGAATTGCGGAAACCGACCCACGCGTGGCGGt                                 <  1:12671786/35‑1 (MQ=255)
   cAGGGAATTGCGGAAACCGACCCACGCGTGGCGGt                                 <  1:7540370/35‑1 (MQ=255)
   cAGGGAATTGCGGAAACCGACCCACGCGCGGCggg                                 >  1:351356/1‑34 (MQ=255)
    aGGGAATTGCGGAAACCGACCCACGCGTGGCgggg                                >  1:7265708/1‑33 (MQ=255)
    aGGGAATTGCGGAAACCGACCCACGCGTGGCgggg                                >  1:8081721/1‑33 (MQ=255)
    aGGGAATTGCGGAAACCGACCCACGCGTGGCgggg                                >  1:959459/1‑33 (MQ=255)
    aGGGAATTGCGGAAACCGACCCACGCGTGGCgggg                                >  1:13597710/1‑33 (MQ=255)
    aGGGAATTGCGGAAACCGACCCACGCGTGGCgggg                                >  1:13915626/1‑33 (MQ=255)
    aGGGAATTGCGGAAACCGACCCACGCGTGGCGGTg                                <  1:2278946/35‑1 (MQ=255)
    aGGGAATTGCGGAAACCGACCCACGCGTGGCGGTg                                <  1:14549745/35‑1 (MQ=255)
    aGGGAATTGCGGAAACCGACCCACGCGTGGCGGTg                                <  1:13012947/35‑1 (MQ=255)
    aGGGAATTGCGGAAACCGACCCACGCGTGGCGGTg                                <  1:3874064/35‑1 (MQ=255)
    aGGGAATTGCGGAAACCGACCCACGCGTGGCGGTg                                <  1:16884120/35‑1 (MQ=255)
     gggAATTGCGGAAACCGACCCACGCGTGGCggggg                               >  1:3327548/1‑32 (MQ=255)
     gggAATTGCGGAAACCGACCCACGCGTGGCggggg                               >  1:3336305/1‑32 (MQ=255)
     gggAATTGCGGAAACCGACCCACGCGTGGCGGTGa                               <  1:3623843/35‑1 (MQ=255)
      ggAATTGCGGAAACCGACCCACGCGTGGCGGTGaa                              <  1:174449/35‑1 (MQ=255)
      ggAATTGCGGAAACCGACCCACGCGTGGCGGTGaa                              <  1:20535381/35‑1 (MQ=255)
      ggAATTGCGGAAACCGACCCACGCGTGGCGGTGaa                              <  1:8360665/35‑1 (MQ=255)
      ggAATTGCGGAAACCGACCCACGCGTGGCGGTGaa                              <  1:9863767/35‑1 (MQ=255)
       gAATTGCGGAAACCGACCCACGCGTGGCGGTGaaa                             <  1:16470070/35‑1 (MQ=255)
       gAATTGCGGAAACCGACCAACGCGTGGCGGTGaaa                             <  1:5716704/35‑1 (MQ=255)
        aaTTGCGGAAACCGACCCACGCGTGGCGGTGaaaa                            <  1:15068954/35‑1 (MQ=255)
        aaTTGCGGAAACCGACCCACGCGTGGCGGTGaaaa                            <  1:11390226/35‑1 (MQ=255)
        aaTTGCGGAAACAGACCCACGCGTGGCGGTGaaaa                            <  1:1343258/35‑1 (MQ=255)
         aTTGCGGAAACCGACCCACGCGTGGCGGTGAAAAt                           <  1:15407618/35‑1 (MQ=255)
            gCGGAAACCGACCCACGCGTGGCGGTGAAAAtttt                        <  1:9260292/35‑1 (MQ=255)
            gCGGAAACCGACCCACGCGTGGCGGTGAAAAtttt                        <  1:13648777/35‑1 (MQ=255)
            gCGGAAACAGACCCACGCGTGGCGGTGAAAAtttt                        <  1:14273763/35‑1 (MQ=255)
             cGGAAACCGACCCACGCGTGGCGGTGAAAATTTTc                       <  1:8878689/35‑1 (MQ=255)
             cGGAAACCGACCCACGCGTGGCGGTGAAAATTTTc                       <  1:18055808/35‑1 (MQ=255)
             cGGAAACCGACCCACGCGTGGCGGTGAAAATTTTc                       <  1:11279882/35‑1 (MQ=255)
              ggAAACCGACCCACGCGTGGCGTTGAAAATTTTCa                      >  1:2795921/1‑35 (MQ=255)
              ggAAACCGACCCACGCGTGGCGGTGAAAATTTTCa                      <  1:107360/35‑1 (MQ=255)
              ggAAACCGACCCACGCGTGGCGGTGAAAATTTTCa                      <  1:1281750/35‑1 (MQ=255)
              ggAAACCGACCCACGCGTGGCGGTGAAAATTTTCa                      <  1:8803032/35‑1 (MQ=255)
               gAAACCGACCCACGCGTGGCGGTGAAAATTTTCaa                     <  1:9709078/35‑1 (MQ=255)
                cccccGACACACGCGTGGCGGTGAAAATTTTCAAc                    <  1:6825395/32‑1 (MQ=37)
                aacccGGCCCACGCGTGGCGGTGAAAATTTTCAAc                    <  1:3791129/32‑1 (MQ=25)
                aaaCCGACCCACGCGTGGCGGTGAAAATTTTCAAc                    <  1:15142444/35‑1 (MQ=255)
                aaaCCGACCCACGCGTGGCGGTGAAAATTTTCAAc                    <  1:6695766/35‑1 (MQ=255)
                aaaCCGACCCACGCGTGGCGGTGAAAATTTTCAAc                    <  1:11636320/35‑1 (MQ=255)
                 agccGACCCACGCGTGGCGGTGAAAATTTTCAACg                   <  1:17499713/33‑1 (MQ=255)
                 aaCCGACCCACGCGTGGCGGTGAAAATTTTCAACg                   <  1:10399015/35‑1 (MQ=255)
                 aaCCGACCCACGCGTGGCGGTGAAAATTTTCAACg                   <  1:17914899/35‑1 (MQ=255)
                 aaCCGACCCACGCGTGGCGGTGAAAATTTTCAACg                   <  1:9505596/35‑1 (MQ=255)
                 aaCCGACCCACGCGTGGCGGTGAAAATTTTCAACg                   <  1:16046117/35‑1 (MQ=255)
                 aaCCGACCCACGCGTGGCGGTGAAAATTTTCAACg                   <  1:19590839/35‑1 (MQ=255)
                 aaCCGACCCACGCGTGGCGGTGAAAATTTTCAACg                   <  1:10488618/35‑1 (MQ=255)
                  aCCGACCCACGCGTGGCGGTGAAAATTTTCAACGt                  <  1:3652555/35‑1 (MQ=255)
                  aCCGACCCACGCGTGGCGGTGAAAATTTTCAACGt                  <  1:19497664/35‑1 (MQ=255)
                   ccGACCCACGCGTGGCGGTGAAAATTTTCAACGTc                 <  1:11020739/35‑1 (MQ=255)
                   ccGACCCACGCGTGGCGGTGAAAATTTTCAACGTc                 <  1:3072320/35‑1 (MQ=255)
                   ccGACCCACGCGTGGCGGTGAAAATTTTCAACGTc                 <  1:2088396/35‑1 (MQ=255)
                    cGACCCACGCGTGGCGGTGAAAATTTTCAACGTCg                <  1:4682113/35‑1 (MQ=255)
                    cGACCCACGCGTGGCGGTGAAAATTTTCAACGTCg                <  1:12133772/35‑1 (MQ=255)
                    cGACCCACGCGTGGCGGTGAAAATTTTCAACGTCg                <  1:17980807/35‑1 (MQ=255)
                     gACCCACGCGTGGCGGTGAAAATTTTCAACGTCGc               <  1:4261304/35‑1 (MQ=255)
                      aCCCACGCGTGGCGGTGAAAATTTTCAACGTCGcc              <  1:7655112/35‑1 (MQ=255)
                      aCCCACGCGTGGCGGTGAAAATTTTCAACATCGcc              <  1:8278549/35‑1 (MQ=255)
                       cccACGCGTGGCGGTGAAAATTTTCAACGTCGccc             <  1:7668815/35‑1 (MQ=255)
                       cccACGCGTGGCGGTGAAAATTTTCAACGTCGccc             <  1:17373626/35‑1 (MQ=255)
                       cccACGCGTGGCGGTGAAAATTTTCAACGTCGccc             <  1:20287739/35‑1 (MQ=255)
                       cccACGCGTGGCGGTGAAAATTTTCAACGTCGccc             <  1:2677662/35‑1 (MQ=255)
                       cccACGCGTCGCGGTGAAAATTTTCAACGTCGccc             >  1:6765423/1‑35 (MQ=255)
                         cACGCGTGGCGGTGAAAATTTTCAACGTCGCCCGa           <  1:6352956/35‑1 (MQ=255)
                          aCGCGTGGCGGTGAAAATTTTCAACGTCGCCCGaa          <  1:6800209/35‑1 (MQ=255)
                          aCGCGTGGCGGTGAAAATTTTCAACGTCGCCCGaa          <  1:11817362/35‑1 (MQ=255)
                          aCGCGTGGCGGTGAAAATTTTCAACGTCGCCCGaa          <  1:8497672/35‑1 (MQ=255)
                           cgcgTGGCGGTGAAAATTTTCAACGTCGCCCGAAg         <  1:16780343/35‑1 (MQ=255)
                           cgcgTGGCGGTGAAAATTTTCAACGTCGCCCGAAg         <  1:10835821/35‑1 (MQ=255)
                            gcgTGGCGGTGAAAATTTTCAACGTCGCCCGAAGc        <  1:1022230/35‑1 (MQ=255)
                            gcgTGGCGGTGAAAATTTTCAACGTCGCCCGAAGc        <  1:19920426/35‑1 (MQ=255)
                            gcgTGGCGGTGAAAATTTTCAACGTCGCCCGAAGc        <  1:1402974/35‑1 (MQ=255)
                            gcgTGGCGGTGAAAATTTTCAACGTCGCCCGAAGc        <  1:13176857/35‑1 (MQ=255)
                             cgTGGCGGTGAAAATTTTCAACGTCGCCCGAAGCg       <  1:3745100/35‑1 (MQ=255)
                             cgTGGCGGTGAAAATTTTCAACGTCGCCCGAAGCg       <  1:5441994/35‑1 (MQ=255)
                               tGGCGGTGAAAATTTTCAACGTCGCCCGAAGCGAt     <  1:7961148/35‑1 (MQ=255)
                               tGGCGGTGAAAATTTTCAACGTCGCCCGAAGCGAt     <  1:17308736/35‑1 (MQ=255)
                               tGGCGGTGAAAATTTTCAACGTCGCCCGAAGCGAt     <  1:6955838/35‑1 (MQ=255)
                                ggCGGTGAAAATTTTCAACGTCGCCCGAAGCGATa    <  1:11276068/35‑1 (MQ=255)
                                ggCGGTGAAAATTTTCAACGTCGCCCGAAGCGATa    <  1:5046556/35‑1 (MQ=255)
                                 gCGGTGAAAATTTTCAACGTCGCCCGAAGCGATaa   <  1:7240676/35‑1 (MQ=255)
                                  cGGTGAAAATTTTCAACGTCGCCCGAAGCGATaaa  <  1:14671242/35‑1 (MQ=255)
                                                                     
GCCCAGGGAATTGCGGAAACCGACCCACGCGTGGCGGTGAAAATTTTCAACGTCGCCCGAAGCGATAAA  >  REL606/2728152‑2728220

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.