Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 578,568 G→A *154* (TAG→TAA rzpD → DLP12 prophage, putative murein endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913578,5680GA100.0% 242.0 / NA 74*154* (TAG→TAArzpDDLP12 prophage, putative murein endopeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/35);  total (39/35)

ATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT  >  NC_000913/578444‑578698
                                                                                                                            |                                                                                                                                  
aTGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCg                                                                                                                      <  2:672815/139‑1 (MQ=255)
      ccTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGc                                                                                                                >  2:747893/1‑139 (MQ=255)
      ccTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGc                                                                                                                >  2:631240/1‑139 (MQ=255)
      ccTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGc                                                                                                                >  2:395489/1‑139 (MQ=255)
                 cTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCaa                                                                                                     <  1:631240/139‑1 (MQ=255)
                      aGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAattatt                                                                                                <  2:1237272/139‑1 (MQ=255)
                            cGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGc                                                                                          <  1:747893/139‑1 (MQ=255)
                              cTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCaa                                                                                        >  2:28984/1‑139 (MQ=255)
                              cTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCaa                                                                                        >  1:817631/1‑139 (MQ=255)
                                               aCCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                       <  1:822139/139‑1 (MQ=255)
                                                cccTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGAACATATCGATGGGCAACTCATGCAATTATTTTGAGCAAGACACACGCGCTTCCAga                                                                      >  1:1355466/1‑138 (MQ=255)
                                                 ccTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                     <  1:395489/139‑1 (MQ=255)
                                                   tCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCaa                                                                                                     >  1:1270157/1‑105 (MQ=255)
                                                   tCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCaa                                                                                                     <  2:1270157/105‑1 (MQ=255)
                                                    cAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAGTACACACGCGCTTCCAGCGGAg                                                                  >  2:1485625/1‑139 (MQ=255)
                                                     agagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                     <  1:803079/135‑1 (MQ=255)
                                                     agagagagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                     >  2:803079/1‑135 (MQ=255)
                                                         agagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCaa                                                                                                     >  1:960171/1‑99 (MQ=255)
                                                         agagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCaa                                                                                                     <  2:960171/99‑1 (MQ=255)
                                                         agagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATaa                                                             <  2:1320467/139‑1 (MQ=255)
                                                         agagGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATaa                                                             <  1:28984/139‑1 (MQ=255)
                                                           agGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATaaag                                                                                                                                  >  2:1005387/1‑68 (MQ=255)
                                                           agGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATaaag                                                                                                                                  <  1:1005387/68‑1 (MQ=255)
                                                           agGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAAt                                                           >  2:45975/1‑139 (MQ=255)
                                                             gCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCg                                                                                                                      <  2:220616/78‑1 (MQ=255)
                                                             gCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCg                                                                                                                      >  1:220616/1‑78 (MQ=255)
                                                             gCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacacac                                                                                 >  1:452559/1‑115 (MQ=255)
                                                             gCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATacacac                                                                                 <  2:452559/115‑1 (MQ=255)
                                                             gCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGc                                                         >  2:331968/1‑139 (MQ=255)
                                                                 aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa                                                     <  2:1420353/139‑1 (MQ=255)
                                                                 aTCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa                                                     >  2:627942/1‑139 (MQ=255)
                                                                      tATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAg                                                   <  2:1792415/136‑1 (MQ=255)
                                                                      tATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAg                                                   >  1:1792415/1‑136 (MQ=255)
                                                                      tATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGtaa                                                >  2:657356/1‑139 (MQ=255)
                                                                      tATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGtaa                                                <  2:817631/139‑1 (MQ=255)
                                                                       aTGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCaa                                                                                        >  1:781618/1‑98 (MQ=255)
                                                                       aTGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCaa                                                                                        <  2:781618/98‑1 (MQ=255)
                                                                        tGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATaa                                                             <  2:59814/124‑1 (MQ=255)
                                                                        tGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATaa                                                             >  1:59814/1‑124 (MQ=255)
                                                                           aaaaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATaaaa                                           >  1:883232/1‑139 (MQ=255)
                                                                             aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa                                                     >  2:119728/1‑127 (MQ=255)
                                                                             aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAAcc                                         >  1:176731/1‑139 (MQ=255)
                                                                             aaaCAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTAACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTaa                                                     <  1:119728/127‑1 (MQ=255)
                                                                                caacTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAg                                      <  2:176731/139‑1 (MQ=255)
                                                                                  acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAattatt                                                                                                >  2:790962/1‑79 (MQ=255)
                                                                                  acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAattatt                                                                                                <  1:790962/79‑1 (MQ=255)
                                                                                  acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAGTACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCa                                    <  1:1485625/139‑1 (MQ=255)
                                                                                  acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAg                                      >  2:36881/1‑137 (MQ=255)
                                                                                  acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAg                                      <  1:36881/137‑1 (MQ=255)
                                                                                  acTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCa                                    >  1:1771526/1‑139 (MQ=255)
                                                                                    tGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAAt                                  <  1:627942/139‑1 (MQ=255)
                                                                                     ggaaggaaCCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATc                                 >  1:157165/1‑139 (MQ=255)
                                                                                            aCCCAGAAGTATATTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTAc                          <  2:1771526/139‑1 (MQ=255)
                                                                                                        atTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAAt                                  <  2:432190/119‑1 (MQ=255)
                                                                                                        atTAATGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAAt                                  >  1:432190/1‑119 (MQ=255)
                                                                                                            aTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                       >  2:1487244/1‑78 (MQ=255)
                                                                                                            aTGAGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                       <  1:1487244/78‑1 (MQ=255)
                                                                                                               aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAg                                                                                           <  1:918784/55‑1 (MQ=255)
                                                                                                               aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAg                                                                                           >  2:918784/1‑55 (MQ=255)
                                                                                                               aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                       <  2:1686334/75‑1 (MQ=255)
                                                                                                               aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                       >  1:1686334/1‑75 (MQ=255)
                                                                                                               aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGcc                                                        >  1:71280/1‑90 (MQ=255)
                                                                                                               aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGcc                                                        <  2:71280/90‑1 (MQ=255)
                                                                                                               aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCt       >  2:817107/1‑139 (MQ=255)
                                                                                                               aGCAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCt       >  1:284628/1‑139 (MQ=255)
                                                                                                                 cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCaa                                                                                                     <  1:486093/43‑1 (MQ=255)
                                                                                                                 cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCaa                                                                                                     >  2:486093/1‑43 (MQ=255)
                                                                                                                 cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                       >  2:689582/1‑73 (MQ=255)
                                                                                                                 cAGTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAg                                                                       <  1:689582/73‑1 (MQ=255)
                                                                                                                   gTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                     >  2:1784640/1‑73 (MQ=255)
                                                                                                                   gTGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCg                                                                     <  1:1784640/73‑1 (MQ=255)
                                                                                                                    tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc                                                                               <  2:754423/62‑1 (MQ=255)
                                                                                                                    tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACAcgc                                                                               >  1:754423/1‑62 (MQ=255)
                                                                                                                    tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGcc                                                        <  2:204807/85‑1 (MQ=255)
                                                                                                                    tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGcc                                                        >  1:204807/1‑85 (MQ=255)
                                                                                                                    tGCAGATAAAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGtttt  >  2:55183/1‑139 (MQ=255)
                                                                                                                            |                                                                                                                                  
ATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTT  >  NC_000913/578444‑578698

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