Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,180,2490TG42.9% 223.9 / 48.8 133G24G (GGT→GGGyjaHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/55);  new base G (57/0);  total (78/55)
Fisher's exact test for biased strand distribution p-value = 3.16e-20
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 262 total aligned reads displayed.

GGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAG  >  REL606/4180174‑4180324
                                                                           |                                                                           
ggCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGt                                                                             <  1:5286569/76‑1 (MQ=255)
   aGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTgcg                                                                          <  2:7128926/76‑1 (MQ=255)
     tGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGGGCt                                                                        >  1:1549878/1‑76 (MQ=255)
       aaCAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGt                                                                      >  1:4419348/1‑76 (MQ=255)
        aCAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGtt                                                                     <  2:5571063/76‑1 (MQ=255)
         cAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGttt                                                                    <  2:204773/76‑1 (MQ=255)
         cAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGttt                                                                    <  1:5846416/76‑1 (MQ=255)
         cAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGgtt                                                                    >  2:2436287/1‑76 (MQ=255)
            ttttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTtctc                                                                 <  1:4519159/76‑1 (MQ=255)
            ttttAACGAAGGGGCGGGTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGGGGTTtttc                                                                 >  2:2825649/1‑76 (MQ=255)
             tttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCa                                                                <  2:5670372/76‑1 (MQ=255)
             tttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGGTtctta                                                                >  2:4170936/1‑74 (MQ=255)
             tttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGGTtctcc                                                                >  1:5610701/1‑75 (MQ=255)
              ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCAc                                                               <  2:3694328/76‑1 (MQ=255)
               tAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACt                                                              <  2:5630691/76‑1 (MQ=255)
               tAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACt                                                              <  1:1422811/76‑1 (MQ=255)
                aaCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTg                                                             <  2:4222131/76‑1 (MQ=255)
                aaCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTg                                                             <  2:340048/76‑1 (MQ=255)
                aaCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCCGTTTCTCACTg                                                             >  1:5765661/1‑76 (MQ=255)
                 aCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGc                                                            <  1:4668685/76‑1 (MQ=255)
                 aCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGc                                                            >  2:7829498/1‑76 (MQ=255)
                  cGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTctctctc                                                           >  2:7013745/1‑73 (MQ=255)
                  cGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCTCTTcc                                                           >  2:5146349/1‑76 (MQ=255)
                   gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCt                                                          <  1:1155115/76‑1 (MQ=255)
                   gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTCACTGCCt                                                          >  2:6627085/1‑76 (MQ=255)
                    aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTg                                                         <  2:5026000/76‑1 (MQ=255)
                    aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTg                                                         >  2:3713729/1‑76 (MQ=255)
                    aaGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTTCCTGCCtc                                                         >  1:2446682/1‑75 (MQ=255)
                     aGGGGTGGTTTCACCCCTTTTGTCTTTCTGGTGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTTACTGCCTGt                                                        >  1:5171166/1‑76 (MQ=255)
                     aGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGt                                                        <  1:2930775/76‑1 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTa                                                       <  2:2139497/76‑1 (MQ=255)
                      ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTTACTGCCCGTa                                                       >  2:3799433/1‑76 (MQ=255)
                       gggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAg                                                      <  2:4130091/76‑1 (MQ=255)
                       gggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAg                                                      <  1:1730321/76‑1 (MQ=255)
                       gggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAg                                                      <  1:5879075/76‑1 (MQ=255)
                       gggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACCGCCTGTAg                                                      >  2:1752078/1‑76 (MQ=255)
                        ggtggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGc                                                     >  1:1931625/1‑76 (MQ=255)
                         gtggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGGTGTTTTTCACTGCCTGTAGcc                                                    >  1:3940783/1‑76 (MQ=255)
                         gtggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGGTGTTTCTTACTGCCCTTAGcc                                                    >  1:5568552/1‑76 (MQ=255)
                         gtggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTAggc                                                    >  1:4830456/1‑74 (MQ=255)
                         gtggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACCGCCTTTAggc                                                    >  2:489486/1‑74 (MQ=255)
                           ggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCCCTGCCTGTAGCcac                                                  >  1:7330220/1‑76 (MQ=255)
                           ggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGcaca                                                  >  2:5908210/1‑73 (MQ=255)
                             tttCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCTCTGCCTGTAGCCCCaa                                                >  1:2156198/1‑76 (MQ=255)
                             tttCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACaa                                                >  1:3066686/1‑76 (MQ=255)
                             tttCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCCGTTTCTTACTGACTGTAGCCACaa                                                >  1:6100865/1‑76 (MQ=255)
                              ttCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTTTCACTTCCTGTAGCCACAAc                                               >  1:6524915/1‑76 (MQ=255)
                               tccccctTTTTGTTTTTCTGGCGTCGATCATTGATGGTGGCGGGTGCGTTGTTTCTCCCTGCTTTTAGCCACAAct                                              <  2:271639/73‑1 (MQ=255)
                               tccccctTTTTGTTTTGCGGGCGTCAATCTTTGATGCTGGCAGGTGCTTTTTTTCTCCCTGCTTTTAGCCACAAct                                              <  2:4540045/73‑1 (MQ=255)
                               ccACCCTTTTTGTTTTTCTGGCGTCTTTCATTGATGCTGGCAGGGCCGCTTTTTCTCCCTCCTTTTAGCCACAAct                                              <  2:6470471/75‑1 (MQ=255)
                                cACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGCCCCCActc                                             >  2:6765298/1‑76 (MQ=255)
                                 aCCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCCGTAGCCACCActct                                            >  2:1611763/1‑76 (MQ=255)
                                 aCCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCCGTTTCTCCCTCCCCGTAGCCACAActcc                                            >  2:4634482/1‑75 (MQ=255)
                                  ccccTTTTGTTTTTGTGGCGTCGATCATTGATGCTGGCAGGTGCGTTTTTTCCCCCTGCCTTTAGCCACAACTCtt                                           <  2:5771553/76‑1 (MQ=255)
                                   cccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTc                                          <  1:1015625/76‑1 (MQ=255)
                                   cccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGCTTCTCACTGGCAACTGCCACAACACttg                                          >  2:3194971/1‑75 (MQ=255)
                                    ccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTTTCTCTGCCTTTAGCCACCACTCTTCa                                         >  1:4723899/1‑76 (MQ=255)
                                    ccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGCCACAACTCTTCa                                         >  2:6258454/1‑76 (MQ=255)
                                    ccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACCACTCTTCa                                         >  2:3027826/1‑76 (MQ=255)
                                     cTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCGGTTTTTCCCTGCCTGTTGCCACAACTCTTcca                                        >  1:3267964/1‑74 (MQ=255)
                                       ttttTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCCCTGCCTTTAGCCACAACTCTTCACtt                                      <  2:4261640/76‑1 (MQ=255)
                                       tttGTCTTTCTGGCGTCGATCATTGATGCTTGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACtt                                      >  1:136785/1‑76 (MQ=255)
                                        ttGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTTTCCCTGCCTGTAGGCACAACTCTTCCCTtc                                     >  1:4179894/1‑76 (MQ=255)
                                         tGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGCCTGTTTTTCCCTGCCTGGAAGCACAACTCTTCACTtcc                                    >  1:3589790/1‑76 (MQ=255)
                                         tGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGCCACAACTCTTCACTtcc                                    >  1:7067446/1‑76 (MQ=255)
                                          gTCTTTCAGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTTTCACTGCCTGTAGACACAACTCTTCACTtcct                                   >  2:858677/1‑76 (MQ=255)
                                           tCTTTCTGGCGTCTATCATTGATGCTGGCAGGTGCGCTGTTTCTCCCTGCCTGTAGCCACAACTCTTCACTtcctc                                  <  1:6864874/76‑1 (MQ=255)
                                           tCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGGTGTTTATCACCCCCTGTAGCCACCACCCCTCACTtcccc                                  >  2:7769144/1‑74 (MQ=255)
                                           tCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAAATCTTCACTttctc                                  >  2:4764490/1‑76 (MQ=255)
                                             tttCTGGCGTCGATCATTGATGCTGGCAGGGGCGCGGTTTGTCTCTGCCCGTAGCCCCTACTTTTTCTTTCCTccc                                >  2:1471827/1‑76 (MQ=255)
                                              ttCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCCCTTCCTCCCt                               >  1:4092003/1‑76 (MQ=255)
                                              ttCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCt                               >  2:2287073/1‑76 (MQ=255)
                                              ttCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGGTTCACCCTGCCTGTAGTCACCACGCATCCCCTCCTCCCt                               >  2:6923688/1‑76 (MQ=255)
                                               tcgggCGTCGATTATTGTTGCTGGCAGGTTCGTGGTTTCCCCCGGCCTTTAGCCACAACTCTTCACTTCCTCCCtt                              <  2:2347210/73‑1 (MQ=255)
                                               tcgggCGTCATTCATTCATGTTGGCAGGTGCCCTGTTCCTCCCTGCCTGTAGCCACAACTCTTCACTTCCTCCCtt                              <  2:1297166/73‑1 (MQ=255)
                                               tCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCCCTGCCTTTAGCCACAACTCTTCACTTCCTCCCtt                              <  1:2662266/76‑1 (MQ=255)
                                               tCTGGCGTCGATCATTGATGCTGGCAGGTGCGCCGTTTCTCACTGCCTGTAGCCCCAAATCTTCACTTTCTTCCtt                              >  2:4060242/1‑76 (MQ=255)
                                                 tGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCCCTGCCTGGAGCCACCAATCCTCACTTCCTCCCTTTa                            >  1:7547441/1‑76 (MQ=255)
                                                 tGGCGTCGATCATTGATGCGGGCTGGTGCGTTTTTTCGCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTa                            <  2:7056287/76‑1 (MQ=255)
                                                  ggCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCttttc                           >  1:4354814/1‑74 (MQ=255)
                                                  ggCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAAATCTTCAATTCCTCCCTTTAc                           >  1:3452071/1‑76 (MQ=255)
                                                  ggCGTCGATCATTGAGGCTGGCAGGTGCGCTTTTTCTCCCTGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTAc                           <  1:6261795/76‑1 (MQ=255)
                                                   gCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTAcc                          >  1:6657431/1‑76 (MQ=255)
                                                      tCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGcc                       <  1:3852001/76‑1 (MQ=255)
                                                       cGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCCCTTCCTCCCTTTTCCGCCa                      >  2:1046745/1‑76 (MQ=255)
                                                         aTCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGt                    >  2:5853695/1‑76 (MQ=255)
                                                           cATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCCCCCGTgg                  >  1:1596408/1‑76 (MQ=255)
                                                           cATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACCCTTCACTTCCTCCCTTTACCGCCAGCgg                  >  1:895416/1‑76 (MQ=255)
                                                             ttGATGCTGGCAGGGGCGCTGTTTCTCACTGGCTGTAGACACAACCCTTCCCTTTCTCCCTTTCATGCCAGTGGAt                >  2:1977863/1‑76 (MQ=255)
                                                                    tGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGa         >  1:6339919/1‑76 (MQ=255)
                                                                       ggggTGCTCTTTTTCTCTCTGCTTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGa      <  2:7253732/74‑1 (MQ=255)
                                                                       cAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGa      >  2:6834776/1‑76 (MQ=255)
                                                                       cAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGa      >  1:2395466/1‑76 (MQ=255)
                                                                        aGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAc     >  1:5469117/1‑76 (MQ=255)
                                                                        aGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAc     >  1:7798520/1‑76 (MQ=255)
                                                                        aGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAc     >  1:1518223/1‑76 (MQ=255)
                                                                        aGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTCGCTTTGCTGAAGAc     >  2:7389263/1‑76 (MQ=255)
                                                                          gTGCGTTTTTTCCCCCTCCTTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCa   <  2:2210572/76‑1 (MQ=255)
                                                                          gTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCa   >  2:4439836/1‑76 (MQ=255)
                                                                           tGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAg  >  1:6919864/1‑76 (MQ=255)
                                                                           |                                                                           
GGCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAG  >  REL606/4180174‑4180324

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