Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606931,3140TC39.4% 112.0 / 50.0 104E84G (GAG→GGG) ybjEpredicted transporter
Reads supporting (aligned to +/- strand):  ref base T (37/26);  new base C (1/40);  total (38/66)
Fisher's exact test for biased strand distribution p-value = 7.68e-10
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 222 total aligned reads displayed.

AGCGACTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCA  >  REL606/931239‑931389
                                                                           |                                                                            
aGCGACTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCt                                                                              >  2:2056238/1‑76 (MQ=255)
  ggACTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTcc                                                                            <  1:724259/75‑1 (MQ=255)
  cGACTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCTCGCTcc                                                                            >  2:2827373/1‑76 (MQ=255)
  cGACTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTcc                                                                            <  2:6738814/76‑1 (MQ=255)
   gACTCCAGCGCCATCGCAATATGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCCCGCTccc                                                                           >  2:2885939/1‑75 (MQ=255)
   gACTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGctccctccc                                                                           >  2:1342538/1‑75 (MQ=255)
   gACTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCa                                                                           <  2:5739497/76‑1 (MQ=255)
   gACTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCa                                                                           <  1:7438386/76‑1 (MQ=255)
    aCTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCTCTCTCCAg                                                                          >  2:5394346/1‑76 (MQ=255)
    aCTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCTCCCCCCAg                                                                          >  2:4555939/1‑76 (MQ=255)
     cTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCTCGCTCCAGc                                                                         >  2:3058514/1‑76 (MQ=255)
      tCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCTCGCTCCAGcc                                                                        >  1:5931451/1‑76 (MQ=255)
       ccAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCCCGCTCCAGccc                                                                       >  1:3387440/1‑75 (MQ=255)
        cAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCCCGCTccacccac                                                                      >  2:2479509/1‑76 (MQ=255)
         aGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCTCGCCCCCGCccca                                                                     >  1:5813650/1‑76 (MQ=255)
         aGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCCCGCTCCCGCcaca                                                                     >  2:3152843/1‑76 (MQ=255)
           cgcCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATc                                                                   <  1:6391867/76‑1 (MQ=255)
            gcCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCTCGCTCCAGCCACCTCa                                                                  >  1:6082156/1‑76 (MQ=255)
             ccATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAg                                                                 <  2:1871630/76‑1 (MQ=255)
             ccATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCGCGCTCCAGCCACATCAg                                                                 >  2:3875640/1‑76 (MQ=255)
             ccATCACAATCCGGGACGGGAGTTTTTCTTGCTGTTGGTGGTCGCGCCACGGCAGGCCCCCCCCCAGCCACATCAg                                                                 <  2:5314113/76‑1 (MQ=255)
              cATCGCAATACGCGACGGGATTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAgg                                                                <  1:5642929/76‑1 (MQ=255)
              cATCCCAATACAGCACGGGAGTTTTTCTTGCTGATGGTGGTTGCCCCACGGCAGACCCCCCTCCAGCCACATCAgg                                                                <  2:2007533/76‑1 (MQ=255)
                 cGCAATACGCGAGGGGATTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGggcg                                                             <  2:3348750/76‑1 (MQ=255)
                   cAATCCGCGCCGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGggcggc                                                           <  2:2595233/76‑1 (MQ=255)
                   cAATACGCGACGGGATTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGggcggc                                                           <  2:6501234/76‑1 (MQ=255)
                   cAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCCCGCTCCCCCCCCATCAAgggggg                                                           >  1:3058393/1‑75 (MQ=255)
                    aaTACGCGGGGGGATTTTTTCTTGCTGGTGGTGGTTGCCCCGCGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCa                                                          <  2:3984251/76‑1 (MQ=255)
                    aaTACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCTCGCTCCAGCCACATCCGGGCGGCa                                                          >  2:5545709/1‑76 (MQ=255)
                    aaTACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCCCGCTCCAGCCACATCCAGGCGGCa                                                          >  2:8283871/1‑76 (MQ=255)
                    aaTACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCa                                                          <  1:7782997/76‑1 (MQ=255)
                     aTACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCaa                                                         <  1:8578545/76‑1 (MQ=255)
                      tACGGGAGGGGTTTTTTTTTTGCGGATGGGGGTTGCGCCACGGCAGGCCCCGCTCCAGCCACATCAGGGCGGCAat                                                       <  2:5318944/76‑1 (MQ=255)
                      tACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAat                                                       <  1:404349/76‑1 (MQ=255)
                       aCGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAata                                                      <  1:2340277/76‑1 (MQ=255)
                          ggACGGGAGTTTTTCTTGCTGATGGGGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTa                                                   <  1:1347049/75‑1 (MQ=255)
                            gcGGGATTTTTTCTTGCTGATGGGGGTTGCCCCCCGGCGGGCCTCCCCCCAGCCACATCAGGGCGGCAATATTaca                                                 <  2:6126013/75‑1 (MQ=255)
                            aggggAGTTTTTCTTGCTGATGGTGGTTGCCCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTaca                                                 <  2:600313/74‑1 (MQ=255)
                               ggAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCTCGCTCCCGCCCCATCAGGGCGGGAATATTAAACAg                                              >  1:7085544/1‑76 (MQ=255)
                               ggAGTTTTTCTTGCTGAGGGTGGTTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAg                                              <  1:4288563/76‑1 (MQ=255)
                                gAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCCCGCTCCCGCCACATCAGGGCGGCACAATTAcacaag                                             >  2:1614636/1‑74 (MQ=255)
                                gAGTTTTTCTTGCGGATGGTGGTTGCGCCACGGCAGCCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGt                                             <  1:2445911/76‑1 (MQ=255)
                                gAGTTTTTCTGGCTGATGGTGGTGGCGCCAGGGCAGCCCCCGCCCCAGCCACATCAGGGCGGCAATATTACACAGt                                             <  2:865705/76‑1 (MQ=255)
                                 aGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGGCTCGCTCCCGCCACATCAGGGCGGCAATATTACACAGTa                                            >  1:512413/1‑76 (MQ=255)
                                 aGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCCATATTACACAGTa                                            >  1:4201887/1‑76 (MQ=255)
                                   tttttCTTGCTGATGGTGTTTGCCCCACGGCAGCCCTCCCCCCAGCCACATCAGGGCGGCAATATTACACAGTaaa                                          <  1:5063885/76‑1 (MQ=255)
                                    ttttCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTaaaa                                         <  2:5896808/76‑1 (MQ=255)
                                      ttCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCCCATCAGGGCGGCAATATTACACAGTAAAATa                                       >  1:3793207/1‑76 (MQ=255)
                                       tCTTGTTGATGGTGGTTCCTCCACGGCAGCCCCCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATaa                                      <  1:8201584/76‑1 (MQ=255)
                                       tCTTGCTGATGGTGGTTGCGCCACGGCAGGGCCCGCTCCAGCCACCTCAGGGCGGCAATAATACACAGTAAAATaa                                      >  2:4883536/1‑76 (MQ=255)
                                       tCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGGAATATTACACAGTAaaaaaa                                      >  1:5346654/1‑76 (MQ=255)
                                        cTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAATAaaaaaaa                                     >  1:2795390/1‑75 (MQ=255)
                                         ttGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCACGCACATTACGGCGGCAATATTACACCGTAAAAAAACg                                    >  2:6200435/1‑76 (MQ=255)
                                         ttGCTGAGGGTGTTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACg                                    <  2:8838014/76‑1 (MQ=255)
                                           gCTGATGGTGGTTGCGCCCCGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGt                                  <  2:146966/76‑1 (MQ=255)
                                           gCTGATGGTGGTTGCGCCCCGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGt                                  <  1:6911368/76‑1 (MQ=255)
                                           gCTGATGGTGGTTGCGCCACGGCAGGCCTTGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACggc                                  >  1:4451549/1‑75 (MQ=255)
                                           gCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGt                                  >  1:3493137/1‑76 (MQ=255)
                                           gCTGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGt                                  <  2:2654364/76‑1 (MQ=255)
                                            cTGATGGTGGTTGCGCCACGGCAGGCCCCCCTCCCGCCACATCAGGGCGGCAATATTACACAGTAAAATAACCGTa                                 >  2:4360164/1‑76 (MQ=255)
                                             tGATGGTGGTTGCGCCACGGCAGGGCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTaa                                >  2:5764228/1‑76 (MQ=255)
                                             tGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCCGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTaa                                >  2:6886793/1‑76 (MQ=255)
                                             tGATGGTGGTTGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTaa                                <  1:1612086/76‑1 (MQ=255)
                                             tGATGGTGGTTGCGCACCGGCAGGCCCCCCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTaa                                <  2:8507831/76‑1 (MQ=255)
                                              gATGGTGGTTGCGCCACGGCAGGGCCCCCCCCAGCCACAACAGGGCGGCAATATTACACAGTAAAAAAAAGGTaaa                               >  1:8884307/1‑75 (MQ=255)
                                              gATGGTGGTGCCGCCACGGCAGGCCCCCCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAAt                               <  1:6510714/76‑1 (MQ=255)
                                               aTGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCCAGGGGGCAATATTACACAGTAAAAAAACGGGAATa                              >  1:151848/1‑76 (MQ=255)
                                               aTGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATa                              >  1:5664899/1‑76 (MQ=255)
                                                 ggtggtTGCCCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACt                            <  2:7429867/76‑1 (MQ=255)
                                                 ggtggtGGCGCCACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACt                            <  2:8372768/76‑1 (MQ=255)
                                                  gtggtTGCGCCACGGCAGGCCTCCCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTg                           <  1:7132590/76‑1 (MQ=255)
                                                   gggtGGCCCCAGGGGAGGCCTCCCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGa                          <  1:5156/75‑1 (MQ=255)
                                                    ggtTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGAc                         >  2:3771556/1‑76 (MQ=255)
                                                     gtTGCGCCACGGCAGGCCTCGCACCTTCCACAGCAGTTCGGCAAGTATTTTTAGTTACATAACGGCAATACTgtac                        >  2:2733651/1‑73 (MQ=255)
                                                       tGCGCCACGGCAGGGCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGc                      >  2:3520771/1‑76 (MQ=255)
                                                       tGCGCCACGGCAGGCCTCGCTCCAGCCACCTCAGGGCGGCAATATTACACAGTAAAAAAACGGTAATACTGACAGc                      >  2:7154765/1‑76 (MQ=255)
                                                         cgcCAGGGCAGCCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAg                    <  1:1394542/76‑1 (MQ=255)
                                                         cgcCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAg                    <  2:9003736/76‑1 (MQ=255)
                                                           ccAGGGCGGCCCCCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaa                  <  1:1245167/76‑1 (MQ=255)
                                                           ccACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAAAATTACACAGTAAAATAACGGGAATACTGACGGCAGaa                  >  1:4580701/1‑76 (MQ=255)
                                                            cgcGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaat                 <  1:6575696/74‑1 (MQ=255)
                                                            cACGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaat                 <  2:8753517/76‑1 (MQ=255)
                                                              cggcggCCCTCCCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaataa               <  1:7983430/76‑1 (MQ=255)
                                                              cGGCAGGCCTCGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaataa               <  2:7241003/76‑1 (MQ=255)
                                                              cGGCAGGCCACCCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaataa               <  1:1029592/76‑1 (MQ=255)
                                                              cGGCAGCCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaataa               <  2:6793687/76‑1 (MQ=255)
                                                               ggCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGaataat              >  2:6633135/1‑76 (MQ=255)
                                                                gCAGGTCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATg             >  2:7829104/1‑76 (MQ=255)
                                                                  aGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCa           >  2:1199576/1‑76 (MQ=255)
                                                                    gCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGa         >  2:566917/1‑76 (MQ=255)
                                                                      cccGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAAt       <  2:3622588/74‑1 (MQ=255)
                                                                       ccGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATc      <  1:2583278/75‑1 (MQ=255)
                                                                       ccGCCCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATc      <  2:4054010/75‑1 (MQ=255)
                                                                        cGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCg     >  1:8983036/1‑76 (MQ=255)
                                                                        cGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCg     >  1:6180436/1‑76 (MQ=255)
                                                                         gCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGc    >  1:7303640/1‑76 (MQ=255)
                                                                          cTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGcc   >  2:731024/1‑76 (MQ=255)
                                                                          cTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGcc   >  2:7954618/1‑76 (MQ=255)
                                                                          cTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGcc   >  1:8901955/1‑76 (MQ=255)
                                                                          cTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGcc   >  2:5083555/1‑76 (MQ=255)
                                                                           |                                                                            
AGCGACTCCAGCGCCATCGCAATACGCGACGGGAGTTTTTCTTGCTGATGGTGGTTGCGCCACGGCAGGCCTCGCTCCAGCCACATCAGGGCGGCAATATTACACAGTAAAATAACGGTAATACTGACGGCAGAATAATGCAGAATCGCCA  >  REL606/931239‑931389

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