Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,854,4010AC50.6% 60.2 / 110.8 81K88N (AAA→AACyeaOhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (12/28);  new base C (0/41);  total (12/69)
Fisher's exact test for biased strand distribution p-value = 7.90e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.14e-07
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 115 total aligned reads displayed.

CTTTCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTA  >  REL606/1854326‑1854473
                                                                             |                                                                        
cTTTCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGGAGAGCCaaaaaa                                                                          >  1:1840138/1‑76 (MQ=255)
cTTTCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGCAACGCCcaaaaa                                                                          >  2:3471090/1‑76 (MQ=255)
 tttCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGGAAACGTCCAAAAAc                                                                         >  2:605080/1‑76 (MQ=255)
  ttCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGCAATCCCCaaaaaaa                                                                         >  1:1904765/1‑75 (MQ=255)
  ttCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGCAACGCCcacaaaaa                                                                        >  1:6365080/1‑74 (MQ=255)
  ttCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACa                                                                        <  1:2316265/76‑1 (MQ=255)
   tCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAg                                                                       <  2:7340922/76‑1 (MQ=255)
   tCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAg                                                                       <  2:1069371/76‑1 (MQ=255)
   tCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAg                                                                       <  1:4100697/76‑1 (MQ=255)
   tCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAg                                                                       <  1:2077896/76‑1 (MQ=255)
   tCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAg                                                                       <  1:3627623/76‑1 (MQ=255)
         gCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTg                                                                 <  1:1437770/76‑1 (MQ=255)
         gCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCAAAAAACCAGCCGCTg                                                                 <  2:3453779/76‑1 (MQ=255)
         gCAATATCTTGCAGAACTGGCGCAACACGAGAAAGAAGGAAAGCGGCGGGCGGACATCGCAAAAAACCAGCCGCTg                                                                 <  2:2206508/76‑1 (MQ=255)
          cAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGCAACGCCCAAAAAAAACCGccgc                                                                >  2:7203079/1‑73 (MQ=255)
           aaTATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGGATAGCCCAACAAAAACCGCTGaa                                                               >  1:2110952/1‑75 (MQ=255)
              atCTTGCAGACCGGGCGCAACACGAGAAAGAAGGAAGGCGGGTGGCGGACATCCCAAAAAACCAGCCGCTGACCCt                                                            <  2:3951984/76‑1 (MQ=255)
              atCTGGCAGAGCGGGCACAACCCGAGCAAAAAGGAAACCGGGTGGCAGACACCCCCAAAAACCAGCCGCTGACCCt                                                            <  2:4424317/76‑1 (MQ=255)
                  tGCAGAACGGGCGCAACACGACCAAGAGGGAAAGCGGCTGGCGGACATCGCCAAAACCAGCCGCTGACCCTGctct                                                       <  2:7334416/76‑1 (MQ=255)
                   gCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGctct                                                       <  1:2588418/76‑1 (MQ=255)
                   gCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGctct                                                       <  1:1131008/76‑1 (MQ=255)
                   gCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCGGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGctct                                                       <  2:698358/76‑1 (MQ=255)
                   gCAGAACTGGCGCAACACGAGCAAGAAGAAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGctct                                                       <  1:1740810/76‑1 (MQ=255)
                    cAGAACTGGCGCAACGCGAGCAAAAAGAAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTa                                                      <  1:652462/76‑1 (MQ=255)
                    cAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTa                                                      <  1:4068113/76‑1 (MQ=255)
                    cAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTa                                                      <  1:5484783/76‑1 (MQ=255)
                     aaaaCGGGCGAAACGCGAGAAAGAAGGAAAGCGGCGGGCGGACTTCGCCAAAAACCAGCCGCTGACCCTGCTCTAc                                                     <  2:2938619/74‑1 (MQ=255)
                     aGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTAc                                                     <  1:7467020/76‑1 (MQ=255)
                     aGAACTGGCGCAACACGAGAAAGAAGGAAGGCGGCGGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTAc                                                     <  2:1622684/76‑1 (MQ=255)
                        aCTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgca                                                  <  1:5305759/76‑1 (MQ=255)
                         cgggCGCAACACGAGCAAGAAGAAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACgcag                                                 <  2:4130920/74‑1 (MQ=255)
                         cTGGCCAAACACGAGAAAGAGGGAAAGCGGCGGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACgcag                                                 <  2:4699950/76‑1 (MQ=255)
                          tGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgcagc                                                <  1:6642700/76‑1 (MQ=255)
                          gggCGCAACACGAGCAAGAAGGAAAGGGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgcagc                                                <  1:1433391/75‑1 (MQ=255)
                          gggCGCAACACGAGAAAGAAGGAAAGCGGCGGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACgcagc                                                <  2:636810/75‑1 (MQ=255)
                          gggCGAAACACGAGCAAGAAGGAAAGCGGCGGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACgcagc                                                <  2:6577214/75‑1 (MQ=255)
                           ggCGCACCACGAGCAAGAAGAAAAGCGGCTGGCGGACAGCGCCAAAAACCAGCCGCTGACCCTGCTCTACgcagca                                               <  2:6451900/76‑1 (MQ=255)
                           ggCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgcagca                                               <  1:2460508/76‑1 (MQ=255)
                           ggCGCAACACGAGAAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACgcagca                                               <  2:3809587/76‑1 (MQ=255)
                               cAACACGAGAAAGAAGGAAGGCGGCTGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCaaaaa                                           <  2:5274366/76‑1 (MQ=255)
                                aaCACAAGCAAGAAGAAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCaaaaaa                                          <  1:81894/76‑1 (MQ=255)
                                  cacGAGAAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAaca                                        <  2:4422626/76‑1 (MQ=255)
                                  aaggagaAAGAAGGAAAGCGGCTGGCGTACTTCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAaca                                        <  1:3479779/73‑1 (MQ=255)
                                   acGAGCAAGAAGGAAAGCGGCTGGGGGGGATCGCCCAAAAAAAAACGCCGGCCCTGCCCTCTCCAGCcaaaaaaaa                                       >  2:5962254/1‑73 (MQ=255)
                                     gAGCAAGAAGGAAAGCGGCTGGGGGGTTTTCTCCAAAAAAAACCCCCGGACCCGCTTCCCCGCACCAAAAAAaacc                                      >  2:5262654/1‑76 (MQ=255)
                                      aGCAAGAGGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAAcaccac                                    <  1:2973865/76‑1 (MQ=255)
                                        aaggaagGAAAGGGGCTGGCGGCCATCGCAAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGc                                  <  1:7080202/73‑1 (MQ=255)
                                          agaagaAAAGCGGCGGGCGAACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAg                                <  2:5741600/76‑1 (MQ=255)
                                          agaagGAAAGCGGCGGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAg                                <  2:7020283/76‑1 (MQ=255)
                                          aaaagGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAg                                <  2:2757288/74‑1 (MQ=255)
                                            aagGAAAGGGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGaa                              <  1:2836197/76‑1 (MQ=255)
                                            aagGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGaa                              <  1:1444594/76‑1 (MQ=255)
                                             agGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAAc                             <  2:1934417/76‑1 (MQ=255)
                                              ggAAAGCGGCTGGCGGACATCGCAAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAAcc                            <  1:2754559/76‑1 (MQ=255)
                                              ggAAAGCGGCGGGGGGACACCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAAcc                            <  2:1970400/76‑1 (MQ=255)
                                               gAAGGCGGCGGGCCGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCa                           <  1:6935915/76‑1 (MQ=255)
                                                aaaGCGGCGGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCAt                          <  2:241698/76‑1 (MQ=255)
                                                 aggCGGCTGGCGGAAATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATg                         <  2:5143063/74‑1 (MQ=255)
                                                 aaGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATg                         <  1:1550402/76‑1 (MQ=255)
                                                  aGCGGTTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATgc                        <  2:6625457/76‑1 (MQ=255)
                                                  aGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATgc                        <  1:1757370/76‑1 (MQ=255)
                                                  aGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATgc                        <  1:5379438/76‑1 (MQ=255)
                                                  aGCGGCGGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATgc                        <  2:7023101/76‑1 (MQ=255)
                                                    gggCTGGCGGACAGCGCCAAAAACCAGCCGCTGTCCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATgcgc                      <  1:3690686/75‑1 (MQ=255)
                                                      gCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCta                    >  2:3865310/1‑75 (MQ=255)
                                                       cTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTgg                   >  2:281789/1‑76 (MQ=255)
                                                       cTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAAACATGCGCTgg                   >  2:4994288/1‑76 (MQ=255)
                                                        tGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACAACGCAGAACCATGCGCTGGt                  >  2:6807884/1‑76 (MQ=255)
                                                        tGGCGGACACCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGt                  <  1:1357182/76‑1 (MQ=255)
                                                        cgggCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTgg                   <  2:6903156/74‑1 (MQ=255)
                                                         gggCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGt                  <  2:1302077/75‑1 (MQ=255)
                                                         gggCGGACACCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGt                  <  1:2788276/75‑1 (MQ=255)
                                                         gggCGAACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGt                  <  1:393505/75‑1 (MQ=255)
                                                           gcgcACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGc                <  2:4420104/76‑1 (MQ=255)
                                                            cgcaCATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt               <  2:404699/73‑1 (MQ=255)
                                                            cGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt               >  1:2671861/1‑76 (MQ=255)
                                                             ggACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAACACACCACGCAGAACCATGCGCTGGTGCTg              >  2:3993064/1‑76 (MQ=255)
                                                             ggACATCCCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTg              <  2:5596573/76‑1 (MQ=255)
                                                               aCATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGc            >  2:3387624/1‑76 (MQ=255)
                                                               aCATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGc            <  2:4819774/76‑1 (MQ=255)
                                                                cATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGcc           >  2:3582535/1‑76 (MQ=255)
                                                                 caccGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGcc           <  2:689944/73‑1 (MQ=255)
                                                                  aTCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCg          >  2:1948147/1‑76 (MQ=255)
                                                                  aTCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCg          >  1:1028633/1‑76 (MQ=255)
                                                                    ccccAAAAACCAGCCGCTGCCCCTTCTCTACGCAGCAAAAACCACCACGCAGAACCATGCGCTGGTGCTGGCCGAc        <  1:577398/74‑1 (MQ=255)
                                                                     gCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACt       >  2:1578696/1‑76 (MQ=255)
                                                                      ccAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTg      <  1:3840661/76‑1 (MQ=255)
                                                                      ccAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTg      >  2:4376376/1‑76 (MQ=255)
                                                                        aaaaaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGc    >  2:4516298/1‑76 (MQ=255)
                                                                        aaaaaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGc    <  1:508085/76‑1 (MQ=255)
                                                                         aaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCt   <  2:5414604/76‑1 (MQ=255)
                                                                         aaaaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCt   >  2:1019923/1‑76 (MQ=255)
                                                                          aaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCCGAACCATGCGCTGGTGCTGGCCGACTGGCTa  <  2:5798657/76‑1 (MQ=255)
                                                                          aaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa  <  1:400389/76‑1 (MQ=255)
                                                                          aaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa  <  2:6658218/76‑1 (MQ=255)
                                                                          aaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa  <  2:3857354/76‑1 (MQ=255)
                                                                          aaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa  <  2:329109/76‑1 (MQ=255)
                                                                          aaaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa  >  2:5540138/1‑76 (MQ=255)
                                                                          aaaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa  >  2:6799147/1‑76 (MQ=255)
                                                                          aaaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa  >  1:6190319/1‑76 (MQ=255)
                                                                             |                                                                        
CTTTCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTA  >  REL606/1854326‑1854473

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: