Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,157,900 +A coding (32/1434 nt) ptsG → fused glucose‑specific PTS enzymes: IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,157,8961.A91.2% 164.8 / 9.2 57coding (28/1434 nt)ptsGfused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base . (5/0);  new base A (21/31);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTCTGTGTAATAAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCC  >  NC_000913/1157755‑1158030
                                                                                                                                                |                                                                                                                                         
tACTCTGTGGAATAAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCa                                                                                                                                          <  1:73925/143‑1 (MQ=255)
 aCTCTGTGTAATAAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCaa                                                                                                                                        <  1:671757/143‑1 (MQ=255)
         taataaATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAAATATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCg                                                                                                                               >  1:43739/1‑143 (MQ=255)
         taataaACAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCCTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCg                                                                                                                               <  1:481778/143‑1 (MQ=255)
          aataaataAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCgg                                                                                                                              <  1:134624/143‑1 (MQ=255)
          aataaataAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATCTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCgg                                                                                                                              <  1:586801/143‑1 (MQ=255)
            taaataaaGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGTGTGAGAACGTAAAAACAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTa                                                                                                                            >  2:300592/1‑143 (MQ=255)
               ataaaGGGCGCTTAGATGCCCTGTACACGGCGAGTCTCTCCCCCCTTGCCAAGCGTGAGCACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTCTGTTTAACAACGCATTTGCTAACCTGCAAAAAGGTCGGTAAAt                                                                                                                         <  2:790821/143‑1 (MQ=255)
               ataaaGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAAt                                                                                                                         <  2:752011/143‑1 (MQ=255)
                      ggccTTGGATGCCCTGCACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGAt                                                                                                                  <  1:508822/140‑1 (MQ=255)
                            aGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGcc                                                                                                            <  1:133339/143‑1 (MQ=255)
                            aGATGCCCTGGACACGGCGAGGCACTCCAACACATGTACTCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGcc                                                                                                            <  1:613628/143‑1 (MQ=255)
                             gATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCg                                                                                                           <  1:38339/143‑1 (MQ=255)
                             gATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCAAGCGTGAGAACGTAAAAAAAGCACCCATACGCAGGGGCAGTCGCAATTAGGTTTAAGAATGCATTTGCTAACCTGCCAAATGGTGCGTGAATCGCGGAtgctgcag                                                                                                           >  2:862659/1‑141 (MQ=255)
                             gATGCCCTGTACACGGCGAGGCGCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCg                                                                                                           <  1:849340/143‑1 (MQ=255)
                                    tGTACACGGCGAGGCACTCCCCCCTTGCCACGCGTGAGAACGTGAAGAGAGCACCCAAACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCg                                                                                                    <  1:658832/143‑1 (MQ=255)
                                      tACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTa                                                                                                  <  1:23880/143‑1 (MQ=255)
                                       acacGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAACAGCACCCATACTCAGGAACAGTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCACAAAGGTCGGTAAAGCGCTGATGCACCAGGTAACCGtcc                                                                                                 >  2:837973/1‑141 (MQ=255)
                                         acGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACACTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGGAAATCGCTGATGCTGGCGGTATCCGTACTg                                                                                               >  1:298478/1‑143 (MQ=255)
                                          cGGCGAGGCGCTCCCCCCTGGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGc                                                                                              <  1:132918/143‑1 (MQ=255)
                                           ggCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGcc                                                                                             >  1:588574/1‑143 (MQ=255)
                                             caacgCTCTCCCCCCTTGCCACGCATGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTa                                                                                           <  2:898461/139‑1 (MQ=255)
                                              gAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTAt                                                                                          >  2:619143/1‑143 (MQ=255)
                                                  ctctCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTCTCGCa                                                                                      >  1:397072/1‑143 (MQ=255)
                                                      ccccccTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCCTTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTa                                                                                  >  1:824575/1‑143 (MQ=255)
                                                      ccccccTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTAGGTAAATCACTGATGCTGCCGGTATCCGTACTGCCTATAGAAGGTa                                                                                  >  1:476741/1‑143 (MQ=255)
                                                      ccccccTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCGCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTa                                                                                  <  2:101012/143‑1 (MQ=255)
                                                         ccccTTGCCACGCGCGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTAtt                                                                                <  2:482299/143‑1 (MQ=255)
                                                          cccTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTc                                                                               <  1:438479/143‑1 (MQ=255)
                                                              tGCCACGCGGGAGAACGTAGAAAAAGCACCCATACTCAGGAGCACTCTCAGTTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTACATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTctgct                                                                           <  1:738100/143‑1 (MQ=255)
                                                                ccACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCCATCATGTTTAAGAATGCATTTGCTCACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCCGGTATTCGGCTgg                                                                         >  2:737980/1‑143 (MQ=255)
                                                                 cACGCGTGAGAACGTAAAAAAAGCACCCATACTCCGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAACTCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCGCGCATTctgctgcg                                                                        >  2:685431/1‑141 (MQ=255)
                                                                    gggTGGGAACGTACAAAAAGCACCCGTACTCAGGAGCACGCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAACCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGt                                                                     <  2:159564/141‑1 (MQ=255)
                                                                     cgTGAGAACGTAGCAAAAGCTCCCTTACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGTTGGGCGTc                                                                    <  2:758551/143‑1 (MQ=255)
                                                                     cgTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGCTATCCGTACTGCCTATGGCTGGTATTGTGCGTGGCGTc                                                                    >  1:245001/1‑143 (MQ=255)
                                                                           caCGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTcc                                                              <  2:508387/142‑1 (MQ=255)
                                                                              gTAAAAAAAGCACCCATACTCAGAAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCCGTTCcgcg                                                           >  2:315482/1‑143 (MQ=255)
                                                                                aagaaaaGCACCCTTGCTCAGGAGCCCTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGaa                                                         <  2:621737/140‑1 (MQ=255)
                                                                                 aaaaaaGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTCACCTGCAAGACGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGGAt                                                        >  1:475187/1‑143 (MQ=255)
                                                                                 aaaaaaGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAAt                                                        <  2:798578/143‑1 (MQ=255)
                                                                                 aaaaaaGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAAt                                                        <  1:132805/143‑1 (MQ=255)
                                                                                  aaaaaGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAAACTGCAAAAAGGTCGGTCAATCGCTGATGCTGCAGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCcgcacatt                                                       >  2:26273/1‑138 (MQ=255)
                                                                                       gCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAgc                                                  <  1:496930/143‑1 (MQ=255)
                                                                                          cccGTACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTACATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAgctgg                                               <  2:514113/143‑1 (MQ=255)
                                                                                           ccatgttCAGGAGCCCCCTCAATTATGTTTAAGATTGCACTTGCGAACCTGCAAAAAGGTCGGTAAATCGCTGATCCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAgctggc                                              <  2:293771/137‑1 (MQ=255)
                                                                                           ccATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTTCTTGGCGTAGGTTCCGCGAATTTCAgctggc                                              >  1:562745/1‑143 (MQ=255)
                                                                                           ccATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGCATTTCAgctggc                                              >  2:412479/1‑143 (MQ=255)
                                                                                           ccATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGCCCCGCGCATTTCAgctggc                                              >  2:685954/1‑143 (MQ=255)
                                                                                                   aGGAGCACTCTCAATTATGTTTAAGAATGCATGTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCccgcc                                      <  2:667617/143‑1 (MQ=255)
                                                                                                    gtaggcctctCAATTATGTGTAAGAATGTATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCGATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCccgccg                                     <  2:539734/137‑1 (MQ=255)
                                                                                                     gAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGGATTTCGGGTGGCTGCGTGCCGt                                    >  1:78599/1‑143 (MQ=255)
                                                                                                        cACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGCACTGCCTATCGCAGGTCTTCTGCTGGGCGTCGGTCCCGCGACTCACAGCTGTCGGCCCGCCGttct                                 >  2:183883/1‑141 (MQ=255)
                                                                                                        cACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAACGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTAGCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGATGCCCGCCGTTGt                                 >  2:853460/1‑143 (MQ=255)
                                                                                                           tctcAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAAGTTCAGCTGGCTGCCCGCCGATATATc                              >  1:167664/1‑143 (MQ=255)
                                                                                                              cAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTGTCGCa                           >  1:620748/1‑143 (MQ=255)
                                                                                                                     gTTTAAGAATGCATTTGCTAGCCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATg                    <  2:360064/143‑1 (MQ=255)
                                                                                                                                 aTTTGCTAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGATGTATCGCATGTTATGGCAGAAGCAGGc        >  2:512084/1‑143 (MQ=255)
                                                                                                                                     tcTAACCTGCAAAAAGCGCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTGTCGCATGTTATGGCAGAAGCAGGCGGtt    <  2:739114/142‑1 (MQ=255)
                                                                                                                                       tAACCTGCAAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTcc  >  1:717650/1‑143 (MQ=255)
                                                                                                                                                |                                                                                                                                         
TACTCTGTGTAATAAATAAAGGGCGCTTAGATGCCCTGTACACGGCGAGGCTCTCCCCCCTTGCCACGCGTGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCC  >  NC_000913/1157755‑1158030

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