Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,428,685 A→C T14T (ACT→ACG)  insH1 ← IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,428,6850AC98.2% 201.8 / ‑4.1 57T14T (ACT→ACG) insH1IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (30/26);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01

GAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGAAAACGCACGACGTGCTGACCGGA  >  NC_000913/1428537‑1428825
                                                                                                                                                    |                                                                                                                                            
gAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGt                                                                                                                                            <  2:606787/151‑1 (MQ=255)
 aaTACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTg                                                                                                                                           >  2:128451/1‑151 (MQ=255)
         gcaTGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTc                                                                                                                                   >  2:710253/1‑151 (MQ=255)
              gTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTAt                                                                                                                              <  1:102704/151‑1 (MQ=255)
                    aGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGa                                                                                                                        <  1:333868/151‑1 (MQ=255)
                       ggATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAgg                                                                                                                     >  1:662426/1‑151 (MQ=255)
                         aTAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTg                                                                                                                   <  1:106010/151‑1 (MQ=255)
                                      ccATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTc                                                                                                       >  2:669543/1‑150 (MQ=255)
                                       cATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCat                                                                                                     >  2:848573/1‑151 (MQ=255)
                                        aTTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCatg                                                                                                    >  2:9973/1‑151 (MQ=255)
                                          tGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCatga                                                                                                   <  2:907631/150‑1 (MQ=255)
                                             ccGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTc                                                                                               >  1:693914/1‑151 (MQ=255)
                                                 cTTGGGATAAAACGGCTCGATGACAGCGGTCATATGCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCctc                                                                                           >  1:663720/1‑151 (MQ=255)
                                                       aTAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGAt                                                                                     >  1:182574/1‑151 (MQ=255)
                                                         aaaaCGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATgc                                                                                   >  2:557696/1‑151 (MQ=255)
                                                         aaaaCGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATgc                                                                                   <  2:630497/151‑1 (MQ=255)
                                                         aaaaCGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATgc                                                                                   <  1:345562/151‑1 (MQ=255)
                                                             cGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACAGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTc                                                                               <  2:542154/151‑1 (MQ=255)
                                                               gCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCg                                                                             >  1:358660/1‑151 (MQ=255)
                                                               gCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCg                                                                             >  2:425656/1‑151 (MQ=255)
                                                                cTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCg                                                                             >  1:730945/1‑150 (MQ=255)
                                                                 tCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCgg                                                                            >  1:517939/1‑150 (MQ=255)
                                                                      gACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGaa                                                                       <  2:894643/150‑1 (MQ=255)
                                                                      gACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAat                                                                      >  2:862458/1‑151 (MQ=255)
                                                                         aGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTTACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATatg                                                                   <  1:670069/151‑1 (MQ=255)
                                                                          gCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATatga                                                                  <  1:139526/151‑1 (MQ=255)
                                                                              tCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATgatgat                                                              >  1:569158/1‑151 (MQ=255)
                                                                              tCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATgatgat                                                              >  1:159018/1‑151 (MQ=255)
                                                                               cATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATgatgat                                                              <  2:588045/150‑1 (MQ=255)
                                                                                atatTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGAtct                                                            <  2:168142/151‑1 (MQ=255)
                                                                                atatTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGAtct                                                            >  1:112224/1‑151 (MQ=255)
                                                                                atatTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGAtct                                                            >  2:507844/1‑151 (MQ=255)
                                                                                       gCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCatat                                                       <  2:897524/149‑1 (MQ=255)
                                                                                       gCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCata                                                        <  1:579710/148‑1 (MQ=255)
                                                                                         cATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAgg                                                   <  2:440858/151‑1 (MQ=255)
                                                                                            ggCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGaa                                                 <  1:757955/150‑1 (MQ=255)
                                                                                              cAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACtt                                              >  2:522343/1‑151 (MQ=255)
                                                                                                  aTCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTc                                          >  2:555832/1‑151 (MQ=255)
                                                                                                  aTCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTc                                          >  1:462462/1‑151 (MQ=255)
                                                                                                    cTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGc                                        >  1:560051/1‑151 (MQ=255)
                                                                                                      gCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTc                                          <  2:728673/147‑1 (MQ=255)
                                                                                                      gCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTc                                          <  2:434530/147‑1 (MQ=255)
                                                                                                          cATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTc                                  >  2:69666/1‑151 (MQ=255)
                                                                                                           aTGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTcc                                 <  2:652636/151‑1 (MQ=255)
                                                                                                            tGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTccc                                >  2:597533/1‑151 (MQ=255)
                                                                                                               gggAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAg                             <  2:427650/151‑1 (MQ=255)
                                                                                                                 gagaGGAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAggg                           >  1:807990/1‑151 (MQ=255)
                                                                                                                     ggAAAATCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGaaa                       <  2:449967/151‑1 (MQ=255)
                                                                                                                       aaaaTCTCTTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGAAAAc                     <  1:592190/151‑1 (MQ=255)
                                                                                                                              ctTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGAAAACGCac                 <  2:26025/148‑1 (MQ=255)
                                                                                                                              ctTTTCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTATTCGCACCTTCCCAAGGGGAAAACGCacgac              >  2:333770/1‑151 (MQ=255)
                                                                                                                                tttCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGAAAACGCacgac              <  1:288100/149‑1 (MQ=255)
                                                                                                                                 ttCGGGTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGAAAACGCACGACGTg           >  1:499952/1‑151 (MQ=255)
                                                                                                                                      gTCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGAAAACGCACGACGTGCTGAc      <  2:877464/151‑1 (MQ=255)
                                                                                                                                       tCTGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGAAAACGCACGACGTGCTGAcc     >  1:112188/1‑151 (MQ=255)
                                                                                                                                         tGACGGCGCTTCGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGAAAACGCACGACGTGCTGACCgg   >  1:351441/1‑151 (MQ=255)
                                                                                                                                          gACGGCGCTTAGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGAAAACGCACGACGTGCTGACCGGa  <  1:444049/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GAATACGCAGCATGGTCTCCAGCGGATAGGGCCGTCGGCCATTGCCCGCCTTGGGATAAAACGGCTCGATGACAGCGGTCATATTCTGCCATGGCAGAATCTGCTCCATGCGGGAGAGGAAAATCTCTTTTCGGGTCTGACGGCGCTTAGTGCTGAATTCACTATCGGCGAAGGTGAGTTGATGGCTCATGATGTCCCTCTGGGATGCGCTCCGGATGAATATGATGATCTCATATCAGGAACTTGTTCGCACCTTCCCAAGGGGAAAACGCACGACGTGCTGACCGGA  >  NC_000913/1428537‑1428825

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