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

CTTTCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTACG  >  REL606/1854326‑1854475
                                                                             |                                                                          
cTTTCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGGAGCGCCacaaaa                                                                            >  1:3046670/1‑76 (MQ=255)
cTTTCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCaaaaaa                                                                            <  2:1361630/76‑1 (MQ=255)
 gttCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAAc                                                                           <  2:1743238/75‑1 (MQ=255)
  ttCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGGGGGGAACGCCacaaaaaa                                                                          >  2:1986975/1‑74 (MQ=255)
   tCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAg                                                                         <  1:3885678/76‑1 (MQ=255)
      cgACCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCg                                                                      <  2:922673/76‑1 (MQ=255)
         tcAATATCTTGAAGAACTGGCGCAACCCGAGCAAGAAGGAAGGCGGCGGGCGGACACCGCAAAAAAACAGCCGCTg                                                                   <  2:2578884/75‑1 (MQ=255)
             tatCTTGCAGAACTGGCGCAACACGAGAAAGAAGGAAGGCGGCTGGCGGCCATCGCCAAAAACCAGCCGCTGAccc                                                               <  2:3255458/76‑1 (MQ=255)
              atCTTGTAGAACGGGCGCAACACGAGAAAGAAGGAAAGCGGCTGGCGGAAACCGCCAAAAACCAGCCGCTGACCCt                                                              <  2:3486185/76‑1 (MQ=255)
               tCTTGCAAAACTGGCGAAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTg                                                             <  1:1807412/76‑1 (MQ=255)
                   gCAGAACGGGCGCAACCGGAAAAAGAAGAAAGGCGGCTGGCGAACATCCCAAAAAACCAGCCGCTGACCCTGctct                                                         <  2:4627933/76‑1 (MQ=255)
                   gCAGAACGGGCGCAACACGAGAAAGAAGGAAAGCGGCGGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGctct                                                         <  2:1505596/76‑1 (MQ=255)
                    cAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTa                                                        <  1:523320/76‑1 (MQ=255)
                    cAGAACTGGCGCAACACGAGAAAGAAGGAAAGCGGCGGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTa                                                        <  2:3843336/76‑1 (MQ=255)
                    cAGAACTGGCGAAACACGAGAAAGAGGGAAGGCGGCGGGCGGACATCCCAAAAAACCAGCCGCTGACCCTGCTCTa                                                        <  2:2833392/76‑1 (MQ=255)
                    aataaCGGGCGCAACACGAGCAAGAGGAAAAGCGGCTGGCGGAAATCGCAAAAAAACAGCCGCTGACCCTGCTCTa                                                        <  2:1080690/73‑1 (MQ=255)
                     aGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTAc                                                       <  1:2552163/76‑1 (MQ=255)
                      gAACTGGCGCAACCCGACCAAGAAGGAAAGCGGCTGGCGAACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACg                                                      <  2:2111396/76‑1 (MQ=255)
                       aaCTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACgc                                                     <  2:4701250/76‑1 (MQ=255)
                       aaCGGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACgc                                                     <  2:2905180/76‑1 (MQ=255)
                          gggCGAACCACGAGAAAGAAGAAAAGCGGCTGGGGGACATCCCAAAAAACCAGCCGCTGACCCTGCTCTACgcagc                                                  <  2:411487/75‑1 (MQ=255)
                             cgcACCACGAGAAAGAAGGAAAGGGGCTGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCaaa                                               <  1:4046540/76‑1 (MQ=255)
                             cgcAACACGGAAAAGAAGATAGCGGTGTGCCGGACTGCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCaaa                                               <  2:1568639/76‑1 (MQ=255)
                             cgcAACACGAGCAGGAAGGAAAGCGGCGGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCaaa                                               <  1:131794/76‑1 (MQ=255)
                               cAACACGAGCAAGAAGGAAAGCGGCTGGGGGGCATCGCCCAAAAAAAACCGCCCACCCCCCCCTCCGCAGCaaaaa                                             >  1:3950943/1‑76 (MQ=255)
                               cAACACGAGCAAGAAGGAAAGCGGCTGGGGGGAATCGACCAAAAAAAAACGGCCGCCCCGCTCCACGCAGCaaaaa                                             >  2:4904270/1‑76 (MQ=255)
                               cAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCaaaaa                                             <  1:4716045/76‑1 (MQ=255)
                               cAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCaaaaa                                             <  1:3475999/76‑1 (MQ=255)
                                 acccGAGCAAGAAGGAAAGCGGCTGGCGGCCATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAac                                           <  1:4639516/73‑1 (MQ=255)
                                    gaaGCAAAAAGAAAGACGGCTGGCGGACACCCCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAAcacc                                        <  2:3967225/74‑1 (MQ=255)
                                    cGAGCAAGAAGGAAGGCGGCTGGCGGACATCGCAAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAAcacc                                        <  1:2018079/76‑1 (MQ=255)
                                    cGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAAcacc                                        <  1:2695950/76‑1 (MQ=255)
                                     gagaAAGAAGGAAAGCGGCGGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAAcacca                                       <  2:1197194/76‑1 (MQ=255)
                                       gCAAGAAGGAAAGGGGCTGGCGGACACCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACg                                     <  1:2795091/76‑1 (MQ=255)
                                           gaagGAAAGCGGCTGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGa                                 <  1:351915/76‑1 (MQ=255)
                                           gaagGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGa                                 <  1:4998108/76‑1 (MQ=255)
                                           gaagGAAAGCGGCTGGCGGACATCGCAAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGa                                 <  1:1755319/76‑1 (MQ=255)
                                            aagGAAAGCGGCTGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGaa                                <  1:4665477/76‑1 (MQ=255)
                                              ggAAGGCGGCTGGCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAAcc                              <  2:1288200/76‑1 (MQ=255)
                                              gaaaaGCGGCTGGGGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAAcc                              <  2:215413/74‑1 (MQ=255)
                                               gAAGGCGCCGGGCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCa                             <  1:377469/76‑1 (MQ=255)
                                               gAAAGCGGCGGCCGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCa                             <  2:2895650/76‑1 (MQ=255)
                                                  aGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATgc                          <  1:97648/76‑1 (MQ=255)
                                                     ggCTGGCGGGCATCGCCAAAAAAAAGCCGCTGACCCCGCTCTTCGCAGCCAAAAAAACCACGCAGAACCATGCGCt                       >  2:971449/1‑76 (MQ=255)
                                                      gCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTg                      >  2:1218452/1‑76 (MQ=255)
                                                      gCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTg                      >  2:1879421/1‑76 (MQ=255)
                                                      gCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAAACATGCGCTg                      >  2:1383146/1‑76 (MQ=255)
                                                      gCTGGCGGACATCCCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTg                      <  1:4552072/76‑1 (MQ=255)
                                                      gCTGGCGGACACCGCAAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTg                      <  2:749801/76‑1 (MQ=255)
                                                       gcgggCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTg                      <  2:34621/73‑1 (MQ=255)
                                                       gcgggCGGACATCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTg                      <  2:4286291/73‑1 (MQ=255)
                                                       gcgggCGAACACCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTg                      <  2:2619858/73‑1 (MQ=255)
                                                        cTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTgg                     >  2:495232/1‑76 (MQ=255)
                                                        cTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTgg                     >  2:1294201/1‑76 (MQ=255)
                                                        cTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTgg                     >  2:1292791/1‑76 (MQ=255)
                                                        cTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTgg                     >  1:4968139/1‑76 (MQ=255)
                                                         tGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGTGCTGGt                    >  1:3185966/1‑76 (MQ=255)
                                                         tGGCGGACACCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGt                    <  2:4241732/76‑1 (MQ=255)
                                                         gggCGGACACCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGt                    <  2:141373/75‑1 (MQ=255)
                                                          ggCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTg                   >  2:2968479/1‑76 (MQ=255)
                                                           gCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGc                  >  1:4872279/1‑76 (MQ=255)
                                                            gggACACCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt                 <  2:3191875/75‑1 (MQ=255)
                                                            cGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt                 <  1:863648/76‑1 (MQ=255)
                                                            cGGACATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt                 <  1:3762383/76‑1 (MQ=255)
                                                            cGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt                 >  2:1024504/1‑76 (MQ=255)
                                                            cGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCt                 >  2:2905020/1‑76 (MQ=255)
                                                              gACACCGCAAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTgg               <  1:155352/76‑1 (MQ=255)
                                                               tcATCGCAAAAACCCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGc              <  2:2628727/75‑1 (MQ=255)
                                                               aaaTCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGc              <  2:2656160/74‑1 (MQ=255)
                                                               aCATCGCCAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGc              <  2:1387937/76‑1 (MQ=255)
                                                                 accGCAAAAACCCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCg            <  2:380010/74‑1 (MQ=255)
                                                                 aTCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCg            >  2:685712/1‑76 (MQ=255)
                                                                 aTCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCg            >  2:1326471/1‑76 (MQ=255)
                                                                  tCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGa           >  2:3181079/1‑76 (MQ=255)
                                                                  tCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGa           >  2:4670302/1‑76 (MQ=255)
                                                                  tCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGa           >  2:167099/1‑76 (MQ=255)
                                                                      ccAAAAACCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTg        <  2:4828988/76‑1 (MQ=255)
                                                                      ccAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTg        >  1:3975901/1‑76 (MQ=255)
                                                                      ccAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTg        >  2:1664465/1‑76 (MQ=255)
                                                                      caaaaaaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTg        <  2:227976/74‑1 (MQ=255)
                                                                       cAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTgg       >  1:1194133/1‑76 (MQ=255)
                                                                         caaacCAGCCGCGAACCCTGCTCTACGCAGCAAAAACCACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCt     <  2:560506/75‑1 (MQ=255)
                                                                         aaacaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCt     <  1:227735/76‑1 (MQ=255)
                                                                         aaaacCAGCCGCTGATCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCt     <  2:1922998/76‑1 (MQ=255)
                                                                         aaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCt     <  2:1447327/76‑1 (MQ=255)
                                                                         aaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCt     <  2:1745754/76‑1 (MQ=255)
                                                                         aaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCt     <  1:1790209/76‑1 (MQ=255)
                                                                         aaaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCt     <  2:560505/76‑1 (MQ=255)
                                                                         aaaaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCt     >  1:1220412/1‑76 (MQ=255)
                                                                          aaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa    <  2:4972543/76‑1 (MQ=255)
                                                                          aaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa    <  2:2692766/76‑1 (MQ=255)
                                                                          aaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa    <  2:2290011/76‑1 (MQ=255)
                                                                          aaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa    <  1:3870420/76‑1 (MQ=255)
                                                                          aaacCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa    <  1:2369462/76‑1 (MQ=255)
                                                                          aaaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa    >  2:3426992/1‑76 (MQ=255)
                                                                          aaaaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTa    >  2:1528587/1‑76 (MQ=255)
                                                                            aaCAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTaag  >  2:1739447/1‑74 (MQ=255)
                                                                             |                                                                          
CTTTCGCGAGCAATATCTTGCAGAACTGGCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTACGCAGCAAAAAACACCACGCAGAACCATGCGCTGGTGCTGGCCGACTGGCTACG  >  REL606/1854326‑1854475

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