Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,419,344 G→A A38V (GCT→GTT)  ydaG ← Rac prophage, uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,419,3440GA100.0% 236.3 / NA 72A38V (GCT→GTT) ydaGRac prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/34);  total (38/34)

TTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGC  >  NC_000913/1419220‑1419476
                                                                                                                            |                                                                                                                                    
ttGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctgcctgc                                                                                                                        >  1:1099147/1‑139 (MQ=255)
  ggTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCTCATCTTGTTGCTACGTATAACGACTTCTGCCTGCCa                                                                                                                      >  1:1482050/1‑139 (MQ=255)
    tataCGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGATGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCaga                                                                                                                    <  2:1099147/139‑1 (MQ=255)
            aTGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAg                                                                                                            >  2:1182717/1‑139 (MQ=255)
             tGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctgcctg                                                                                                                         <  2:606991/125‑1 (MQ=255)
             tGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctgcctg                                                                                                                         >  1:606991/1‑125 (MQ=255)
              ggCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGTGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTc                                                                                                          >  1:1325564/1‑139 (MQ=255)
              ggCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTc                                                                                                          >  2:540442/1‑139 (MQ=255)
              ggCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTc                                                                                                          >  2:633982/1‑139 (MQ=255)
                   tGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                     <  1:677618/139‑1 (MQ=255)
                   tGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                     <  2:1436325/139‑1 (MQ=255)
                        gTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCTCATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                <  2:1482050/139‑1 (MQ=255)
                        gTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                <  1:981211/139‑1 (MQ=255)
                        gTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                <  1:1182717/139‑1 (MQ=255)
                               ccTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCt                                                                                         >  2:685522/1‑139 (MQ=255)
                                   gTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                     >  2:59714/1‑123 (MQ=255)
                                   gTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                     <  1:59714/123‑1 (MQ=255)
                                   gTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                <  2:148910/128‑1 (MQ=255)
                                   gTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                >  1:148910/1‑128 (MQ=255)
                                      ggggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                  <  1:633982/139‑1 (MQ=255)
                                        ggAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCt                                                                                >  2:1312199/1‑139 (MQ=255)
                                             acGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctgcct                                                                                                                          <  2:714695/92‑1 (MQ=255)
                                             acGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctgcct                                                                                                                          >  1:714695/1‑92 (MQ=255)
                                             acGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCaga                                                                                                                    <  1:26907/98‑1 (MQ=255)
                                             acGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCaga                                                                                                                    >  2:26907/1‑98 (MQ=255)
                                               gccgcACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                >  1:330228/1‑116 (MQ=255)
                                               gccgcACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                <  2:330228/116‑1 (MQ=255)
                                                 cgcACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                  <  1:340240/128‑1 (MQ=255)
                                                 cgcACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                  >  2:340240/1‑128 (MQ=255)
                                                 cgcACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAg                                                                       >  2:1237071/1‑139 (MQ=255)
                                                          gggATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGc                                                              <  1:870122/139‑1 (MQ=255)
                                                            gATTCGTTGTTGCCGAGCTCGATTTTTTGCAATTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGACAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGcaa                                                            <  1:801676/139‑1 (MQ=255)
                                                                cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                <  1:1144042/99‑1 (MQ=255)
                                                                cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                >  2:1144042/1‑99 (MQ=255)
                                                                cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGca                                                                                      <  1:662854/109‑1 (MQ=255)
                                                                cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGca                                                                                      >  2:662854/1‑109 (MQ=255)
                                                                cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                  <  2:349431/113‑1 (MQ=255)
                                                                cGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                  >  1:349431/1‑113 (MQ=255)
                                                                  ttgttgCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCg                                                                                                         <  1:844522/88‑1 (MQ=255)
                                                                  ttgttgCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCg                                                                                                         >  2:844522/1‑88 (MQ=255)
                                                                     ttgCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCt                                                                                >  1:1341757/1‑110 (MQ=255)
                                                                     ttgCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCt                                                                                <  2:1341757/110‑1 (MQ=255)
                                                                     ttgCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGcaaca                                                          >  2:769951/1‑132 (MQ=255)
                                                                     ttgCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGcaaca                                                          <  1:769951/132‑1 (MQ=255)
                                                                      tgCCGAGGTCGATTTTTTGTATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                <  2:1165626/93‑1 (MQ=255)
                                                                      tgCCGAGGTCGATTTTTTGTATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGt                                                                                                >  1:1165626/1‑93 (MQ=255)
                                                                      tgCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                  >  1:337828/1‑107 (MQ=255)
                                                                      tgCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGc                                                                                  <  2:337828/107‑1 (MQ=255)
                                                                            ggTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACt                                            >  1:207198/1‑139 (MQ=255)
                                                                             gTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctg                                                                                                                             >  2:1514483/1‑57 (MQ=255)
                                                                             gTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAACGACTTctg                                                                                                                             <  1:1514483/57‑1 (MQ=255)
                                                                                            tttGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATg                            >  2:1137863/1‑139 (MQ=255)
                                                                                              tGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCa                          <  2:935056/139‑1 (MQ=255)
                                                                                              tGCGAATCTCACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCa                          <  1:1312199/139‑1 (MQ=255)
                                                                                                    tctcACATCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCt                    <  1:1137863/139‑1 (MQ=255)
                                                                                                        acaTCTTGTTGCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTAtc                >  1:1192809/1‑139 (MQ=255)
                                                                                                             ttgttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTcc           >  2:1357932/1‑139 (MQ=255)
                                                                                                               gttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                     >  2:149143/1‑47 (MQ=255)
                                                                                                               gttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                     <  1:149143/47‑1 (MQ=255)
                                                                                                               gttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTcccc         >  2:825578/1‑139 (MQ=255)
                                                                                                                ttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                     <  1:658907/46‑1 (MQ=255)
                                                                                                                ttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTtgc                                                                                                     >  2:658907/1‑46 (MQ=255)
                                                                                                                ttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGca                                                                                      <  2:847368/61‑1 (MQ=255)
                                                                                                                ttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGca                                                                                      >  1:847368/1‑61 (MQ=255)
                                                                                                                ttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCt        <  2:1492092/139‑1 (MQ=255)
                                                                                                                ttgCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCt        >  2:1110778/1‑139 (MQ=255)
                                                                                                                  gCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCg                                                                                                         <  1:1060624/40‑1 (MQ=38)
                                                                                                                  gCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCg                                                                                                         >  2:1060624/1‑40 (MQ=38)
                                                                                                                  gCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCAc                         >  1:247492/1‑120 (MQ=255)
                                                                                                                  gCTACGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCAAAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCAc                         <  2:247492/120‑1 (MQ=255)
                                                                                                                      cGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGc  >  1:1380492/1‑139 (MQ=255)
                                                                                                                      cGTATAACGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGc  >  1:4357/1‑139 (MQ=255)
                                                                                                                            |                                                                                                                                    
TTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACACCAGGGATTCGTTGTTGCCGAGGTCGATTTTTTGCATTTTGCGAATCTCACATCTTGTTGCTACGTATAGCGACTTCTGCCTGCCAGAGATCCCAGTCGTTGCTGCGTAAAGCCTGCACAGCCTGGTTGTAAGTGATACCGCAACAATCCATCAAATACTGAACTACTTCGTAATGCACCATCTTATCTCTCCCCTTAACGC  >  NC_000913/1419220‑1419476

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: ‑