Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,535,0090AC36.3% 288.1 / 24.4 146T86P (ACC→CCC) focBpredicted formate transporter
Reads supporting (aligned to +/- strand):  ref base A (71/22);  new base C (0/53);  total (71/75)
Fisher's exact test for biased strand distribution p-value = 2.19e-22
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 211 total aligned reads displayed.

TCTTCACAGGCATTAACCCATCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTGGCGAACTTGGCTGA  >  REL606/2534934‑2535084
                                                                           |                                                                           
 cTTCACAGGCATTAACCCATCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCAc                                                                            <  1:4978503/76‑1 (MQ=255)
 cTTCACAGGCATTAACCCATCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCAc                                                                            >  1:2074163/1‑76 (MQ=255)
   tCACAGGCATTAACCCATCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCa                                                                          <  2:4152721/76‑1 (MQ=255)
    cacaGGCATTAACCCATCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCAcccg                                                                         >  1:1881027/1‑74 (MQ=255)
    cacaGGCATTAACCCATCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAg                                                                         <  2:1148521/76‑1 (MQ=255)
     acaGGCATTAACCCATCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGGTTGCGGGACCAGc                                                                        >  2:7500641/1‑76 (MQ=255)
      caGGCATTAACCCATCTGGTGGGCGGCCTTTGCTTTGCCCTCGGCTTTATTTTGCTGGCGGTTTGCGGCCCCAGcc                                                                       <  2:5658429/76‑1 (MQ=255)
      caGGCATTAACCCATCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGcc                                                                       >  1:353473/1‑76 (MQ=255)
        ggCATTAACCCATCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTg                                                                     <  1:6257831/76‑1 (MQ=255)
             tAACCCATCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCAc                                                                >  1:7324281/1‑76 (MQ=255)
              aaCCCATCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTTCGGCACCAGCCTGTTCAcc                                                               >  2:8389629/1‑76 (MQ=255)
               aCCCATCTGGGGGGCGCCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCt                                                              <  2:4811191/76‑1 (MQ=255)
                cccATCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCtc                                                             >  1:1315883/1‑76 (MQ=255)
                 ccATCTGGTGGGCGGCCTTTGCTTTGCCCTCGGCTTTATTTGGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCtcg                                                            <  1:6084792/76‑1 (MQ=255)
                 ccATCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCtcg                                                            >  1:4196787/1‑76 (MQ=255)
                  cATCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCtcgt                                                           <  1:401573/76‑1 (MQ=255)
                   agcTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCGGGCGGTTTGCGGCCCCAGCCTGTTCACCtcgtc                                                          <  1:7858015/74‑1 (MQ=255)
                   aTCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCtcgtc                                                          >  2:8505186/1‑76 (MQ=255)
                     cTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGGACCAGCCTGTTCACCTCGTCgg                                                        >  1:239378/1‑76 (MQ=255)
                     cTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCGCCAGCCTGCTGACCTTGCCgg                                                        >  1:5980078/1‑76 (MQ=255)
                     cTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCgg                                                        >  1:6653532/1‑76 (MQ=255)
                      tgggggGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGt                                                       <  1:5872625/76‑1 (MQ=255)
                       ggtggGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTa                                                      <  1:8636620/76‑1 (MQ=255)
                         tggGCGTCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAAt                                                    <  1:305846/76‑1 (MQ=255)
                         tggGCGTCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAAt                                                    <  2:6789012/76‑1 (MQ=255)
                         ggggCGTCCTTTGCTTTGCACTCGGTTTTATTTTGCGGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAAt                                                    <  2:405179/75‑1 (MQ=255)
                         ggggCGGCCTTTGCTTTGCCCTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAAt                                                    <  1:683468/75‑1 (MQ=255)
                          gggCGGCCTTTGCTTTGCCCTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATg                                                   <  2:5459966/76‑1 (MQ=255)
                             cggcCTTTGCTTTGCCCTCGGCTTTATTTTGCGGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACg                                                <  1:3842873/76‑1 (MQ=255)
                             cggcCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACg                                                <  1:8068609/76‑1 (MQ=255)
                             cggcCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACg                                                <  1:3329416/76‑1 (MQ=255)
                             cggcCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACg                                                <  1:5315435/76‑1 (MQ=255)
                             cggcCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACg                                                >  1:6553655/1‑76 (MQ=255)
                             cggcCTTTGCTTTCCCCTCGGCTTTATTTGGGTGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACg                                                <  1:7450899/76‑1 (MQ=255)
                               ccttttgttttccACTGGGTTTTTTTTTGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGt                                              <  1:7927500/73‑1 (MQ=255)
                                cttttgcttcacCCTCGGTTTTATTTTGATGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTg                                             <  1:312205/74‑1 (MQ=255)
                                 ctttgctttgCTCTCGGTTTTATTTTGCTGGCGATTGGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGa                                            <  2:2185244/76‑1 (MQ=255)
                                  tttgctttgcCCTCGGCTTTTTTTTGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGAt                                           <  1:4873184/76‑1 (MQ=255)
                                  tttgctttgcACTCGGGTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGAt                                           <  1:6619768/76‑1 (MQ=255)
                                   ttgttttgcACTCGGCTTTATTTTGCTGGCGTTTTGCGCCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATg                                          <  1:5089955/76‑1 (MQ=255)
                                   ttgctttgctCTCGGTTTTTTTTGGCGGGGGGTTGGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATg                                          <  1:6779768/76‑1 (MQ=255)
                                   ttgctttgcCCTCGTCTTTATTTTGCGGGCGTTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATg                                          <  1:715941/76‑1 (MQ=255)
                                   ttgctttgcACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCCGCCTGTTCACCTCGTCGGTAATGACGGTGATg                                          >  2:1438498/1‑76 (MQ=255)
                                   ttgctttgcACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATg                                          >  1:2652316/1‑76 (MQ=255)
                                   ttgctttgcACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCAACTCGTCGGTAATGACGGGGATg                                          >  2:8113398/1‑76 (MQ=255)
                                    tgctttgcACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGGGATgg                                         >  2:3722852/1‑76 (MQ=255)
                                     gctttgcACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGGGATGGc                                        >  2:8810310/1‑76 (MQ=255)
                                      cttttcACTCGGCTTTTTTTTGCGGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCa                                       <  1:6335245/76‑1 (MQ=255)
                                      ctttccACTCGGCTTTATTTTGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCa                                       <  2:5071550/76‑1 (MQ=255)
                                       tttgcACTGGGCTTTATTTTGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCaa                                      <  2:4454811/76‑1 (MQ=255)
                                       tttgcACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCaa                                      >  1:8977581/1‑76 (MQ=255)
                                        ttgcACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCaaa                                     >  1:6966991/1‑76 (MQ=255)
                                          gcACTCGGCTTTATTTTGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCaaaaa                                   <  2:3317238/76‑1 (MQ=255)
                                           cACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCCCGTCGGTAATGACGGTGATAGCAAAAAg                                  >  1:184947/1‑76 (MQ=255)
                                             cTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTTACCTCGTCGGTAATTACGGAGATGGGAAAAAgac                                >  1:4504194/1‑74 (MQ=255)
                                             cTCGGCTTTATTTTGCGGGGGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTc                                <  1:7663360/76‑1 (MQ=255)
                                               cGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGGTGGCGAAAAGACgg                              >  2:255342/1‑76 (MQ=255)
                                                ggCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGGGATGGCAAAAAGTCggg                             >  1:2920192/1‑76 (MQ=255)
                                                ggCTTTATTGGTCGGGGGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCggg                             <  1:2955162/76‑1 (MQ=255)
                                                  cTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGc                           >  1:8888062/1‑76 (MQ=255)
                                                  cTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAAGCGGGGc                           >  2:1845689/1‑76 (MQ=255)
                                                  cTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGGAATGACGGGGATGGCAAAAAGGCGGGGc                           >  2:3219324/1‑76 (MQ=255)
                                                   tttATTTTGCTGGCGGGTTGCGGCACCAGCCTGTTCAACTCGTCGGTAATGACGGTGATGGCAACAAGTCgggggg                          >  2:7583201/1‑74 (MQ=255)
                                                    ttATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGt                         >  1:4042722/1‑76 (MQ=255)
                                                     tATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGtt                        >  1:6386916/1‑76 (MQ=255)
                                                      aTTTTGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGtta                       <  2:5492867/76‑1 (MQ=255)
                                                      aTTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGttt                       >  1:5398169/1‑75 (MQ=255)
                                                        tttGGTGGGGTTTCGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGttatt                     <  1:4922664/76‑1 (MQ=255)
                                                        tttGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGttatt                     >  1:1627532/1‑76 (MQ=255)
                                                        tttGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGttatt                     <  2:2319057/76‑1 (MQ=255)
                                                         ttGCTGGCGTTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGttatta                    <  1:1960422/76‑1 (MQ=255)
                                                         ttGCTGGCGGTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGttatta                    <  1:467234/76‑1 (MQ=255)
                                                         ttGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGGGATGGCAAAAAGTCGGGGGGttatta                    >  1:321205/1‑76 (MQ=255)
                                                         ttGCGGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGttatta                    <  2:733835/76‑1 (MQ=255)
                                                         ttGCGGGCGGTTTGCCGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGttatta                    <  2:8538657/76‑1 (MQ=255)
                                                         tggCTGGGGTTTTGCGCCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGttatta                    <  1:5728852/74‑1 (MQ=255)
                                                          tGCTGGGGGTTTGGGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAg                   <  2:4373143/76‑1 (MQ=255)
                                                          tGCTGGCGTTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAg                   <  1:8605875/76‑1 (MQ=255)
                                                          tGCGGGCGTTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAg                   <  1:6033527/76‑1 (MQ=255)
                                                           gCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGGGATGGCAAAAAGTCGGGGCGGTTTTAGt                  >  2:3448333/1‑76 (MQ=255)
                                                            cTGGCGTTTTGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGtt                 <  2:5533593/76‑1 (MQ=255)
                                                            cTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGtt                 >  1:4323434/1‑76 (MQ=255)
                                                            cTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGGTATTAGtt                 >  2:2436166/1‑76 (MQ=255)
                                                             tGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGGTATTAGTTg                >  2:694977/1‑76 (MQ=255)
                                                             tGGCGGATTGCGGCGCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTg                <  2:1515670/76‑1 (MQ=255)
                                                              ggCGTTTGGCGGCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTgg               <  1:5424341/76‑1 (MQ=255)
                                                               gCGGTGTGCGACCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATAGCAAAAAGTCGGGGCGTTATTAGTTGGc              <  2:6175950/76‑1 (MQ=255)
                                                                gggTTTGCGACCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTGGCg             <  1:3977007/75‑1 (MQ=255)
                                                                cgttttGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTCATTAGTTGGCg             <  2:229817/73‑1 (MQ=255)
                                                                cGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTGGCg             >  1:2537208/1‑76 (MQ=255)
                                                                cGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTGGCg             >  1:2124802/1‑76 (MQ=255)
                                                                 ggTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTGGCGa            >  1:3297455/1‑76 (MQ=255)
                                                                 ggTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTGGCGa            >  1:3435376/1‑76 (MQ=255)
                                                                  gTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTGGCGaa           >  1:6181643/1‑76 (MQ=255)
                                                                  gTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTGGCGaa           >  2:818669/1‑76 (MQ=255)
                                                                   tttGCGCCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTGGCGAAc          <  1:2732207/76‑1 (MQ=255)
                                                                     tGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTGGCGAACtt        >  1:6593450/1‑76 (MQ=255)
                                                                     tGCGCCCCCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTGGCGAACtt        <  1:4799787/76‑1 (MQ=255)
                                                                     ggcgcCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTGGCGAACtt        <  2:576749/75‑1 (MQ=255)
                                                                      gCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTGGCGAACTTg       >  1:763866/1‑76 (MQ=255)
                                                                           |                                                                           
TCTTCACAGGCATTAACCCATCTGGTGGGCGGCCTTTGCTTTGCACTCGGCTTTATTTTGCTGGCGGTTTGCGGCACCAGCCTGTTCACCTCGTCGGTAATGACGGTGATGGCAAAAAGTCGGGGCGTTATTAGTTGGCGAACTTGGCTGA  >  REL606/2534934‑2535084

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: