Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,854,4010AC48.6% 111.5 / 54.7 109K88N (AAA→AACyeaOhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (18/38);  new base C (0/53);  total (18/91)
Fisher's exact test for biased strand distribution p-value = 1.66e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 145 total aligned reads displayed.

CTTTCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTACGT  >  REL606/1854326‑1854476
                                                                            |                                                                           
cTTTCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCaaaaaa                                                                             <  1:7750424/76‑1 (MQ=255)
   tCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAg                                                                          <  1:7202467/76‑1 (MQ=255)
    cgcgAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGc                                                                         <  1:4131809/76‑1 (MQ=255)
    cgcgAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGc                                                                         <  2:6369234/76‑1 (MQ=255)
    cgcgAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGc                                                                         <  1:2142315/76‑1 (MQ=255)
    cgcgAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGc                                                                         <  2:737773/76‑1 (MQ=255)
     gcgAGCAAGATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGcc                                                                        <  1:1199118/76‑1 (MQ=255)
      cgAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCg                                                                       <  2:5609737/76‑1 (MQ=255)
         gaaaTATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCGGGCGGACATCGCCAAAAACCAGCCGCTg                                                                    <  2:1639122/74‑1 (MQ=255)
           aaTATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGAc                                                                  <  1:124133/76‑1 (MQ=255)
            atatCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGCCTGGCGGACATCGCCAAAAAACAGCCGCTGAcc                                                                 <  1:1156521/76‑1 (MQ=255)
             tatCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGCGACGGCCAAAAAAAACCCGCGAccc                                                                >  2:1304114/1‑76 (MQ=255)
              atCTTGCAGAACGGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCt                                                               <  1:7052305/76‑1 (MQ=255)
               tCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTg                                                              <  1:1050260/76‑1 (MQ=255)
               tCTGGCAGACCGGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACACCGCCAAAAACCAGCCGCTGACCCTg                                                              <  2:2747430/76‑1 (MQ=255)
                 ttGCAGACCGGGCGCAACAGGAGCAAGAGGGAAAGCGGCGGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGct                                                            <  2:6682834/76‑1 (MQ=255)
                   gaaGAACTGGCGAAACACGAAAAAGAAGGAAAGGGGCTGGCGGACATCGCAAAAAAACAGCCGCTGACCCTGctct                                                          <  2:1203480/74‑1 (MQ=255)
                   gaaGAACGGCCGCAACACGAGAAAGAAGGAAGGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGctct                                                          <  2:5085647/74‑1 (MQ=255)
                   gCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGGAACGCCCAAAAACAACCGCCGATCACGctct                                                          >  1:1510893/1‑76 (MQ=255)
                   gCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGctct                                                          <  1:5309419/76‑1 (MQ=255)
                   gCAGAACGGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGctct                                                          <  1:311716/76‑1 (MQ=255)
                   gCAAAACTGGCGCAACGCGAGCAAGAGGGAAAGCGGCTGGCGGACATCGCAAAAAAACAGCCGCTGACCCTGctct                                                          <  1:3490345/76‑1 (MQ=255)
                    cAGACCTGGCGCAACACGAGAAAGAAGGAAAGCGGCGGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTa                                                         <  2:285834/76‑1 (MQ=255)
                     aaaaCGGGCGCAACACGAGAAAGAAGAAAAGCGGCGGGCGGACACCGCAAAAAACCAGCCGCTGACCCTGCTCTAc                                                        <  2:7752806/74‑1 (MQ=255)
                      gAACTGGCGCAACACGAGAAAGAAGGAAAGCGGCGGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACg                                                       <  2:5776410/76‑1 (MQ=255)
                       agcTGGCGCCACCCGAGAAAAAAGAAAGGGGGGTGGCCGCCATCGCAAAAAACCAGCCGCTGACCCTGCTCTACgc                                                      <  2:3000187/74‑1 (MQ=255)
                       aaCTGGCGCAACAGGAGCAAGAAGGAAAGCGCCTGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACgc                                                      <  1:6326757/76‑1 (MQ=255)
                       aaCTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgc                                                      <  1:1657328/76‑1 (MQ=255)
                       aaCTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgc                                                      <  1:7469120/76‑1 (MQ=255)
                         cgggCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgcag                                                    <  1:2816847/74‑1 (MQ=255)
                          tGGCGCAACACGAGCAAGAAGAAAAGCGGCGGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACgcagc                                                   <  2:5658659/76‑1 (MQ=255)
                          tGGCCCAACACGAGAAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgcagc                                                   <  1:3810073/76‑1 (MQ=255)
                          gggCGCAACGCGAGAAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgcagc                                                   <  1:5476475/75‑1 (MQ=255)
                           ggCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgcagca                                                  <  1:2079171/76‑1 (MQ=255)
                            gcgcAACACGAGAAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCaa                                                 <  1:6729351/76‑1 (MQ=255)
                             cgcAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCaaa                                                <  1:3576024/76‑1 (MQ=255)
                              gcAACAGGAGAAAGAAGGAAAGCGGCGGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCaaaa                                               <  2:23758/76‑1 (MQ=255)
                               cAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCaaaaa                                              <  1:5469768/76‑1 (MQ=255)
                               cAACACGAGAAAGAAGGAAAGCGGCGGGCGAACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCaaaaa                                              <  2:2111783/76‑1 (MQ=255)
                               aaaCACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCaaaaa                                              <  1:5106136/75‑1 (MQ=255)
                                aaCAGGAGCAGAAAGGAAGGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCaaaaaa                                             <  2:2509700/76‑1 (MQ=255)
                                 acacGAGAAAGAAGGAAAGGGGCTGGCGGACACCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAac                                            <  1:5927221/76‑1 (MQ=255)
                                 acacGAGAAAAAAGGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAac                                            <  2:3063922/76‑1 (MQ=255)
                                  cacGAGCAAGAAGGAAAGCGGCTGGGGGGTATTGCCCAAAAAAAGCCGCCGGCCCCCCTCTCCCCCGCcaaaaaaa                                           >  1:7242004/1‑74 (MQ=255)
                                  cacGAGAAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAaca                                           <  1:4313582/76‑1 (MQ=255)
                                  cacGACAAAGAAGAAAAGCGGCTGGCGGCCACCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAaca                                           <  1:4460469/76‑1 (MQ=255)
                                   gagAGCAAGAAGGAAAGCGGCTGGGGAACACCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAacac                                          <  2:389754/74‑1 (MQ=255)
                                    cGAGCAAGAAGAAAGGCGGCTGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAAcacc                                         <  2:2757459/76‑1 (MQ=255)
                                    cGAGAAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAAcacc                                         <  1:999516/76‑1 (MQ=255)
                                       gCAAGAAGGAAAGCGGCTGGCGGCCATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACg                                      <  1:5764476/76‑1 (MQ=255)
                                       gCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACg                                      <  2:2235548/76‑1 (MQ=255)
                                        cAAGAAGGAAAGCGGCTGGGGGGGAACGCCCAAAAAAAACCGCTGAACCCGCTCTACCCCGCAAAAAAaaacacac                                     >  1:5087258/1‑74 (MQ=255)
                                         aagaagAAAAGCGGCTGGCGGACATCCCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCa                                    <  1:5035650/76‑1 (MQ=255)
                                          agatgGAAAGGTGCTGGCGGACGTCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAg                                   <  1:4387404/76‑1 (MQ=255)
                                           gaagaaAAGCGGGTGGCGAACACCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGa                                  <  2:4890829/76‑1 (MQ=255)
                                            agggAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGaa                                 <  2:7781694/74‑1 (MQ=255)
                                            aagGAAACCGGCGGGCGGACATCGCAAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGaa                                 <  1:5330502/76‑1 (MQ=255)
                                             agGAAAGCGGCTGGGGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAAc                                <  2:488623/76‑1 (MQ=255)
                                                aaaGGGGCGGGCGAACACCCCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCAt                             <  2:2371626/76‑1 (MQ=255)
                                                aaaGCGGCGGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCAt                             <  1:6856248/76‑1 (MQ=255)
                                                  aGCGGCTGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATgc                           <  2:6507687/76‑1 (MQ=255)
                                                   gcggcgGGCGCACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATgcg                          <  2:245834/76‑1 (MQ=255)
                                                      gcgggCGCACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTg                       <  2:5311414/73‑1 (MQ=255)
                                                       cTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTgg                      >  1:7561404/1‑76 (MQ=255)
                                                       cTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTgg                      >  1:7289819/1‑76 (MQ=255)
                                                        tGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGt                     >  2:7554755/1‑76 (MQ=255)
                                                        tGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAACAAAAAACACCACGCAGAACCATGCGCTGGt                     >  2:6252722/1‑76 (MQ=255)
                                                        tGGCGGACACCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGt                     <  1:2138753/76‑1 (MQ=255)
                                                         ggCGGACCCCACCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTg                    <  2:6602335/76‑1 (MQ=255)
                                                         ggCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTg                    <  2:4351923/76‑1 (MQ=255)
                                                         ggCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTg                    <  2:3632128/76‑1 (MQ=255)
                                                           gagACTTCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt                  <  2:7680179/74‑1 (MQ=255)
                                                           cGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt                  >  1:6645250/1‑76 (MQ=255)
                                                           cGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt                  >  2:402689/1‑76 (MQ=255)
                                                           cGGACATCGCCAAAAAACAGCCGCTGACACTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGAGCt                  >  2:5941982/1‑76 (MQ=255)
                                                           cGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt                  <  2:3246265/76‑1 (MQ=255)
                                                            ggACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGTAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTg                 <  2:2939303/76‑1 (MQ=255)
                                                             gACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTgg                <  1:4288492/76‑1 (MQ=255)
                                                               cATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGcc              >  2:119209/1‑76 (MQ=255)
                                                               cATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGcc              >  2:4291225/1‑76 (MQ=255)
                                                                 accGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCg             <  1:2045987/74‑1 (MQ=255)
                                                                  tCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGa            <  1:3078184/76‑1 (MQ=255)
                                                                   cGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGAc           <  1:2842733/76‑1 (MQ=255)
                                                                    gcaaaaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACt          <  1:5993989/73‑1 (MQ=255)
                                                                    gcaaaaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACt          <  2:2926827/73‑1 (MQ=255)
                                                                    gaaaaaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACt          <  1:6535675/73‑1 (MQ=255)
                                                                    gCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACt          <  1:6285402/76‑1 (MQ=255)
                                                                    gCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACt          >  2:3570750/1‑76 (MQ=255)
                                                                    ccaaaaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACt          <  1:5302606/73‑1 (MQ=255)
                                                                    ccaaaaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACt          <  1:4664089/73‑1 (MQ=255)
                                                                     ccAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTg         <  1:7072892/76‑1 (MQ=255)
                                                                     ccAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTg         >  2:7736600/1‑76 (MQ=255)
                                                                     ccAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTg         >  2:5659738/1‑76 (MQ=255)
                                                                      cAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTgg        <  1:4389059/76‑1 (MQ=255)
                                                                      cAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTgg        >  1:2383409/1‑76 (MQ=255)
                                                                       aaaaaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGc       >  1:940322/1‑76 (MQ=255)
                                                                         aaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa     <  1:5636300/76‑1 (MQ=255)
                                                                          aaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTAc    >  1:172480/1‑76 (MQ=255)
                                                                           aaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGATGCTGGCCGACTTGCTaag   >  2:1937359/1‑74 (MQ=255)
                                                                            aCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTAAGt  >  2:2824474/1‑76 (MQ=255)
                                                                            |                                                                           
CTTTCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTACGT  >  REL606/1854326‑1854476

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