Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,034,536 G→A pseudogene (486/663 nt) yedS → pseudogene, outer membrane protein homology, putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,034,5360GA100.0% 244.8 / NA 76pseudogene (486/663 nt)yedSpseudogene, outer membrane protein homology, putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/38);  total (38/38)

GTCGTATATGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAGGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTA  >  NC_000913/2034402‑2034663
                                                                                                                                      |                                                                                                                               
gtCGTATATGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                             <  2:726829/139‑1 (MQ=255)
     atatGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGc                                                                                                                        <  1:1615333/139‑1 (MQ=255)
          aCGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGaa                                                                                                                   <  2:1620284/139‑1 (MQ=255)
           cGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAAc                                                                                                                  >  1:781130/1‑139 (MQ=255)
                     tGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACtttgctttgc                                                                                                        >  2:1460173/1‑139 (MQ=255)
                            atatGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCa                                                                                                 <  1:1460173/139‑1 (MQ=255)
                            atatGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAg                                                                                                >  2:573372/1‑139 (MQ=255)
                            atatGCTGCCAGAATTTGGTGGCGATACCTATCCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAatat                                                                                                   >  1:993298/1‑137 (MQ=255)
                            atatGCAGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTTCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAatat                                                                                                   <  2:993298/137‑1 (MQ=255)
                                   gcCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACtttgctttgc                                                                                                        >  1:951777/1‑125 (MQ=255)
                                   gcCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACtttgctttgc                                                                                                        <  2:951777/125‑1 (MQ=255)
                                   gcCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGtaat                                                                                          <  2:781130/139‑1 (MQ=255)
                                      aGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTaataat                                                                                       >  2:901590/1‑139 (MQ=255)
                                                    aTACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                             <  2:271596/87‑1 (MQ=255)
                                                    aTACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCg                                                                                                                             >  1:271596/1‑87 (MQ=255)
                                                    aTACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAatat                                                                                                   <  2:737979/113‑1 (MQ=255)
                                                    aTACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAatat                                                                                                   >  1:737979/1‑113 (MQ=255)
                                                           tGCTGGCGCTGACAACTTCATGAACGTTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGc                                                                                                                        <  1:1060589/85‑1 (MQ=255)
                                                           tGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGc                                                                                                                        >  2:1060589/1‑85 (MQ=255)
                                                                gcgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGAt                                                                                                                           <  2:385814/77‑1 (MQ=255)
                                                                gcgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGAt                                                                                                                           >  1:385814/1‑77 (MQ=255)
                                                                 cgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGTTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGaa                                                            >  1:149893/1‑139 (MQ=255)
                                                                 cgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctt                                                                                                           >  1:346787/1‑92 (MQ=255)
                                                                 cgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctt                                                                                                           <  2:346787/92‑1 (MQ=255)
                                                                 cgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgc                                                                                                             >  2:812631/1‑90 (MQ=255)
                                                                 cgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgc                                                                                                             <  1:812631/90‑1 (MQ=255)
                                                                 cgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGaa                                                            >  2:182929/1‑139 (MQ=255)
                                                                 cgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGaa                                                            >  1:1350567/1‑139 (MQ=255)
                                                                 cgcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGACGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGaa                                                            >  1:958638/1‑139 (MQ=255)
                                                                  gcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgc                                                                                                             <  2:1268416/89‑1 (MQ=255)
                                                                  gcTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgc                                                                                                             >  1:1268416/1‑89 (MQ=255)
                                                                     gACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCa                                                                                                 <  2:1329894/98‑1 (MQ=255)
                                                                     gACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCa                                                                                                 >  1:1329894/1‑98 (MQ=255)
                                                                      aCAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTc                                                       >  1:262816/1‑139 (MQ=255)
                                                                       cAACTTCATGAACGGTCGTGCTATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCg                                                     <  1:573372/139‑1 (MQ=255)
                                                                              aTGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGa                                                                                     >  1:962729/1‑101 (MQ=255)
                                                                              aTGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGa                                                                                     <  2:962729/101‑1 (MQ=255)
                                                                                  aCGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGc                                                                                                                        <  1:737910/62‑1 (MQ=255)
                                                                                  aCGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGc                                                                                                                        >  2:737910/1‑62 (MQ=255)
                                                                                  aCGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGa                                           <  1:182929/139‑1 (MQ=255)
                                                                                  aCGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGa                                           >  1:669001/1‑139 (MQ=255)
                                                                                  aCGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGa                                           <  2:175540/139‑1 (MQ=255)
                                                                                                ggTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGtggtgg                                                                            <  2:1664606/92‑1 (MQ=255)
                                                                                                ggTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGtggtgg                                                                            >  1:1664606/1‑92 (MQ=255)
                                                                                                ggTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAgtggtg                                            <  1:1515583/124‑1 (MQ=255)
                                                                                                ggTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAgtggtg                                            >  2:1515583/1‑124 (MQ=255)
                                                                                                   gtTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACtttgctttgc                                                                                                        >  2:1158352/1‑61 (MQ=255)
                                                                                                   gtTGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACtttgctttgc                                                                                                        <  1:1158352/61‑1 (MQ=255)
                                                                                                    ttGCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCaa                         <  2:958638/139‑1 (MQ=255)
                                                                                                      gCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctttg                                                                                                         >  1:1664762/1‑57 (MQ=255)
                                                                                                      gCAACCTACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctttg                                                                                                         <  2:1664762/57‑1 (MQ=255)
                                                                                                            tACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGa                                                                           >  2:413113/1‑81 (MQ=255)
                                                                                                            tACCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGa                                                                           <  1:413113/81‑1 (MQ=255)
                                                                                                             aCCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAAct                                                                                                                 >  2:1378551/1‑42 (MQ=255)
                                                                                                             aCCGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAAct                                                                                                                 <  1:1378551/42‑1 (MQ=255)
                                                                                                              ccGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctttg                                                                                                         >  2:949608/1‑49 (MQ=255)
                                                                                                              ccGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctttg                                                                                                         <  1:949608/49‑1 (MQ=255)
                                                                                                               cGTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTg              >  2:1012409/1‑139 (MQ=255)
                                                                                                                gTAACAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGa             <  1:901590/139‑1 (MQ=255)
                                                                                                                    cAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGa                                           <  2:1365928/105‑1 (MQ=255)
                                                                                                                    cAATGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGa                                           >  1:1365928/1‑105 (MQ=255)
                                                                                                                      aTGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGttt       >  1:698768/1‑139 (MQ=255)
                                                                                                                       tGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctttg                                                                                                         <  2:1064648/40‑1 (MQ=255)
                                                                                                                       tGGCTTCTTTGGTTAAGTCGATGGCTTGAActttgctttg                                                                                                         >  1:1064648/1‑40 (MQ=255)
                                                                                                                       tGGCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTc      >  1:136940/1‑139 (MQ=255)
                                                                                                                        ggCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCAt                                                                       <  2:410699/73‑1 (MQ=255)
                                                                                                                        ggCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCAt                                                                       >  1:410699/1‑73 (MQ=255)
                                                                                                                         gCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCa                                                                                                       <  2:708924/40‑1 (MQ=255)
                                                                                                                         gCTTCTTTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCa                                                                                                       >  1:708924/1‑40 (MQ=255)
                                                                                                                          cttcttTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGt   <  2:698768/139‑1 (MQ=255)
                                                                                                                           ttcttTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTa  <  1:1012409/139‑1 (MQ=255)
                                                                                                                           ttcttTGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTa  <  2:149893/139‑1 (MQ=255)
                                                                                                                           ttcttTGGTTAAGATGATGGCTTGAACTTTGCTTTGCAATACCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTGCCAAAGAAAATGGTGATGGTTTCGGTa  <  1:458919/139‑1 (MQ=255)
                                                                                                                              tttGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTc                                                       >  1:1504907/1‑83 (MQ=255)
                                                                                                                              tttGGTTAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTc                                                       <  2:1504907/83‑1 (MQ=255)
                                                                                                                                   ttAAGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTa  >  1:1529229/1‑131 (MQ=255)
                                                                                                                                      |                                                                                                                               
GTCGTATATGACGTAGAAGCATGGACCGATATGCTGCCAGAATTTGGTGGCGATACCTATGCTGGCGCTGACAACTTCATGAACGGTCGTGCTAATGGTGTTGCAACCTACCGTAACAATGGCTTCTTTGGTTAGGTCGATGGCTTGAACTTTGCTTTGCAATATCAGGGTAATAATGAAAGTGGTGGACCATTCGGGCAGGAAGGTTCCGGTAGTGGTGATGGCCGTAGCCTTTCCAAAGAAAATGGTGATGGTTTCGGTA  >  NC_000913/2034402‑2034663

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