Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 293,940 C→T 25.0% intergenic (+57/‑471) USA300HOU_0261 → / → ispD2 zinc (Zn2+)‑binding dehydrogenase/2‑C‑methyl‑D‑erythritol 4‑phosphate cytidylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730293,9400CT25.0% 96.4 / 32.4 60intergenic (+57/‑471)USA300HOU_0261/ispD2zinc (Zn2+)‑binding dehydrogenase/2‑C‑methyl‑D‑erythritol 4‑phosphate cytidylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (20/25);  new base T (10/5);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

CCGTTAGAAAAAGGCCCGGAGACATTTGATAAATTAGTTAACAAGAAAGAACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCATTTTGTATACTTAAAATGAGATGGCAATATT  >  CP000730/293800‑294074
                                                                                                                                             |                                                                                                                                      
ccGTTAGAAAAAGGCCCGGAGACATTTGATAAATTAGTTAACAAGAAAGAACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGc                                                                                                                                        <  1:139310/141‑1 (MQ=255)
 cgtgcaaaaaaGGCTCGGAGACATTTGATAAATGAGTTCACAAGAATGAGACGTTTGTTTAAGTCATGTTTACGATTTATTAGGATTGTCCTTTTAGGACGAAAACGTTGGTATAGTTATAGCTATGAAAGTGCGAATGcc                                                                                                                                       <  2:161575/135‑1 (MQ=255)
        aaaaGGCCCGGAGACATTTGATAAATTAGTTAACAAGAAAGAACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGt                                                                                                                                >  1:623513/1‑141 (MQ=255)
           aGGCCCGGAGACATTTGATAAATTAGTTAACAAGAAAGAACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCta                                                                                                                             >  1:138942/1‑141 (MQ=255)
           aGGCCCGGAGACATTTGATAAATTAGTTAACAAGAAAGAACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGGCTGGTCta                                                                                                                             >  1:159618/1‑141 (MQ=255)
                  gagaTATTTGTTAAATTAGTTAACAAGAAAGAACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGATAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATAc                                                                                                                      <  2:576725/141‑1 (MQ=255)
                   agaCATTTGATAAATTAGTTAACAAGAAAGAACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACt                                                                                                                     >  1:584100/1‑141 (MQ=255)
                    gaCATTTGATAAATTAGTTAACAAGAAAGAACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTa                                                                                                                    <  1:720522/141‑1 (MQ=255)
                    gaCATTTGATAAATTAGTTAACAAGAAAGAACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTa                                                                                                                    <  1:119902/141‑1 (MQ=255)
                     aCATTTGATAAATTAGTTAACAAGAAAGAACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGATATGAAAGTGTGAATGCCGTCTGGACTATAGAtattat                                                                                                                   >  1:560773/1‑141 (MQ=255)
                              aaaTTAGTTAACAAGAAAGAACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGCGAACCCAAACGCTGGTATATGCCTTGCTAGGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAaataat                                                                                                          <  1:764461/141‑1 (MQ=255)
                                 ttAGTTAACAAGAAAGAACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTCTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCa                                                                                                       >  1:595410/1‑141 (MQ=255)
                                    attAACAAGAAAGAACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCt                                                                                                    <  2:582003/140‑1 (MQ=255)
                                           agaaagaaCGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGGCGAAAAGGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATATTATCGTAATAATTCATCTTCGAatt                                                                                             >  1:572968/1‑140 (MQ=255)
                                               agaaCGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACg                                                                                         >  1:333971/1‑141 (MQ=255)
                                                   aCGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTg                                                                                      <  1:427791/141‑1 (MQ=255)
                                                    cGATTTGATAAAGTCATGTTCACGACTTATTAGTATGCACCTTTGAGGACGAAAACGCTGCTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCCAATCTCCGTTGa                                                                                     >  2:164206/1‑141 (MQ=255)
                                                     aaTTTGATAAAGTCATGTTTACGATTTATTAGTACGCACCTTTGAGGACGAAAACGCTGGTATAGTTATATCTATGAAAGTGCGAATGCCGTCTGGCCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGAt                                                                                    <  1:76427/140‑1 (MQ=255)
                                                         tGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAAt                                                                                <  2:311918/141‑1 (MQ=255)
                                                          gATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATa                                                                               >  2:384378/1‑141 (MQ=255)
                                                                gTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGt                                                                         >  1:52089/1‑141 (MQ=255)
                                                                gTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGt                                                                         >  1:35068/1‑141 (MQ=255)
                                                                     gTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACt                                                                    <  2:244839/141‑1 (MQ=255)
                                                                      tttACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACtt                                                                   >  1:64480/1‑141 (MQ=255)
                                                                                      aTGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTtg                                                   >  2:90258/1‑141 (MQ=255)
                                                                                        gCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTtgtg                                                 >  2:30449/1‑141 (MQ=255)
                                                                                         cACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGa                                                >  2:706131/1‑141 (MQ=255)
                                                                                             tttGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCg                                            <  1:5696/141‑1 (MQ=255)
                                                                                             tttGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCg                                            <  1:195360/141‑1 (MQ=255)
                                                                                                 aGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGgttg                                        >  2:469200/1‑141 (MQ=255)
                                                                                                 aGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTACCACTAAAGTAGACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGgttg                                        >  1:177723/1‑141 (MQ=255)
                                                                                                 aGGACCAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGgttg                                        <  2:133447/141‑1 (MQ=255)
                                                                                                  ggACGAAAACGCTGGTATAGTTATAGCTATGAGAGAGCGAAAACCTTCTGTTCTAGAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGgttgt                                       <  2:159618/141‑1 (MQ=255)
                                                                                                  ggACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGgttgt                                       <  2:138942/141‑1 (MQ=255)
                                                                                                  ggACGAAAACGCTGGTATAGGTATAGCTATGAAAGTGCGAATGCCGTCTGGGCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAATTATGTTTATGAGTCGgtttt                                       >  2:656688/1‑139 (MQ=255)
                                                                                                   gACGAAAACGCTGGTACAGTTATAGCTATGAAAGTGCGAATGCCGTCTCGTCTATAGATACTATCGAAATAATTCATCTTCGAATATATGTTTATAAATCGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGgttgtt                                      <  1:566882/141‑1 (MQ=255)
                                                                                                    aCGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGttt                                     <  2:273799/141‑1 (MQ=255)
                                                                                                    aCGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGttt                                     <  2:622567/141‑1 (MQ=255)
                                                                                                    aCGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGttt                                     <  2:22832/141‑1 (MQ=255)
                                                                                                    aCGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCTTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGttt                                     <  1:210610/141‑1 (MQ=255)
                                                                                                       aaaaCGCTTGTATAGTTATAGCGATGAAAGTGCGAATGCCGTCTGGTCTAGAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGc                                  <  1:608508/141‑1 (MQ=255)
                                                                                                         aaCGCTGGTGTAGTTATAGCTATGAAAGTGCGAGTGCCGTCTGGTCTATAGATGCTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCAt                                <  2:232011/141‑1 (MQ=255)
                                                                                                          aCGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCAtt                               >  1:64683/1‑141 (MQ=255)
                                                                                                           cGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCAttt                              <  1:638939/141‑1 (MQ=255)
                                                                                                           cGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTAACTTGTGTGAAATATGCTGGTGAATCGGTTGTTTTGCAttt                              >  2:122759/1‑141 (MQ=255)
                                                                                                            gCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCAtttt                             >  2:328914/1‑141 (MQ=255)
                                                                                                             cTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCCTTTTg                            >  1:519152/1‑141 (MQ=255)
                                                                                                              tGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCATTTTGt                           >  1:62589/1‑141 (MQ=255)
                                                                                                              tGGTATAGTTATAGCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGACGGTTTACTTGTGTGAAATATGATTGTGAATCGGTTGTTTTGGATTTTGt                           >  2:711840/1‑141 (MQ=255)
                                                                                                               ggTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCATTTTGta                          <  2:627120/141‑1 (MQ=255)
                                                                                                                 tataGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCATTTTGtata                        <  2:333971/141‑1 (MQ=255)
                                                                                                                  ataGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCATTTTGTATAc                       <  1:30449/141‑1 (MQ=255)
                                                                                                                         taGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCATTTTGTATACTTAAAAt                <  1:770144/141‑1 (MQ=255)
                                                                                                                           gCTATGAAAGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCATTTTGTATACTTAAAATga              >  1:772483/1‑141 (MQ=255)
                                                                                                                           gCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCATTTTGTATACTTAAAATga              <  1:378517/141‑1 (MQ=255)
                                                                                                                           gCTATGAAAGTGCAAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCACCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCATTTTGTATACTTAAATTga              <  1:556122/141‑1 (MQ=255)
                                                                                                                             tATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCATTTTGTATACTTAAAATgaga            >  2:397373/1‑141 (MQ=255)
                                                                                                                               tGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGCGAAATATGCTTGTGAATCGGTTGTTTTGCATTTTGTATACTTAAAATGAGATg          <  2:473442/141‑1 (MQ=255)
                                                                                                                                   aGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCATTTTGTATACTTAGAATGAGATGGCaa      >  2:556833/1‑141 (MQ=255)
                                                                                                                                   aGTGCGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCATTTTGTATACTTAAAATGAGATGGCaa      >  2:719698/1‑141 (MQ=255)
                                                                                                                                       cGAATGTCGTCTGGTCTATAGATACTATCGAAATAATTCATCATCGAATATACGCTGATAAATAGCCGGTTGACTTGTGTGAAATATGCTTGTGCATCGGCTGTTTTTCATTTTGTATACTTAAAATGAGATCGCAATCtt  >  1:739389/1‑141 (MQ=255)
                                                                                                                                       cGAACGCCGTCTGGTCGATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCATTTTGTATACTTAAAATGAGATGGCAATAtt  <  1:164206/141‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                      
CCGTTAGAAAAAGGCCCGGAGACATTTGATAAATTAGTTAACAAGAAAGAACGATTTGATAAAGTCATGTTTACGATTTATTAGTATGCACCTTTGAGGACGAAAACGCTGGTATAGTTATAGCTATGAAAGTGCGAATGCCGTCTGGTCTATAGATACTATCGAAATAATTCATCTTCGAATATACGTTGATAAATAGCCGGTTTACTTGTGTGAAATATGCTTGTGAATCGGTTGTTTTGCATTTTGTATACTTAAAATGAGATGGCAATATT  >  CP000730/293800‑294074

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