Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,991,256 A→T 6.7% Y35F (TAC→TTC)  cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,991,2560AT 6.7% 183.7 / 3.2 60Y35F (TAC→TTC) cyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (31/25);  new base T (2/2);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTC  >  NC_000913/3991109‑3991403
                                                                                                                                                   |                                                                                                                                                   
ggCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTa                                                                                                                                                     <  4:219127/148‑1 (MQ=255)
 gCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACa                                                                                                                                                   <  7:37207/149‑1 (MQ=255)
      aTCACAGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTa                                                                                                                                              <  2:115365/149‑1 (MQ=255)
          caGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGc                                                                                                                                          <  5:21937/149‑1 (MQ=255)
          caGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCGTGCTGCTATGGTGCCTGAATTCCAACAGGTCTACAGTCTACTGc                                                                                                                                          >  8:122220/1‑149 (MQ=255)
           aGTCATGACGGGGAGCAAATCAGTCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTccc                                                                                                                                         <  6:42651/149‑1 (MQ=255)
            gTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGATACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTCCTGCTATCCGGAATGCATTCCAAAAGGACTACAGTCTACTGCCg                                                                                                                                        >  6:20490/1‑149 (MQ=255)
               aTGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACa                                                                                                                                     >  3:26267/1‑149 (MQ=255)
                    gggaaccAAAGCCGGCGCTACGCCATGCACCTCTATATTGCGACTCGGCAACAGAGACTGGAGGCCATAACTCAATTGCGGCTGGCCCGCGCGCGTGCTGCTACGGGGCCTGCATTCCACCCGGGCTACAGTCTACTGCCGACAttgtt                                                                                                                                <  8:236906/143‑1 (MQ=255)
                    gggTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACAttgtt                                                                                                                                >  3:76972/1‑149 (MQ=255)
                       tAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCa                                                                                                                             >  6:129332/1‑149 (MQ=255)
                         gCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCAc                                                                                                                            <  7:192242/148‑1 (MQ=255)
                              tCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCt                                                                                                                                               <  1:71574/124‑1 (MQ=255)
                              tCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCt                                                                                                                                               >  2:71574/1‑124 (MQ=255)
                               cAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCAc                                                                                                                      >  3:157327/1‑148 (MQ=255)
                                         gTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAAt                                                                                                           >  1:112257/1‑149 (MQ=255)
                                                aCCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGt                                                                                                    >  7:258/1‑149 (MQ=255)
                                                       tatTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCt                                                                                               <  7:35059/147‑1 (MQ=255)
                                                       tatTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTAGATCTCGCTCTTCCACCTATCATTAAAAAAAAACAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCtt                                                                                              >  2:34948/1‑148 (MQ=255)
                                                        atTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTTCAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGa                                                                                            >  4:179401/1‑149 (MQ=255)
                                                        atTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTTCAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGa                                                                                            >  4:182661/1‑149 (MQ=255)
                                                        atTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGa                                                                                            >  4:54339/1‑149 (MQ=255)
                                                          tGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATAGTTGCACTATCACCATCCGCTAATGCCGGGTTACAatgatg                                                                                          >  5:143694/1‑149 (MQ=255)
                                                           gagaCTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATg                                                                                          <  3:79456/148‑1 (MQ=255)
                                                           gagaCTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATg                                                                                          <  5:129332/148‑1 (MQ=255)
                                                            agaCTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATgg                                                                                         >  7:53043/1‑148 (MQ=255)
                                                                  tGAAACAGAGACTGGATGCCATAAATCAATTGCGTCTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATACCAACAGTTCTACATTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGCGTTACATTGAAGGTAACGtt                                                                                  >  2:137548/1‑149 (MQ=255)
                                                                        agagACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCa                                                                              >  1:79122/1‑147 (MQ=255)
                                                                         gagaCTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGc                                                                                                                              <  6:223552/98‑1 (MQ=255)
                                                                         gagaCTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGc                                                                                                                              >  5:223552/1‑98 (MQ=255)
                                                                         gagaCTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAg                                                                           <  4:26267/149‑1 (MQ=255)
                                                                                gATGCCATAAATGTATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCAATCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTTATGCCGGGTTACCTTGATGGTAACGTTCCGAAAGGCAttt                                                                     <  1:137548/148‑1 (MQ=255)
                                                                                 aTGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGc                                                                   <  2:79122/149‑1 (MQ=255)
                                                                                     cATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATgg                                                                                         >  4:7100/1‑123 (MQ=255)
                                                                                     cATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATgg                                                                                         <  3:7100/123‑1 (MQ=255)
                                                                                      aTAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTa                                                              >  5:68226/1‑149 (MQ=255)
                                                                                         aaTCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTa                                                              <  4:161192/146‑1 (MQ=255)
                                                                                         aaTCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTa                                                              >  3:161192/1‑146 (MQ=255)
                                                                                                        gATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGc                                                         >  8:129830/1‑136 (MQ=255)
                                                                                                        gATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGc                                                         <  7:129830/136‑1 (MQ=255)
                                                                                                         aTCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCCTCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCaa                                           >  1:84806/1‑149 (MQ=255)
                                                                                                          tCGCGCGCTTGCTGCTATGTGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAAc                                          <  4:106342/149‑1 (MQ=255)
                                                                                                          tCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTTCAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAAc                                          <  3:182661/149‑1 (MQ=255)
                                                                                                          tCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTTCAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAAc                                          <  3:179401/149‑1 (MQ=255)
                                                                                                              gcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAAc                                          <  1:183500/145‑1 (MQ=255)
                                                                                                              gcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAAc                                          >  2:183500/1‑145 (MQ=255)
                                                                                                                    gctgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCt                                >  7:9252/1‑149 (MQ=255)
                                                                                                                      tgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACAATGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTg                               >  8:130732/1‑148 (MQ=255)
                                                                                                                         tATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGaaa                                                >  8:140304/1‑128 (MQ=255)
                                                                                                                         tATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGaaa                                                <  7:140304/128‑1 (MQ=255)
                                                                                                                               gCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCCCCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCtt                      >  4:221773/1‑148 (MQ=255)
                                                                                                                               gCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCtt                      >  7:15192/1‑148 (MQ=255)
                                                                                                                                   gCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACt                 >  7:120371/1‑149 (MQ=255)
                                                                                                                                    cATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACt                 <  4:157327/148‑1 (MQ=255)
                                                                                                                                      ttCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTa              <  8:53043/149‑1 (MQ=255)
                                                                                                                                       tCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTcc                                                                                >  6:156047/1‑82 (MQ=255)
                                                                                                                                       tCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTcc                                                                                <  5:156047/82‑1 (MQ=255)
                                                                                                                                          aaCAGGTCTACAGTCTACTGCCGACATTGTTGCACAATCACCATCCGCTAATTCCGGGTTACCTTGATGGTAACGTTCCCAAACGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGt          >  4:197100/1‑149 (MQ=255)
                                                                                                                                           aCAAGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTg         >  8:5265/1‑149 (MQ=255)
                                                                                                                                             aGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGa       <  1:34948/149‑1 (MQ=255)
                                                                                                                                                  tACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTc  <  2:112257/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCGGAATCACAGTCATGACGGGTAGCAAATCAGGCGATACGTCTTGTACCTCTATATTGAGACTCTGAAACAGAGACTGGATGCCATAAATCAATTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGTGGAATGTC  >  NC_000913/3991109‑3991403

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