Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,258,251 T→A 8.5% intergenic (‑241/+34) yhaO ← / ← tdcG putative transporter/L‑serine dehydratase 3, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,258,2510TA 8.5% 201.6 / 8.5 59intergenic (‑241/+34)yhaO/tdcGputative transporter/L‑serine dehydratase 3, anaerobic
Reads supporting (aligned to +/- strand):  ref base T (30/24);  new base A (3/2);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

ATAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAG  >  NC_000913/3258105‑3258389
                                                                                                                                                  |                                                                                                                                          
aTAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACAttt                                                                                                                                          <  4:660656/149‑1 (MQ=255)
aTAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACAttt                                                                                                                                          <  1:572533/149‑1 (MQ=255)
   aaGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTtgtg                                                                                                                                       >  1:581197/1‑149 (MQ=255)
     gAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCa                                                                                                                                     >  4:547154/1‑149 (MQ=255)
     gAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCa                                                                                                                                     >  4:268530/1‑149 (MQ=255)
           aTACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAg                                                                                                                               >  1:603871/1‑149 (MQ=255)
             aCGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGa                                                                                                                             <  4:18880/149‑1 (MQ=255)
             aCGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGa                                                                                                                             <  2:390347/149‑1 (MQ=255)
              cgcgcgAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGAt                                                                                                                            <  3:268530/149‑1 (MQ=255)
                         agagAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACa                                                                                                                 <  2:603871/149‑1 (MQ=255)
                                   ttCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAg                                                                                                       >  2:266572/1‑149 (MQ=255)
                                      acaATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCg                                                                                                    >  2:498972/1‑149 (MQ=255)
                                        aaTCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCa                                                                                                  >  2:575373/1‑149 (MQ=255)
                                          tCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGa                                                                                                >  4:151393/1‑149 (MQ=255)
                                            aaGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGAcc                                                                                              <  2:581197/149‑1 (MQ=255)
                                             aGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCa                                                                                             >  2:44820/1‑149 (MQ=255)
                                                     aaaTAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATg                                                                                     <  3:151393/149‑1 (MQ=255)
                                                              agaTGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGtcct                                                                            >  2:681718/1‑149 (MQ=255)
                                                                aTGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGtcctcc                                                                          >  3:287638/1‑149 (MQ=255)
                                                                    aTAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGt                                                                      <  2:624101/149‑1 (MQ=255)
                                                                        gAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATg                                                                  >  2:301876/1‑149 (MQ=255)
                                                                        gAATTTTCTACGAAATTAATCACAGGAGAATATTCATTCTCTATTTGCTGCAGATCACAAAAAAAAGGTGCACATTTGTTCTCCCAAAGATGAACGCTGACAACAATGTCAGCAGAAAACCACTTTAATGGACAGTCCTCCGCGTGagg                                                                  >  2:340986/1‑147 (MQ=255)
                                                                                 actaAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGt                                                         >  4:551634/1‑149 (MQ=255)
                                                                                  ctaAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGtcctc                                                                           >  1:632150/1‑130 (MQ=255)
                                                                                  ctaAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGtcctc                                                                           <  2:632150/130‑1 (MQ=255)
                                                                                  ctaAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTa                                                        >  3:567642/1‑149 (MQ=255)
                                                                                   taAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGc                                                                                   <  1:159568/121‑1 (MQ=255)
                                                                                   taAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGc                                                                                   >  2:159568/1‑121 (MQ=255)
                                                                                   taAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTAt                                                       <  1:498972/149‑1 (MQ=255)
                                                                                   taAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTAt                                                       >  2:218014/1‑149 (MQ=255)
                                                                                     aaTTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTa                                                        <  3:243765/146‑1 (MQ=255)
                                                                                     aaTTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTa                                                        >  4:243765/1‑146 (MQ=255)
                                                                                           tCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGtt                                               <  4:567642/149‑1 (MQ=255)
                                                                                            cGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCATCCGCGTGATGTTTCGCGGTATTTATCGTTc                                              >  1:197580/1‑149 (MQ=255)
                                                                                                   gaATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCtt                                       <  1:44820/149‑1 (MQ=255)
                                                                                                      tatTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGc                                    >  1:428088/1‑149 (MQ=255)
                                                                                                        ttcattcTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGa                                                                                                                             <  4:483540/58‑1 (MQ=255)
                                                                                                        ttcattcTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGa                                                                                                                             >  3:483540/1‑58 (MQ=255)
                                                                                                               tctATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCAt                           <  3:6038/149‑1 (MQ=255)
                                                                                                                 tATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATAc                         >  3:356893/1‑149 (MQ=255)
                                                                                                                   tttGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACAt                       >  1:42247/1‑149 (MQ=255)
                                                                                                                    ttGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCAGTATTTATCGTTCATATCTTTGCCGGTTTCATACATc                      <  3:551634/149‑1 (MQ=255)
                                                                                                                         gCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGtctc                 <  1:575373/149‑1 (MQ=255)
                                                                                                                              tCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTc                                              >  2:354693/1‑115 (MQ=255)
                                                                                                                              tCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTc                                              <  1:354693/115‑1 (MQ=255)
                                                                                                                                acaAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACt          <  1:266572/149‑1 (MQ=255)
                                                                                                                                acaAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACt          >  4:276926/1‑149 (MQ=255)
                                                                                                                                  aaaaaaaaGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGtctc                 >  3:254665/1‑140 (MQ=255)
                                                                                                                                  aaaaaaaaGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGtctc                 <  4:254665/140‑1 (MQ=255)
                                                                                                                                  aaaaaaaaGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACttt        >  1:240548/1‑149 (MQ=255)
                                                                                                                                    aaaaaaGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTAt      <  1:681718/149‑1 (MQ=255)
                                                                                                                                     ccttgggTGCACAAATGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACttt        <  3:154083/141‑1 (MQ=255)
                                                                                                                                     ccttgggTGCACAAATGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACttt        >  4:154083/6‑146 (MQ=255)
                                                                                                                                     ccttgggTGCACAAATGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATc     >  1:466116/6‑149 (MQ=255)
                                                                                                                                     aaaaaGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATc     <  1:301876/149‑1 (MQ=255)
                                                                                                                                      aaaaGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCTGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCg    >  1:512235/1‑149 (MQ=255)
                                                                                                                                       ttgggTGCACAAATGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCa                                                                                 <  2:116914/68‑1 (MQ=38)
                                                                                                                                       ttgggTGCACAAATGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCa                                                                                 >  1:116914/4‑71 (MQ=38)
                                                                                                                                        aaGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAg  <  3:276926/149‑1 (MQ=255)
                                                                                                                                           gTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAg  <  2:66112/146‑1 (MQ=255)
                                                                                                                                           gTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAg  >  1:66112/1‑146 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
ATAAAGAAATGATACGCGCGAGGTCAGAGAAAATCTTCACAATCAAGACTTTGAAATAACGGAGATGGATAAGAATTTTCTACTAAATTAATCGCAGGAGAATATTCATTCTCTATTTGCGGCAGATCACAAAAAAAAGGTGCACATTTGTGCACCCAAGGATGAAAGCTGACAGCAATGTCAGCCGCAGACCACTTTAATGGCCAGTCCTCCGCGTGATGTTTCGCGGTATTTATCGTTCATATCTTTGCCGGTTTCATACATCGTCTCGATCACTTTATCGAG  >  NC_000913/3258105‑3258389

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