Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 649,079 T→C T575A (ACC→GCC)  mrdA ← transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606649,0790TC96.7% 671.1 / 7.3 180T575A (ACC→GCC) mrdAtranspeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  ref base T (1/5);  new base C (86/87);  total (87/93)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

Only 100 of 220 total aligned reads displayed.

CGGACCCGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGACCGAA  >  REL606/649008‑649151
                                                                       |                                                                        
cgGACCCGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGTGGCAAAGGCGGccatc                                                                      >  2:3796172/1‑76 (MQ=255)
cgGACCCGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCCAAAGGGGccatc                                                                      >  2:1048540/1‑76 (MQ=255)
cgGACCCGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAATCGGgcata                                                                      >  2:1338935/1‑75 (MQ=255)
     ccgccccgccTTCCTCCAGAATCATGGCGCCAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGttt                                                                 <  2:7323705/76‑1 (MQ=255)
     ccGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTTGATGGCGCAAAAGCGGGCCTCCGttt                                                                 >  2:123065/1‑76 (MQ=255)
     ccGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGGGCAAAGGCGGCCATCAGttt                                                                 >  2:8555591/1‑76 (MQ=255)
     ccGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCCTCAGgtt                                                                 >  2:8403570/1‑76 (MQ=255)
     ccGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGttt                                                                 >  2:9258769/1‑76 (MQ=255)
     ccGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGttt                                                                 >  2:9034024/1‑76 (MQ=255)
     ccGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGttt                                                                 >  2:7103039/1‑76 (MQ=255)
     ccGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGttt                                                                 >  2:2693992/1‑76 (MQ=255)
        ccccgccgTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTtgtg                                                              <  2:1196497/74‑1 (MQ=255)
         gccgccGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTTTgg                                                             >  2:4456475/1‑76 (MQ=255)
         gccgccGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTgg                                                             >  2:3124353/1‑76 (MQ=255)
         gccgccGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAAGCGGCCATCAGTTTGTgg                                                             >  2:2577504/1‑76 (MQ=255)
         cccgcctTTCTCCAGAATCATGGCCCCAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTgg                                                             <  2:9361586/75‑1 (MQ=255)
         cccgccgTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTgg                                                             <  1:2039459/75‑1 (MQ=255)
          ccgccgTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGGTTGTGGt                                                            >  1:8300180/1‑76 (MQ=255)
          ccgccgTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGGTTGTGGt                                                            >  1:6151546/1‑76 (MQ=255)
           cgccgTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTc                                                           >  2:9568487/1‑76 (MQ=255)
           cgccgTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTc                                                           <  1:3013504/76‑1 (MQ=255)
           cgccgTTCCCCAGAATCAGGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTc                                                           <  1:4903655/76‑1 (MQ=255)
              cgTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGGCATCAGTTTGTGGTCACg                                                        >  1:8251200/1‑76 (MQ=255)
              cgTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGGCACg                                                        >  1:4019746/1‑76 (MQ=255)
               gtactcCAGAAAAAGGGGGGCAGCTACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGt                                                       <  1:7599312/73‑1 (MQ=255)
               gTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCCGTTTGTGGTCACGt                                                       >  2:3373856/1‑76 (MQ=255)
               gTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGt                                                       <  1:1363017/76‑1 (MQ=255)
                ttCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGGCACGTa                                                      >  2:2496968/1‑76 (MQ=255)
                 tctcCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTaa                                                     <  1:492122/76‑1 (MQ=255)
                  ctcCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTaaa                                                    <  2:2348152/76‑1 (MQ=255)
                   tcCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAAc                                                   <  1:1048540/76‑1 (MQ=255)
                   tcCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAAc                                                   <  1:1338935/76‑1 (MQ=255)
                   cccAGAATCTGGGCCACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAAc                                                   <  1:9208914/75‑1 (MQ=255)
                   cccAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAAc                                                   <  1:2358531/75‑1 (MQ=255)
                   cccAGAATCATGCCGCCAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAAc                                                   <  1:3796172/75‑1 (MQ=255)
                      aGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCt                                                >  2:6974889/1‑76 (MQ=255)
                      aGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCt                                                <  1:186332/76‑1 (MQ=255)
                      aGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGGAAACGCt                                                >  2:8306940/1‑76 (MQ=255)
                        aaTCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCg                                              <  1:4701069/76‑1 (MQ=255)
                          tCATGGCGACAGCCACTTGCGGGTTGTTGTATGGCGCAAAGGCGGCCATCAGCTTGTGGTCACGTAAAAGCTCggg                                            >  2:7776323/1‑75 (MQ=255)
                          tCATGGCAACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGc                                            <  2:9524900/76‑1 (MQ=255)
                           cATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCa                                           >  1:2231114/1‑76 (MQ=255)
                           cATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCAGCa                                           >  1:2835902/1‑76 (MQ=255)
                              ggCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCCTCAGTTTGTGGTCACGTAAACGCTCGGGAAtt                                        >  2:5302298/1‑76 (MQ=255)
                              ggCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCCTCAGTTTGTGGTCACGTAAACGCTCGGCAAtt                                        >  1:5489093/1‑76 (MQ=255)
                              ggCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAAtt                                        >  1:2580011/1‑76 (MQ=255)
                               gcgcCAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAAttt                                       <  1:1842821/76‑1 (MQ=255)
                               gcgcCAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAAttt                                       <  2:6923387/76‑1 (MQ=255)
                                  aCAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTtgt                                    >  2:5755278/1‑76 (MQ=255)
                                    aGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCCTCAGTTTGTGGTCACGGAAACGCTCGGCAATTTTGTgc                                  >  2:6482063/1‑76 (MQ=255)
                                     gCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTgcg                                 <  2:5886454/76‑1 (MQ=255)
                                     gCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTgcg                                 <  1:5767731/76‑1 (MQ=255)
                                         cTTGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCAtt                             <  1:3452167/76‑1 (MQ=255)
                                          ttGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATta                            >  2:2088898/1‑76 (MQ=255)
                                          ttGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATta                            >  2:1500420/1‑76 (MQ=255)
                                          ttGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATta                            >  2:2946679/1‑76 (MQ=255)
                                          ttGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATta                            <  1:932641/76‑1 (MQ=255)
                                          ttGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATta                            <  1:8916690/76‑1 (MQ=255)
                                           tGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATtat                           <  1:8623020/76‑1 (MQ=255)
                                           tGCGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATtat                           <  2:6194791/76‑1 (MQ=255)
                                             cGGATTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGGGCATTATAg                         <  2:406898/76‑1 (MQ=255)
                                                aTTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGtt                      <  1:288670/76‑1 (MQ=255)
                                                aTTGTTGTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGtt                      >  1:2173779/1‑76 (MQ=255)
                                                 ttgttgTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGttt                     <  1:1971200/76‑1 (MQ=255)
                                                 ttgttgTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGttt                     <  1:3124050/76‑1 (MQ=255)
                                                 ttgttgTATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGttt                     >  2:6571632/1‑76 (MQ=255)
                                                  tgttgtATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTttc                    <  2:5807809/76‑1 (MQ=255)
                                                   gttgtATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTttcg                   >  1:1670716/1‑76 (MQ=255)
                                                   gttgtATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTttcg                   >  1:6039371/1‑76 (MQ=255)
                                                   gttgtATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTttcg                   >  1:9140066/1‑76 (MQ=255)
                                                   gttgtATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTttcg                   <  2:3695909/76‑1 (MQ=255)
                                                    ttgtATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTttcgt                  >  2:2607594/1‑76 (MQ=255)
                                                      gtATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTttcgttc                >  2:6102365/1‑76 (MQ=255)
                                                       tATGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTttcgttcg               >  2:1858534/1‑76 (MQ=255)
                                                        aTGGCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGc              <  2:7640993/76‑1 (MQ=255)
                                                          ggCGCAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCtt            <  2:5119089/76‑1 (MQ=255)
                                                            cgcAAAGGCGTCCATCAGGTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTc          <  2:3907001/76‑1 (MQ=255)
                                                            cgcAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTc          <  2:1670716/76‑1 (MQ=255)
                                                            cgcAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTc          <  1:7309090/76‑1 (MQ=255)
                                                              cAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAg        >  1:4584860/1‑76 (MQ=255)
                                                              cAAAGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAg        <  1:1889168/76‑1 (MQ=255)
                                                               aaaGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGa       >  2:4717955/1‑76 (MQ=255)
                                                               aaaGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGa       >  2:1766029/1‑76 (MQ=255)
                                                               aaaGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGa       >  2:4206278/1‑76 (MQ=255)
                                                                aaGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTGGTTCGCTTTCAGAc      >  1:6571246/1‑76 (MQ=255)
                                                                aaGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGAc      <  1:3918606/76‑1 (MQ=255)
                                                                aaGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGAc      >  1:5584187/1‑76 (MQ=255)
                                                                aaGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGAc      <  1:4119866/76‑1 (MQ=255)
                                                                aaGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGAc      <  1:8462539/76‑1 (MQ=255)
                                                                 aGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGAcc     >  1:4734662/1‑76 (MQ=255)
                                                                 aGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGAcc     >  1:9525750/1‑76 (MQ=255)
                                                                 aGGCGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGAcc     >  2:6443314/1‑76 (MQ=255)
                                                                  ggcggcCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGACCg    <  1:9566323/76‑1 (MQ=255)
                                                                  ggcggcCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGACCg    >  2:6785833/1‑76 (MQ=255)
                                                                  ggcggcCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGACCg    <  1:155378/76‑1 (MQ=255)
                                                                  ggcggcCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGACCg    >  2:2840814/1‑76 (MQ=255)
                                                                    cGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGACCGaa  >  1:7460332/1‑76 (MQ=255)
                                                                    cGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGACCGaa  >  1:6969842/1‑76 (MQ=255)
                                                                    cGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGACCGaa  >  1:1284156/1‑76 (MQ=255)
                                                                    cGGCCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGACCGaa  >  1:3660884/1‑76 (MQ=255)
                                                                       |                                                                        
CGGACCCGCGCCGCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAATTTTGTGCGCATTATAGGTTTCGTTCGCTTTCAGACCGAA  >  REL606/649008‑649151

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