Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,938,441 Δ1 bp coding (215/993 nt) rbsR → DNA‑binding transcriptional dual regulator RbsR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,938,4410C.100.0% 167.5 / NA 36coding (215/993 nt)rbsRDNA‑binding transcriptional dual regulator RbsR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (21/15);  total (21/15)

TTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCG  >  NC_000913/3938204‑3938683
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
ttttAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAAc                                                                                                                                                                                                                                      >  1:249844/1‑251 (MQ=255)
      aCAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGc                                                                                                                                                                                                                                <  1:287856/251‑1 (MQ=255)
                        tgggtACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCACAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTGGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAg                                                                                                                                                                                                                <  2:251131/245‑1 (MQ=255)
                               aaTGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAACCACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCg                                                                                                                                                                                                        >  1:120171/1‑250 (MQ=255)
                                       aTGTTGCCCGCATGGCAAGAGTTTCTACCACAACAGTTTCTCACGATATCACTAAAGATCGCTTCGTCAAGGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTCAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAACCTCAATCCAACACATACCATTGGCATGCTGATCGCTGCCAGAACCAATCTTTCTATTCCGAACTGGTGCGTGGCGTTGACAGCCGCTGCTTCGTACGCGGtt                                                                                                                                                                                               >  2:265175/1‑251 (MQ=255)
                                                  cgggcgggcgTGTCTATCTCAACAGTTTCGCACGGGATCAATAAAGATCGCCTCGTCAGTGAAGCGATTACCGCCAAAGGTGAAGCGGCGATTAAAGAACTCAATTACGCGCCAGCAGCTCTGGCGCGTAGCCTAAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCGGTACCAATCTTTCTATTCAGATCTGGTGCGTGGCGTTGAACGCAGCTGTTTCGAACGCGGTTATAGTCTCGTc                                                                                                                                                                                    <  2:264786/249‑1 (MQ=255)
                                                  cTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGACCGCGGTTATAGTCTCGTc                                                                                                                                                                                    >  1:5439/1‑251 (MQ=255)
                                                    ggcgggcgTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATGAAAGAACTCAATGACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAAGCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGGTGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCt                                                                                                                                                                                  <  1:259457/251‑1 (MQ=255)
                                                       gggcgTTGCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTg                                                                                                                                                                               <  2:1608/251‑1 (MQ=255)
                                                       gggcgATTCTACCTCAACAGTTTCTCACGTTCTCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAACGTTGAACCGGCGATTAAAGAACTCAATTCCGCGCCATCAGCTCTGGCGCGTAGCCTCAACCTCAATCACACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGCACTGGTGCGTGGCGTTGAACGCCACTGCTTCGAACGCGGTTCTAGTCTCGTCCTTTg                                                                                                                                                                               >  2:163188/1‑251 (MQ=255)
                                                                                  cGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAACCTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGACCAGCGGAt                                                                                                                                                    >  2:313874/1‑251 (MQ=255)
                                                                                              aGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCt                                                                                                                                        <  1:601/251‑1 (MQ=255)
                                                                                                   gCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGaa                                                                                                                                    >  1:77331/1‑250 (MQ=255)
                                                                                                       cGTCAGTGAAGCTATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCt                                                                                                                               <  1:375799/251‑1 (MQ=255)
                                                                                                              gAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTCGCCTCAAACTCAATCAAACACATAACATTGGCATGTTGATCACTGCCAGTCCCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCa                                                                                                                         >  1:437904/1‑250 (MQ=255)
                                                                                                                          gCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGAt                                                                                                            >  1:127604/1‑251 (MQ=255)
                                                                                                                             aaaGGTGAAGCGGCGGTTAAAGTTCTCAATTACGCGCCATCACCTCTGGCGCTTAGCCTTCAACTCAATCAAACACATACCATTGGAATGTTGATCACTGCCAGTGCCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCGGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGCAGGCGATGAACAGCGGATGAATCGCAATGTGGAAACGTTGATGCAAAAACGCGCTGATGGc                                                                                                         >  2:248576/1‑251 (MQ=255)
                                                                                                                                           gATTAAAGAACTCAGTTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACAGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACtgtg                                                                                           >  1:144854/1‑251 (MQ=255)
                                                                                                                                                    aCTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAAc                                                                                  >  2:243338/1‑251 (MQ=255)
                                                                                                                                                       cAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACg                                                                                 <  1:257864/249‑1 (MQ=255)
                                                                                                                                                       cAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGc                                                                                <  1:19181/250‑1 (MQ=255)
                                                                                                                                                           tACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCa                                                                            >  2:26386/1‑250 (MQ=255)
                                                                                                                                                                  cATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTcgc                                                                    >  1:84382/1‑251 (MQ=255)
                                                                                                                                                                      aGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGa                                                                <  1:304813/251‑1 (MQ=255)
                                                                                                                                                                            ggCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCa                                                           <  1:73671/250‑1 (MQ=255)
                                                                                                                                                                                        cAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATc                                               <  2:310310/250‑1 (MQ=255)
                                                                                                                                                                                               aatcaaACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTg                                       >  2:367878/1‑251 (MQ=255)
                                                                                                                                                                                                   aaaCACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGa                                                                                     >  1:348273/1‑201 (MQ=255)
                                                                                                                                                                                                   aaaCACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGa                                                                                     <  2:348273/201‑1 (MQ=255)
                                                                                                                                                                                                   aaaCACATACCATTGGCATGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCt                                    >  1:422418/1‑250 (MQ=255)
                                                                                                                                                                                                                   aaTGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGAGATGATGGACTg                   >  2:290716/2‑251 (MQ=255)
                                                                                                                                                                                                                     tGTTGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTggg                 >  1:170599/1‑251 (MQ=255)
                                                                                                                                                                                                                        tGATCACTGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCt               <  1:355709/250‑1 (MQ=255)
                                                                                                                                                                                                                               tGCCAGTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCg        >  2:178807/1‑250 (MQ=255)
                                                                                                                                                                                                                                  cAGTACCAATCTTTCTATTCAGAACTGGTGCGGGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATgg    <  2:437904/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gTACCAATCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCg  >  2:387475/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
TTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCG  >  NC_000913/3938204‑3938683

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