Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,739,5410C→A 2.5% 0.9 ‑1.3 ‑0.2246A24A (GCG→GCThycFformate hydrogenlyase complex iron‑sulfur protein
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (2/5):  ref base (122/116):  total (124/122)
Fisher's exact test strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.28e-01

Only 100 of 284 total aligned reads displayed.

TCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCG  >  REL606/2739507‑2739575
                                  |                                  
tCTGCTCTGGCTTACCACGGAAGGTTTTATCAAcc                                    >  1:3534878/1‑35 (MQ=255)
tCTGCTCTGGCTTACCACGGAAGGTTTTATCAAcc                                    >  1:15837740/1‑35 (MQ=255)
 cTGCTCTGGCTTACCACGGAAGTTTTTATCAAccc                                   >  1:5390782/1‑34 (MQ=255)
 cTGCTCTGGCTTACCACGGAAGTTTTTATCAACCg                                   >  1:4846247/1‑35 (MQ=255)
 cTGCTCTGGCTTACCACGGAAGTTTTTATCAACCg                                   <  1:16029968/35‑1 (MQ=255)
 cTGCTCTGGCTTACCACGGAAGTTTTTATCAACCg                                   >  1:15157434/1‑35 (MQ=255)
 cTGCTCTGGCTTACCACGGAAGGTTTTATCAACCg                                   >  1:18428400/1‑35 (MQ=255)
  ttctctGGCTTCCCACGGAAGTTTTTATCAACCGc                                  <  1:16862682/33‑1 (MQ=37)
  tGCTCTGGCTTCCCACGGATGTTTTTATCAACCGc                                  <  1:2966951/35‑1 (MQ=37)
  tGCTCTGGCTTACCACGGAAGGTTTTATCAACCGc                                  >  1:8671453/1‑35 (MQ=255)
   gCTCTGGCTTACCACGGAAGTTTTTATCAACCGCa                                 <  1:14759138/35‑1 (MQ=255)
   gCTCTGGCTTACCACGGAAGTTTTTATCAACCGCa                                 >  1:15615493/1‑35 (MQ=255)
   gCTCTGGCTTACCACGGAAGTTTTTATCAACCGCa                                 <  1:16763320/35‑1 (MQ=255)
   gCTCTGGCTTACCACGGAAGTTTTTATCAAACGCa                                 <  1:7551489/35‑1 (MQ=255)
    ctctGGCTTACCGCGGAAGGTTTTATCAAACGCaa                                >  1:5849814/1‑35 (MQ=25)
    ctctGGCTTACCACGGAAGTTTTTATCAACCGCaa                                >  1:11917611/1‑35 (MQ=255)
    ctctGGCTTACCACGGAAGTTTTTATCAACCGCaa                                >  1:11464264/1‑35 (MQ=255)
    ctctGGCTTACCACGGAAGTTTTTATAAACCGCaa                                <  1:15337233/35‑1 (MQ=255)
     tctGGCTTACCACGTAAGTTTTTATCAACCGCAAt                               >  1:385425/1‑35 (MQ=255)
     tctGGCTTACCACGGAAGTTTTTATCAACCGCAAt                               <  1:2750811/35‑1 (MQ=255)
     tctGGCTTACCACGGAAGTTTTTATCAACCGCAAt                               <  1:16677355/35‑1 (MQ=255)
      ctGGCTTTCCACGGAAGTTTTTATCAACCGCAATc                              <  1:14864426/35‑1 (MQ=255)
      ctGGCTTACCACGGAAGTTTTTATCAACCGCAATc                              <  1:5648678/35‑1 (MQ=255)
      ctGGCTTACCACGGAAGTTTTTATCAACCGCAATc                              >  1:16516936/1‑35 (MQ=255)
      ctGGCTTACCACGGAAGGTTTTATCCACCGCAATc                              >  1:17115612/1‑35 (MQ=37)
       tGGCTTACCACGTAAGTTTTTATCAACCGCAATCg                             >  1:16212971/1‑35 (MQ=255)
       tGGCTTACCACGGAAGTTTTTATCAACCGCAATCg                             >  1:15477794/1‑35 (MQ=255)
       tGGCTTACCACGGAAGTTTTTATCAACCGCAATCg                             <  1:10595945/35‑1 (MQ=255)
       tGGCTTACCACGGAAGTTTTTATCAACCGCAATCg                             >  1:6490826/1‑35 (MQ=255)
       tGGCTTACCACGGAAGTTTTTATCAACCGCAATCg                             <  1:6594664/35‑1 (MQ=255)
       tGGCTTACCACGGAAGTTTTTATCAACCGCAATCg                             >  1:6650899/1‑35 (MQ=255)
       tGGCTTACCACGGAAGTTTTTATCAACCGCAATCg                             >  1:5203582/1‑35 (MQ=255)
       tGGCTTACCACGGAAGTTTTTATCAACCGCAATCg                             >  1:19261670/1‑35 (MQ=255)
        ggCTTACCACGGAAGTTTTTATCAACCGCAATCgg                            >  1:13400427/1‑35 (MQ=255)
        ggCTTACCACGGAAGTTTTTATCAACCGCAATCgg                            <  1:10232446/35‑1 (MQ=255)
        ggCTTACCACGGAAGTTTTTATCAACCGCAATCgg                            >  1:12649448/1‑35 (MQ=255)
        ggCTTACCACGGAAGTTTTTATCAACCGCAATCgg                            >  1:19928079/1‑35 (MQ=255)
         gCTTACCACGGAAGTTTTTATCAACCGCAATCGGc                           >  1:18282553/1‑35 (MQ=255)
         gCTTACCACGGAAGTTTTTATCAACCGCAATCGGc                           >  1:2058133/1‑35 (MQ=255)
         gCTTACCACGGAAGTTTTTATCAACCGCAATCGGc                           >  1:14921799/1‑35 (MQ=255)
         gCTTACCACGGAAGTTTTTATCAACCGCAATCGGc                           >  1:14178066/1‑35 (MQ=255)
         gCTTACCACGGAAGTTTTTATCAACCGCAATCGGc                           >  1:11737127/1‑35 (MQ=255)
         gCTTACCACGGAAGTTTTTATCAACCGCAATCGGc                           <  1:12120950/35‑1 (MQ=255)
          cTTACCACGGAAGTTTTTATCAACCGCAATCGGCt                          <  1:20147474/35‑1 (MQ=255)
          cTTACCACGGAAGTTTTTATCAACCGCAATCGGCt                          <  1:10110963/35‑1 (MQ=255)
          cTTACCACGGAAGTTTTTATCAACCGCAATCGGCt                          >  1:18115059/1‑35 (MQ=255)
           ttACCACGGAAGTTTTTATCAACCGCAATCGGCTc                         >  1:12382698/1‑35 (MQ=255)
            tACCACGGAAGTTTTTATCAACCGCCATCGGGTcc                        >  1:13866487/1‑35 (MQ=25)
            tACCACGGAAGTTTTTATCAACCGCAATCGGCTcc                        >  1:14338314/1‑35 (MQ=255)
            tACCACGGAAGTTTTTATCAACCGCAATCGGCTcc                        <  1:904319/35‑1 (MQ=255)
             aCCACGGAAGTTTTTATCAACCGCAATCGGCTCCa                       <  1:17333511/35‑1 (MQ=255)
             aCCACGGAAGTTTTTATCAACCGCAATCGGCTCCa                       <  1:4937875/35‑1 (MQ=255)
             aCCACGGAAGTTTTTATCAACCGCAATCGGCTCCa                       <  1:8486352/35‑1 (MQ=255)
             aCCACGGAAGTTTTTATCAACCGCAATCGGCTCCa                       <  1:3367445/35‑1 (MQ=255)
              ccACGGAAGTTTTTATCAACCGCAATCGGCTcccg                      >  1:3203843/1‑33 (MQ=255)
              ccACGGAAGTTTTTATCAACCGCAATCGGCTCCAg                      <  1:1269351/35‑1 (MQ=255)
              ccACGGAAGTTTTTATCAACCGCAATCGGCTCCAg                      <  1:5809602/35‑1 (MQ=255)
              ccACGGAAGGTTTTATCAACCGCAATCGGCTcccc                      >  1:5100195/1‑33 (MQ=37)
               cACGGAAGTTTTTATCAACCGCAATCGGCTCCAGc                     <  1:10783231/35‑1 (MQ=255)
               cACGGAAGTTTTTATCAACAGCAATCGGCTCCAGc                     <  1:4139531/35‑1 (MQ=255)
                 cGGGAGTTGTTATCAACCGCAATCGGCTCCAGCgg                   <  1:3576005/35‑1 (MQ=37)
                 cGGAAGTTTTTATCAACCGCAATCGGCTCCAGCgg                   <  1:3628587/35‑1 (MQ=255)
                 cGGAAGTTTTTATCAACCGCAATCGGCTCCAGCgg                   <  1:291777/35‑1 (MQ=255)
                  ggAAGTTTTTATCAACCGCAATCGGCTCCCGCggg                  >  1:6774777/1‑34 (MQ=255)
                  ggAAGTTTTTATCAACCGCAATCGGCTCCAGCGGa                  <  1:13306552/35‑1 (MQ=255)
                   gAAGTTTTTATCACCCGCAATCGGCTCCAGCGGAt                 <  1:16228760/35‑1 (MQ=255)
                   gAAGTTTTTATCAACCGCAATCGGCTCCAGCGGAt                 >  1:630699/1‑35 (MQ=255)
                   gAAGTTTTTATCAACCGCAATCGGCTCCAGCGGAt                 >  1:9839363/1‑35 (MQ=255)
                   gAAGTTTTTATCAACCGCAATCGGCTCCAGCGGAt                 >  1:2476444/1‑35 (MQ=255)
                    aaGTTTTTATCAACCGCAATCGGCTCCAGCGGATa                >  1:11170534/1‑35 (MQ=255)
                     aGTTTTTATCAACCGCAATCGGCTCCCGCGGGTaa               >  1:7427172/1‑35 (MQ=37)
                     aGTTTTTATCAACCGCAATCGGCTCCAGCGGATaa               >  1:17756270/1‑35 (MQ=255)
                     aGTTTTTATCAACCGCAATCGGCTCCAGCGGATaa               <  1:17034624/35‑1 (MQ=255)
                      gTTTTTATCAACCGCAATCGGCTCCAGCGGGTAAg              >  1:13544645/1‑35 (MQ=255)
                      gTTTTTATCAACCGCAATCGGCTCCAGCGGATAAg              >  1:15737358/1‑35 (MQ=255)
                      gTTTTTATCAACCGCAATCGGCTCCAGCGGATAAg              <  1:18790602/35‑1 (MQ=255)
                      gTTTTTATCAACCGCAATCGGCTCCAGCGGATAAg              >  1:8304482/1‑35 (MQ=255)
                       tttttATCAACCGCAATCGGCTCCAGCGGATAAGa             >  1:20249683/1‑35 (MQ=255)
                       tttttATCAACCGCAATCGGCTCCAGCGGATAAGa             <  1:19530743/35‑1 (MQ=255)
                       tttttATCAACCGCAATCGGCTCCAGCGGATAAGa             <  1:5973958/35‑1 (MQ=255)
                        ttttATCAACCGCAATCGGCTCCAGCGGATAAGAc            >  1:7231650/1‑35 (MQ=255)
                         tttATCAACCGCAATCGGCTCCAGCGGATAAGACg           >  1:7922054/1‑35 (MQ=255)
                         tttATCAACCGCAATCGGCTCCAGCGGATAAGACg           >  1:10381830/1‑35 (MQ=255)
                         tttATCAACCGCAATCGGCTCCAGCGGATAAGACg           <  1:9291085/35‑1 (MQ=255)
                         tttATCAACCGCAATCGGCTCCAGCGGATAAGACg           <  1:10033720/35‑1 (MQ=255)
                          ttATCAACCGCAATCGGCTCCAGCGGATAAGACGa          >  1:17212367/1‑35 (MQ=255)
                          ttATCAACCGCAATCGGCTCCAGCGGATAAGACGa          <  1:17839137/35‑1 (MQ=255)
                          ttATCAACCGCAATCGGCTCCAGCGGATAAGACGa          <  1:4357858/35‑1 (MQ=255)
                            aTCAACCGCAATCGGCTCCAGCGGATAAGACGAgg        >  1:12535550/1‑35 (MQ=255)
                            aTCAACCGCAATCGGCTCCAGCGGATAAGACGAgg        <  1:10348707/35‑1 (MQ=255)
                            aTCAACCGCAATCGGCTCCAGCGGATAAGACGAgg        <  1:9579657/35‑1 (MQ=255)
                              cAACAGCAATCGGCTCCAGCGGATAAGACGAGGTc      <  1:12972887/35‑1 (MQ=255)
                               aaCCGCAATCGGCTCCAGCGGATAAGACGAGGTCg     <  1:10993032/35‑1 (MQ=255)
                                aCCGCAATCGGCTCCAGCGGATAAGACGAGGTCGc    <  1:14002210/35‑1 (MQ=255)
                                aCCGCAATCGGCTCCAGCGGATAAGACGAGGTCGc    <  1:8718628/35‑1 (MQ=255)
                                aCCGCAATCGGCTCCAGCGGATAAGACGAGGTCGc    <  1:2296249/35‑1 (MQ=255)
                                aCCGCAATCGGCTCCAGCGGATAAGACGAGGGCGc    >  1:3464561/1‑35 (MQ=255)
                                 ccGCAATCGTCTCCAGCGGATAAGACGAGGTCGcc   >  1:18516840/1‑35 (MQ=255)
                                 ccGAAATCGGCTCCAGCGGATAAGACGAGGTCGcc   <  1:7073443/35‑1 (MQ=255)
                                 ccGAAATCGGCTCCAGCGGATAAGACGAGGTCGcc   <  1:17871952/35‑1 (MQ=255)
                                  |                                  
TCTGCTCTGGCTTACCACGGAAGTTTTTATCAACCGCAATCGGCTCCAGCGGATAAGACGAGGTCGCCG  >  REL606/2739507‑2739575

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ 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.