Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,726,9991.→C 1.9% 7.0 ‑8.5 ‑0.0215coding (151/1515 nt)ygaAanaerobic nitric oxide reductase transcription regulator
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (188/23):  total (188/27)
Fisher's exact test strand distribution p-value = 2.03e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.71e-01

Only 100 of 268 total aligned reads displayed.

TCCTTTGCCAGACCGTCGATGGCAAGCGGAATAAACTGCCGCGAATCGTAACGTAGCAACGCAGACGCA  >  REL606/2726965‑2727033
                                   |                                  
tCCTTTGCCAGCCCGTCGATGGCAAGCGGAATaaa                                     <  1:11558911/35‑1 (MQ=255)
tCCTTTGCCAGCCCGTCGATGGCAAGCGGAATaaa                                     <  1:14132144/35‑1 (MQ=255)
 cttttGCCAGACCGTCGATGGCAAGCGGAATAAAc                                   <  1:11712315/33‑1 (MQ=255)
 ccTTTGCCAGACCGTCGATGGCAAGCGGAATAAAc                                   >  1:10041563/1‑35 (MQ=255)
 ccTTTGCCAGACCGTCGATGGCAAGCGGAATAAAc                                   >  1:7148114/1‑35 (MQ=255)
  cTTTGCCAGCCCGTCGATGGCAAGCGGAATAAACt                                  <  1:3920267/35‑1 (MQ=255)
  cTTTGCCAGCCCGTCGATGGCAAGCGGAATAAACt                                  <  1:5773567/35‑1 (MQ=255)
  cTTTGCCAGACCGTCGATGGCAAGCGGAATAAACt                                  >  1:1077816/1‑35 (MQ=255)
  cTTTGCCAGACCGTCGATGGCAAGCGGAATAAACt                                  >  1:15671159/1‑35 (MQ=255)
  cTTTGCCAGACCGTCGATGGCAAGCGGAATAAACt                                  <  1:14480824/35‑1 (MQ=255)
   tttGCCAGGCCGTCGATGGCAAGCGGAATAAACTg                                 <  1:928550/35‑1 (MQ=255)
   tttGCCAGACCGTCGATGGCAAGCGGAATAAACTg                                 >  1:9049378/1‑35 (MQ=255)
   tttGCCAGACCGTCGATGGCAAGCGGAATAAACTg                                 >  1:19229998/1‑35 (MQ=255)
   tttGCCAGACCGTCGATGGCAAGCGGAATAAACTg                                 <  1:13819416/35‑1 (MQ=255)
   tttGCCAGACAGTCGATGGCAAGCGGAATAAACTg                                 <  1:14851853/35‑1 (MQ=255)
   tttCCCGGACCGTCGATGGCAAGCGGAATAAACTg                                 <  1:1306795/35‑1 (MQ=37)
   tgtGCCAGCCCGTCGATGGCAAGCGGAATAAACTg                                 <  1:16911813/33‑1 (MQ=37)
   gttGCCAGCCCGTCGATGGCAAGCGGAATAAACTg                                 <  1:1122332/34‑1 (MQ=255)
    ttGCCAGACCGTCGATGGCAAGCGGAATAAACTGc                                >  1:5591777/1‑35 (MQ=255)
    gccccAGACCGTCGATGGCAAGCGGAATAAACTGc                                <  1:14780985/32‑1 (MQ=255)
      gCCAGACCGTCGATGGCAAGCGGAATAAACTGCcg                              >  1:18577520/1‑35 (MQ=255)
       ccAGATCGTCGATGGCAAGCGGAATAAACTGCcgc                             >  1:6556635/1‑35 (MQ=255)
       ccAGACCGTCGATGGCCAGCGGAATAAACTGCcgc                             >  1:4471251/1‑35 (MQ=255)
       ccAGACCGTCGATGGCAAGCGGAATAAACTGcccc                             >  1:3957450/1‑33 (MQ=255)
       ccAGACCGTCGATGGCAAGCGGAATAAACTGCcgc                             >  1:20247813/1‑35 (MQ=255)
       ccAGACCGTCGATGGCAAGCGGAATAAACTGCcgc                             >  1:5540112/1‑35 (MQ=255)
       ccAGACCGTCGATGGCAAGCGGAATAAACTGCcgc                             >  1:14280321/1‑35 (MQ=255)
        cAGACCGTCGATGGCAAGCGGAATAAACTGCcgcg                            >  1:7215798/1‑35 (MQ=255)
         aGACCGTCGATTGCAAGCGGAATAAACTGCCGCGa                           >  1:7811694/1‑35 (MQ=255)
         aGACCGTCGATGGCAAGCGGAATAAACTGCCGCGa                           >  1:3639517/1‑35 (MQ=255)
         aGACCGTCGATGGCAAGCGGAATAAACTGCCGCGa                           >  1:16265483/1‑35 (MQ=255)
         aGACCGTCGATGGCAAGCGGAATAAACTGCCGCGa                           >  1:17352255/1‑35 (MQ=255)
          gACCGTCGATGGCAAGCGGAATAAACTGCCGCGaa                          >  1:1707731/1‑35 (MQ=255)
           aCCGTCGATGGCAAGCGGAATAAACTGCCGCGAAt                         >  1:8999199/1‑35 (MQ=255)
           aCCGTCGATGGCAAGCGGAATAAACTGCCGCGAAt                         >  1:8335666/1‑35 (MQ=255)
           aCCGTCGATGGCAAGCGGAATAAACTGCCGCGAAt                         >  1:2665439/1‑35 (MQ=255)
           aCCGTCGATGGCAAGCGGAATAAACTGCCGCGAAt                         >  1:14163540/1‑35 (MQ=255)
           aCCGTCGATGGCAAGCGGAATAAACTGCCGCGAAt                         >  1:13197979/1‑35 (MQ=255)
            ccGTCGATGGCAAGCGGAATAAACTGCCGCGaacc                        >  1:756468/1‑33 (MQ=255)
            ccGTCGATGGCAAGCGGAATAAACTGCCGCGAATc                        >  1:18021350/1‑35 (MQ=255)
             cGTCGATGGCAAGCGGAATAAACTGCCGCGAATCg                       >  1:19375553/1‑35 (MQ=255)
             cGTCGATGGCAAGCGGAATAAACTGCCGCGAATCg                       >  1:12926417/1‑35 (MQ=255)
             cGTCGATGGCAAGCGGAATAAACTGCCGCGAATCg                       >  1:13105113/1‑35 (MQ=255)
              gTCGATGGCAAGCGTAATAAACTGCCGCGAATCGt                      >  1:15401877/1‑35 (MQ=255)
              gTCGATGGCAAGCGGAATAAACTGCCGCGAATCGt                      >  1:16729276/1‑35 (MQ=255)
              gTCGATGGCAAGCGGAATAAACTGCCGCGAATCGt                      >  1:8126249/1‑35 (MQ=255)
              gTCGATGGCAAGCGGAATAAACTGCCGCGAATCGt                      >  1:368668/1‑35 (MQ=255)
               tCGATGGCAAGCGGAATAAACTGCCGCGAATCGTa                     >  1:8677402/1‑35 (MQ=255)
                cGATGGCAAGCGGAATAAACTGCCGCGAATCGTaa                    >  1:4408662/1‑35 (MQ=255)
                cGATGGCAAGCGGAATAAACTGCCGCGAATCGTaa                    >  1:6748061/1‑35 (MQ=255)
                 gATGGCAAGCGGAATAAACTGCCGCGAATCGTAAc                   >  1:4269868/1‑35 (MQ=255)
                 gATGGCAAGCGGAATAAACTGCCGCGAATCGTAAc                   >  1:5439271/1‑35 (MQ=255)
                 gATGGCAAGCGGAATAAACTGCCGCGAATCGTAAc                   >  1:17571505/1‑35 (MQ=255)
                 gATGGCAAGCGGAATAAACTGCCGCGAATCGTAAc                   >  1:20066426/1‑35 (MQ=255)
                 gATGGCAAGCGGAATAAACTGCCGCGAATCGTAAc                   >  1:1116624/1‑35 (MQ=255)
                 gATGGCAAGCGGAATAAACTGCCGCGAATCGTAAc                   >  1:19691403/1‑35 (MQ=255)
                 gATGGCAAGCGGAATAAACTGCCGCGAATCGCAAc                   >  1:89140/1‑35 (MQ=255)
                  aTGGCAAGCGGAATAAACTGCCGCGAATCGTAACg                  >  1:6325751/1‑35 (MQ=255)
                  aTGGCAAGCGGAATAAACTGCCGCGAATCGTAACg                  >  1:18824070/1‑35 (MQ=255)
                   tGGCAAGCGGAATAAACTGCCGCGAATCGTAACGt                 >  1:10552635/1‑35 (MQ=255)
                   tGGCAAGCGGAATAAACTGCCGCGAATCGTAACGt                 >  1:8241328/1‑35 (MQ=255)
                   tGGCAAGCGGAATAAACTGCCGCGAATCGTAACGt                 >  1:19935366/1‑35 (MQ=255)
                   tGGCAAGCGGAATAAACTGCCGCAAATCGTAACGt                 >  1:12628044/1‑35 (MQ=255)
                    ggCAAGCGGAATAAACTGCCGCGAATCGTAACGTa                >  1:2150149/1‑35 (MQ=255)
                    ggCAAGCGGAATAAACTGCCGCGAATCGTAACGTa                >  1:3814791/1‑35 (MQ=255)
                      cAAGCGGAATAAACTGCCGCGAATCGTAACGTAGc              >  1:4894318/1‑35 (MQ=255)
                      cAAGCGGAATAAACTGCCGCGAATCGTAACGTAGc              >  1:2887409/1‑35 (MQ=255)
                      cAAGCGAAATAACCTGCCGCGAATCGTAACGTAGc              <  1:20150834/35‑1 (MQ=37)
                       aaGCGGAATAAACTGCCGCGAATCGTAACGTAGCa             >  1:11401636/1‑35 (MQ=255)
                       aaGCGGAATAAACTGCCGCGAATCGTAACGTAGCa             >  1:599905/1‑35 (MQ=255)
                       aaGCGGAATAAACTGCCGCGAATCGTAACGTAGCa             >  1:9518799/1‑35 (MQ=255)
                       aaGCGGAATAAACTGCCGCGAATCGTAACGCAGCa             >  1:20605601/1‑35 (MQ=255)
                        ggCGAAAAAACCTGCCGCGAATCGTAACGTAGCaa            <  1:6880808/34‑1 (MQ=25)
                        aGCGGAATAAACTGCCGCGAATCGTAACGTAGCaa            >  1:240240/1‑35 (MQ=255)
                         gCGGAATAAACTGCCGCGAATCGTAACGTAGCAAc           >  1:2575551/1‑35 (MQ=255)
                         gCGGAATAAACTGCCGCGAATCGTAACGTAGCAAc           >  1:2789434/1‑35 (MQ=255)
                         gCGGAATAAACTGCCGCGAATCGTAACGTAGCAAc           >  1:3502757/1‑35 (MQ=255)
                         gCGGAATAAACTGCCGCGAATCGTAACGTAGCAAc           >  1:13434028/1‑35 (MQ=255)
                         gCGGAATAAACTGCCGCGAATCGTAACGTAGCAAc           >  1:16226569/1‑35 (MQ=255)
                          gggAATCCACTGCCGCGAATCGTAACGTAGCAACg          <  1:12893120/34‑1 (MQ=37)
                          cGGAATAAACTGCCGCGAATCGTAACGTAGCAACg          >  1:1530242/1‑35 (MQ=255)
                          cGGAATAAACTGCCGCGAATCGTAACGTAGCAACg          >  1:16859658/1‑35 (MQ=255)
                           ggAATAAACTGCCGCGAATCGTAACGTAGCAACGc         >  1:18148460/1‑35 (MQ=255)
                           ggAATAAACTGCCGCGAATCGTAACGTAGCAACGc         >  1:11000683/1‑35 (MQ=255)
                           ggAAAAAACCTGCCGCGAATCGTAACGTAGCAACg          <  1:3017536/35‑1 (MQ=25)
                            gAATAAACTGCCGCGAATCGTAACGTAGCAACGCa        >  1:6271301/1‑35 (MQ=255)
                            gAATAAACTGCCGCGAATCGTAACGTAGCAACGCa        >  1:16065127/1‑35 (MQ=255)
                            gAATAAACTGCCGCGAATCGTAACGTAGCAACGCa        >  1:6980780/1‑35 (MQ=255)
                            gAATAAACTGCCGCGAATCGTAACGTAGCAACGCa        >  1:13382827/1‑35 (MQ=255)
                            gAATAAACTGCCGCGAATCGTAACGTAGCAACGCa        >  1:8851867/1‑35 (MQ=255)
                             aaTAAACTGCCGCGAATCGTAACGTAGCAACGCAg       >  1:10103713/1‑35 (MQ=255)
                              aTAAACTGCCGCGAATCGTAACGTAGCAACGCAGa      >  1:7668983/1‑35 (MQ=255)
                               tAACCTGCCGCGAATCGTAACGTAGCAACGCAGAc     <  1:18331639/35‑1 (MQ=255)
                               tAAACTGCCGCGAATCGTAACGTAGCAACGCAGAc     >  1:17677633/1‑35 (MQ=255)
                                aaaCTGCCGCGAATCGTAACGTAGCAACGCAGACg    >  1:13636676/1‑35 (MQ=255)
                                aaaCTGCCGCGAATCGTAACGTAGCAACGCAGACg    >  1:436968/1‑35 (MQ=255)
                                 aaCTGCCGCGAATCGTAACGTAGCAACGCAGACGc   >  1:4813915/1‑35 (MQ=255)
                                 aaCTGCCGCGAATCGTAACGTAGCAACGCAGACGc   >  1:10843835/1‑35 (MQ=255)
                                 aaCTGCCGCGAATCGTAACGTAGCAACGCAGACGc   >  1:5126562/1‑35 (MQ=255)
                                 aaCTGCCGCGAATCGTAACGTAGCAACGCAGACGc   >  1:19880962/1‑35 (MQ=255)
                                   |                                  
TCCTTTGCCAGACCGTCGATGGCAAGCGGAATAAACTGCCGCGAATCGTAACGTAGCAACGCAGACGCA  >  REL606/2726965‑2727033

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