Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,646,189 G→C A155G (GCA→GGA)  pepP ← M24B family Xaa‑Pro dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,646,1890GC100.0% 229.5 / NA 64A155G (GCA→GGA) pepPM24B family Xaa‑Pro dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (38/26);  total (38/26)

ATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATAT  >  CP000730/1646056‑1646322
                                                                                                                                     |                                                                                                                                     
aTCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTac                                                                                                                                >  1:211889/1‑141 (MQ=255)
 tCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTGATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTaca                                                                                                                               >  2:372233/1‑141 (MQ=255)
  cACTTGCAACACGATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaa                                                                                                                              <  1:163522/141‑1 (MQ=255)
   aCTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaac                                                                                                                             >  2:524731/1‑141 (MQ=255)
    cTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaaca                                                                                                                            >  1:157366/1‑141 (MQ=255)
     ttGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAg                                                                                                                           >  2:519498/1‑141 (MQ=255)
      tGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGt                                                                                                                          <  2:396353/141‑1 (MQ=255)
      tGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGt                                                                                                                          <  2:583441/141‑1 (MQ=255)
         aaCACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTaa                                                                                                                       <  2:13960/141‑1 (MQ=255)
           cacCATGTGGTAATGCACCTCTATGACCTGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTaaaa                                                                                                                     <  2:510918/141‑1 (MQ=255)
            acCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTACATCTGCTCATAATTCTAGAATTTTGCTTCCTAATATACTCTTTATTTCTTTTTTAGTCATACTTCCTTTTACATCAGTTTAAat                                                                                                                    >  2:643765/1‑141 (MQ=255)
            acCATGTGGTAATGCAACTCTATGACCAGATGCTACAATCGTATCGAATGATGGTAAATCTGCAAATAAATCTAGAATTTTGCTTTCTAATATTGCCTTTAATTCTATTTCAGTCATGCCTCCTTTTACAACAGATaaaaa                                                                                                                    <  2:674929/141‑2 (MQ=255)
             ccATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAata                                                                                                                   <  1:648230/141‑1 (MQ=255)
             ccATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAata                                                                                                                   >  2:451518/1‑141 (MQ=255)
              cATGTGGTAATGCACCTCTATGACCAGATTCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatat                                                                                                                  >  2:210696/1‑141 (MQ=255)
                     tAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTTATTCTCGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCata                                                                                                           <  1:317738/141‑1 (MQ=255)
                      aaTGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCatat                                                                                                          <  2:193111/141‑1 (MQ=255)
                            ccTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtca                                                                                                    >  1:194600/1‑141 (MQ=255)
                               ctATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtcatca                                                                                                 >  2:346922/1‑141 (MQ=255)
                                 aTGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaac                                                                                               <  2:194600/141‑1 (MQ=255)
                                  tGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaaca                                                                                              >  1:412259/1‑141 (MQ=255)
                                        gATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTa                                                                                        >  1:416975/1‑141 (MQ=255)
                                         aTGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAg                                                                                       >  2:413356/1‑141 (MQ=255)
                                            cTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTACCACTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctg                                                                                    <  1:12515/141‑1 (MQ=255)
                                                         gAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTaa                                                                       <  2:602513/141‑1 (MQ=255)
                                                         gAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTaa                                                                       <  1:524731/141‑1 (MQ=255)
                                                           atgatgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAg                                                                     >  2:141329/1‑141 (MQ=255)
                                                             gatgGTCCATCTGCTAATCCTTCGTTCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCa                                                                   <  1:361967/141‑1 (MQ=255)
                                                               tgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAAt                                                                 >  2:616585/1‑141 (MQ=255)
                                                               tgGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAAt                                                                 >  2:616616/1‑141 (MQ=255)
                                                                   ccATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCg                                                             >  2:482727/1‑141 (MQ=255)
                                                                   ccATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCg                                                             <  2:682705/141‑1 (MQ=255)
                                                                   ccATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCg                                                             >  2:630955/1‑141 (MQ=255)
                                                                    cATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGt                                                            >  2:35982/1‑141 (MQ=255)
                                                                     aTCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTc                                                           <  1:164143/141‑1 (MQ=255)
                                                                          cTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCAt                                                      >  2:392095/1‑141 (MQ=255)
                                                                            ccTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCt                                                    >  2:288141/1‑141 (MQ=255)
                                                                                aTTCTAGCATTGTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACATCAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGa                                                <  2:372967/141‑1 (MQ=255)
                                                                                 ttCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGAc                                               >  1:522813/1‑141 (MQ=255)
                                                                                   cTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGt                                             >  1:192516/1‑141 (MQ=255)
                                                                                    tAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCt                                                    >  2:618232/1‑133 (MQ=255)
                                                                                    tAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCCTCAGCATCt                                                    <  1:618232/133‑1 (MQ=255)
                                                                                           ttGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAAtt                                     >  2:496721/1‑141 (MQ=255)
                                                                                             gCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCttt                                                                                >  2:485138/1‑96 (MQ=255)
                                                                                             gCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCttt                                                                                <  1:485138/96‑1 (MQ=255)
                                                                                              cTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTa                                  >  2:417878/1‑141 (MQ=255)
                                                                                               tttCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTAt                                 <  1:616616/141‑1 (MQ=255)
                                                                                               tttCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTAt                                 <  1:616585/141‑1 (MQ=255)
                                                                                                        ttGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTAt                        >  1:483632/1‑141 (MQ=255)
                                                                                                        ttGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTAt                        <  2:483632/141‑1 (MQ=255)
                                                                                                        ttGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTAt                        <  2:109356/141‑1 (MQ=255)
                                                                                                        ttGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTAt                        >  1:453264/1‑141 (MQ=255)
                                                                                                        ttGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTAGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTAt                        >  2:559875/1‑141 (MQ=255)
                                                                                                          gCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTTAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTa                      >  2:577912/1‑141 (MQ=255)
                                                                                                           ccTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAg                     <  1:259231/141‑1 (MQ=255)
                                                                                                                aaTTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATg                <  2:102215/141‑1 (MQ=255)
                                                                                                                    cTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCtta            >  1:134956/1‑141 (MQ=255)
                                                                                                                    cTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCtta            >  1:271465/1‑141 (MQ=255)
                                                                                                                       tttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCttatta         >  2:252279/1‑141 (MQ=255)
                                                                                                                          cAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATg      >  2:412738/1‑141 (MQ=255)
                                                                                                                           aGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGa     >  2:210201/1‑141 (MQ=255)
                                                                                                                           aGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGa     >  1:184053/1‑141 (MQ=255)
                                                                                                                            gTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGat    >  1:285192/1‑141 (MQ=255)
                                                                                                                              cATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTATCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGatat  <  2:663940/141‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
ATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATAT  >  CP000730/1646056‑1646322

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