Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,180,2490TG45.8% 172.9 / 54.0 131G24G (GGT→GGGyjaHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (22/49);  new base G (60/0);  total (82/49)
Fisher's exact test for biased strand distribution p-value = 7.67e-19
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 220 total aligned reads displayed.

GCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAG  >  REL606/4180175‑4180324
                                                                           |                                                                            
gCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTg                                                                             <  2:3823986/76‑1 (MQ=255)
gCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTg                                                                             <  1:6352957/76‑1 (MQ=255)
  aGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTgcg                                                                           <  1:6548920/76‑1 (MQ=255)
  aGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTgcg                                                                           <  1:371439/76‑1 (MQ=255)
    tGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCt                                                                         <  2:2346412/76‑1 (MQ=255)
    tGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCt                                                                         <  1:9274876/76‑1 (MQ=255)
    tGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCt                                                                         <  2:6113252/76‑1 (MQ=255)
    tGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCt                                                                         <  2:8700261/76‑1 (MQ=255)
    tGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCt                                                                         <  1:4828868/76‑1 (MQ=255)
    tGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCt                                                                         <  1:2648951/76‑1 (MQ=255)
    tGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCt                                                                         >  2:9073765/1‑76 (MQ=255)
     gAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTgcgtgt                                                                        >  2:7172136/1‑73 (MQ=255)
     gAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTg                                                                       >  1:7044962/1‑76 (MQ=255)
      aaCAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGt                                                                      <  1:8797562/76‑1 (MQ=255)
      aaCAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGt                                                                      <  1:2080532/76‑1 (MQ=255)
      aaCAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGt                                                                      <  1:3304208/76‑1 (MQ=255)
       aCAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGtt                                                                     <  1:4030009/76‑1 (MQ=255)
          gTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTtct                                                                  <  1:8857705/76‑1 (MQ=255)
            tttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCGGTTTTTCa                                                                >  1:2121315/1‑76 (MQ=255)
             ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCAc                                                               <  1:5334428/76‑1 (MQ=255)
             ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCAc                                                               <  2:1863303/76‑1 (MQ=255)
             ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCAc                                                               <  1:2983595/76‑1 (MQ=255)
             ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCAc                                                               <  1:2287350/76‑1 (MQ=255)
             ttAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTtctctc                                                               >  1:6762379/1‑74 (MQ=255)
               aaCGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGGTGTTTCTCACtt                                                             >  2:6412601/1‑75 (MQ=255)
                 cGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGcc                                                           <  2:3239544/76‑1 (MQ=255)
                 cGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGGTGGTTCTTCCTGcc                                                           >  2:3122110/1‑76 (MQ=255)
                  gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCt                                                          <  1:6696649/76‑1 (MQ=255)
                  gAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGGTGTTTTTCACTGCCt                                                          >  1:3027890/1‑76 (MQ=255)
                     ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTCCCTGTCTgga                                                       >  2:601726/1‑74 (MQ=255)
                     ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTa                                                       >  2:6445722/1‑76 (MQ=255)
                     ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTa                                                       >  2:2968615/1‑76 (MQ=255)
                     ggggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTa                                                       >  2:5188657/1‑76 (MQ=255)
                      gggTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAg                                                      <  2:8118878/76‑1 (MQ=255)
                         tggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCca                                                   <  1:4213059/76‑1 (MQ=255)
                         tggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCca                                                   <  2:4665962/76‑1 (MQ=255)
                         tggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCca                                                   >  2:6508561/1‑76 (MQ=255)
                          ggtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTTGCcac                                                  >  1:8976437/1‑76 (MQ=255)
                           gtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGGTGTTTTTTACTGCCTGTAGGcaca                                                 >  1:690243/1‑76 (MQ=255)
                           gtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTTACTGCCTGTAGCcaca                                                 >  1:8209668/1‑76 (MQ=255)
                           gtTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTCCCTGACTTTTGGcaca                                                 >  1:7582449/1‑76 (MQ=255)
                               cACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGGTGTTTCTTCCTGCCTTTTGCCACaaacc                                             >  1:7224004/1‑73 (MQ=255)
                               cACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGTTTCTTACCTCCTGTAGCCACAAccc                                             >  1:556651/1‑74 (MQ=255)
                               cACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTTGCCACAActc                                             >  1:2209029/1‑76 (MQ=255)
                               cACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGGTTCTTACTGCCCGTAGCCACAActc                                             >  1:2803984/1‑76 (MQ=255)
                                aCCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTTACTGCCTGTAGCCCCAActct                                            >  2:2148843/1‑76 (MQ=255)
                                aCCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCCCAActct                                            >  2:257127/1‑76 (MQ=255)
                                 ccccTTTGTGTCTTTCTGGCGTCGATCATTGATGCTGGCGGTTGGCTTCTTTCCCCTTGCCGTTAGCCACAActct                                            <  2:1390424/76‑1 (MQ=255)
                                 ccccTTTTTGTTTTCTGGCGTCGATATTTGTTGCTGCCGGGGGCGCGGTTCCCCGCTGCCTTTAGCCACAACTCtt                                           <  2:1785745/76‑1 (MQ=255)
                                 ccccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCCCAACTCtt                                           >  2:820068/1‑76 (MQ=255)
                                 ccccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGAGTGCTCTTTCTCACTCGCTGTAGCCGCAAATCtt                                           >  2:4549269/1‑76 (MQ=255)
                                  cctttttGTCTTTCTGGCGTCGATCATTGATGCGGGCGGGTGCGCTTTTTCTCCCTGCCTTTAGCCACAACTCTTc                                          <  2:8976437/73‑1 (MQ=255)
                                  cctttttGTCTTTCTGGCGTCGATCATTGATCCTGGGGGGGGCGCTTTTTCCCCCTGCCTTTAGCCACAACTCTTc                                          <  2:120501/73‑1 (MQ=255)
                                  cctttttGGTTTCCTGCCGTCGATCATTGATGCGGGCAGGTGCGCTTTTTCCCCCTGCCTTTAGCCACAACTCTTc                                          <  2:1083870/73‑1 (MQ=255)
                                  cccTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGGTTCTCACCGCCTGTGGCCACCACTCTTc                                          >  1:1710233/1‑76 (MQ=255)
                                   ctttttGTCTTTTTGGCGTTAATCGTTGGTGCTGGGAGGTGCGTTTTTTCTCTCTGCTTGTAGCCACAACTCTTCa                                         <  2:8014375/74‑1 (MQ=255)
                                    cTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTTGCCCCAACTCTTCAc                                        >  1:8327073/1‑76 (MQ=255)
                                      tttGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTTGCCCCCACTCTTCACtt                                      >  1:1588531/1‑76 (MQ=255)
                                         tGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGGTTTTTCCTTCCTTTTGCCCCCACTCTTCACTtcc                                    >  1:8398116/1‑76 (MQ=255)
                                          gtttttCGGGGGTCGTTCTTTGATGCTGTCAGGTGCGTTGTTTCTCCCTGCCTTTAGCCACAACTCTTCACTtcct                                   <  2:8227003/73‑1 (MQ=255)
                                           tttttGTGGCGTGGATCTTTGATGCTGGCGGGTGCGTTGTTTCCCCCGGCCTTTAGCCACAACTCTTCACTtcctc                                  <  2:7665353/74‑1 (MQ=255)
                                             tttCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCCCTGCCTGTTGCCACAACACTTTACCTcccccc                                >  1:121785/1‑76 (MQ=255)
                                             tttCTGGCGTCGATCATTGATGCTGGCAGGGGCGCTGGTTCTCCCTGCCCGTTGACACCACTCTTTCCTTcccccc                                >  1:4944767/1‑76 (MQ=255)
                                              ttCTGGCGTCGATCATTGATGCTGGCAGGGGGGCTGGTTCTCACTGGCTGTAACCCCCAATCTTCACTTTCTCCCt                               >  1:1078750/1‑76 (MQ=255)
                                              ttCGGGCGTCGATCATGGATGCGGGCGGGTGGGTTGTTTCTCCCTGCCTGTAGCCACAACTCTTCACTTCCTCCCt                               <  2:8832147/76‑1 (MQ=255)
                                                 tGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTTACTGCCTGTAGCCCCAACTCTTCACTTCCTCCCtttt                            >  2:521544/1‑75 (MQ=255)
                                                 tGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTTTCACTGCCTGTAGCCACCACTCCTCCCTTCCTCCCCTTa                            >  2:274838/1‑76 (MQ=255)
                                                 tGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACCACTCTTCACTTTCTCCCTTTa                            >  2:3466948/1‑76 (MQ=255)
                                                 tGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTa                            >  2:1046345/1‑76 (MQ=255)
                                                 tGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCCTTa                            >  2:862163/1‑76 (MQ=255)
                                                 tGGCGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACCCTTCCCCTTCTCCCCTTa                            >  2:480781/1‑76 (MQ=255)
                                                    ggTCGATCATTGATGCTGGCAGGTGCGCTTTTTCTCCCTGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTACCg                         <  1:7540783/75‑1 (MQ=255)
                                                    ggTCGATCATGGTTGCGGAAAGTTGCGTTGTTCCTCCCTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCg                         <  1:1883947/75‑1 (MQ=255)
                                                    cGTCGATCATTGATGCTGGCGGGTGCGCTGTTTCTCCCTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCg                         <  1:5040929/76‑1 (MQ=255)
                                                    cGTCGATCATTGATGCTGGCGGGGGCGTTGTTTCTCCCTGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTACCg                         <  1:3247010/76‑1 (MQ=255)
                                                    cGTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACCGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCg                         >  2:8290284/1‑76 (MQ=255)
                                                     gTCGATCATTGATGCTGGCAGGGGGGGTGGTTTTCACTGACTGGAGGGACAAACATTTACTTTATCCCTTTACTGc                        >  2:4258553/1‑76 (MQ=255)
                                                     gTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACCACTCTTCACCTCCCCCCCTTTCCGc                        >  2:5344323/1‑76 (MQ=255)
                                                     gTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTTCCGc                        >  1:5192281/1‑76 (MQ=255)
                                                     gTCGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCCTTTCCGc                        >  2:4877545/1‑76 (MQ=255)
                                                     gTCGATCATTGATGCTGGCAGGGCCGCTGGTTTTTCCTCCCCTGTGCCACAAAACTTCTCCTCCTCCCCTTACCGc                        >  2:745208/1‑76 (MQ=255)
                                                       cGATCATTGATGCTGGCAGGTGCGCCGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTATCCCCa                      >  2:2800462/1‑76 (MQ=255)
                                                       cGATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCa                      >  1:8582997/1‑76 (MQ=255)
                                                        gATCATTGATGCTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAg                     >  1:9344525/1‑76 (MQ=255)
                                                        gATCATTGAGGCTGGCAGGTGCGCTGTTTCTCCCTGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAg                     <  1:3584788/76‑1 (MQ=255)
                                                               gATGCTGGCGGGTCCGCTTTTTCTCCCTCCCTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGAttt              <  1:2748728/76‑1 (MQ=255)
                                                                 ttgtGGCGGGTGCGCTGTTCCCCCCTCCTTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGc            <  1:9962/73‑1 (MQ=255)
                                                                 tgctgcCGGGTGCTCTGTTCCCCCCCGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGc            <  1:1391659/76‑1 (MQ=255)
                                                                 tGCTGGCGGTGGCTCTGTTTCTCCCTGCCTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGc            <  1:3440546/76‑1 (MQ=255)
                                                                 tGCTGGCGGGTGCGCTGTTTCCCCCTGCTTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGc            <  1:4295931/76‑1 (MQ=255)
                                                                 tGCGGGCGGGTGCTCTGTTTCTCCCTGCTTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGc            <  1:8239591/76‑1 (MQ=255)
                                                                   cTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTg          >  2:2399255/1‑76 (MQ=255)
                                                                   cTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTg          >  2:4104422/1‑76 (MQ=255)
                                                                   cTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTg          >  2:4512423/1‑76 (MQ=255)
                                                                   cTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTg          >  2:7726887/1‑76 (MQ=255)
                                                                   cTGGCAGGGGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTg          >  2:6881714/1‑76 (MQ=255)
                                                                      gcgggTGCCCTTTTCCCCCCTGCTTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAg       <  1:4900573/73‑1 (MQ=255)
                                                                       cAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGa      >  2:1650121/1‑76 (MQ=255)
                                                                         ggTGCGTTTTTTCCCCCTCCTTTTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGAcc    <  1:7922262/76‑1 (MQ=255)
                                                                           tGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAg  >  1:7776330/1‑76 (MQ=255)
                                                                           |                                                                            
GCAGTGAACAGTTTTAACGAAGGGGTGGTTTCACCCCTTTTGTCTTTCTGGCGTCGATCATTGATGCTGGCAGGTGCGCTGTTTCTCACTGCCTGTAGCCACAACTCTTCACTTCCTCCCTTTACCGCCAGTGGATTTGCTGAAGACCAG  >  REL606/4180175‑4180324

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