Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606931,3140TC40.3% 117.6 / 59.2 119E84G (GAG→GGG) ybjEpredicted transporter
Reads supporting (aligned to +/- strand):  ref base T (35/36);  new base C (1/47);  total (36/83)
Fisher's exact test for biased strand distribution p-value = 3.15e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 225 total aligned reads displayed.

GCGACTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCA  >  REL606/931240‑931389
                                                                          |                                                                           
 cGACTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTcc                                                                           <  2:1176590/76‑1 (MQ=255)
  gACTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCTCCCTCCa                                                                          >  2:5156371/1‑76 (MQ=255)
  gACTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCa                                                                          <  1:7193866/76‑1 (MQ=255)
  gACTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCCCGCTccc                                                                          >  1:1049261/1‑75 (MQ=255)
   accccAGCGCCATCGCAATACGCGACGGGAGTTGTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAg                                                                         <  2:1894098/73‑1 (MQ=255)
   aCTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCGGCGCTcccg                                                                         >  2:7320805/1‑74 (MQ=255)
     tCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCTCGCTCCAGcc                                                                       >  1:2539764/1‑76 (MQ=255)
       cAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCCGCcac                                                                     >  1:238815/1‑76 (MQ=255)
       cAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCcac                                                                     <  1:2979132/76‑1 (MQ=255)
        aGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCTCCCTCCCTCcaca                                                                    >  2:1527590/1‑76 (MQ=255)
         gcgcCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACAt                                                                   <  2:4950445/76‑1 (MQ=255)
          ccccATCGCAATGGGCGACGGGAGTTTTTCTTGCTGAGGGTGGTTGCCCCACGGCAGGCCTCCCCCCAGCCACATc                                                                  <  2:696363/74‑1 (MQ=255)
           gcCATCGCAATCCGCGACGGGATTTTTTTTTGCTGATGGTGGTTGCGCCACGGCAGGCCCCGCCCCAGCCACATCa                                                                 <  2:6144443/76‑1 (MQ=255)
           gcCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCCCGCTCCCGCCACATCa                                                                 >  1:6632532/1‑76 (MQ=255)
           gcCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCa                                                                 <  1:3997938/76‑1 (MQ=255)
           gcCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCCCGCTCCCCCCCCATCa                                                                 >  1:593504/1‑76 (MQ=255)
            ccATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCTCGCTCCACCCACCTCAg                                                                >  1:6687558/1‑76 (MQ=255)
                  cAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCCCGCCCCAGCCACCTCAAggcggc                                                          >  2:5975571/1‑76 (MQ=255)
                   aaTGCGCGACGGTAGTTTTTCTGGCTGGTGGTGTGTGCTCCACGGCGGCCCCCGCCCCAGCCACATCAGGGCGGCa                                                         <  2:6836658/76‑1 (MQ=255)
                   aaTACGCGACGGGATTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCa                                                         <  2:846596/76‑1 (MQ=255)
                   aaTACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCTCGCTCCAGCCACATCAGGGGGGCa                                                         >  2:2711735/1‑76 (MQ=255)
                   aaTACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCCCGCTCCCGCCACATCCGggcgggc                                                         >  1:3097994/1‑74 (MQ=255)
                     tccgcgCCGGGATTTTTTCTTGCGGATGGTGTTGGCCCCAGGGCGGGCCTCGCCCCAGCCACATCAGGGCGGCAat                                                       <  2:1284172/74‑1 (MQ=255)
                      tcgcgGCGGGAGTTTTTCTTGCTGATGGTGGTTCCGCCACGGCAGGCCTCCCCCCAGCCACATCAGGGCGGCAata                                                      <  1:6011139/75‑1 (MQ=255)
                      aCGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCCAATCACGGCCGCactt                                                      >  2:3331951/1‑73 (MQ=255)
                      aCGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCCCGCTCCCGCCACATCAGGGCGGGAata                                                      >  1:6494298/1‑76 (MQ=255)
                        gcgACGGGATTTTTTCTTGCGGATGGTGGTTGCCCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATAtt                                                    <  1:7125911/76‑1 (MQ=255)
                        gcgACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCCCGGCAGGCCCCGCCCCAGCCACATCAGGGCGGCAATAtt                                                    <  2:3122348/76‑1 (MQ=255)
                        gcgACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATAtt                                                    <  1:1833159/76‑1 (MQ=255)
                        gcgACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGCCCTCGCCCCAGCCACATCAGGGCGGCAATAtt                                                    <  1:5811711/76‑1 (MQ=255)
                         cgAGGGGATTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTa                                                   <  2:2924356/76‑1 (MQ=255)
                         cgACGGTATTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTa                                                   <  2:7397068/76‑1 (MQ=255)
                         cgACGGTAGTTTTTCTTGCTGATGGTGGTGGCGCCACGGCAGGCCTCCCCCCAGCCACATCAGGGCGGCAATATTa                                                   <  2:6524159/76‑1 (MQ=255)
                         cgACGGGAGTTTTTCTTGCTGATGGGGGTTCCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTa                                                   <  2:3785262/76‑1 (MQ=255)
                          ggcGGGAGTTTTTCTTGCTGATGGTGGTTCCCCCGCGGCGCGCCTCGCTCCAGCCACATCAGGGCGGCAATATTac                                                  <  2:4811701/74‑1 (MQ=255)
                            cGGGATTTTTTCTGGCTGATGTGGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTacac                                                <  1:4345009/76‑1 (MQ=255)
                            cGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGCCCTCGCTCCAGCCACATCAGGGCGGCAATATTacac                                                <  2:1610607/76‑1 (MQ=255)
                             gggTGTTTTTCTTGCTGATGGTGGTTGCCCCACGGCAGGCCGCCCTCCAGCCACATCAGGGCGGCAATATTAcaca                                               <  2:5349536/76‑1 (MQ=255)
                              ggATTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAg                                              <  2:4622820/76‑1 (MQ=255)
                               gAGTTTTTCTTGCTGATGGTGGTTGCCCCCCGGCAGCCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGt                                             <  1:6838149/76‑1 (MQ=255)
                                attttttCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTa                                            <  1:6889053/74‑1 (MQ=255)
                                aGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCCCGCTCCAGCCCCATCAGGGCGGGAATATTACACAGTa                                            >  1:7327280/1‑76 (MQ=255)
                                  tttttCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTaaa                                          <  1:1093331/76‑1 (MQ=255)
                                    tttCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCCGCCACATCCAGGCGGCGATATTACACAGTaaaaa                                        >  1:5529488/1‑75 (MQ=255)
                                     ttCTTGCTGATGGTGGTGGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATa                                       <  2:4679105/76‑1 (MQ=255)
                                     ttCTTGCGGAGGGTGGTTGCGCCACGGCGGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATa                                       <  2:7615238/76‑1 (MQ=255)
                                      tCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATaa                                      <  1:7440564/76‑1 (MQ=255)
                                      tCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATaa                                      <  2:5684241/76‑1 (MQ=255)
                                      tCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCCGCCACAACCGGGCGGCAATACTACACCGTAaaaaaa                                      >  2:5515417/1‑76 (MQ=255)
                                       tttGCTGATGGTGGTTGCGCCGCGGCGGGCCTCCCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAAc                                     <  2:1040066/75‑1 (MQ=255)
                                         tGCTGTTGGTGGTGGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACgg                                   <  2:6062392/76‑1 (MQ=255)
                                           cTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTa                                 >  1:7617173/1‑76 (MQ=255)
                                           cTGATGGTGGTTCCCCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTa                                 <  2:5142885/76‑1 (MQ=255)
                                           cTGATGGTGGTGGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTa                                 <  2:339503/76‑1 (MQ=255)
                                            tgttggtggTTGCCCCCCGGCAGGCCGCTCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTaa                                <  1:5560364/73‑1 (MQ=255)
                                            tGATGGTGGTTGCGCCACGGCAGGGCTCTCTCCCGCCCCATCCAGGCGGCAATATTAAACACTAAAATAAACGTaa                                >  2:2301665/1‑76 (MQ=255)
                                            tGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTaa                                <  2:5207182/76‑1 (MQ=255)
                                            tGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTaa                                <  2:3011949/76‑1 (MQ=255)
                                            tGATGGGGGCTCCGCCACGGCAGCCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTaa                                <  2:6363281/76‑1 (MQ=255)
                                             gagggtggtGGCCCCACGGCAGCCCTCCCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAAt                               <  1:137124/73‑1 (MQ=255)
                                              gtggtggTTGCCCCACGGCGGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATa                              <  1:3382060/75‑1 (MQ=255)
                                              aTGGTGGTTGCGCCCCGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATa                              <  2:837465/76‑1 (MQ=255)
                                              aTGGTGGTTGCGCCACGGCAGGGCTCGCTCCCGCCACATCAGGGCGGCCATATTACACAGTAAAATAAAGGGAATa                              >  2:1972367/1‑76 (MQ=255)
                                              aTGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCCATATTACACAGTAAAAAAAAGGGAATa                              >  2:4158932/1‑76 (MQ=255)
                                                ggtggtTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACt                            <  1:185930/76‑1 (MQ=255)
                                                gggggtTGCGCCACCGCAGGCCCCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACt                            <  1:6034174/73‑1 (MQ=255)
                                                 tgtgtTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTg                           <  2:1783498/73‑1 (MQ=255)
                                                 gtggtTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTg                           <  1:804049/76‑1 (MQ=255)
                                                 gtggtTGCGCCACGGCAGGCCCCGCTCCAACCACATCCAGGCGGCACTATTTAACAGCAAAATAACGGGAAAACTg                           >  2:6727587/1‑76 (MQ=255)
                                                  tggtggCGCCACGGCAGGCCTCCCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGa                          <  1:4870860/76‑1 (MQ=255)
                                                  tggtTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGa                          >  2:4466172/1‑76 (MQ=255)
                                                  gggtTGCGCCGCGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGa                          <  1:5425294/75‑1 (MQ=255)
                                                      tGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGc                      >  1:2921135/1‑76 (MQ=255)
                                                      tGCCCCACGGCAGCCCTCCCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGc                      <  1:5279889/76‑1 (MQ=255)
                                                        cgccgcGGCGGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAg                    <  2:4151001/76‑1 (MQ=255)
                                                        cgcCACGGCAGTCCTGGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAg                    <  2:4023470/76‑1 (MQ=255)
                                                        cgcCACGGCAGGGCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATAATGACGGCAg                    >  2:2580430/1‑76 (MQ=255)
                                                         gcCCCGGGAGCCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGa                   <  1:6066674/76‑1 (MQ=255)
                                                         gcCAGGGCAGACCCCCCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGa                   <  1:5792203/76‑1 (MQ=255)
                                                         gcCACGGCAGGGCTCCCTCCAGCCACACCCGGGCGGCAATATTACACACTAAAATAACGGGAATACTTATGGCAGa                   >  2:4242638/1‑76 (MQ=255)
                                                           cccGGCCGACCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaat                 <  1:4490346/74‑1 (MQ=255)
                                                           caaggCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaat                 <  1:1942653/73‑1 (MQ=255)
                                                            cgggCAGCCCCCCCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaata                <  2:2128792/74‑1 (MQ=255)
                                                             cGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaaaaa               >  2:5485894/1‑76 (MQ=255)
                                                             cGGCAGCCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaataa               <  1:3471483/76‑1 (MQ=255)
                                                                gaGGCCCCCCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGc            <  1:3686416/75‑1 (MQ=255)
                                                                cAGGCCTCACTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGc            <  2:5757566/76‑1 (MQ=255)
                                                                cAGCCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGc            <  1:2606918/76‑1 (MQ=255)
                                                                 aGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCa           >  1:4134842/1‑76 (MQ=255)
                                                                 aGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCa           >  1:3878709/1‑76 (MQ=255)
                                                                 aGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCa           >  2:7430084/1‑76 (MQ=255)
                                                                 aGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCa           >  1:5149217/1‑76 (MQ=255)
                                                                 aGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAAAATAATgaa           >  2:1145100/1‑74 (MQ=255)
                                                                 aGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCa           <  1:2142760/76‑1 (MQ=255)
                                                                    ccTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGaa        >  1:851668/1‑76 (MQ=255)
                                                                     cTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAAt       >  2:1363393/1‑76 (MQ=255)
                                                                      ccGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATc      <  1:4778764/75‑1 (MQ=255)
                                                                       cGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAAtag     <  2:6546805/76‑3 (MQ=255)
                                                                         cTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGcc   >  2:3811509/1‑76 (MQ=255)
                                                                          tCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCa  >  1:1581803/1‑76 (MQ=255)
                                                                          |                                                                           
GCGACTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCA  >  REL606/931240‑931389

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: