Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,180,2490TG46.1% 195.4 / 65.2 154G24G (GGT→GGGyjaHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/63);  new base G (71/0);  total (91/63)
Fisher's exact test for biased strand distribution p-value = 7.30e-26
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 303 total aligned reads displayed.

GGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAG  >  REL606/4180174‑4180324
                                                                           |                                                                           
ggCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGt                                                                             <  2:1206246/76‑1 (MQ=255)
 gCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTg                                                                            <  1:6615980/76‑1 (MQ=255)
 gCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTg                                                                            <  2:4712796/76‑1 (MQ=255)
  cAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTgc                                                                           <  2:10166038/76‑1 (MQ=255)
   aGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTgcg                                                                          <  2:6832514/76‑1 (MQ=255)
   aGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGgcg                                                                          >  2:9680204/1‑76 (MQ=255)
       aaCAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGt                                                                      <  2:5576898/76‑1 (MQ=255)
          aGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTtc                                                                   <  1:4704082/76‑1 (MQ=255)
            ttttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTtctc                                                                 <  2:1492093/76‑1 (MQ=255)
             tttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGGTGTTtctta                                                                >  1:8294462/1‑74 (MQ=255)
             tttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTtttcc                                                                >  1:2549489/1‑75 (MQ=255)
              ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCAc                                                               <  1:5733551/76‑1 (MQ=255)
              ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCCGTTTTTCAc                                                               >  1:4382970/1‑76 (MQ=255)
               tAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACt                                                              <  1:4051242/76‑1 (MQ=255)
               tAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACt                                                              <  1:4669857/76‑1 (MQ=255)
               tAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACt                                                              <  1:4637491/76‑1 (MQ=255)
                aaCGAAGGGGTTGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTg                                                             <  2:7750650/76‑1 (MQ=255)
                 aCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGc                                                            <  2:420169/76‑1 (MQ=255)
                 aCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTCTCTGc                                                            >  2:6995309/1‑76 (MQ=255)
                  cGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGcc                                                           >  2:9818106/1‑76 (MQ=255)
                   gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTCCCTGCCt                                                          >  2:10210914/1‑76 (MQ=255)
                   gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCGGGTTTTCACTGCCt                                                          >  1:1314723/1‑76 (MQ=255)
                   gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCCGGTTTTCACTTccc                                                          >  1:6246645/1‑75 (MQ=255)
                   gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCCGTTTCTCTCTGCCt                                                          >  2:9071260/1‑76 (MQ=255)
                   gAAGGCGTGCTTTCTCCCCTTTTGTGTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCt                                                          <  1:198076/76‑1 (MQ=255)
                    aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCCCTTccct                                                         >  2:6166094/1‑74 (MQ=255)
                    aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGGTTTTTcctgcctg                                                         >  2:8311281/1‑76 (MQ=255)
                    aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGCGCGCTGTTTCTCTCTGCCTg                                                         >  2:7396433/1‑76 (MQ=255)
                    aaGGGGTGGTTTAACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCtt                                                         >  2:4321525/1‑75 (MQ=255)
                     aGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGt                                                        <  1:9439464/76‑1 (MQ=255)
                     aGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCCGTTTCTccctccctct                                                        >  2:3555232/1‑74 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGGCCTGGTTCTCGCTTCCTGTa                                                       >  2:9587107/1‑76 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTTTCCCTGCCTGTa                                                       >  2:2283915/1‑76 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTa                                                       <  1:5297782/76‑1 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTa                                                       <  1:2170891/76‑1 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGGTTCTTCCTGCCCGtt                                                       >  1:1016971/1‑75 (MQ=255)
                       gggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAg                                                      <  1:9870689/76‑1 (MQ=255)
                       gggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGGTTTTCACTGCCTTTAg                                                      >  2:9388443/1‑76 (MQ=255)
                        ggtggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGc                                                     <  2:10043932/76‑1 (MQ=255)
                          tggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCca                                                   <  1:8957955/76‑1 (MQ=255)
                           ggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCcac                                                  >  1:9318640/1‑76 (MQ=255)
                             tttCCCCCTTTTTGTCTTTCGGGCGTCGATCATTGATGCTGGCGGGTGCTCTTTTTCTCCCTGCCTTTAGCCACaa                                                <  2:2883925/76‑1 (MQ=255)
                               tCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGAGGCTTTTCACTGCCTGTAAACACAAca                                              >  2:8057813/1‑75 (MQ=255)
                               tCACCCCTTTTGTCTTTCTGACGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACCGCCTGCTGCCACaaat                                              >  2:6606395/1‑74 (MQ=255)
                                cgccccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCGGTTTTTCACTCGCCGTAGGCACAActc                                             >  1:4248890/3‑76 (MQ=255)
                                cACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTTTCACTGCCTGTAGGCACCActc                                             >  2:6040678/1‑76 (MQ=255)
                                  ccccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTTTCACTGCCTGTAGCCACAACTCtt                                           >  1:161630/1‑76 (MQ=255)
                                   cccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTCCCTGCCCGTCGCCACAACTCTTc                                          >  1:4207467/1‑76 (MQ=255)
                                    ccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTTTCACTGCCTGTAGCCACAAATCTTCa                                         >  2:1650843/1‑76 (MQ=255)
                                    ccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGGTGTTTTTCCCTGCCTGGAGCCACAACTTTTcc                                         >  1:6835451/1‑75 (MQ=255)
                                     cTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCCGTTGCCACCAATCtttcc                                        >  1:5851439/1‑73 (MQ=255)
                                      ttttGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTTTCCCTGCCCGCCGCCACCACCCTTTACt                                       >  1:365133/1‑76 (MQ=255)
                                      ttttGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCGGTTTTTCTCGTCCTATAAGCACAACTCTTCCCt                                       >  1:3942424/1‑76 (MQ=255)
                                       tttGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCCGTTTCTCACTGCCTGTAGCCACCACTCTTCACtt                                      >  1:4614340/1‑76 (MQ=255)
                                        ttGTCTTTCTGGCGTCGATCATTGATGCTGGCGGGTGCGCTGTTTCTCCCTGCCTGTAGCCACAACTCTTCACTtc                                     <  2:58868/76‑1 (MQ=255)
                                         tGTCTTTCTGGGGTCGATCATTGATGCTGGCAGGGGCGCTTTTTTTCCCCGCCCTTAGCCACAACTCTTCACCtcc                                    >  2:6242203/1‑76 (MQ=255)
                                         tGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCCCTGCCTTTAGCCCCAACACTTCACTtcc                                    >  1:8431477/1‑76 (MQ=255)
                                         tGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTtcc                                    >  1:5969575/1‑76 (MQ=255)
                                         tGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCCGGTTTTCTCTGCCTGTAGCCACAAATCTCCACtttc                                    >  1:970149/1‑74 (MQ=255)
                                          gTCTTTCGGGCGGCGATCATTGATGCGGGCGGGTGCGCTGTTTCCCCCTGCCTTTAGCCACAACTCTTCACTtcct                                   <  2:1928224/76‑1 (MQ=255)
                                           tcgttCTGGCTGCGGTCTTAGATGCTGCCAGGTGGGCTGTTTCTCACCGCTTTTAGCCACAACTCTTCACTtcctc                                  <  2:2468725/73‑1 (MQ=255)
                                           tCTTTCTGGCGTCTATCATTGATGCTGGCGGGTGCGCTTTTTCTCCCTGCCTTTAGCCACAACTCTTCACTtcctc                                  <  2:5256549/76‑1 (MQ=255)
                                           tCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCCCTGCCTGTAGCCACAACTCTTCACTtcctc                                  >  1:9622325/1‑76 (MQ=255)
                                           tCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGCCACAACTCTTTACTtcctc                                  >  1:7245200/1‑76 (MQ=255)
                                            cTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTttctcc                                 >  2:8676577/1‑76 (MQ=255)
                                             tttCTGGCGTCGATCTTTGATGGTGGGGGGTGCCCTTTTTCTCCCTGCTTTTAGCCACAACTCTTCACTTCCTccc                                <  2:1387454/76‑1 (MQ=255)
                                             tttCTGGCGTCGATCATTGATGCTGGCAGGTGCGCGGTTTCTCACTGCCTGTTGCCACCACTCTTCACTTcccccc                                >  1:4849654/1‑76 (MQ=255)
                                             tttCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTTTCCCCGCCTGTAGCCCCAACTCTTCACTTCCTccc                                >  1:6440044/1‑76 (MQ=255)
                                             tttCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGCCACAACTCTTCACTTCCTccc                                >  2:3249321/1‑76 (MQ=255)
                                             tttCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCCGTAGCCCAAAATCTTCACTTCCTccc                                >  1:7343960/1‑76 (MQ=255)
                                              ttCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCt                               >  1:2730903/1‑76 (MQ=255)
                                              ttCTGGCGTCGATCATTGATGCTGGCAGGGGCGCCGTTTCTCACACCCCGTAGCCACAAATCTTCACTTTCTcccc                               >  1:8719624/1‑75 (MQ=255)
                                               tCTGGCGTCGATCATTGATGCTGGCAGGTGCCCTGTTTCTCCCTGCCTGTAGCCACAACTCTTCACTTCCTCCCtt                              <  1:239749/76‑1 (MQ=255)
                                                cTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTTTCACTGCCTGTAGCCACAACTCTTCACCTCCTCCCctt                             >  2:4509520/1‑76 (MQ=255)
                                                 tGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTa                            >  2:3774960/1‑76 (MQ=255)
                                                  ggggTCGATCATGGATGCTGGCGGGTGCGCTTTTTCTCCCTGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTAc                           <  2:4072658/73‑1 (MQ=255)
                                                  ggCGTCGGTCTTTGGTGTGGGCAGGTGCGTTTTTTCTCACTGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTAc                           <  1:3451322/76‑1 (MQ=255)
                                                  ggCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCCCTTCCTTTAGCCACAAATCTTCACTTCCTCTCttttt                           >  2:1547618/1‑74 (MQ=255)
                                                  ggCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCAATGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTAc                           >  2:5106525/1‑76 (MQ=255)
                                                   gCGTCGATCTTTGATGCTGGCGGGGGCGCTGTTTCTCACTCCCTGTAGCCACAACTCTTCACTTCCTCCCTTTAcc                          <  2:2548552/76‑1 (MQ=255)
                                                      tCGATCTTGGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGca                       <  1:3106260/76‑2 (MQ=255)
                                                      tCGATCATTGATGCTGGCGGGTGCGCTTTTTCTCACTGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGcc                       <  1:6071180/76‑1 (MQ=255)
                                                       cGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCa                      >  1:7726003/1‑76 (MQ=255)
                                                        gATTATTGATGGTGGCAGTGGCCTTTTTTCTCCCTGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAg                     <  1:2908736/76‑1 (MQ=255)
                                                        gATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTATCCCCAg                     >  2:6410092/1‑76 (MQ=255)
                                                          tCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGCCACAACTCTTCACTTTCTCCCTTTTCCGCCAGTg                   >  2:5359274/1‑76 (MQ=255)
                                                          tCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCCGTg                   >  1:1045659/1‑76 (MQ=255)
                                                             ttgttgCTGGCGGGTGCCTTTTTCCCCCCTGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGAt                <  1:7056598/76‑1 (MQ=255)
                                                              tgttGCGGGCAGGGGCGCGGTTTCTCCCTCCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGAtt               <  1:5146649/73‑1 (MQ=255)
                                                              tGATGCTGGGAGGTGCGTTGTTTCTCCCTGCTTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGAtt               <  1:765928/76‑1 (MQ=255)
                                                              tGATGCTGGCGGGGGCCCTTTTTCCCCCTGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGAtt               <  1:5624165/76‑1 (MQ=255)
                                                              tGATGCGGGCGGTTGCTCTTTTTCTCCCTGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGAtt               <  1:6865792/76‑1 (MQ=255)
                                                               gATGCTGGCAGGTGCGCTGTTTCTCCCCGCCTGGAGCCACAACCCTTCTCTTCCTCCCTTTACCGCCAGTGGAttt              >  2:890418/1‑76 (MQ=255)
                                                               gATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCCGTGGAttt              >  2:10325561/1‑76 (MQ=255)
                                                                   gtGGCGGGTGCGTTTTTCCCCCCTCCTTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTg          <  2:5667688/75‑1 (MQ=255)
                                                                    tGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGa         >  1:7575341/1‑76 (MQ=255)
                                                                         ggTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAcc    >  2:8593541/1‑76 (MQ=255)
                                                                         ggTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAcc    >  2:3026560/1‑76 (MQ=255)
                                                                          gTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCa   >  1:1870369/1‑76 (MQ=255)
                                                                           tGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAg  >  1:8668613/1‑76 (MQ=255)
                                                                           |                                                                           
GGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAG  >  REL606/4180174‑4180324

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: