Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 511,219 T→C E54G (GAA→GGA) 
E54G (GAA→GGA) 		 copA ←
copA ←		 Cu(+) exporting P‑type ATPase
Cu(+) exporting P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913511,2190TC100.0% 226.2 / NA 61E54G (GAA→GGA) 
E54G (GAA→GGA) 		copA
copA		Cu(+) exporting P‑type ATPase
Cu(+) exporting P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/30);  total (31/30)

TCGTCATCGGCGGTCGCTGCCGGAAGCGCCTCAGAAACCGCTGTCAGTGCTTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGACAGGCCGTC  >  NC_000913/511087‑511352
                                                                                                                                    |                                                                                                                                     
tcgtcatcgGCGGTCGCTGCCGGAAGCGCCTCAGAAACCGCTGTCAGTGCTTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAg                                                                                                                                <  2:256026/140‑1 (MQ=255)
      tcgGCGGTCGCTGCCGGAAGCGCCTCAGAAACCGCTGTCAGTGCTTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTc                                                                                                                          <  1:1254420/140‑1 (MQ=255)
              cGCTGCCGGAAGCGCCTCAGAAACCGCTGTCAGTGCTTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTgg                                                                                                                  <  1:384268/140‑1 (MQ=255)
                  gCCGGAAGCGCCTCAGAAACCGCTGTCAGTGCTTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGt                                                                                                              <  1:83717/140‑1 (MQ=255)
                             cTCAGAAACCGCTGTCAGTGCTTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGt                                                                                                   >  1:1330564/1‑140 (MQ=255)
                              tCAGAAACCGCTGTCAGTGCTTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTg                                                                                                  <  1:993983/140‑1 (MQ=255)
                              tCAGAAACCGCTGTCAGTGCTTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTg                                                                                                  <  1:147399/140‑1 (MQ=255)
                                aGAAACCGCTGTCAGTGCTTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTgcg                                                                                                >  2:478355/1‑140 (MQ=255)
                                   aaCCGCTGTCAGTGCTTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCtt                                                                                             <  1:1063459/140‑1 (MQ=255)
                                     ccGCTGTCAGTGCTTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCa                                                                                           <  1:182100/140‑1 (MQ=255)
                                         tGTCAGTGCTTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCtt                                                                                             <  1:146406/134‑1 (MQ=255)
                                         tGTCAGTGCTTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCtt                                                                                             >  2:146406/1‑134 (MQ=255)
                                             aGTGCTTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGa                                                                                   >  1:1209048/1‑140 (MQ=255)
                                                   tCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATc                                                                             >  1:484290/1‑140 (MQ=255)
                                                    ccGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAg                                                                                    >  1:1115125/1‑132 (MQ=255)
                                                    ccGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAg                                                                                    <  2:1115125/132‑1 (MQ=255)
                                                     cGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCg                                                                           >  2:280868/1‑140 (MQ=255)
                                                      gACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGc                                                                          <  1:522193/140‑1 (MQ=255)
                                                       aCGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATcc                                                                            <  1:1294824/137‑1 (MQ=255)
                                                       aCGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATcc                                                                            >  2:1294824/1‑137 (MQ=255)
                                                           gATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCt                                                                     <  1:561058/140‑1 (MQ=255)
                                                           gATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCt                                                                     <  1:577444/140‑1 (MQ=255)
                                                            aTTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTc                                                                    >  2:704186/1‑140 (MQ=255)
                                                                   aCTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCaa                                                                  <  1:116043/135‑1 (MQ=255)
                                                                   aCTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCaa                                                                  >  2:116043/1‑135 (MQ=255)
                                                                     tCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCgg                                                           >  2:75811/1‑140 (MQ=255)
                                                                     tCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCgg                                                           >  1:801630/1‑140 (MQ=255)
                                                                     tCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCgg                                                           >  2:135488/1‑140 (MQ=255)
                                                                     tCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCgg                                                           >  1:555238/1‑140 (MQ=255)
                                                                      ccGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCaa                                                                  <  1:1032243/132‑1 (MQ=255)
                                                                      ccGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCaa                                                                  >  2:1032243/1‑132 (MQ=255)
                                                                       cGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGAc                                                         >  1:1195551/1‑140 (MQ=255)
                                                                       cGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGAc                                                         >  1:536427/1‑140 (MQ=255)
                                                                             cGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGa                                                          <  2:1355089/133‑1 (MQ=255)
                                                                             cGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGa                                                          >  1:1355089/1‑133 (MQ=255)
                                                                             cGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATAAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGtt                                                   >  1:716184/1‑140 (MQ=255)
                                                                               gTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCa                                                 <  2:1195551/140‑1 (MQ=255)
                                                                                  tAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGa                                              >  1:881681/1‑140 (MQ=255)
                                                                                  tAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGa                                              >  1:695368/1‑140 (MQ=255)
                                                                                       tttGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGActttc                                         >  2:153928/1‑140 (MQ=255)
                                                                                            gTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCa                                    <  2:536427/140‑1 (MQ=255)
                                                                                               gCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCAcgc                                 >  2:673523/1‑140 (MQ=255)
                                                                                                         gCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCa                       >  1:985805/1‑140 (MQ=255)
                                                                                                         gCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCa                       >  2:22201/1‑140 (MQ=255)
                                                                                                         gCGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCa                       >  2:261764/1‑140 (MQ=255)
                                                                                                          cGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAg                      >  2:81785/1‑140 (MQ=255)
                                                                                                          cGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGAAACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCAc                                   <  1:917034/127‑1 (MQ=255)
                                                                                                          cGTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGAAACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCAc                                   >  2:917034/1‑127 (MQ=255)
                                                                                                           gTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGt                     <  2:881681/140‑1 (MQ=255)
                                                                                                           gTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGt                     <  1:22201/140‑1 (MQ=255)
                                                                                                           gTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGt                     <  2:695368/140‑1 (MQ=255)
                                                                                                           gTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGt                     <  2:985805/140‑1 (MQ=255)
                                                                                                           gTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGt                     <  1:261764/140‑1 (MQ=255)
                                                                                                           gTCATAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGt                     <  1:75811/140‑1 (MQ=255)
                                                                                                              aTAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGt                     >  2:364562/1‑137 (MQ=255)
                                                                                                              aTAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGt                     <  1:364562/137‑1 (MQ=255)
                                                                                                               tAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGt                     <  2:853511/136‑1 (MQ=255)
                                                                                                               tAACCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGt                     >  1:853511/1‑136 (MQ=255)
                                                                                                                aaCCCGCTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCg                <  1:81785/140‑1 (MQ=255)
                                                                                                                      cTTGTTTGATGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGAc          >  1:466714/1‑140 (MQ=255)
                                                                                                                              aTGGTTCCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGACAGGCCGTc  <  1:280868/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                     
TCGTCATCGGCGGTCGCTGCCGGAAGCGCCTCAGAAACCGCTGTCAGTGCTTCCGACGGGATTGATGACTCCGCCAGCGGTTTAGCCTTTGGGTGGCTTACAGATGCGTCATAACCCGCTTGTTTGATGGTTTCGATCAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGACAGGCCGTC  >  NC_000913/511087‑511352

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑