Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,897,109 C→T H60Y (CAC→TAC)  sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,897,1090CT100.0% 219.5 / NA 62H60Y (CAC→TAC) sdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/36);  total (26/36)

ATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGCACAGGGA  >  NC_000913/1896976‑1897234
                                                                                                                                     |                                                                                                                             
aTCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACcaca                                                                                                                         <  1:218204/140‑1 (MQ=255)
  cTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACAcc                                                                                                                       <  1:874476/140‑1 (MQ=255)
      ccATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGata                                                                                                                   <  1:746945/140‑1 (MQ=255)
       cATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGatat                                                                                                                  <  2:660622/140‑1 (MQ=255)
       cATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGatat                                                                                                                  <  2:745516/140‑1 (MQ=255)
         tACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCg                                                                                                                >  2:253773/1‑140 (MQ=255)
          aCCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGc                                                                                                               <  2:1176764/140‑1 (MQ=255)
          aCCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGc                                                                                                               <  1:802336/140‑1 (MQ=255)
           ccGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGcc                                                                                                              <  1:483505/140‑1 (MQ=255)
              tAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCatt                                                                                                           >  2:1129764/1‑140 (MQ=255)
              tAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCatt                                                                                                           >  1:462278/1‑140 (MQ=255)
               aGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCatta                                                                                                          >  2:665632/1‑140 (MQ=255)
                gggCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCattat                                                                                                         <  2:462278/140‑1 (MQ=255)
                 ggCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCattatt                                                                                                        >  1:646643/1‑140 (MQ=255)
                   ccTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATtattat                                                                                                      <  1:665632/140‑1 (MQ=255)
                     tATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATgg                                                                                                    >  1:1124691/1‑140 (MQ=255)
                       tGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATg                                                                                                     <  2:706038/137‑1 (MQ=255)
                       tGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATg                                                                                                     >  1:706038/1‑137 (MQ=255)
                       tGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGt                                                                                                  <  2:1124691/140‑1 (MQ=255)
                          aGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCtt                                                                                               <  2:874343/140‑1 (MQ=255)
                           ggCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTg                                                                                              <  1:956538/140‑1 (MQ=255)
                           ggCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTg                                                                                              <  1:978668/140‑1 (MQ=255)
                            gCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGc                                                                                             >  1:453213/1‑140 (MQ=255)
                                gTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTg                                                                                              >  2:113129/1‑135 (MQ=255)
                                gTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTg                                                                                              <  1:113129/135‑1 (MQ=255)
                                gTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGt                                                                                         >  2:633822/1‑140 (MQ=255)
                                gTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGt                                                                                         >  1:1083526/1‑140 (MQ=255)
                                gTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGt                                                                                         >  1:988622/1‑140 (MQ=255)
                                        ttCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc                                                                                 >  2:373474/1‑140 (MQ=255)
                                        ttCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAAcc                                                                                 >  2:1127199/1‑140 (MQ=255)
                                          cgtcgATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTg                                                                               >  1:262473/1‑140 (MQ=255)
                                           gtcgATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc                                                                              <  2:969052/140‑1 (MQ=255)
                                           gtcgATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGc                                                                              <  2:453213/140‑1 (MQ=255)
                                             cgATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCa                                                                            <  1:554437/140‑1 (MQ=255)
                                                       gTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGatat                                                                  >  2:557046/1‑140 (MQ=255)
                                                       gTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGatat                                                                  <  1:190883/140‑1 (MQ=255)
                                                       gTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGatat                                                                  >  1:608004/1‑140 (MQ=255)
                                                         cGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCg                                                                <  1:1019551/140‑1 (MQ=255)
                                                          gAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGa                                                               <  1:193223/140‑1 (MQ=255)
                                                            aaaaGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc                                                                           >  2:935841/1‑126 (MQ=255)
                                                            aaaaGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCAc                                                                           <  1:935841/126‑1 (MQ=255)
                                                            aaaaGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACa                                                             >  2:704670/1‑140 (MQ=255)
                                                                               gTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcg                                          <  1:1127199/140‑1 (MQ=255)
                                                                               gTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcg                                          <  1:318237/140‑1 (MQ=255)
                                                                               gTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcg                                          <  1:633822/140‑1 (MQ=255)
                                                                               gTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATGcg                                          >  2:744988/1‑140 (MQ=255)
                                                                                ttACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTc                                                       >  2:603152/1‑126 (MQ=255)
                                                                                ttACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTc                                                       <  1:603152/126‑1 (MQ=255)
                                                                                ttACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgc                                         >  1:737810/1‑140 (MQ=255)
                                                                                 tACTCGTGTTGCCGTGGGCGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgcg                                        <  2:102532/140‑1 (MQ=255)
                                                                                 tACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTcgcg                                        <  2:1107850/140‑1 (MQ=255)
                                                                                   cTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTcc                                                      <  1:875364/124‑1 (MQ=255)
                                                                                   cTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTcc                                                      >  2:875364/1‑124 (MQ=255)
                                                                                                           gTTCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgctg              <  2:988622/140‑1 (MQ=255)
                                                                                                            ttCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgctgc             <  1:932868/140‑1 (MQ=255)
                                                                                                             tCACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgctgct            <  1:1015730/140‑1 (MQ=255)
                                                                                                              cACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgctgctg           >  1:817032/1‑140 (MQ=255)
                                                                                                              cACTGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTctgctgctg           >  1:925585/1‑140 (MQ=255)
                                                                                                                 tGTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGca        >  2:912016/1‑140 (MQ=255)
                                                                                                                  gTCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGcac       <  1:557046/140‑1 (MQ=255)
                                                                                                                   tCGCTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGcaca      >  2:707276/1‑140 (MQ=255)
                                                                                                                      cTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGCACAggg   <  1:704670/140‑1 (MQ=255)
                                                                                                                      cTGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGCACAggg   <  1:707276/140‑1 (MQ=255)
                                                                                                                       tGACGGGTAAAGGCTACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGCACAGGGa  <  2:609498/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                             
ATCTTCCCATACCGTAGGGCCTATGAAGGCAGGTAAACAGTTCGTCGATGATCTGGTCGAAAAAGGCTTACTGGATAGCGTTACTCGCGTTGCCGTGGACGTTTATGGTTCACTGTCGCTGACGGGTAAAGGCCACCACACCGATATCGCCATTATTATGGGTCTTGCAGGTAACGAACCTGCCACCGTGGATATCGACAGTATTCCCGGTTTTATTCGCGACGTAGAAGAGCGCGAACGTCTGCTGCTGGCACAGGGA  >  NC_000913/1896976‑1897234

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