Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,889,5840CA37.4% 155.9 / 82.6 91T114N (ACC→AAC) xdhCxanthine dehydrogenase, Fe‑S binding subunit
Reads supporting (aligned to +/- strand):  ref base C (15/42);  new base A (34/0);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 1.41e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 175 total aligned reads displayed.

CTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATC  >  REL606/2889509‑2889658
                                                                            |                                                                            
  cATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCa                                                                            <  1:2201228/76‑1 (MQ=255)
  cATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCa                                                                            <  1:1968239/76‑1 (MQ=255)
  cATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACa                                                                            >  1:3133285/1‑76 (MQ=255)
   aTGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAcacc                                                                           >  1:1110000/1‑73 (MQ=255)
   aTGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAacac                                                                           >  1:5280681/1‑76 (MQ=255)
   aTGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCAc                                                                           <  1:5631773/76‑1 (MQ=255)
   aTGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCAc                                                                           <  2:5943563/76‑1 (MQ=255)
   aTGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCAc                                                                           <  2:6334853/76‑1 (MQ=255)
   aTGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACAAc                                                                           >  1:5203479/1‑76 (MQ=255)
    tGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAacacc                                                                          >  2:2592416/1‑75 (MQ=255)
    tGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACg                                                                          >  1:7706946/1‑76 (MQ=255)
     gTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAacacag                                                                         >  1:2869390/1‑74 (MQ=255)
     gTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACgg                                                                         >  1:803335/1‑76 (MQ=255)
     gTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACgg                                                                         >  1:4522472/1‑76 (MQ=255)
        cagcagGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCCCAGCaa                                                                      >  1:2316581/1‑76 (MQ=255)
        cagcagGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACCCGGCaa                                                                      >  1:4067091/1‑76 (MQ=255)
        cagcagGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACGGCaa                                                                      >  1:565715/1‑76 (MQ=255)
        cagcagGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACGGCaa                                                                      >  2:5500383/1‑76 (MQ=255)
         agcagGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAAt                                                                     <  1:3173554/76‑1 (MQ=255)
          gcagGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACCGCCATg                                                                    >  1:623597/1‑76 (MQ=255)
          gcagGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAAAACGGCAATg                                                                    >  1:691513/1‑76 (MQ=255)
             ggCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACCGCAAAGCTg                                                                 >  1:4842590/1‑76 (MQ=255)
             ggCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACCCCTGGCCTGATTATGGCTACCACGGCAATGCTg                                                                 <  2:4856216/76‑1 (MQ=255)
              gCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACCGACATGCtcg                                                                >  2:2334961/1‑74 (MQ=255)
                ttATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACGGCAAtggtgggg                                                              >  1:2541822/1‑74 (MQ=255)
                ttATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAAACCAGCAATGCTgggg                                                              >  2:7592135/1‑74 (MQ=255)
                  aTGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTTGCGaa                                                            <  2:6706391/76‑1 (MQ=255)
                  aTGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGaa                                                            <  2:6792638/76‑1 (MQ=255)
                  aTGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAACGCCGGCGaa                                                            >  1:435592/1‑76 (MQ=255)
                  aTGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACGGCAATGCTGGCgga                                                            >  1:3035013/1‑74 (MQ=255)
                   tGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACGGCCAAGCCGGCGgaa                                                           >  1:3527659/1‑76 (MQ=255)
                    gCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAAc                                                          <  1:4533565/76‑1 (MQ=255)
                    gCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACAGCAAGGTCGGAGAAAc                                                          >  1:1262990/1‑76 (MQ=255)
                    gCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAAAACCGCAATGCTGGCGAAAc                                                          >  1:1531002/1‑76 (MQ=255)
                     cGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAAcc                                                         <  2:6097662/76‑1 (MQ=255)
                      gAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACGGCAATGCTGGCGaaaaac                                                        >  1:310618/1‑73 (MQ=255)
                       aaaTCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCAc                                                       <  2:1039904/76‑1 (MQ=255)
                       aaaTCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCAc                                                       <  2:6957274/76‑1 (MQ=255)
                       aaaTCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCAc                                                       <  2:920664/76‑1 (MQ=255)
                        aaTCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACAGCCATGATGGCGGAAacacc                                                      >  2:4227552/1‑75 (MQ=255)
                          tCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCAcgcg                                                    <  2:1762159/76‑1 (MQ=255)
                          tCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCAcgcg                                                    <  2:5025808/76‑1 (MQ=255)
                          tCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCAcgcg                                                    <  2:6369847/76‑1 (MQ=255)
                          tCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCAcgcg                                                    <  2:7491730/76‑1 (MQ=255)
                          tCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACAACGGCAATGCTGGGCAAAACAcgcg                                                    >  2:3252628/1‑76 (MQ=255)
                            cGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCgag                                                  <  1:556315/76‑1 (MQ=255)
                            cGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCgag                                                  <  1:6596887/76‑1 (MQ=255)
                             ggCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCgaga                                                 <  1:5497890/76‑1 (MQ=255)
                             ggCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCgaga                                                 <  1:1960563/76‑1 (MQ=255)
                              gcgcAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGaa                                                <  1:965242/76‑1 (MQ=255)
                              gcgcAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACGGGCATGCTGGCGCAACCACGCGAGaa                                                >  2:3553015/1‑76 (MQ=255)
                               cgcAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAg                                               <  1:4457821/76‑1 (MQ=255)
                                gcagtgcagtGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGc                                             <  2:3938398/76‑1 (MQ=255)
                                gcagtgcagtGCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACGGCAATGCTGGCGAAACCACGCGAGAAGc                                             >  1:1209533/1‑76 (MQ=255)
                                 cagtgcagtgCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACCGCCATGCTGGCGGAACCACGCGAGAAGcc                                            >  1:4240828/1‑76 (MQ=255)
                                 cagtgcagtgCGGGTTTTGTACGCCTGGCCTGATTATGGCTAACAAGGCAATGCTGGCGCAAACACGAgaggagac                                             >  1:1782607/1‑74 (MQ=255)
                                    tgcagtgcGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACCGCAATGCTGGGGAAACCACCCGAGGAGCCAtt                                         >  1:4295398/1‑76 (MQ=255)
                                    tgcagggcGGTTTTGTACCGCCTGCCTTGATTATGGCTACCACGGCAATGCTGGCAAAACCCCGCGAAAAGCCAtt                                         <  2:7696898/76‑1 (MQ=255)
                                      cagtggGGGTTGTATACGTGAGGCCTGATTATGCCTCCCACGGAAATGCTGGAAAACCCCCGCAAGAAGCCATTaa                                       <  2:3930732/76‑1 (MQ=255)
                                      cagtgcGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACGGCCAAGCTGGCGGAACCCAGCGAGAAACCATTaa                                       >  1:3274531/1‑76 (MQ=255)
                                      cagtgcGGGTTTTGTACGCCTGGCCTGATTATGGCTAACACCGCAATGCTGGCGGAACCCCGCCAGAAGGCAATaa                                       >  2:1095205/1‑76 (MQ=255)
                                       ggtgcGGTTTTTGTCCGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAAc                                     <  2:5406474/75‑1 (MQ=255)
                                         tgcGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCa                                   <  1:874737/76‑1 (MQ=255)
                                         ggcGGGTTTTGTACGCCGGGCCTGATTTTGGCTACCACGGAAATGCTGGCGAAACCACGCGAGAAGCCATTAACCa                                   <  2:3505491/75‑1 (MQ=255)
                                          gcGGGTTTTGTACGCCTGGCCTGATTATGGCTAACAACGCAAAGCTGGCGGAAACACACGCGGAGCCCTTTAACCa                                   >  1:4994063/1‑76 (MQ=255)
                                            gggTTTTGTCCCCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTa                                <  1:4237055/76‑1 (MQ=255)
                                             ggTTTTGTCCGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTAc                               <  1:705114/76‑1 (MQ=255)
                                              gTTTTGTACGCCTGGCCTGATTATGGCTAACACGGCAAAGCTGGCGGAACCACGCGAGAAGCCATTAAACATTACg                              >  1:643158/1‑76 (MQ=255)
                                               ttttGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACgg                             <  1:243781/76‑1 (MQ=255)
                                               ttttGTACGCCTGGCCTGATTATGGCTAACACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAAACATTACgg                             >  1:5811961/1‑76 (MQ=255)
                                                 ttGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCGGGCGAAACCACGCGAGAAGCCATTAACCATTACGGaa                           <  1:5409595/76‑1 (MQ=255)
                                                  tGTACGCCTGGCCTGATTATGGCTAACACCGGAATGCTGGCGGAAACACGCGAGAAACCATTAACCCATTCCGGaa                           >  2:2412947/1‑76 (MQ=255)
                                                  tGTACGCCTGGCCTGATTATGGCTAACACCGCAAAGCTGGCGGAAACACGCGAGAAGCCATTAAACAATACAGaaa                          >  2:4499767/1‑76 (MQ=255)
                                                    tACGCCTGGCCTGATTATGGCTAACACCGCAAAGCTGGGGAAACCCCGCGAGAAGCCAATAAACATTACGGaaaat                        >  2:1453247/1‑74 (MQ=255)
                                                    tACGCCTGGCCTGATTATGGCTAAAAAGCCACAGTTGGGGACAAAACCCGAGAAGCCATTAACCCATACTGAAAtt                        >  2:1533281/1‑76 (MQ=255)
                                                    gtccGCCTGCCCTGATTATGGCTACCACGGAAATGCTGGCGAAACCACGCGAAAAGCCATTAACCATTACGGAAAt                         <  1:3364183/73‑1 (MQ=255)
                                                         cctggcctgATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATtcgtc                    <  1:3872197/76‑1 (MQ=255)
                                                         cctggcctgATTAGGGCTACCACGGCAATGCTGGCGAACCCACGCGAGAAGCCATTAACCATTACGGAAATtcgtc                    <  1:6057219/76‑1 (MQ=255)
                                                           tggcctgATTATGGCTAACACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTcgc                  >  1:560150/1‑76 (MQ=255)
                                                            ggcctgATTATGGCTAACACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTcgcg                 >  1:7568560/1‑76 (MQ=255)
                                                            ggcctgATTATGGCTAACACGGCAAAGCTGGCGAAACCACGCGAGAAGCCATTAAACATTACGGAAAATTGTcgcg                 >  2:1961631/1‑76 (MQ=255)
                                                            ggcctgATTATGCCTACCACGGCAATGCTGGCAAAACCACGGAAGAAGCCATTAACCATTACGGAAATTCGTcgcg                 <  2:7316847/76‑1 (MQ=255)
                                                              ctggATTAGGGCTCCCACGGCAATGCTGGCGAAACCACGCGAAAAGCCATTAACCATTACGGAAATTCGTCGCGGa               <  2:576124/73‑1 (MQ=255)
                                                               ctgTTTAGGGCTACCACGGAAATGCTGGCGAAACCCCGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGAc              <  1:6987602/76‑1 (MQ=255)
                                                               ctgATTATGGCTACCACGGCAATGCTGCCAAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGAc              <  1:6529405/76‑1 (MQ=255)
                                                                ttaTTATGCCTACCACGGCAATGCTGGCGAACCCCCGCGAAAAGCCATTAACCATTACGGAAATTCGTCGCGGACt             <  1:7209010/74‑1 (MQ=255)
                                                                tgATTATGCCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACt             <  1:3925552/76‑1 (MQ=255)
                                                                 gATTATGGCTACCACGGGAATGCTGGCGAAACCACGCGAGAAGCCCTTAACCATTACGGAAATTCGTCGCGGACTg            >  1:7804868/1‑76 (MQ=255)
                                                                 gATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTg            <  1:3589400/76‑1 (MQ=255)
                                                                  aTTATGGCTAACACCGCAATGCTGGCGGAAACCCCCGAGAAGGCATTAACCATTACGGAAAATCGTCGCGGGCTgg           >  2:6463922/1‑76 (MQ=255)
                                                                   ttATGGCTAACACGGCAATGCTGGCGGAAACACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGGCTGGc          >  2:3654003/1‑76 (MQ=255)
                                                                   ttATGGCTAACACGGCAAAGCTGGCGAAAACACACGAGAAACCACTAACCATTACGGAAAATCGTCCCGGGCTggg          >  2:4708298/1‑75 (MQ=255)
                                                                     aTGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCgg        >  2:6321318/1‑76 (MQ=255)
                                                                      tGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCggg       >  1:7698410/1‑76 (MQ=255)
                                                                       ggCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCgggg      >  1:374628/1‑75 (MQ=255)
                                                                       ggCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCgggg      >  2:3397365/1‑75 (MQ=255)
                                                                        ggtACCCCGGCAATGCTGGCGAATCCACGTGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGaa     <  1:1610325/74‑1 (MQ=255)
                                                                         ctcccACGGCAATGCGGGCAAAACCACGGGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGaaa    <  1:6720640/73‑1 (MQ=255)
                                                                         cTACCACGGCAATGCGGGCAAAACCACGCAAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGaaa    <  1:1980585/76‑1 (MQ=255)
                                                                           aCCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGgaaga  >  1:1041155/1‑74 (MQ=255)
                                                                            |                                                                            
CTCATGTTCAGCAGGCTTATGCGAAATCCGGCGCAGTGCAGTGCGGGTTTTGTACGCCTGGCCTGATTATGGCTACCACGGCAATGCTGGCGAAACCACGCGAGAAGCCATTAACCATTACGGAAATTCGTCGCGGACTGGCGGGAAATC  >  REL606/2889509‑2889658

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