Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,854,4010AC44.0% 87.2 / 27.1 75K88N (AAA→AACyeaOhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (15/27);  new base C (0/33);  total (15/60)
Fisher's exact test for biased strand distribution p-value = 5.88e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 115 total aligned reads displayed.

CTTTCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTA  >  REL606/1854326‑1854473
                                                                            |                                                                         
ctgtCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCaaaaaa                                                                           <  2:475233/73‑1 (MQ=255)
 tttCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGGGTGGCCaaaaaaa                                                                          >  1:5700950/1‑75 (MQ=255)
 tttCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAAc                                                                         <  1:993996/76‑1 (MQ=255)
 tttCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAAc                                                                         <  1:2900167/76‑1 (MQ=255)
  ttCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACa                                                                        <  1:7996846/76‑1 (MQ=255)
    cgcgAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGc                                                                      <  1:1698280/76‑1 (MQ=255)
         gCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTg                                                                 <  1:3574431/76‑1 (MQ=255)
           aaTATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGCACACCCAAAAAAAAACCGCCGAc                                                               >  1:8079821/1‑76 (MQ=255)
           aaTATCTTGCAGAACGGGCGAAACACGAGCAAGAAGAAAGGGGTCGGGCGGACATCGCAAAAAACCAGCCGCTGAc                                                               <  2:2865468/76‑1 (MQ=255)
           aaTATCTTGAAGACCTGGCGCAACACGAGAAAGAAGGAAAGGGGCTGGCGGACATCGCCAAAAACCAGCCGCTGAc                                                               <  2:2689205/76‑1 (MQ=255)
             tatCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGAccc                                                             <  1:832573/76‑1 (MQ=255)
              tccTTGCAGCCCGGGCGGAAAACGACCAAGAAGGAAAGCGGTGGGCGGAAAGCGCAAAAAACCAGCCGCTGACCCt                                                            <  2:5861525/74‑1 (MQ=255)
              atCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCt                                                            <  1:1490487/76‑1 (MQ=255)
              atCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCt                                                            <  1:3911003/76‑1 (MQ=255)
              atCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCt                                                            <  1:7957931/76‑1 (MQ=255)
               tCTTGCAGAACGGGCGCAACACGAGCAAGAAGAAAAGCGGCGGGCGGACATCGCAAAAAACCAGCCGCTGACCCTg                                                           <  2:2936530/76‑1 (MQ=255)
                   gCAGAACTGGCGCAACACGAGCAAGAAGAAAGGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGctct                                                       <  2:6736090/76‑1 (MQ=255)
                    cAGACCTGGCGCAACAGGAGAAAGAAGGAAGGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTa                                                      <  1:7613149/76‑1 (MQ=255)
                     aGAACGGGCGCAACACGAGCAAAAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTAc                                                     <  1:1009726/76‑1 (MQ=255)
                      gAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACg                                                    <  2:7370990/76‑1 (MQ=255)
                      gAACTGGCGAAACACGAGCAAGAAGGAAGGCGGCGGGCGGACCTCGCAAAAAACCAGCCGCTGACCCTGCTCTACg                                                    <  2:155535/76‑1 (MQ=255)
                      aaaCGGGCGCAACACGAGAAAGAAGGAAAGCGGCTGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACg                                                    <  2:214739/75‑1 (MQ=255)
                       aaCGGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgc                                                   <  1:2933049/76‑1 (MQ=255)
                         caggCGCAAGACGAGCAAGAAGGAATGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgcag                                                 <  1:7278023/74‑1 (MQ=255)
                         cTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGCATCGCCCAAAAACAAACCCCGCGCCTGCGCTACgcag                                                 >  1:8250232/1‑76 (MQ=255)
                          tGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgcagc                                                <  1:8107696/76‑1 (MQ=255)
                           ggCGCAACACGAGCAAGAAGAAAAGCGGCTGGCGAACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgcagca                                               <  1:5250415/76‑1 (MQ=255)
                             cgcAACAGGAGCAAGAAGGAAGGCGGCTGGCGGACACCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCaaa                                             <  2:655166/76‑1 (MQ=255)
                             cgcAACACGAGCAAGAAGGAAAGCGGCGGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCaaa                                             <  1:5449943/76‑1 (MQ=255)
                             cgcAACACGAGAAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCaaa                                             <  2:5172817/76‑1 (MQ=255)
                             cacAACACGAGCAAGAAGGAAAGGGGCGGGCGGACATCACAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCaaa                                             <  2:3120724/74‑1 (MQ=255)
                              gcAACACGAGCAAGAAGGAAAGCGGCTGGGGGGGAACGCCCAAAAACAAACGCTGCCCCCCCCCTCCGCCGCaaaa                                            >  1:2686482/1‑76 (MQ=255)
                               catcacGAGCAAGAAGGAAAGCGGCTGGCGAACTTCGCAAAAAAACAGCCGCTGACCCTGCTCTACGCAGCaaaaa                                           <  2:7288578/73‑1 (MQ=255)
                                  cccGAGAAAAAGGGAAAGGGGCGGGCGGACATCCCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAaca                                        <  1:4925349/74‑1 (MQ=255)
                                  cacGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAaca                                        <  1:2674913/76‑1 (MQ=255)
                                   acGAGCAAGAAGGAAAGCGGCTGGGGGGCATCGCCAAAAAACAACCGCCCACCCCCCTCTGCCCACCAAAAAaaac                                       >  2:4923207/1‑76 (MQ=255)
                                   acGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAacac                                       <  1:5394183/76‑1 (MQ=255)
                                    cgcgCAAGAAGGAAAGCGGCGGGCGGACATGCCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAAcacc                                      <  2:1070333/73‑1 (MQ=255)
                                     gagaAAGAAGGAAAGCGGCGGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAAcacca                                     <  2:1306139/76‑1 (MQ=255)
                                     gAGCAAGAAGGAAGGCGGCTGGCGGACATCGCAAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAAcacca                                     <  2:1039345/76‑1 (MQ=255)
                                      ggaaagaagGAACGCGGCTGGCGGAACCCGAAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAAcaccac                                    <  2:454344/73‑1 (MQ=255)
                                      gacAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAAcaccac                                    <  1:5924980/74‑1 (MQ=255)
                                      aggaagaagGGAGGGGGGTGGCCGACGTGGCAGAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAAcaccac                                    <  2:6011758/73‑1 (MQ=255)
                                       gaaagaagGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACg                                   <  1:1160820/74‑1 (MQ=255)
                                       gCAAGAAGGAAAGCGGCTGGGGGGCCACGCCCAAAAAAAACCCCCGGCCCCGCCCTACCGCGGAAAAAAaaccaca                                   >  1:5628817/1‑75 (MQ=255)
                                       gCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACg                                   <  1:1454721/76‑1 (MQ=255)
                                        cAAGAAGAAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGc                                  <  1:1678892/76‑1 (MQ=255)
                                        aaagaggGAAAGCGCCGGGCGGACATCCCAAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGc                                  <  2:7469100/75‑1 (MQ=255)
                                        aaagaagAAAGGGGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGc                                  <  2:1499938/75‑1 (MQ=255)
                                         aagaagGAAAGCGGCTGGGGGGCAGCGCCAAAAAAAAGCCGCCGACCCCGCTCTTCCCCGCCAAAAAcaacaccga                                 >  1:6214555/1‑73 (MQ=255)
                                         aagaagGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCa                                 <  1:5861210/76‑1 (MQ=255)
                                         aagaagGAAAGCGGCGGCCGAACACCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCa                                 <  1:579434/76‑1 (MQ=255)
                                           gaagaaAGGCGGCTGCCGGCAGACGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGa                               <  2:1567258/76‑1 (MQ=255)
                                              ggAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAAcc                            <  1:3400463/76‑1 (MQ=255)
                                               aaaaGCGGCTGGCGGACATCGCAAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCa                           <  2:3098303/75‑1 (MQ=255)
                                                aaaGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCAt                          <  1:2918591/76‑1 (MQ=255)
                                                aaaGCGGCGGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCAt                          <  1:8101516/76‑1 (MQ=255)
                                                aaaGCGGCGGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCAt                          <  1:8091319/76‑1 (MQ=255)
                                                 acgCGGCTGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATg                         <  2:478762/74‑1 (MQ=255)
                                                  ggcggcgGGCGGACATCGCAAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATgc                        <  2:7273242/75‑1 (MQ=255)
                                                   gcggcgGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATgcg                       <  2:3472772/76‑1 (MQ=255)
                                                   gcggcgGCCGAACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATgcg                       <  2:3097083/76‑1 (MQ=255)
                                                   gCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATgcg                       <  2:3886996/76‑1 (MQ=255)
                                                     ggCTGGCGGACATCGCCCAAAAACACCCGCTTAACCTTTTCTTGCCACACAAAAACACCACGCAGAAACATgcccc                     >  2:6076041/1‑73 (MQ=255)
                                                      gcgggCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTg                    <  1:2670956/73‑1 (MQ=255)
                                                       cTGGGGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTgg                   <  1:5838668/76‑1 (MQ=255)
                                                        tGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGt                  >  2:5688287/1‑76 (MQ=255)
                                                        tGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGt                  >  2:4525720/1‑76 (MQ=255)
                                                        gggCGAACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGt                  <  2:2676928/75‑1 (MQ=255)
                                                         ggCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTg                 >  2:2685668/1‑76 (MQ=255)
                                                         ggCGGACATCGCCAAAAAACAGACGCTGACCCTGTTCTACGCAGTAAAAAACACCACGCAGAACCATGCGCTGGTg                 >  2:1482133/1‑76 (MQ=255)
                                                         ggCGGACATCCCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTg                 <  1:2448570/76‑1 (MQ=255)
                                                          ggggACCTCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGc                <  2:6567163/74‑1 (MQ=255)
                                                          gCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGc                <  1:4550899/76‑1 (MQ=255)
                                                          gCGGACATCCCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGc                <  2:4888939/76‑1 (MQ=255)
                                                           cGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt               >  1:3117228/1‑76 (MQ=255)
                                                           cGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt               >  2:56619/1‑76 (MQ=255)
                                                           cGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt               <  1:7259852/76‑1 (MQ=255)
                                                            ggACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTg              <  2:2128770/76‑1 (MQ=255)
                                                            gcacacCGCAAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTg              <  1:4671967/74‑1 (MQ=255)
                                                              ccATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGc            <  2:7600960/75‑1 (MQ=255)
                                                               cATCGCCAAAAACCCGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGcc           <  1:2080313/76‑1 (MQ=255)
                                                               cATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGcc           >  1:1254937/1‑76 (MQ=255)
                                                                accGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCg          <  1:6795871/74‑1 (MQ=255)
                                                                aTCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCCGAACCATGCGCTGGTGCTGGCCg          >  1:2240594/1‑76 (MQ=255)
                                                                aTCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCg          >  1:451307/1‑76 (MQ=255)
                                                                  cGCCAAAACATCAGCCGCTGACCCTGCTCTACGCAGCAAAAACCGCCACGCAGAACCATGCGCTGGTGCTGGCCGa         <  2:469124/76‑1 (MQ=255)
                                                                   gcaaaaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACt       <  2:4007538/73‑1 (MQ=255)
                                                                   ccaaaaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACt       <  2:5702860/73‑1 (MQ=255)
                                                                    ccAAAACCCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCGGGTGCTGGCCGACTg      <  2:7466432/76‑1 (MQ=255)
                                                                    ccAAAAACCCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACGCCACGCAGAACCATGCGCTGGTGCTGGCCGACt       <  1:6819979/76‑1 (MQ=255)
                                                                    ccAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTg      <  2:360412/76‑1 (MQ=255)
                                                                    ccAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTg      >  2:4147566/1‑76 (MQ=255)
                                                                     cAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCTGACTgg     >  1:2000172/1‑76 (MQ=255)
                                                                      aaaaaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGc    >  2:357965/1‑76 (MQ=255)
                                                                       aaaacCAGCCGCTTACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCt   <  2:6895467/76‑1 (MQ=255)
                                                                       aaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCt   <  1:6251652/76‑1 (MQ=255)
                                                                       aaaaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCt   >  2:4747752/1‑76 (MQ=255)
                                                                        aaacCAGCCGCTGGCCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa  <  2:1306871/76‑1 (MQ=255)
                                                                        aaaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa  >  2:6330769/1‑76 (MQ=255)
                                                                            |                                                                         
CTTTCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTA  >  REL606/1854326‑1854473

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