Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,072,818 Δ1 bp 14.8% coding (689/1308 nt) mdtA ← multidrug transporter membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,072,8180A.14.8% 135.8 / 23.0 54coding (689/1308 nt)mdtAmultidrug transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base A (20/25);  new base . (4/4);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

CACCGGTGGATTGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAGCCCCT  >  NC_002947/4072669‑4072965
                                                                                                                                                     |                                                                                                                                                    
cacCGGTGGATTGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCAccacca                                                                                                                                                      <  1:403264/150‑1 (MQ=255)
    ggTGGATTGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAgcgc                                                                                                                                                  <  1:791173/150‑1 (MQ=255)
    ggTGGATTGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAgcgc                                                                                                                                                  >  1:348937/1‑150 (MQ=255)
        gATTGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCggtg                                                                                                                                             <  2:348937/150‑1 (MQ=255)
        gATTGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCggtg                                                                                                                                             <  1:658314/150‑1 (MQ=255)
          ttGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCggtggt                                                                                                                                           >  1:112108/1‑150 (MQ=255)
           tGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAACGTGACGTTGATCGGCTCGGTCTGGGTGATCACCACCGCGCTGgtggtg                                                                                                                                          <  2:120398/150‑1 (MQ=255)
                     tGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTggcgg                                                                                                                                >  2:823457/1‑150 (MQ=255)
                         aCCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGccac                                                                                                                            >  1:85748/1‑150 (MQ=255)
                                  ccAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCGCGCTGgtggtgat                                                                                                                                        >  1:58286/1‑127 (MQ=255)
                                  ccAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCGCGCTGgtggtgat                                                                                                                                        <  2:58286/129‑3 (MQ=255)
                                   cAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCGCGCTGgtggtga                                                                                                                                         >  1:598753/1‑126 (MQ=255)
                                   cAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCGCGCTGgtggtga                                                                                                                                         <  2:598753/127‑2 (MQ=255)
                                    aGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCa                                                                                                                 <  2:85748/150‑1 (MQ=255)
                                              cGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGtt                                                                                                                                     >  1:771909/1‑120 (MQ=255)
                                              cGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGtt                                                                                                                                     <  2:771909/120‑1 (MQ=255)
                                                gCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGgc                                                                                                     >  1:521708/1‑150 (MQ=255)
                                                    gCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAg                                                                                                 <  2:568571/150‑1 (MQ=255)
                                                     cTGCCGCCACTGCGGTAGCGCTCCAGCCCAGTGCTCAGTGGGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGTTGGCGTTGCCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAgg                                                                                                <  2:274455/150‑1 (MQ=255)
                                                       ggCGCCACTGCGGTGGCGCGCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCCCCAGGTTGCCCAGATCGAGCTGGCGCAGGcc                                                                                              <  1:84780/150‑1 (MQ=255)
                                                       ggCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGcc                                                                                              <  2:521708/150‑1 (MQ=255)
                                                             aCTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACg                                                                                        >  2:799850/1‑150 (MQ=255)
                                                                  ggTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCGGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCAc                                                                                   <  1:823457/150‑1 (MQ=255)
                                                                      gcgcTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGAt                                                                               <  1:550410/150‑1 (MQ=255)
                                                                            cAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCa                                                                                    <  2:740608/139‑1 (MQ=255)
                                                                            cAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCa                                                                                    >  1:740608/1‑139 (MQ=255)
                                                                                  aGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGAt                                                                   <  1:799850/150‑1 (MQ=255)
                                                                                       tCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCCGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCAAGAACGAGCTAGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCAggg                                                              >  2:623636/1‑150 (MQ=255)
                                                                                       tCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTggg                                                              >  1:105951/1‑150 (MQ=255)
                                                                                       tCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTggg                                                              >  1:62726/1‑150 (MQ=255)
                                                                                           ttCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGaa                                                          >  1:288542/1‑150 (MQ=255)
                                                                                             cGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAAt                                                        <  2:773673/150‑1 (MQ=255)
                                                                                                tttCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATggg                                                                            >  2:166345/1‑127 (MQ=255)
                                                                                                tttCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATggg                                                                            <  1:166345/127‑1 (MQ=255)
                                                                                                tttCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTgg                                                               <  2:654930/140‑1 (MQ=255)
                                                                                                tttCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTgg                                                               >  1:654930/1‑140 (MQ=255)
                                                                                                tttCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGa                                                     <  1:623636/150‑1 (MQ=255)
                                                                                                 ttCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGc                                                                                                                    >  2:212302/1‑86 (MQ=255)
                                                                                                 ttCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGc                                                                                                                    <  1:212302/86‑1 (MQ=255)
                                                                                                 ttCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTCGTGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAg                                                    <  1:794479/150‑1 (MQ=255)
                                                                                                    cGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAACCGAGGtt                                                 <  2:105951/150‑1 (MQ=255)
                                                                                                     ggCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAg                                                    >  1:242260/1‑146 (MQ=255)
                                                                                                     ggCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAg                                                    <  2:242260/146‑1 (MQ=255)
                                                                                                        aGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAgg                                             >  2:31422/1‑150 (MQ=255)
                                                                                                                ggCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGcgtc                                     >  1:664049/1‑150 (MQ=255)
                                                                                                                      gTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGAc                               >  2:304625/1‑150 (MQ=255)
                                                                                                                      gTTGATCGGCCAGGTCTGGGTACTCACCACCAGAGCGGTGGTGTCGTTGGCGGCCACCAGCTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCAGTACATTGGCACACGGTTCTGCGTGAAAGCGAGGTTCAGGCGGGCGTGGTTGAc                               <  2:120647/150‑1 (MQ=255)
                                                                                                                       ttGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGAcc                              >  1:456390/1‑150 (MQ=255)
                                                                                                                             ggCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGAcagcaagc                        <  2:771574/150‑8 (MQ=255)
                                                                                                                                       gggTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCt              <  2:664049/150‑1 (MQ=255)
                                                                                                                                            gacaccacCGCGCTGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGc                                                                                               <  2:751702/62‑1 (MQ=37)
                                                                                                                                            gacaccacCGCGCTGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGc                                                                                               >  1:751702/3‑64 (MQ=37)
                                                                                                                                             tCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCa        <  2:288542/150‑1 (MQ=255)
                                                                                                                                              caccacCGCGCTGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGACTGATTGGATTTGACCAg       >  2:832316/1‑150 (MQ=255)
                                                                                                                                               accaccAGCGCGGTGGCGTAGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAGc      >  2:780120/1‑150 (MQ=255)
                                                                                                                                                ccaccaGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCa        <  2:563618/147‑1 (MQ=255)
                                                                                                                                                ccaccaGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCa        >  1:563618/1‑147 (MQ=255)
                                                                                                                                                   ccaGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAGCCCCt  >  2:414531/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CACCGGTGGATTGCAGCTTGCTGTCACCGCGGTCCCAGGCCTCGACCGGCAGGCTGGCGCCACTGCGGTAGCGCTCCAGCACAGTGCTCAGTTCGGTTTCCGGCAGGGTGAAGGCGACGTTGATCGGCTCGGTCTGGGTGATCACCACCAGCGCGGTGGTGTCGTTGGCGGCCACCAGGTTGCCCAGATCGAGCTGGCGCAGGCCCACACGGCCACTGATGGGCGCACGGATCTGGGTGAAATCGAGGTTCAGGCGGGCGTCGTTGACCTGTGCCTGATTGGTCTTGACCAGCCCCT  >  NC_002947/4072669‑4072965

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