Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,854,4010AC45.8% 76.3 / 47.0 83K88N (AAA→AACyeaOhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/26);  new base C (0/38);  total (19/64)
Fisher's exact test for biased strand distribution p-value = 1.56e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 102 total aligned reads displayed.

TTTCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTAC  >  REL606/1854327‑1854474
                                                                            |                                                                         
tttCGCGAGCAATATCTTGCAGAACTGGTGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAAc                                                                          <  1:2080986/76‑1 (MQ=255)
tttCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGCAACGCCccaaaaa                                                                          >  1:521684/1‑75 (MQ=255)
 ttCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGCAACGCGCaaaaaaa                                                                          >  2:3578878/1‑75 (MQ=255)
 ttCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACa                                                                         <  1:6194016/76‑1 (MQ=255)
 ttCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACa                                                                         <  2:5502654/76‑1 (MQ=255)
 ttCGCGAGCAATATCTTGCAGAACTGACGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACa                                                                         <  2:215091/76‑1 (MQ=255)
   cgcgAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGc                                                                       <  2:3945296/76‑1 (MQ=255)
            tatCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATAGCCAAAAAACAGCCGCTGAccc                                                              <  1:600699/76‑1 (MQ=255)
             atCTTGCAGACCGGGCGCAACACGAGCAAGAAGGAAAGCGTCGGGCGGACACCGCCAAAAACCAGCCGCTGACCCt                                                             <  2:4723708/76‑1 (MQ=255)
             atCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCt                                                             <  1:4451688/76‑1 (MQ=255)
             atCTTGCAGAACTGGCGCAACACGAGCAAGAAGAAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCt                                                             <  2:6139352/76‑1 (MQ=255)
             accTTGAAGAACGGGCGCACCAGGAGCAAAAAGAAAGGCGGCGGGCGGACATCCCAAAAAACCAGCCGCTGACCCt                                                             <  2:5200859/74‑1 (MQ=255)
              tCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGGATCGCCCAAAAAAAACCGCCGAccccc                                                            >  1:3569694/1‑74 (MQ=255)
              tCTGGCAGAACTGGCTCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTg                                                            <  1:1483510/76‑1 (MQ=255)
              tCTGGCAGAACGGGCGCAACACGAAAAAGAAGAAAAGCGGCGGGGGGACATCGCCAAAAACCAGCCGCTGACCCTg                                                            <  2:3222307/76‑1 (MQ=255)
                ttGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGct                                                          <  1:4135917/76‑1 (MQ=255)
                  gCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGctct                                                        <  2:3219222/76‑1 (MQ=255)
                   cAGAGCTGCCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTa                                                       <  1:5343026/76‑1 (MQ=255)
                   cAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTa                                                       <  1:5198298/76‑1 (MQ=255)
                   cAGAACGGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTa                                                       <  1:3697506/76‑1 (MQ=255)
                   aagagCTGGCGCAACACGAACAAGAGGGAAGGCGGCGGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTa                                                       <  2:2444942/75‑1 (MQ=255)
                    gaaaCTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTAc                                                      <  1:2804163/74‑1 (MQ=255)
                    aGAACTGGGCCACACGGAGAAAGAAGGAAGACGGCTGGCGACCACCCCCAAAAAACAGCCGCTGACCCTGCTCTAc                                                      <  2:299451/76‑1 (MQ=255)
                    aGAACGGGCGCAACACGAGAAAGAAGGAAAGGGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTAc                                                      <  1:5395268/76‑1 (MQ=255)
                    aGAACGGGCGCAACACGAGAAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTAc                                                      <  1:193250/76‑1 (MQ=255)
                    aGAACGGCCGCAACCCAAAAAAGAAGAAAAGCGGCTGGCGGACATCCCCAAAAACCAGCCGCTGACCCTGCTCTAc                                                      <  2:3229656/76‑1 (MQ=255)
                     gagcTGGAGCACCACGGACAAGAGGAAAGGCGGCTGGCGAACATCCCAAAAAACCAGCCGCTGACCCTGCTCTACg                                                     <  2:4551263/73‑1 (MQ=255)
                     gAACGGGCGCAACACGAGCAAGAGGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACg                                                     <  1:5846977/76‑1 (MQ=255)
                        cgggCCCAACACGAGCAAGAAGGAAGGCGGCGGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACgcag                                                  <  1:3643637/74‑1 (MQ=255)
                        cgggCACAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACgcag                                                  <  2:2183884/74‑1 (MQ=255)
                        cTGGCCCAACACAAGCAAAAAGGAACGCGGTTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACgcag                                                  <  1:3579856/76‑1 (MQ=255)
                          ggggcAACACGAGCAAGAAGGAAAGGGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgcagca                                                <  1:2303948/73‑1 (MQ=255)
                          ggCGCAACACGAGAAAGAAGGAAAGCGGCGGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACgcagca                                                <  2:474903/76‑1 (MQ=255)
                           gcgcAACAAAAGAAAGAAGGAAGGCGGCGGGCGGCCACCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCaa                                               <  2:3119445/76‑1 (MQ=255)
                            cgcAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCaaa                                              <  1:4362249/76‑1 (MQ=255)
                            cgcAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCaaa                                              <  1:5005916/76‑1 (MQ=255)
                            ccaaaCACGAGAAAGAAGGAAAGCGGCGGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCaaa                                              <  2:6020337/73‑1 (MQ=255)
                             gcAACACGAGAAAGAGGGAAAGCGGCTGGCGGACACCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCaaaa                                             <  2:4167500/76‑1 (MQ=255)
                              cAACACGAGCAAGAAGGAAAGCGGCGGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCaaaaa                                            <  2:2459666/76‑1 (MQ=255)
                               aaCACGAGAAAGAAGGAAGGCGGCGGGCGGACACCCCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCaaaaaa                                           <  1:4959839/76‑1 (MQ=255)
                                acacGAGAAAAAAGGAAAGCGGCTGGCGGAAATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAac                                          <  1:3622496/76‑1 (MQ=255)
                                   cggagaAAAAAGAAAAGTGGTTGGCGGACACCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAacac                                        <  2:2568468/74‑1 (MQ=255)
                                   cGAGAAAGAAGGAAAGCGGCGGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAAcacc                                       <  2:475907/76‑1 (MQ=255)
                                     gAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAAcacca                                      <  1:840988/76‑1 (MQ=255)
                                       gCAAAAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACg                                    <  1:3902718/76‑1 (MQ=255)
                                         tctaagGAAAGCGGCTGGCGGCCACCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCa                                  <  1:2683783/73‑1 (MQ=255)
                                         aagaagGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCa                                  <  1:4200987/76‑1 (MQ=255)
                                         aagaagGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCa                                  <  1:643042/76‑1 (MQ=255)
                                         aagaagAAAAGCGGCTGGCGGACACCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCa                                  <  2:3011637/76‑1 (MQ=255)
                                          agaagGAAGGGGGCGGGCGAACAGCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAg                                 <  2:2667227/76‑1 (MQ=255)
                                          agaagGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAg                                 <  1:1814641/76‑1 (MQ=255)
                                            aagAAAAGGGGCTGCGGAACAACGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGaa                               <  2:1088817/76‑1 (MQ=255)
                                            aaagAAAGAGCCTGCGGGACACCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGaa                               <  2:3559283/73‑1 (MQ=255)
                                             agGAAAGCGGCTGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAAc                              <  1:270027/76‑1 (MQ=255)
                                              ggAAGGCGCCGGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAAcc                             <  2:5692298/76‑1 (MQ=255)
                                              ggAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAAcc                             <  1:5584161/76‑1 (MQ=255)
                                              ggAAAGCGGCTGGCGGACATCCCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAAcc                             <  1:3616715/76‑1 (MQ=255)
                                              ggAAAGCGGCTGGCGAACAGCCCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAAcc                             <  2:550893/76‑1 (MQ=255)
                                                aaggCGGCTGGCGGCCATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCAt                           <  2:1652877/73‑1 (MQ=255)
                                                aaaGCGGCTGGCGGACATCGCCAAAAATCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCAt                           <  1:1588423/76‑1 (MQ=255)
                                                aaaGCGGCTGGCGGACATCGCAAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCAt                           <  1:3386516/76‑1 (MQ=255)
                                                aaaGCGCCGGGCGGACACCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCAt                           <  2:3054216/76‑1 (MQ=255)
                                                 aggCGGTTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATg                          <  2:856018/74‑1 (MQ=255)
                                                 aaccGGCTGGCGGACAGCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATg                          <  1:4983710/73‑1 (MQ=255)
                                                 aaGCGGCGGGCGGACATCGCAAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATg                          <  1:2765777/76‑1 (MQ=255)
                                                   gCGGGTGGCGGACACCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATgcg                        <  2:1097607/76‑1 (MQ=255)
                                                   gCGGCTGGCGGACACCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATgcg                        <  1:6169275/76‑1 (MQ=255)
                                                   gCGGCTGGCGGACACCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATgcg                        <  1:3994951/76‑1 (MQ=255)
                                                     ggCTGGCGGGCATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCt                      >  2:5541224/1‑76 (MQ=255)
                                                     ggCTGGCGGGCATCGCCAAAAAACAGCCGCTGACCCTGCTCCAAGCAACACAAAACACCACGCCGAACCATGAGCt                      >  2:3187053/1‑76 (MQ=255)
                                                      gCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTg                     >  2:4075936/1‑76 (MQ=255)
                                                      gCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTg                     >  1:325001/1‑76 (MQ=255)
                                                       cTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTgg                    >  2:191714/1‑76 (MQ=255)
                                                       cTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTgg                    >  1:2865213/1‑76 (MQ=255)
                                                        tGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGt                   >  1:2995371/1‑76 (MQ=255)
                                                         ggCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTg                  <  2:5468844/76‑1 (MQ=255)
                                                         ggCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTg                  >  2:4415053/1‑76 (MQ=255)
                                                          gCGGACGCCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACGCCACGCAGAACCATGGGCTGGTGc                 <  2:5869873/76‑1 (MQ=255)
                                                           cgaaCATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt                <  2:1811003/73‑1 (MQ=255)
                                                           cgaaCATCCCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt                <  2:2751669/73‑1 (MQ=255)
                                                           cGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt                >  1:3465319/1‑76 (MQ=255)
                                                           cGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt                >  2:1537452/1‑76 (MQ=255)
                                                           cGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt                >  2:2599227/1‑76 (MQ=255)
                                                           cGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt                >  2:1002589/1‑76 (MQ=255)
                                                            ggACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTg               >  1:4732444/1‑76 (MQ=255)
                                                            ggACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTg               >  2:2923897/1‑76 (MQ=255)
                                                            ggACATCCCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTg               <  2:638107/76‑1 (MQ=255)
                                                            gaaCATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTg               <  2:4230404/74‑1 (MQ=255)
                                                             gACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTgg              >  2:2707306/1‑76 (MQ=255)
                                                             gACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTgg              >  2:2556000/1‑76 (MQ=255)
                                                               cATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGcc            >  2:2080260/1‑76 (MQ=255)
                                                                    ccAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTg       <  1:1767647/76‑1 (MQ=255)
                                                                     caaaaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTg       <  2:6307099/74‑1 (MQ=255)
                                                                      cAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTgg      >  1:2118610/1‑76 (MQ=255)
                                                                       aaaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGc     <  1:1797053/76‑1 (MQ=255)
                                                                        aaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCt    <  1:5254120/76‑1 (MQ=255)
                                                                         aaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa   <  1:5694155/76‑1 (MQ=255)
                                                                         aaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa   <  2:1500713/76‑1 (MQ=255)
                                                                          aaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTaa  >  1:2256550/1‑75 (MQ=255)
                                                                          aaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTAc  <  1:2798617/76‑1 (MQ=255)
                                                                            |                                                                         
TTTCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTAC  >  REL606/1854327‑1854474

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