Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,412,000 C→T *412* (TAG→TAA intR ← Rac prophage, integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,412,0000CT100.0% 245.0 / NA 75*412* (TAG→TAAintRRac prophage, integrase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/40);  total (35/40)

CCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTC  >  NC_000913/1411869‑1412130
                                                                                                                                   |                                                                                                                                  
cccACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATa                                                                                                                             <  2:200935/139‑1 (MQ=255)
    cGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                         >  1:225747/1‑139 (MQ=255)
     gTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcgc                                                                                                                        <  2:1255041/139‑1 (MQ=255)
            cgCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCa                                                                                                                 <  2:939716/139‑1 (MQ=255)
              cAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                               <  2:307257/139‑1 (MQ=255)
                                     ggTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCAt                                                                                        <  1:1515915/139‑1 (MQ=255)
                                         gtagtttacttcttttGTTATGTGTGGATTATATTGAATTATTTCAGTACTTTGCTACTAAATGGGCCAAATCGGGCGCAAACTTTGCAATAACGATAACCGCGCATTCAACATCGATATCTGTTCGTCGTTCTTGtca                                                                                    <  1:1490138/126‑1 (MQ=255)
                                         gTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                                    <  2:1581156/139‑1 (MQ=255)
                                               ttctttctTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCaa                                                                                                                <  2:1188275/105‑1 (MQ=255)
                                               ttctttctTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCaa                                                                                                                >  1:1188275/1‑105 (MQ=255)
                                                tctttcttTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                                     <  1:287417/99‑1 (MQ=255)
                                                tctttcttTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                                     >  2:287417/1‑99 (MQ=255)
                                                tctttcttTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                               >  2:1133623/1‑105 (MQ=255)
                                                tctttcttTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                               <  1:1133623/105‑1 (MQ=255)
                                                tctttcttTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATcc                                                                             >  2:7156/1‑139 (MQ=255)
                                                  tttctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCcac                                                                           <  1:649153/139‑1 (MQ=255)
                                                    tctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATa                                                                                                                             <  2:4188/87‑1 (MQ=255)
                                                    tctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATa                                                                                                                             >  1:4188/1‑87 (MQ=255)
                                                    tctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATa                                                                                                                             >  1:1534802/1‑87 (MQ=255)
                                                    tctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATa                                                                                                                             <  2:1534802/87‑1 (MQ=255)
                                                    tctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                                    <  1:1261135/128‑1 (MQ=255)
                                                    tctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                                    >  2:1261135/1‑128 (MQ=255)
                                                    tctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACAt                                                                         <  2:1235608/139‑1 (MQ=255)
                                                         gTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTGCGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGt                                                                    >  1:1151806/1‑139 (MQ=255)
                                                         gTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGtt                                                                                                  >  2:174575/1‑109 (MQ=255)
                                                         gTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGtt                                                                                                  <  1:174575/109‑1 (MQ=255)
                                                         gTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGt                                                                    >  1:453248/1‑139 (MQ=255)
                                                           aaTTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCTTTCAACATGTCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                                  <  1:991085/139‑1 (MQ=255)
                                                            aTTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaaa                                                                 >  1:361086/1‑139 (MQ=255)
                                                               tgttgATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAAttt                                                              <  2:1571740/139‑1 (MQ=255)
                                                                    aTTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATac                                                         <  1:7156/139‑1 (MQ=255)
                                                                    aTTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATac                                                         >  2:571019/1‑139 (MQ=255)
                                                                         cTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCaa                                                    <  1:707060/139‑2 (MQ=255)
                                                                                     cAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTc                                        <  2:453248/139‑1 (MQ=255)
                                                                                        ttCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCAt                                                                                        <  1:323428/88‑1 (MQ=255)
                                                                                        ttCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCAt                                                                                        >  2:323428/1‑88 (MQ=255)
                                                                                        ttCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATg                                     >  1:119999/1‑139 (MQ=255)
                                                                                        ttCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATg                                     <  2:119999/139‑1 (MQ=255)
                                                                                        ttCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATg                                     >  2:692065/1‑139 (MQ=255)
                                                                                          ctctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGcc                                   >  1:1238870/1‑139 (MQ=255)
                                                                                            ctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATg                                     <  2:236252/135‑1 (MQ=255)
                                                                                            ctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATg                                     >  1:236252/1‑135 (MQ=255)
                                                                                            ctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGcccc                                 >  2:1228323/1‑139 (MQ=255)
                                                                                             tGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGtt                                                                                                  >  1:80768/1‑73 (MQ=255)
                                                                                             tGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGtt                                                                                                  <  2:80768/73‑1 (MQ=255)
                                                                                             tGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCaa                                                                                <  1:649272/91‑1 (MQ=255)
                                                                                             tGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCaa                                                                                >  2:649272/1‑91 (MQ=255)
                                                                                                  cTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                           <  2:225747/139‑1 (MQ=255)
                                                                                                     aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtggctg                        >  2:1114253/1‑139 (MQ=255)
                                                                                                     aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtggctg                        >  2:1337190/1‑139 (MQ=255)
                                                                                                     aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                           <  2:1451322/136‑1 (MQ=255)
                                                                                                     aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                           >  1:1451322/1‑136 (MQ=255)
                                                                                                         gggCAAATTGGGGGCAAACTTTGCAATTACGATAAcc                                                                                                                          >  1:426286/1‑37 (MQ=255)
                                                                                                         gggCAAATTGGGGGCAAACTTTGCAATTACGATAAcc                                                                                                                          <  2:426286/37‑1 (MQ=255)
                                                                                                         gggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCta                    >  1:243386/1‑139 (MQ=255)
                                                                                                          ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                                    <  2:74884/74‑1 (MQ=255)
                                                                                                          ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                                    >  1:74884/1‑74 (MQ=255)
                                                                                                          ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                                    <  2:499467/74‑1 (MQ=255)
                                                                                                          ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                                    >  1:499467/1‑74 (MQ=255)
                                                                                                          ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCtat                   >  1:1321195/1‑139 (MQ=255)
                                                                                                              aaTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAAt               <  1:571019/139‑1 (MQ=255)
                                                                                                               aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCAt                                                                                        >  2:826905/1‑65 (MQ=255)
                                                                                                               aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCAt                                                                                        <  1:826905/65‑1 (MQ=255)
                                                                                                               aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                                  >  2:98617/1‑87 (MQ=255)
                                                                                                               aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                                  <  1:98617/87‑1 (MQ=255)
                                                                                                                    gggCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATg                                     >  2:1330166/1‑111 (MQ=255)
                                                                                                                    gggCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATg                                     <  1:1330166/111‑1 (MQ=255)
                                                                                                                     ggCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAAt                                                                <  1:137927/83‑1 (MQ=255)
                                                                                                                     ggCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAAt                                                                <  1:137926/83‑1 (MQ=255)
                                                                                                                     ggCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAAt                                                                >  2:137926/1‑83 (MQ=255)
                                                                                                                     ggCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAAt                                                                >  2:137927/1‑83 (MQ=255)
                                                                                                                         aaCTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                                  <  1:149631/77‑1 (MQ=255)
                                                                                                                         aaCTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                                  >  2:149631/1‑77 (MQ=255)
                                                                                                                           cTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTc  <  2:243386/139‑1 (MQ=255)
                                                                                                                           cTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTc  <  1:1337190/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                  
CCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGGTTCGCTC  >  NC_000913/1411869‑1412130

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