Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606402,0560T→G33.3% 12.566180ribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  new base (21/1):  ref base (17/27):  total (38/28)

CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACT  >  REL606/402021‑402091
                                   |                                   
ccACATGGCGAAGTACACGCGTTGCGTATGGCGGGt                                     <  2:3433825/36‑1 (MQ=255)
ccACATGCCGAAGTACACGCGTTGCGTATGGCGGGt                                     <  2:1181283/36‑1 (MQ=255)
ccACATGCCGAAGTACACGCGTTGCGTATGGCGGGt                                     <  2:2268032/36‑1 (MQ=255)
ccACATGCCGAACTACACGCGTTGCGTATGGCGGGt                                     <  1:3469740/36‑1 (MQ=255)
 cacaTGCCGAAGTACACGCGTTGCGTATGGCGGGTg                                    <  1:1115001/36‑1 (MQ=255)
 cacaTGCCGAAGTACACGCGTTGCGTATGGCGGGTg                                    <  2:2566973/36‑1 (MQ=255)
  ccaTGCCGAAGTACACGCGTTGCGTATGGCGGGTGa                                   <  2:376601/35‑1 (MQ=255)
  acaTGCCGAAGTACACGCGTTGCGTATGGCGGGTgg                                   >  1:3638531/1‑35 (MQ=255)
   caTGCCGAAGTACACGCGTTGCGTATGGCGGGTGaa                                  <  2:3473716/36‑1 (MQ=255)
    aTGCCGTAGTACACGGGTTGCGTATGGCGGGTGaaa                                 <  2:2483726/36‑1 (MQ=38)
    aTGCCGAAGTACACGCGTTGCGTATGGCGGGTGgaa                                 >  2:518885/1‑36 (MQ=255)
    aTGCCGAAGTACACGCGGTGCGTATGGCGGGTGaaa                                 <  2:1516661/36‑1 (MQ=255)
     tGCCGAAGTACACGCGTTGCGTATGGCGGGTGaaaa                                <  1:3261147/36‑1 (MQ=255)
     tGCCGAAGTACACGCGTTGCGTATGGCGGGGGgaaa                                >  2:1793494/1‑36 (MQ=38)
      gCCGAAGTACACGCGTTGCGTATGGCGGGGGaaaaa                               >  1:2290762/1‑36 (MQ=255)
       ccGAAGTACACGCTTTGCGTATGGCGGGTGAAAAAg                              <  1:2466786/36‑1 (MQ=255)
       ccGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAg                              <  2:79044/36‑1 (MQ=255)
       ccGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAg                              >  2:1054424/1‑36 (MQ=38)
       ccGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAg                              >  2:3450056/1‑36 (MQ=38)
        cGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGc                             <  2:1190673/36‑1 (MQ=255)
         gAAGTACACTCGTTGCGTATGGCGGGTGAAAAAGcc                            <  2:2234058/36‑1 (MQ=255)
         gAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGcc                            <  2:2565542/36‑1 (MQ=255)
         gAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGcc                            <  1:3542241/36‑1 (MQ=255)
         gAAGTACACGCGTTGCGTATGGCGGGGGaaaaaaac                            >  2:2122459/1‑33 (MQ=38)
         gAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGcc                            >  1:3359844/1‑36 (MQ=38)
         gAAGTACACGCGTTGCGTATGGCGGGGGAAAAAgac                            >  2:800619/1‑34 (MQ=38)
          aaGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCa                           <  1:3221843/36‑1 (MQ=255)
           aGTACACGCGTTGCGTATGGCGGGGGAAAAAACCaa                          >  1:862659/1‑36 (MQ=38)
           aGTACACGCGTTACGTATGGCGGGGGGAAAAGCCaa                          >  2:3606596/1‑36 (MQ=37)
            gttcacGCGTTGCGTATGGCGGGTGAAAAAGCCaaa                         <  1:3021180/33‑1 (MQ=255)
            gTACACGCGTTGCGTATGGCGGGGGGAAAAGCCcaa                         >  2:3666872/1‑36 (MQ=37)
            gTACACGCGTTGCGTATGGCGGGGGGAAAAGCCaaa                         >  2:568007/1‑36 (MQ=38)
            gTACACGCGTTGCGTATGGCGGGGGAAAAAGCCaaa                         <  2:1145639/36‑1 (MQ=255)
             tACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAg                        <  2:3613038/36‑1 (MQ=255)
             tACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAg                        <  1:1478127/36‑1 (MQ=255)
             tACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAg                        <  1:2018526/36‑1 (MQ=255)
             tACACGCGTTGCGTATGGCGGGGGAAAAAGGCAAAg                        >  2:3067614/1‑36 (MQ=37)
              acacGCGTTGCGTATGGCGGGGGGAAAAGTCAAAgg                       >  2:3095640/1‑36 (MQ=37)
              acacGCGTTGCGTATGGCGGGGGAAAAAGACaaacg                       >  1:749264/1‑34 (MQ=37)
               cacGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGt                      >  1:2465596/1‑36 (MQ=255)
               cacGCGTAGCGTATGGCGGGTGAAAAAGCCAAAGGt                      <  2:2949930/36‑1 (MQ=255)
                acGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTg                     <  2:3484349/36‑1 (MQ=255)
                acGCGTTGCGTATGGCGGGTGAAAAAGACAAAgggg                     >  2:2685094/1‑34 (MQ=38)
                acGCGTTGCGTATGGCGGGGGGAAAAGGCAAAgggg                     >  2:3022357/1‑34 (MQ=21)
                acGCGTTGCGTATGGCGGGGGAAAAAACCCACGGTg                     >  2:12183/1‑36 (MQ=19)
                   cgTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGa                  <  1:3702732/36‑1 (MQ=255)
                   cgTTGCGTATGGCGGGGGGAAAAGGCAAAGGGGCgg                  >  2:751820/1‑35 (MQ=16)
                   cgTTGCGTATGGCGGGGGAAAAAGCCAAAGGCGCgg                  >  1:469555/1‑35 (MQ=37)
                    gTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAc                 >  1:2556950/1‑36 (MQ=255)
                    gTTGCGTATGGCGGGGGGGAAAGCCAAAGGTGGGAc                 >  2:2336276/1‑36 (MQ=16)
                     ttGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAcc                <  1:2078410/36‑1 (MQ=255)
                     ttGCGTATGGCGGGTGAAAAAGCCAAAGGGGGGAcc                >  1:2272928/1‑36 (MQ=37)
                     tggCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAcc                <  1:407756/34‑1 (MQ=255)
                      tGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCg               >  1:2637011/1‑36 (MQ=255)
                      tGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCg               <  2:1643371/36‑1 (MQ=255)
                      tGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGAACg               >  2:1109222/1‑36 (MQ=37)
                        cGTATGGCGGGTGAAAAAACCAAAGGTGCGACCGcc             >  2:1167814/1‑36 (MQ=255)
                         gTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCt            >  1:2276194/1‑36 (MQ=255)
                         gTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCt            <  2:1275498/36‑1 (MQ=255)
                         gTATGGCGGGTGAAAAAGCCAAAAGTGCGAACGCCt            >  1:3210910/1‑36 (MQ=37)
                          tagggCGGGTGAAAAAGCCAAAGGTGCGACCGCCTa           <  2:1707569/33‑1 (MQ=255)
                          tATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTa           <  2:890670/36‑1 (MQ=255)
                          tATGGCGGGGGGAAAAACAAAAGGTGGGACCGGCTa           >  2:790008/1‑36 (MQ=12)
                           aTGGCGGGTGAAAAAGCCAAAAGTGCGACCGCCTac          >  1:223237/1‑35 (MQ=255)
                            tGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATg         >  1:2715832/1‑36 (MQ=255)
                            tGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATg         >  1:2386002/1‑36 (MQ=255)
                            tGGCGGGTGAAAAAGCCAAAGGAGGCACCCCCTATg         >  1:981093/1‑36 (MQ=16)
                               cGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTca      >  2:839249/1‑36 (MQ=255)
                                 ggTGAAAAAGCCAAAGGTGCGACCGCCCATGCcaca    >  1:2274198/1‑36 (MQ=37)
                                  gTGAAAAAGCCAAAGGTGCGACCGCCTATGTCacac   >  1:915935/1‑36 (MQ=255)
                                   tGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACt  >  1:107980/1‑36 (MQ=255)
                                   |                                   
CCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACT  >  REL606/402021‑402091

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 12 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG