Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,485,4500AG34.3% 238.7 / 91.7 137G193G (GGA→GGGydcVpredicted spermidine/putrescine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (0/90);  new base G (47/0);  total (47/90)
Fisher's exact test for biased strand distribution p-value = 7.67e-38
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 160 total aligned reads displayed.

ATCTTGGGGCCAATGGCTGGCAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCATGAACGAACGTTA  >  REL606/1485377‑1485525
                                                                         |                                                                            
aTCTTGGGGCCAATGGCTGGCAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCaggggg                                                                            >  2:4527707/1‑76 (MQ=255)
aTCTTGGGGCCAATGGCTGGCAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGCTCGGCGTTACTGGCaggggg                                                                            >  2:4661221/1‑76 (MQ=255)
 tCTTGGGGCCAATGGCTGGCAACCCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAggagga                                                                           <  1:439198/76‑1 (MQ=255)
  ctgggggCAAGTGGCGGGAAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGaa                                                                          <  1:719810/73‑1 (MQ=255)
  cTTGGGGCCAATGGCTGGCAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGaa                                                                          <  2:3813704/76‑1 (MQ=255)
   ttGGGGCCAATGGCTGACAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCaggggggaa                                                                         >  1:1334921/1‑73 (MQ=255)
   tgggggCCAATGGCTGGCAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAAt                                                                         <  1:3771725/74‑1 (MQ=255)
    tgtggCCAAGGCCTGGCAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATg                                                                       <  2:2543808/73‑1 (MQ=255)
    tGGGGCCAATGGCTGGCAACCCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATg                                                                       <  2:4452886/76‑1 (MQ=255)
    tGGGGCCAATGGCTGGCAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATg                                                                       <  2:1584907/76‑1 (MQ=255)
    tGGGCCAAATGGCTGGCAACCCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATg                                                                       <  2:2663224/76‑1 (MQ=255)
        gCCAATGCCTGGCAATCCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTgg                                                                   <  1:4878184/76‑1 (MQ=255)
          cAATGGCTGGCAACCCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCg                                                                 <  2:297262/76‑1 (MQ=255)
          cAATGGCTGGCAACCCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCg                                                                 <  1:4704672/76‑1 (MQ=255)
           aaTGCCTGGCAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGt                                                                <  1:1643184/76‑1 (MQ=255)
            aTGGCTGGCAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGtt                                                               <  2:593472/76‑1 (MQ=255)
             tggctggcAAATCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGAAGCAGGCGctt                                                              >  1:1389560/1‑76 (MQ=255)
             tgcctggcAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGttt                                                              <  2:3165504/73‑1 (MQ=255)
              tgctggaAAGCCTTCCGCTACGTATTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTg                                                             <  2:264092/75‑1 (MQ=255)
              ggctggcAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGGTGGGGGTTg                                                             >  2:1706330/1‑76 (MQ=255)
               gctggcAACCCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGc                                                            <  1:1428311/76‑1 (MQ=255)
                  ggcAACCCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCtt                                                         <  1:2702362/76‑1 (MQ=255)
                  ggcAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCtt                                                         <  2:2313745/76‑1 (MQ=255)
                  ggcAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCtt                                                         <  2:724795/76‑1 (MQ=255)
                    cAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGAGGCCGGGGTGTGCCtttt                                                       >  1:875143/1‑74 (MQ=255)
                     aacccTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTc                                                      <  1:25552/73‑1 (MQ=255)
                      aaCCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGAAGCTGGCGTTTGGCTTGttg                                                     >  2:3259666/1‑74 (MQ=255)
                        ccTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGGCGGGGGTTTGCTTGTCGtt                                                   >  2:1099187/1‑76 (MQ=255)
                        ccTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGGCGGGGGTTTCCTTGTCGtt                                                   >  1:623994/1‑76 (MQ=255)
                        ccTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGAAGCTGGCGGTTGCCTTTTCGtt                                                   >  2:4323431/1‑76 (MQ=255)
                         cTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGGTGGCGTTTGCCTCTTCGttg                                                  >  1:2208255/1‑75 (MQ=255)
                         cTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTc                                                  <  1:2079610/76‑1 (MQ=255)
                           tCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGAATATGGCGGTTGCCCTGTCGCTCgt                                                >  1:1087335/1‑75 (MQ=255)
                           tCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGa                                                <  2:2374286/76‑1 (MQ=255)
                             cGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGGTGCTGGGGGTTGCCTTGTCGTTCGATg                                              >  1:4328335/1‑76 (MQ=255)
                             cGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATg                                              <  2:1607817/76‑1 (MQ=255)
                             cGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATg                                              <  2:1676669/76‑1 (MQ=255)
                             cGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATg                                              <  1:86680/76‑1 (MQ=255)
                             cGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATg                                              <  2:2926306/76‑1 (MQ=255)
                              gCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGa                                             <  1:1016970/76‑1 (MQ=255)
                              gCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGa                                             <  2:4487418/76‑1 (MQ=255)
                               cTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGAGGGGGGGGGTTTCCTTGTCGTTTTATGaa                                            >  1:2368749/1‑76 (MQ=255)
                               cTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGaa                                            <  1:2907689/76‑1 (MQ=255)
                                tACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGGTGGCGGTTGCCTTGTCGTTTGATGaaa                                           >  1:1040386/1‑76 (MQ=255)
                                 aCGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGGTGGCGTTTTCCCTGTCGTTTGGTGAAat                                          >  2:1890903/1‑76 (MQ=255)
                                 aCGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGCTGGCGGTTTTCCTGTCGTTTGGTGAAat                                          >  1:487075/1‑76 (MQ=255)
                                 aCGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGCTGGCGGTTGCCCTGTCGTTTGGTGAAat                                          >  1:1571325/1‑76 (MQ=255)
                                 aCGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGAGGCTGGCGGTTGCCTTGTTCTTTGATGgaaa                                          >  2:2958248/1‑75 (MQ=255)
                                 aCGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAat                                          <  1:2549822/76‑1 (MQ=255)
                                  cGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGCTTGCGTTTGCCTTTTCGTTGTCTGAAatc                                         >  2:1452816/1‑76 (MQ=255)
                                  cGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGAAGCTGGCGTTTTGCTTGTCGTTTGATGcaaaa                                         >  1:2332437/1‑74 (MQ=255)
                                  cGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAatc                                         <  2:5074112/76‑1 (MQ=255)
                                  cGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAatc                                         <  2:4293273/76‑1 (MQ=255)
                                  cGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAatc                                         <  1:4560617/76‑1 (MQ=255)
                                   ttAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAatca                                        <  1:633308/75‑1 (MQ=255)
                                   gTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGCTGGGGGTTGCCTTGTCGTTCGCTGAAataa                                        >  2:854467/1‑74 (MQ=255)
                                   gTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAatca                                        <  1:2742423/76‑1 (MQ=255)
                                    tAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGCTGGCGTTTTCCTTGTCGTTCGATGAAatcat                                       >  1:3666698/1‑76 (MQ=255)
                                     aGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGCTGGCGTTTGCCTTTTCGTTCGATGAAaacatc                                      >  1:957294/1‑76 (MQ=255)
                                     aGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAatcatc                                      <  1:404305/76‑1 (MQ=255)
                                       tgtTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGCTGGGGTTTGCCCTGTCGTTCGATGAAAATATCGt                                    >  1:4689041/1‑76 (MQ=255)
                                        gtTGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGGTGGCGTTTTCCCTTTCGTTCGGTGAAaacatctct                                   >  1:1789753/1‑73 (MQ=255)
                                        gtTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGtt                                   <  1:4527707/76‑1 (MQ=255)
                                        gtTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGtt                                   <  1:765704/76‑1 (MQ=255)
                                          tGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGCTGGCGTTTGCCTTGTCGTTCGATGAAAACATCGTTAc                                 >  1:2254684/1‑76 (MQ=255)
                                          tGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGAGGCTGGCGTTTGCCTTGTTGTTCGCTGAAAACATCGTTAc                                 >  2:1873381/1‑76 (MQ=255)
                                          tGCCGAATCTCAGTTCGGCGTTACTGGCAGGGGGGAAGGTGGGGGTTGGCTTGTCGTTTGCTGAAATAATCGTTAc                                 >  2:3901157/1‑76 (MQ=255)
                                           gCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTAca                                <  2:5048325/76‑2 (MQ=255)
                                           gCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACg                                <  1:3881169/76‑1 (MQ=255)
                                           gCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACg                                <  1:2954614/76‑1 (MQ=255)
                                            ccGAATCTCAGTTCGGCGTTACTGGCAGGGGGGATGCTGGCGTTTGTCTTGTTGTTCGATGAAATCATCGTTTCGa                               >  1:2830778/1‑76 (MQ=255)
                                             cGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGAc                              <  2:3705965/76‑1 (MQ=255)
                                              gAATCCCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGAcc                             <  1:4939771/76‑1 (MQ=255)
                                               aaTCTCAGTTCGGCGTTACTGGCAGGGGGGATGCTGGCGGTTGGCTTGTCGTTCGCTGAAATCATCGTTTCGACCt                            >  2:4189249/1‑76 (MQ=255)
                                               aaTCTCAGTTCGGCGTTACTGGCAGGGGGGATGCTGGCGGTTGCCTTGTCGTTCGATGAAAACATCATTACGACCt                            >  2:240957/1‑76 (MQ=255)
                                               aaTCTCAGTTCGGCGTTACTGGCAGGGGGGAAGGTGGGGGTTTTCTTGTCGTTCGGTGAAATAATCGTTACGACCt                            >  2:3403530/1‑76 (MQ=255)
                                               aaTCTCAGTTCGGCGTTACTGGCAGGGGGAAATCTGGGGGGTTCCTTCGCGGTCGATGAAATCAATGGTACGATCt                            >  2:4552720/1‑76 (MQ=255)
                                               aaTCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCt                            <  1:4853311/76‑1 (MQ=255)
                                                 tctcAGTTCGGCGTTACTGGCAGGGGGGATGGTGGGGGTTTGCTTGTTGTTCGATGAAAACCTCGTTACGACCttt                          >  2:4598355/1‑76 (MQ=255)
                                                  ctcAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTa                         <  2:1185575/76‑1 (MQ=255)
                                                    cAGTTCGGCGTTACTGGCAGGGGGGATGCTGGGGTTTGCCTTGTCGTTCGCTGAAATCATCGTTTCGACCTTTACg                       >  1:1011175/1‑76 (MQ=255)
                                                      gtgcGGCGTTACTGGCCGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGc                     <  1:5157229/73‑1 (MQ=255)
                                                      gTTCGGCGTTACTGGCAGGAGGAATGCTGGCTTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGc                     <  2:122323/76‑1 (MQ=255)
                                                          ggCGTTACTGGCAGGGGGGAAGCTGGCGGTTTCCTTGTAGTTCGGTGAAAACATCGTTACGACCTTTTCGGGAggg                 >  2:631621/1‑75 (MQ=255)
                                                           gCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTc                <  1:2299092/76‑1 (MQ=255)
                                                             gTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCAt              <  1:4552490/76‑1 (MQ=255)
                                                             gTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCAt              <  2:3277300/76‑1 (MQ=255)
                                                                aCTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCATgaa           <  1:654932/76‑1 (MQ=255)
                                                                 cTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCATgaac          <  1:4493307/76‑1 (MQ=255)
                                                                    gCAGGAGGAATGCTGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCATgaacgaac      <  1:4023363/76‑1 (MQ=255)
                                                                     cAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCATgaacgaac      <  2:767159/76‑1 (MQ=255)
                                                                     cAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCATgaacgaac      <  2:2880853/76‑1 (MQ=255)
                                                                     cAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCATgaacgaac      <  1:5058608/76‑1 (MQ=255)
                                                                     cAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCATgaacgaac      <  1:97798/76‑1 (MQ=255)
                                                                     cAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCATgaacgaac      <  1:1960163/76‑1 (MQ=255)
                                                                     cAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCATgaacgaac      <  2:3025253/76‑1 (MQ=255)
                                                                      aggaggAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCATGaacgaacg     <  2:4974410/76‑1 (MQ=255)
                                                                      aggaggAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCATGaacgaacg     <  1:4179254/76‑1 (MQ=255)
                                                                      aggaggAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTGACGACCTGTACGGCAGGGCGTggacgacac     >  2:2774340/1‑73 (MQ=255)
                                                                         aggaATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCATGAACGAACGTTa  <  1:3779331/76‑1 (MQ=255)
                                                                         |                                                                            
ATCTTGGGGCCAATGGCTGGCAAACCTTCCGCTACGTAGTGTTGCCGAATCTCAGTTCGGCGTTACTGGCAGGAGGAATGCTGGCGTTTGCCTTGTCGTTCGATGAAATCATCGTTACGACCTTTACGGCAGGTCATGAACGAACGTTA  >  REL606/1485377‑1485525

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