Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 66,630 A→G L652P (CTA→CCA)  polB ← DNA polymerase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL60666,6300AG98.9% 349.9 / ‑4.6 88L652P (CTA→CCA) polBDNA polymerase II
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (16/71);  total (16/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 101 total aligned reads displayed.

ACCCGCCATCAGTTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATAGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCTTTAAACACCATCC  >  REL606/66557‑66702
                                                                         |                                                                          
ccccGCCATCAGTTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGc                                                                          <  2:4438687/75‑1 (MQ=255)
  ccGCCATCAGTTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTc                                                                        <  1:77667/76‑1 (MQ=255)
   cGCCATCAGTTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCCCTCCGGGAGGTGCGCAGGGATGGCTcc                                                                       >  1:1411361/1‑76 (MQ=255)
    gCCATCAGTTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCct                                                                      <  2:4729896/76‑1 (MQ=255)
    gCCATCAGTTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCct                                                                      <  2:4346222/76‑1 (MQ=255)
    gCCATCAGTTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCct                                                                      <  1:5077903/76‑1 (MQ=255)
     ccATCAGTTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCctg                                                                     <  1:5542267/76‑1 (MQ=255)
      cATCAGTTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCctgc                                                                    <  2:4732790/76‑1 (MQ=255)
       aTCAGTTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCctgct                                                                   <  1:3793809/76‑1 (MQ=255)
        tCAGTTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCctgctg                                                                  <  2:5045179/76‑1 (MQ=255)
        tCAGGTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCctgctg                                                                  <  2:705070/76‑1 (MQ=255)
          aGTTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGaa                                                                <  1:1704558/76‑1 (MQ=255)
           gTTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGaaa                                                               <  2:1601618/76‑1 (MQ=255)
           gTTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGaaa                                                               <  1:3079785/76‑1 (MQ=255)
             ttGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAAct                                                             <  1:4671262/76‑1 (MQ=255)
              tGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAActg                                                            <  2:2750596/76‑1 (MQ=255)
               gTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAActgc                                                           <  1:857143/76‑1 (MQ=255)
                   aTGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTggg                                                       <  2:5530141/76‑1 (MQ=255)
                     ggTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGcc                                                     <  1:843870/76‑1 (MQ=255)
                     ggTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGcc                                                     <  2:1412686/76‑1 (MQ=255)
                     ggTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGcc                                                     <  2:4079151/76‑1 (MQ=255)
                      gTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCa                                                    <  2:3533646/76‑1 (MQ=255)
                      gTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCa                                                    <  1:2053337/76‑1 (MQ=255)
                      gTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCa                                                    <  2:5130689/76‑1 (MQ=255)
                       tttCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAg                                                   <  2:4033427/76‑1 (MQ=255)
                       tttCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAg                                                   <  1:613177/76‑1 (MQ=255)
                       tttCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAg                                                   <  1:5167260/76‑1 (MQ=255)
                        ttCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAgc                                                  <  2:1029715/76‑1 (MQ=255)
                        ttCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAgc                                                  <  2:3658353/76‑1 (MQ=255)
                        ttCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAgc                                                  <  2:2232528/76‑1 (MQ=255)
                         tCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAgcg                                                 <  1:3121674/76‑1 (MQ=255)
                         tCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAgcg                                                 <  1:3722133/76‑1 (MQ=255)
                          cgcgTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAgcgg                                                <  2:6196673/76‑1 (MQ=255)
                             gTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGt                                             <  1:4563516/76‑1 (MQ=255)
                              tACATATTCCTGATATGGCTCCCTGCGGGAGATGGGGAGGGATAGATCCCTCCCGAACAGCTGGGGCCgcggccgc                                            >  2:2466588/1‑73 (MQ=255)
                               tcATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTcc                                           <  2:1113896/75‑1 (MQ=255)
                               aCATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTcc                                           <  2:5220056/76‑1 (MQ=255)
                               aCATATTCCTGATATGGCTCGTGGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACGGCTGGGCCAGCGGCGTcc                                           <  2:4300658/76‑1 (MQ=255)
                               aCATATTCCTGATATGGCTCCCTGCGGGAGATGGGCAGGGATGGGTCCCCCCGCAACGGCTGGGGCGGCGGCGGcc                                           >  1:3512636/1‑76 (MQ=255)
                                cATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCa                                          <  2:6316731/76‑1 (MQ=255)
                                 atatTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAg                                         <  1:3861935/76‑1 (MQ=255)
                                 atatTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAg                                         <  2:6082808/76‑1 (MQ=255)
                                 atatTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAg                                         <  1:6317300/76‑1 (MQ=255)
                                  tatTCCTGATATGGCTCGTGGCGGAAGATGCTCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGt                                        <  1:1253217/76‑1 (MQ=255)
                                   atTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTc                                       <  1:5669052/76‑1 (MQ=255)
                                   atTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTc                                       <  1:5560355/76‑1 (MQ=255)
                                    ttCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATAGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCg                                      <  1:2795540/76‑1 (MQ=255)
                                     tCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCgg                                     <  1:5176947/76‑1 (MQ=255)
                                      ccTGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGt                                    <  2:28622/76‑1 (MQ=255)
                                        tGATATGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTgc                                  <  2:2958521/76‑1 (MQ=255)
                                          atgtGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTgcgc                                <  2:4929595/73‑1 (MQ=255)
                                           tatGGCTCGTTGGGGAGTATGCGCTGGTAGCGCTCCTACTGAAACGGCTGGCCCAGCGGCGTCCAGTCGGTGCGCa                               <  2:4321410/76‑1 (MQ=255)
                                           tatGGCTCGTGGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCa                               <  1:5238580/76‑1 (MQ=255)
                                            atGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCAc                             <  1:4326323/76‑1 (MQ=255)
                                             tGGCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACg                           <  1:3642131/76‑1 (MQ=255)
                                              ggCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACgg                          <  1:3765540/76‑1 (MQ=255)
                                              ggCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACgg                          <  1:2421557/76‑1 (MQ=255)
                                               gCTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGt                         <  2:514060/76‑1 (MQ=255)
                                                ctctTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGtt                        <  1:299348/76‑1 (MQ=255)
                                                ctctTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGtt                        <  1:2585175/76‑1 (MQ=255)
                                                cTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGGCAGCGGGGGCCAGTCGGCGCGCGCgggt                        >  1:1070425/1‑74 (MQ=255)
                                                cTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGtt                        <  2:5237909/76‑1 (MQ=255)
                                                cTCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGGCGGTGCGCACgggt                        >  1:1040724/1‑74 (MQ=255)
                                                 tCGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGttt                       <  1:4925744/76‑1 (MQ=255)
                                                 tCGGTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAACCTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGttt                       <  2:5471523/76‑1 (MQ=255)
                                                  cGTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTc                      <  1:3546739/76‑1 (MQ=255)
                                                  cGTTGCGGAAGATGCGCAGGTATGGCTCCCGCTGAAAATGCTGGGCCAGCGGCGGCCCGTCCGGGCGCGCGGGTTc                      >  2:4503144/1‑76 (MQ=255)
                                                   gTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGGCAGCGGCGTCCCGTCGGTGCGCCCGGGTTcc                     >  1:4416009/1‑76 (MQ=255)
                                                   gTTGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTcc                     <  1:3583807/76‑1 (MQ=255)
                                                   gTTGCGGAAGATGCGCAGGTATGGCTCCCCCCGAAACTGCTGGGGCCGCGGCGGCCAGTGCGTGCGCCCGGgtttc                     >  2:3286944/1‑74 (MQ=255)
                                                    ttGTGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCa                    <  1:752928/76‑1 (MQ=255)
                                                    ttGCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCa                    <  1:3323400/76‑1 (MQ=255)
                                                     tGCGGAAGATGCGCAGGTATGGCTCCTGCTTTAACTGCTGGGCGCGCGGGGGCCAGGCGGGGCGGACCGGTTTctg                    >  1:2038686/1‑74 (MQ=255)
                                                      gCGGAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGc                  <  1:5917208/76‑1 (MQ=255)
                                                      gCGGAAGATGCGCAGGTATGGCTCCTGCGGAAACTTCTGGGCCGGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGc                  <  2:2040469/76‑1 (MQ=255)
                                                      gCGGAAGATGCGCAGGTATGGCTCCCGCCGAAAACGCCGGGCCCGCGGCGCCCAGCCGGGGGGCCCGGGTTTCCac                   >  2:1651040/1‑75 (MQ=255)
                                                        ggAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCCGTCGGTGGGCACGGGTTCCCGccc                >  1:724750/1‑76 (MQ=255)
                                                        ggAAGATGCGCAGGTATGGCTCCCGCTGAAACTGCTGGGCCAGCGGCGTCCCGTCGGTGCGCCCGGGTTTCAGccc                >  2:3613718/1‑76 (MQ=255)
                                                         gAAGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGGTTCCAGCCCt               >  1:1335608/1‑76 (MQ=255)
                                                          aaGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCtt              >  1:6246227/1‑76 (MQ=255)
                                                          aaGATGCGCAGGTATGGCTCCCGCTGAAACTGCTGGGGCCGCGGCGTCCAGGCGGTGGGCGCGGTTTCCAGCCCtt              >  2:246842/1‑76 (MQ=255)
                                                           aGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCttt             <  2:3343329/76‑1 (MQ=255)
                                                           aGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCttt             <  2:342761/76‑1 (MQ=255)
                                                           aGATGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCttt             <  1:3706997/76‑1 (MQ=255)
                                                             aTGCGCAGTTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCTTTaa           <  1:4506723/76‑1 (MQ=255)
                                                             aTGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGGTTTCAGCCCcttta           >  1:1063666/1‑74 (MQ=255)
                                                              tGCGCAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCTTTaaa          <  2:1222365/76‑1 (MQ=255)
                                                              tGCGCAGGTATGGCTCCCGCTGAAACTGCTGGGCCAGCGGCGGCCCAGCGGTGCTCACGGTTTCCACGCCTTTaaa          >  2:6156436/1‑76 (MQ=255)
                                                              tGCGCAGGTATGCCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCTTTaaa          <  1:2213736/76‑1 (MQ=255)
                                                                cgcAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCTTTAAaca        <  1:2394759/76‑1 (MQ=255)
                                                                cgcAGGTATGGCTCCCGCTGAAACTGCTGGGCCCGCGGCGTCCCGTCGGGGGGCCCGGGTTTCAGCCCCTTTAaca        >  2:4534182/1‑76 (MQ=255)
                                                                 gcAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCTTTAAacac       <  1:4089982/76‑1 (MQ=255)
                                                                 gcAGGTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCTTTAAacac       <  2:1040724/76‑1 (MQ=255)
                                                                   gggTATGGCTCTTGCGGAAGCTGATGGGCCTGTTGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCTTTAAACACCa     <  2:3173122/75‑1 (MQ=255)
                                                                   aGGTATAGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCTTTAAACACCa     <  1:841101/76‑1 (MQ=255)
                                                                   aGGTAGGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCTTTAAACACCa     <  1:1458514/76‑1 (MQ=255)
                                                                    ggTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCTTTAAACACCAt    <  1:3902466/76‑1 (MQ=255)
                                                                    ggTATGGCTCCCGCTGAAACTGCTGGGGCAGCGGCGTCCCGTCGGTGCGCACGGGTTCCCGCCCTTTTAACACCAt    >  2:4047592/1‑76 (MQ=255)
                                                                     gTATGGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCTTTAAACACCATc   <  1:3152394/76‑1 (MQ=255)
                                                                      tatgGCTCCCGCTGAAAATGCTGGGCCAGCGGCGTCCCGTCGGTGGGGACGGGTTCCAGCCCCTTAAACACCAccc  >  2:6121938/1‑76 (MQ=255)
                                                                         |                                                                          
ACCCGCCATCAGTTTGTCGATGGTTTCGCGTACATATTCCTGATATGGCTCGTTGCGGAAGATGCGCAGGTATAGCTCCTGCTGAAACTGCTGGGCCAGCGGCGTCCAGTCGGTGCGCACGGTTTCCAGCCCTTTAAACACCATCC  >  REL606/66557‑66702

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