Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,767,8830G→T 6.5% 20.7 ‑30.3 0.0242L106M (CTG→ATG) ygbEconserved inner membrane protein
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/16):  ref base (199/27):  total (199/43)
Fisher's exact test strand distribution p-value = 6.66e-14
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 298 total aligned reads displayed.

GAACAGAGGAGACAAGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCT  >  REL606/2767849‑2767917
                                  |                                  
gAACAGAGGAGACAAGAAAGTACAAATTAGCCCAg                                    >  1:6685169/1‑35 (MQ=255)
gAACAGAGGAGACAAGAAAGTACAAATTAGCCCAg                                    >  1:3850510/1‑35 (MQ=255)
gAACAGAGGAGACAAGAAAGTACAAATTAGCCCAg                                    >  1:8885956/1‑35 (MQ=255)
gAACAGAGGAGACAAGAAAGTACAAATTAGCCCAg                                    >  1:11554805/1‑35 (MQ=255)
 aaCAGAGGAGACAAGAAAGTACAAATTAGCCCAGt                                   >  1:10113505/1‑35 (MQ=255)
 aaCAGAGGAGACAAGAAAGTACAAATTAGCCCAGt                                   >  1:17357875/1‑35 (MQ=255)
 aaCAGAGGAGACAAGAAAGTACAAATTAGCCCAGt                                   >  1:16316467/1‑35 (MQ=255)
  aCAGAGGAGACAAGAAAGTACAAATTAGCCCAGTa                                  >  1:7962821/1‑35 (MQ=255)
  aCAGAGGAGACAAGAAAGTACAAATTAGCCCAGTa                                  <  1:13685673/35‑1 (MQ=255)
  aCAGAGGAGACAAGAAAGTACAAATTAGCCCAGTa                                  >  1:15514304/1‑35 (MQ=255)
  aCAGAGGAGACAAGAAAGTACAAATTAGCCCAGTa                                  >  1:8602726/1‑35 (MQ=255)
  aCAGAGGAGACAAGAAAGTACAAATTAGCCCAGTa                                  >  1:7705485/1‑35 (MQ=255)
   cAGAGGAGACAAGAAAGTACAAATTAGCCCAGTAg                                 >  1:4599382/1‑35 (MQ=255)
   cAGAGGAGACAAGAAAGTACAAATTAGCCCAGTAg                                 >  1:552971/1‑35 (MQ=255)
   cAGAGGAGACAAGAAAGTACAAATTAGCCCAGTAg                                 >  1:4192942/1‑35 (MQ=255)
   cAGAGGAGACAAGAAAGTACAAATTAGCCCAGTAg                                 >  1:1340492/1‑35 (MQ=255)
   cAGAGGAGACAAGAAAGTACAAATTAGCCCAGTAg                                 <  1:2893026/35‑1 (MQ=255)
    agagGAGACAAGAAAGTACAAATTAGCCCAGTAGc                                >  1:580501/1‑35 (MQ=255)
    agagGAGACAAGAAAGTACAAATTAGCCCAGTAGc                                <  1:16500379/35‑1 (MQ=255)
    agagGAGACAAGAAAGTACAAATTAGCCCAGTAGc                                >  1:18999800/1‑35 (MQ=255)
     gaggagACAAGAAAGTACAAATTAGCCCAGTAGcc                               >  1:4469134/1‑35 (MQ=255)
     gaggagACAAGAAAGTACAAATTAGCCCAGTAGcc                               <  1:10956169/35‑1 (MQ=255)
      agtagaCAAGAAAGTACAAATTAGCCCAGTAGCca                              >  1:9254041/4‑35 (MQ=255)
      aggagACAAGAAAGTACAAATTAGCCCAGTAGCca                              >  1:563373/1‑35 (MQ=255)
      aggagACAAGAAAGTACAAATTAGCCCAGTAGCca                              <  1:1951284/35‑1 (MQ=255)
      aggagACAAGAAAGTACAAATTAGCCCAGTAGCca                              <  1:14549755/35‑1 (MQ=255)
      aggagACAAGAAAGTACAAATTAGACCAGTAGCca                              <  1:11281854/35‑1 (MQ=255)
        gagaCAAGAAAGTACAAATTAGCCCAGTAGCcaca                            >  1:19932222/1‑35 (MQ=255)
        gagaCAAGAAAGTACAAATTAGCCCAGTAGCcaca                            >  1:18018284/1‑35 (MQ=255)
         agaCAAGAAAGTACAAATTAGCCCAGTAGCCACAt                           >  1:7488390/1‑35 (MQ=255)
         agaCAAGAAAGTACAAATTAGCCCAGTAGCCACAt                           >  1:983755/1‑35 (MQ=255)
         agaCAAGAAAGTACAAATTAGCCCAGTAGCCACAt                           <  1:18152709/35‑1 (MQ=255)
         agaCAAGAAAGTACAAATTAGCCCAGTAGCCACAt                           >  1:1058544/1‑35 (MQ=255)
          gaCAAGAAAGTACAAATTAGCCCAGTATCCACATa                          >  1:18698736/1‑35 (MQ=255)
          gaCAAGAAAGTACAAATTAGCCCAGTAGCCACATa                          >  1:4512041/1‑35 (MQ=255)
           aCAAGAAAGTACAAATTAGCCCAGTAGCCACATaa                         >  1:16682748/1‑35 (MQ=255)
            cAAGAAAGTACAAATTAGCCCAGTAGCCACATaaa                        >  1:19188399/1‑35 (MQ=255)
             aaGAAAGTACAAATTAGCCCAGTAGCCACATAAAc                       >  1:12687806/1‑35 (MQ=255)
              aGAAAGTACAAATTAGCCCAGTAGCCACATAAACa                      >  1:3300224/1‑35 (MQ=255)
              aGAAAGTACAAATTAGCCCAGTAGCCACATAAACa                      >  1:17172202/1‑35 (MQ=255)
              aGAAAGTACAAATTAGCCCAGTAGCCACATAAACa                      >  1:11765550/1‑35 (MQ=255)
               gAAAGTACAAATTAGCCCAGTAGCCACATAAACAg                     >  1:15807045/1‑35 (MQ=255)
                aaaGTACAAATTAGCCCAGTAGCCACATAAACAGt                    >  1:17877530/1‑35 (MQ=255)
                aaaGTACAAATTAGCCCAGTAGCCACATAAACAGt                    >  1:12004815/1‑35 (MQ=255)
                aaaGTACAAATTAGCCCAGTAGCCACATAAACAGt                    >  1:6569415/1‑35 (MQ=255)
                 aaGTACAAATTAGCCCAGTAGCCACATAAACATTg                   >  1:8378580/1‑35 (MQ=255)
                 aaGTACAAATTAGCCCAGTAGCCACATAAACATTg                   >  1:7819488/1‑35 (MQ=255)
                 aaGTACAAATTAGCCCAGTAGCCACATAAACAGTg                   >  1:3493029/1‑35 (MQ=255)
                  aGTACAAATTAGCCCAGTAGCCACATAAACAGTgc                  >  1:16053606/1‑35 (MQ=255)
                  aGTACAAATTAGCCCAGTAGCCACATAAACAGTgc                  >  1:10473074/1‑35 (MQ=255)
                  aGTACAAATTAGCCCAGTAGCCACATAAACAGTgc                  >  1:5427321/1‑35 (MQ=255)
                   gTACAAATTAGCCCAGTAGCCACATAAACAGTgcg                 >  1:1147387/1‑35 (MQ=255)
                    taaaaaTTACCCCTTTAGCCACATAAACAGTgcgc                <  1:18475089/32‑1 (MQ=16)
                    tACAAATTACCCCATTAGCCACATAAACAGTgcgc                <  1:11202975/35‑1 (MQ=37)
                     aCAAATTAGCCCAGTAGCCACATAAACAGTGCGcc               >  1:9875537/1‑35 (MQ=255)
                     aCAAATTAGCCCAGTAGCCACATAAACAGTGCGcc               >  1:13519170/1‑35 (MQ=255)
                       aaaTTAGCCTGGTAGCCACATAAACAGTGCGCCaa             <  1:5943052/35‑1 (MQ=37)
                       aaaTTAGCCCAGTAGCCACATAAACAGTGCGCCaa             >  1:5094951/1‑35 (MQ=255)
                       aaaTTAGCCCAGTAGCCACATAAACAGTGCGCCaa             >  1:17716308/1‑35 (MQ=255)
                       aaaTTACCCCATTAGCCACATAAACAGTGCGCCaa             <  1:10040698/35‑1 (MQ=37)
                       aaaTTACCCCATTAGCCACATAAACAGTGCGCCaa             <  1:1032273/35‑1 (MQ=37)
                        aaTTAGCCCATTAGCCACATAAACAGTGCGCCaaa            <  1:16788812/35‑1 (MQ=255)
                        aaTTAGCCCAGTAGCCACATAAACAGTGCGCCaaa            >  1:16078584/1‑35 (MQ=255)
                         aTTAGCCCATTAGCCACATAAACAGTGCGCCAAAc           <  1:4103995/35‑1 (MQ=255)
                         aTTAGCCCAGTAGCCACATAAACAGTGCGCCAAAc           >  1:1313836/1‑35 (MQ=255)
                         aTTAGCCAAGTAGCCACATAAACAGTGCGCCAAAc           <  1:13369509/35‑1 (MQ=255)
                         aTTACCCCATTAGCCACATAAACAGTGCGCCAAAc           <  1:17490940/35‑1 (MQ=37)
                         aTTACCCCAGTAGCCACATAAACAGTGCGCCAAAc           <  1:18210587/35‑1 (MQ=255)
                          ttAGCCCATTAGCCACATAAACAGTGCGCCAAACa          <  1:12772930/35‑1 (MQ=255)
                          ttAGCCCAGTAGCCACATAAACAGTGCGCCAAACa          >  1:9540316/1‑35 (MQ=255)
                          ttAGCCCAGTAGCCACATAAACAGTGCGCCAAACa          <  1:10689521/35‑1 (MQ=255)
                           tAGCCCATTAGCCACATAAACAGTGCGCCAAACAt         <  1:2389524/35‑1 (MQ=255)
                           tAGCCCAGTAGCCACATAAACAGTGCGCCAAACAt         >  1:14255330/1‑35 (MQ=255)
                           aaGCCCATTAGCCACATAAACAGTGAGCCAAACAt         <  1:1472984/34‑1 (MQ=37)
                            aGCCCAGTAGCCACATAAACAGTGCGCCAAACATa        >  1:6902303/1‑35 (MQ=255)
                            aGCCCAGTAGCCACATAAACAGTGCGCCAAACATa        >  1:4916155/1‑35 (MQ=255)
                            aGCCCAGTAGCCACATAAACAGTGCGCCAAACATa        >  1:20501943/1‑35 (MQ=255)
                             gCCCAGTAGCCACATAAACAGTGCGCCAAACATaa       >  1:3374377/1‑35 (MQ=255)
                             ccccATTAGCCACATAAACAGTGCGCCAAACATaa       <  1:13090568/34‑1 (MQ=255)
                              cccATTAGCCACATAAACAGTGCGCCAAACATAAt      <  1:4399644/35‑1 (MQ=255)
                              cccAGTAGCCACATAAACAGTGCGCCAAACATAAt      >  1:6976202/1‑35 (MQ=255)
                              cccAGTAGCCACATAAACAGTGCGCCAAACATAAt      >  1:12256571/1‑35 (MQ=255)
                              cccAGTAGCCACATAAACAGTGCGCCAAACATAAt      >  1:16571716/1‑35 (MQ=255)
                              cccAGTAGCCAAATAAACAGTGTGCCAAACATAAt      <  1:1850144/35‑1 (MQ=37)
                              cccAGTAGCAACATAAACAGTGCGCCAAACATAAt      <  1:7360432/35‑1 (MQ=255)
                               ccATTAGCCACATAAACAGTGCGCCAAACATAATg     <  1:15244096/35‑1 (MQ=255)
                               ccAGTAGCCACATAAACAGTGCGCCAAACATAATg     >  1:8096664/1‑35 (MQ=255)
                               ccAGTAGCCACATAAACAGTGCGCCAAACATAATg     >  1:2088564/1‑35 (MQ=255)
                               ccAGTAGCCACATAAACAGTGCGCCAAACATAATg     >  1:18832731/1‑35 (MQ=255)
                               ccAGTAGCCACATAAACAGTGCGCCAAACATAATg     >  1:309486/1‑35 (MQ=255)
                               ccAGTAGCCACATAAACAGTGCGCAAACATAATGc    <  1:10025851/35‑1 (MQ=37)
                                cAGTAGCCACATAAACAGTGCGCCAAACATAATGc    >  1:14623231/1‑35 (MQ=255)
                                cAGTAGCCACATAAACAGTGCGCCAAACATAATGc    >  1:2666476/1‑35 (MQ=255)
                                cAGTAGCCACATAAACAGTGCGCCAAACATAATGc    >  1:19599034/1‑35 (MQ=255)
                                aaGTAGCCACATAAACAGTGCGCCAAACATAATGc    <  1:14463687/34‑1 (MQ=255)
                                 aGTAGCCACATAAACAGTGCGCCAAACATAATGcc   >  1:3944413/1‑35 (MQ=255)
                                  gTAGCCACATAAACAGTGCGCCAAACATAATGCCt  >  1:625763/1‑35 (MQ=255)
                                  gTAGCCACATAAACAGTGCGCCAAACATAATGCCt  >  1:15083888/1‑35 (MQ=255)
                                  gTAGCCACATAAACAGTGCGCCAAACATAATGCCt  >  1:15311097/1‑35 (MQ=255)
                                  gTACCCACATAAACAGTGCGCCAAACATAATGCCt  >  1:20102210/1‑35 (MQ=255)
                                  |                                  
GAACAGAGGAGACAAGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCT  >  REL606/2767849‑2767917

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