Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,757,2130A→T 4.9% 4.1 ‑17.3 ‑0.0241I7F (ATC→TTC) ygbKhypothetical protein
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Rejected: Polymorphic base substitution creates a homopolymer stretch.
Reads supporting (aligned to +/- strand):  new base (0/13):  ref base (173/55):  total (173/68)
Fisher's exact test strand distribution p-value = 2.94e-08
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 271 total aligned reads displayed.

TACCGGGAGCGAAATCATGATCAAGATTGGCGTTATCGCCGATGATTTTACCGGCGCGACGGATATCGC  >  REL606/2757179‑2757247
                                  |                                  
tACCGGGAGCGAAATCATGATCAAGATTGGCGTTa                                    <  1:2368582/35‑1 (MQ=255)
 aCCGGGAGCGAAATCATGATCAAGATTGGCGTTAt                                   >  1:7029158/1‑35 (MQ=255)
 aCCGGGAGCGAAATCATGATCAAGATTGGCGTTAt                                   <  1:1257009/35‑1 (MQ=255)
 aCCGGGAGCGAAATCATGATCAAGATTGGCGTTAt                                   >  1:2696640/1‑35 (MQ=255)
 aCCGGGAGCGAAATCATGATCAAGATTGGCGTTAt                                   <  1:10470521/35‑1 (MQ=255)
  ccGGGAGCGGAATCATGATCAAGATTGGCGTTATc                                  <  1:15330558/35‑1 (MQ=255)
  ccGGGAGCGAAATCATGATCAAGATTGGCGTTATc                                  >  1:16777833/1‑35 (MQ=255)
  ccGGGAGCGAAATCATGATCAAGATTGGCGTTATc                                  <  1:4049962/35‑1 (MQ=255)
  ccGGGAGCGAAATCATGATCAAGATTGGCGTTATc                                  >  1:3300480/1‑35 (MQ=255)
   cGGGAGCGAAATCATGATCAAGATTGGCGTTATCg                                 <  1:11457460/35‑1 (MQ=255)
   cGGGAGCGAAATCATGATCAAGATTGGCGTTATCg                                 >  1:1334919/1‑35 (MQ=255)
    gggAGCGAAATCATGATCAAGATTGGCGTTATCGc                                <  1:14641673/35‑1 (MQ=255)
    gggAGCGAAATCATGATCAAGATTGGCGTTATCGc                                >  1:14724872/1‑35 (MQ=255)
    gggAGCGAAATCATGATCAAGATTGGCGTTATCGc                                <  1:11794328/35‑1 (MQ=255)
    gggAGCGAAATCATGATCAAGATTGGCGTTATCGc                                <  1:13862240/35‑1 (MQ=255)
    gggAGCGAAATCATGATCAAGATTGGCGTTATCGc                                >  1:16128305/1‑35 (MQ=255)
     ggAGCGAAATCATGATCAAGATTGGCGTTATCGcc                               >  1:1589769/1‑35 (MQ=255)
      gAGCGAAATCATGATCAAGATTGGCGTTATCGCCg                              >  1:870311/1‑35 (MQ=255)
      gAGCGAAATCATGATCAAGATTGGCGTTATCGCCg                              >  1:6949836/1‑35 (MQ=255)
      gAGCGAAATCATGATCAAGATTGGCGTTATCGCCg                              <  1:19384309/35‑1 (MQ=255)
       aGCGAAATCATGATCAAGATTGGCGTTATCGCCga                             >  1:12109353/1‑35 (MQ=255)
       aGCGAAATCATGATCAAGATTGGCGTTATCGCCga                             <  1:11014811/35‑1 (MQ=255)
       aGCGAAATCATGATCAAGATTGGCGTTATCGCCga                             >  1:5851663/1‑35 (MQ=255)
        gCGAAATCATGATCAAGATTGGCGTTATCGCCgat                            >  1:5662145/1‑35 (MQ=255)
        gCGAAATCATGATCAAGATTGGCGTTATCGCCgat                            >  1:20160110/1‑35 (MQ=255)
        gCGAAATCATGATCAAGATTGGCGTTATCGCCgat                            >  1:9678387/1‑35 (MQ=255)
        gCGAAATCATGATCAAGATTGGCGTTATCGCCgat                            >  1:218511/1‑35 (MQ=255)
         cGAAATCATGATCAAGATTGGCGTTATCGCCgatg                           >  1:18286697/1‑35 (MQ=255)
          gAAATCATGATCAAGATTGGCGTTATCGCCgatga                          >  1:2584322/1‑35 (MQ=255)
          gAAATCATGATCAAGATTGGCGTTATCGCCgatga                          >  1:10866409/1‑35 (MQ=255)
          gAAATCATGATCAAGATTGGCGTTATCGCCgatga                          >  1:1036683/1‑35 (MQ=255)
           aaaTCATGATCAAGATTGGCGTTATCGCCgatgat                         >  1:20380232/1‑35 (MQ=255)
           aaaTCATGATCAAGATTGGCGTTATCGCCgatgat                         >  1:137844/1‑35 (MQ=255)
           aaaTCATGATCAAGATTGGCGTTATCGCCgatgat                         >  1:2832132/1‑35 (MQ=255)
            aaTCATGATCAAGATTGGCGTTATCGCCGATGAtt                        <  1:1672349/35‑1 (MQ=255)
             aTCATGATCAAGATTGGCGTTATCGCCGATGAttt                       >  1:17875081/1‑35 (MQ=255)
              tCATGATCAAGATTGGCGTTATCGCCGATGAtttt                      >  1:5358952/1‑35 (MQ=255)
              tCATGATCAAGATTGGCGTTATCGCCGATGAtttt                      >  1:1485042/1‑35 (MQ=255)
              tCATGATCAAGATTGGCGTTATCGCCGATGAtttt                      >  1:10214416/1‑35 (MQ=255)
               cATGATCAAGATTGGCGTTATCGCCGATGATTTTa                     >  1:1715190/1‑35 (MQ=255)
                aTGATCAAGATTGGCGTTATCGCCGATGATTTTAc                    >  1:9011892/1‑35 (MQ=255)
                 tGATCAAGATTGGCGTTATCGCCGATGATTTTAcc                   >  1:5093565/1‑35 (MQ=255)
                 tGATCAAGATTGGCGTTATCGCCGATGATTTTAcc                   >  1:3199695/1‑35 (MQ=255)
                  gATCAAGATTGGCTTTTTCGCCGATGATTTTACCg                  <  1:1491561/35‑1 (MQ=37)
                  gATCAAGATTGGCGTTATCGCCGATGATTTTACCg                  >  1:13656541/1‑35 (MQ=255)
                  gATCAAGATTGGCGTTATCGCCGATGATTTTACCg                  <  1:12308987/35‑1 (MQ=255)
                  gATCAAGATTGGCGTTATCGCCGATGATTTTACCg                  <  1:19547427/35‑1 (MQ=255)
                   ataaaGATTGGCTTTTTCGCCGATGATTTTACCgg                 <  1:14464233/32‑1 (MQ=25)
                   agcAAGATTGGCTTTTTCGCCGATGATTTTACCgg                 <  1:8182469/33‑1 (MQ=25)
                   aTCAAGATTGGCGTTATCGCCGATGATTTTACCgg                 >  1:8210529/1‑35 (MQ=255)
                    tCAAGATTGGCGTTATCGCCGATGATTTCACCGgc                >  1:11842855/1‑35 (MQ=255)
                     cAAGATTGGCGTTATCGCCGATGATTTTACCGgcg               >  1:5581738/1‑35 (MQ=255)
                     cAAGATTGGCGTTACCGCCGATGATTTTACCGgcg               >  1:7612663/1‑35 (MQ=255)
                      aaGATTGGGTTTATCGCCGATGATTTTACCGgcgc              <  1:17262105/35‑1 (MQ=25)
                      aaGATTGGCGTTATCGCCGATGATTTTACCGgcgc              >  1:3401703/1‑35 (MQ=255)
                      aaGATTGGCGTTATCGCCGATGATTTTACCGgcgc              >  1:5182908/1‑35 (MQ=255)
                      aaGATGGGCTTTTTCGCCGATGATTTTACCGgcgc              <  1:10723828/35‑1 (MQ=25)
                        gATTGGCGTTATCGCCGATGATTTTACCGGCGCGa            >  1:7741093/1‑35 (MQ=255)
                        gATTGGCGTTATCGCCGATGATTTTACCGGCGCGa            >  1:8473316/1‑35 (MQ=255)
                        gATTGGCGTTATCGCCGATGATTTTACCGGCGCGa            >  1:951077/1‑35 (MQ=255)
                        gATGGGCTTTTTCGCCGATGATTTTACCGGCGCGa            <  1:3584314/35‑1 (MQ=25)
                         aTTGTCGTTATCGCCGATGATTTTACCGGCGCGAc           >  1:6886694/1‑35 (MQ=255)
                         aTTGGCGTTATCGCCGATGATTTTACCGGCGCGAc           >  1:5170128/1‑35 (MQ=255)
                         aTTGGCGTTATCGCCGATGATTTTACCGGCGCGAc           >  1:4811432/1‑35 (MQ=255)
                         aTTGGCGTTATCGCCGATGATTTTACCGGCGCGAc           >  1:4852836/1‑35 (MQ=255)
                         aTTGGCGTTATCGCCGATGATTTTACCGGCGCGAc           >  1:9427089/1‑35 (MQ=255)
                          ttGGCTTTATCGCCGATGATTTTACCGGCGCGACg          <  1:15991060/35‑1 (MQ=255)
                          ttGGCGTTATCGCCGATGATTTTACCGGCGCGACg          >  1:20035782/1‑35 (MQ=255)
                          tgggCTTTTTCGCCGATGATTTTACCGGCGCGACg          <  1:2416009/33‑1 (MQ=25)
                          tgggCGTTATCGCCGATGATTTTACCGGCGCGACg          <  1:4592624/33‑1 (MQ=255)
                           tGGCGTTATCGCCGATGATTTTACCGGCGCGACgg         >  1:11605629/1‑35 (MQ=255)
                           gggCGTTATCGCCGATGATTTTACCGGCGCGACgg         <  1:15094285/34‑1 (MQ=255)
                            ggCTTTTTCGCCGATGATTTTACCGGCGCGACGGa        <  1:10023762/35‑1 (MQ=37)
                            ggCGTTTTCGCCGATGATTTTACCGGCGCGACGGa        <  1:8633056/35‑1 (MQ=255)
                            ggCGTTATCGCCGATGATTTTACCGGCGCGACGGa        >  1:19871086/1‑35 (MQ=255)
                             gctttATCGCCGATGATTTTACCGGCGCGACGGat       <  1:9897551/32‑1 (MQ=255)
                             gCGTTATCGCCGATGATTTTACCGGCGCGACGGat       >  1:6222833/1‑35 (MQ=255)
                             gCGTTATCGCCGATGATTTTACCGGCGCGACGGat       >  1:4406800/1‑35 (MQ=255)
                              ctttATCGCCGATGATTTTACCGGCGCGACGGata      <  1:8586011/33‑1 (MQ=255)
                              ctttATCGCCGATGATTTTACCGGCGCGACGGata      <  1:17522879/33‑1 (MQ=255)
                              cGTTATCGCCGATGATTTTACCGGCGCGACGGata      >  1:15783579/1‑35 (MQ=255)
                              cGTTATCGCCGATGATTTTACCGGCGCGACGGata      >  1:796661/1‑35 (MQ=255)
                              cGTTATCGCCGATGATTTTACCGGCGCGACGGata      >  1:3962649/1‑35 (MQ=255)
                               gTTTTCGCCGATGATTTTACCGGCGCGACGGatat     <  1:13972007/35‑1 (MQ=255)
                               gTTTTCGCCGATGATTTTACCGGCGCGACGGatat     <  1:14325449/35‑1 (MQ=255)
                               gTTATCGCCGATGATTTTACCGGCGCGACGGatat     >  1:10153589/1‑35 (MQ=255)
                               gTTATCGCCGATGATTTTACCGGCGCGACGGatat     >  1:7301625/1‑35 (MQ=255)
                               gTTATCGCCGATGATTTTACCGGCGCGACGGatat     >  1:20038793/1‑35 (MQ=255)
                                ttATCGCCGATGATTTTACCGGCGCGACGGATATc    >  1:19161305/1‑35 (MQ=255)
                                ttATCGCCGATGATTTTACCGGCGCGACGGATATc    >  1:12658618/1‑35 (MQ=255)
                                ttATCGCCGATGATTTTACCGGCGCGACGGATATc    >  1:12962321/1‑35 (MQ=255)
                                ttATCGCCGATGATTTTACCGGCGCGACGGATATc    >  1:20296852/1‑35 (MQ=255)
                                 tATCGCCGATGATTTTACCGGCGCGACGGATATCg   >  1:13542792/1‑35 (MQ=255)
                                 tATCGCCGATGATTTTACCGGCGCGACGGATATCg   >  1:14979972/1‑35 (MQ=255)
                                 tATCGCCGATGATTTTACCGGCGCGACGGATATCg   >  1:16917316/1‑35 (MQ=255)
                                 tATCGCCGATGATTTTACCGGCGCGACGGATATCg   >  1:18057810/1‑35 (MQ=255)
                                 tATCGCCGATGATTTTACCGGCGCGACGAATATCg   <  1:1890960/35‑1 (MQ=255)
                                  aTCGCCGATGATTTTACCGGCGCGACGGATATCGc  >  1:16442457/1‑35 (MQ=255)
                                  aTCGCCGATGATTTTACCGGCGCGACGGATATCGc  <  1:305223/35‑1 (MQ=255)
                                  aTCGCCGATGATTTTACCGGCGCGACGGATATCGc  >  1:20577907/1‑35 (MQ=255)
                                  |                                  
TACCGGGAGCGAAATCATGATCAAGATTGGCGTTATCGCCGATGATTTTACCGGCGCGACGGATATCGC  >  REL606/2757179‑2757247

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