Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 720,946 C→T intergenic (+106/+110) ybfK → / ← kdpE uncharacterized protein/response regulator in two‑component regulatory system with KdpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913720,9460CT100.0% 301.9 / NA 85intergenic (+106/+110)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (42/43);  total (42/43)

GGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTA  >  NC_000913/720851‑721045
                                                                                               |                                                                                                    
ggTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTa                                                                                                  >  2:1384126/1‑100 (MQ=255)
ggTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTa                                                                                                  >  2:442240/1‑100 (MQ=255)
      gATATTTTATGAAGGAAATATGATTGAAATGTGAAGGTGAGCAGTAAAAAAAGATCAATCTTAATTTATTTAATGATGAGCTTTTTACGTAGTAATATaa                                                                                            >  2:584072/1‑100 (MQ=255)
        tatTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATaaaa                                                                                          >  1:225692/1‑100 (MQ=255)
         atTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAat                                                                                         <  1:833929/100‑1 (MQ=255)
           tCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAat                                                                                         >  1:908221/1‑98 (MQ=255)
            cTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATAtt                                                                                      >  2:500484/1‑100 (MQ=255)
            cTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATAtt                                                                                      >  1:1376634/1‑100 (MQ=255)
              aTCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATAtt                                                                                      >  2:102360/1‑98 (MQ=255)
              aTCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGa                                                                                    <  2:790600/100‑1 (MQ=255)
              aTCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGa                                                                                    >  2:1048664/1‑100 (MQ=255)
                ccAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAt                                                                                  <  2:1314528/100‑1 (MQ=255)
                 cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAtt                                                                                 >  2:1000068/1‑100 (MQ=255)
                  aGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTg                                                                                >  1:910359/1‑100 (MQ=255)
                      aaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTac                                                                            <  1:104518/100‑2 (MQ=255)
                      aaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAt                                                                           >  1:698470/1‑100 (MQ=255)
                        tACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttt                                                                         >  2:360531/1‑100 (MQ=255)
                         aCGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAtttt                                                                        <  1:587733/100‑1 (MQ=255)
                         aCGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAtttt                                                                        <  2:562189/100‑1 (MQ=255)
                          cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttttt                                                                       <  1:530060/100‑1 (MQ=255)
                          cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttttt                                                                       >  1:557110/1‑100 (MQ=255)
                           ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                      >  1:1308224/1‑100 (MQ=255)
                           ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                      <  1:1219849/100‑1 (MQ=255)
                           ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                      <  2:1386669/100‑1 (MQ=255)
                           ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                      <  2:181700/100‑1 (MQ=255)
                            gTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgt                                                                     >  2:248838/1‑100 (MQ=255)
                            gTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                      <  2:430917/99‑1 (MQ=255)
                             ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                    >  2:1502633/1‑100 (MQ=255)
                              tGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtgt                                                                   <  1:391633/100‑1 (MQ=255)
                                aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtgtg                                                                  <  2:124801/99‑1 (MQ=255)
                                 aCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgtt                                                                <  2:376386/100‑1 (MQ=255)
                                  cTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgttg                                                               >  1:1522829/1‑100 (MQ=255)
                                    gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTgttgtt                                                             <  1:907754/100‑1 (MQ=255)
                                     tgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGttt                                                            <  2:738752/100‑1 (MQ=255)
                                     tgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGttt                                                            <  1:140159/100‑1 (MQ=255)
                                      gAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTa                                                           <  1:1252524/100‑1 (MQ=255)
                                       aaGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTaa                                                          <  2:1536236/100‑1 (MQ=255)
                                        aGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAg                                                         <  2:1126558/100‑1 (MQ=255)
                                         ggTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGa                                                        <  1:169610/100‑1 (MQ=255)
                                          gTGATCAATAAAAAATGATCAATCTTAATTTATTTACTGATGAGCATTTTACTTAGTAATATAAAATATTGACTTGTTATTTTTGTGGGGTGTTTAAGAt                                                       >  2:1346669/1‑100 (MQ=255)
                                           tGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGa                                                        <  2:86980/98‑1 (MQ=255)
                                           tGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATa                                                      >  1:810993/1‑100 (MQ=255)
                                           tGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATa                                                      >  1:769543/1‑100 (MQ=255)
                                            gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaa                                                     >  1:675546/1‑100 (MQ=255)
                                               cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                                  >  2:1012147/1‑100 (MQ=255)
                                               cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                                  >  1:315279/1‑100 (MQ=255)
                                                aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                                  <  2:56043/99‑1 (MQ=255)
                                                aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAg                                                 >  1:1409632/1‑100 (MQ=255)
                                                aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAg                                                 <  2:820571/100‑1 (MQ=255)
                                                aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAg                                                 >  1:1175029/1‑100 (MQ=255)
                                                aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAg                                                 >  2:568071/1‑100 (MQ=255)
                                                 aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAg                                                 >  2:592065/1‑99 (MQ=255)
                                                  tAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGcc                                               >  2:615304/1‑100 (MQ=255)
                                                   aaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCg                                              <  1:275641/100‑1 (MQ=255)
                                                    aaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGt                                             <  2:1248073/100‑1 (MQ=255)
                                                       aaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatt                                          <  2:48553/100‑1 (MQ=255)
                                                         tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattta                                        <  1:102360/100‑1 (MQ=255)
                                                          gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttat                                       >  2:1475415/1‑100 (MQ=255)
                                                          gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttat                                       >  2:448183/1‑100 (MQ=255)
                                                             cAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                    <  1:1144276/100‑1 (MQ=255)
                                                             cAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                    <  2:54465/100‑1 (MQ=255)
                                                               aTCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatta                                  <  2:193065/100‑1 (MQ=255)
                                                               aTCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatta                                  <  1:1005583/100‑1 (MQ=255)
                                                                   tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                    <  2:55703/94‑1 (MQ=255)
                                                                   tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                              <  2:822865/100‑1 (MQ=255)
                                                                     atttatttAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                            >  1:1400978/1‑100 (MQ=255)
                                                                       ttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTtta                          <  2:769539/100‑1 (MQ=255)
                                                                         atttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                        >  1:1066488/1‑100 (MQ=255)
                                                                          tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                       <  2:889832/100‑1 (MQ=255)
                                                                          tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                       <  1:1300029/100‑1 (MQ=255)
                                                                           atcttgatgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTACTTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                       <  1:395770/96‑1 (MQ=255)
                                                                            taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                     >  2:383059/1‑100 (MQ=255)
                                                                            taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                     <  2:698467/100‑1 (MQ=255)
                                                                               tgatgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTACTTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattagat                   >  2:395769/1‑97 (MQ=255)
                                                                                gatgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaa                 >  1:311832/1‑100 (MQ=255)
                                                                                 atgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATCTATTATTACGGCTTTAATTAATaaaa                <  1:956950/100‑1 (MQ=255)
                                                                                  tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag               >  2:228927/1‑100 (MQ=255)
                                                                                  tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag               >  1:336240/1‑100 (MQ=255)
                                                                                  tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag               >  2:228906/1‑100 (MQ=255)
                                                                                  tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGACTTTAATTAATaaaa                >  1:1035895/1‑99 (MQ=255)
                                                                                      ttttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAAATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag           <  2:827750/99‑1 (MQ=255)
                                                                                      cTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag           <  2:16220/100‑1 (MQ=255)
                                                                                      cTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca            >  2:344022/1‑99 (MQ=255)
                                                                                       tttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTATGGCTTTAATTAATAAAaggcag           <  2:910351/99‑1 (MQ=255)
                                                                                           tACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt      >  1:584565/1‑100 (MQ=255)
                                                                                               cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa  >  1:666574/1‑100 (MQ=255)
                                                                                               |                                                                                                    
GGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTA  >  NC_000913/720851‑721045

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