Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,140,9260AC25.5% 101.0 / 10.9 51G229G (GGT→GGGasmAputative assembly protein AsmA
Reads supporting (aligned to +/- strand):  ref base A (36/2);  new base C (0/13);  total (36/15)
Fisher's exact test for biased strand distribution p-value = 2.20e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCTGTGCTTGCCCGCTTAGCGTACTATCATTGGCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCAAA  >  NC_000913/2140793‑2141057
                                                                                                                                      |                                                                                                                                   
ccTGTGCTTGCCCGCTTAGCGTACTATCATTGGCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCa                                                                                                                                     >  2:1275251/1‑134 (MQ=255)
ccTGTGCTTGCCCGCTTAGCGTACTATCATTGGCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCa                                                                                                                                     <  1:1275251/134‑1 (MQ=255)
     gCTTGCCCGCTTAGCGTACTATCATTGGCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTg                                                                                                                          >  1:410486/1‑140 (MQ=255)
     gCTTGCCCGCTTAGCGTACTATCATTGGCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTg                                                                                                                          >  2:1041309/1‑140 (MQ=255)
           ccGCTTAGCGTACTATCATTGGCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGtt                                                                                                                    >  1:873902/1‑140 (MQ=255)
              cTTAGCGTACTATCATTGGCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAAt                                                                                                                 >  1:546435/1‑140 (MQ=255)
                         aTCATTGGCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAg                                                                                                      >  1:15824/1‑140 (MQ=255)
                          tCATTGGCGGTCAAACTATTTTGTTAAAATAAAAGGCTTTTATGTGCCCCCTGCCACTGGGCCTGGAAACCCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAgc                                                                                                     <  2:1429411/140‑1 (MQ=255)
                            tttGGGGGTAAAACTAATTTGGTAAAATAAAGGGCTTTTATGTGCCCCCTGCCACTGGGCCTGGAACCCCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTAAATAgccg                                                                                                   <  2:1110080/139‑1 (MQ=255)
                            aTTGGCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAgccc                                                                                                   >  1:883403/1‑139 (MQ=255)
                             ttGGCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAgccgc                                                                                                  >  1:593829/1‑140 (MQ=255)
                             ttGGCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAgccgc                                                                                                  >  1:858967/1‑140 (MQ=255)
                             ttGGCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAgccgc                                                                                                  >  1:990254/1‑140 (MQ=255)
                               ggCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGccgccg                                                                                                >  1:1037029/1‑140 (MQ=255)
                               ggCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGccgccg                                                                                                >  2:1268826/1‑140 (MQ=255)
                               ggCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGccgccg                                                                                                >  2:1369321/1‑140 (MQ=255)
                                  ggTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTa                                                                                             >  1:1514396/1‑140 (MQ=255)
                                  ggTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTa                                                                                             >  2:1438205/1‑140 (MQ=255)
                                       aaCAAATTGGTTAAAATAAAAGGCTTTTAGGTGCCTCCTGCCACTGGGCCGGGAAACCCCCCTGCCCCTGAATCCCTTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTAAATAGCCGCCGTTAAATCa                                                                                        <  2:1514396/140‑1 (MQ=255)
                                         atattttGGTTAAATGAAAGGGTTTTTTGGTCCTCCTGCCACTGGCCCTGGAAACTCCCCTGCCCCTAAATCCTTTGTTTTGGTAAATCGGCACCCTGCAACTGCCGGTTATTTTGTTCAATAGCCGCCGTTAAATCATg                                                                                      <  2:63533/136‑1 (MQ=255)
                                                  ttAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATc                                                                             >  1:63533/1‑140 (MQ=255)
                                                   tAAATGAAGGGCGTTTATGTGCCCCCTGCCACTGGGCCTGGAACCTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCa                                                                            <  1:1397245/140‑1 (MQ=255)
                                                   aaaaTAAAGGGCTTTTGTGTGCCCCCTCCCTCGGGGCCGGAAACCCCCCCTGCCCCTGAATCCCTTGTTTTGGTAAATCGGCCCCCTGCACCTGCCAGTTATTTTGTTAAATAGCCGCCGTTAAATCATGCGGATAATCa                                                                            <  2:885187/139‑1 (MQ=255)
                                                      aTGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGaa                                                                         >  1:1548752/1‑140 (MQ=255)
                                                        gAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGc                                                                       <  1:1369321/140‑1 (MQ=255)
                                                        aaaaGGGTTTTATGTGCCTCCTGCCACTGGGCCTGGAACCCCCCCTGCCCCTGAATACCTTGTTTTGGAAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGc                                                                       <  2:53507/139‑1 (MQ=255)
                                                           aGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAgctgtctc                                                                    >  1:105491/1‑133 (MQ=255)
                                                             gctttttTGTGCCCCCTGCCACTGGGCCTGGAAACCCCCCTGTCCCTGAATCCCTTGTTTTGGTAAACCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCta                                                                  <  2:802293/137‑1 (MQ=255)
                                                                tttATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccg                                                               <  1:1268826/140‑1 (MQ=255)
                                                                tttATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccg                                                               >  1:238259/1‑140 (MQ=255)
                                                                 ttATGTGCCTCCTGCCACTGGGCCTGAAACCCCCCCTGCCCCTAAATCCCTTGTTTTGGTAAATCGCCCCCCTGCAACTGCCATTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgt                                                              <  2:1037029/140‑1 (MQ=255)
                                                                  tttgtgCCCCCTGCCACGGGGCCGGAAACCTCCCCTCCCCCTGAATCCCTTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgta                                                             <  2:773631/138‑1 (MQ=255)
                                                                   ttgtgCCTCCTGCCACTGGGCCTGGAACCCCCCCTGCCCCTGAATCCCTTTTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTATTTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgtac                                                            <  1:133410/139‑1 (MQ=255)
                                                                     gtgCCTCCTGCCACGGGGCCTGGAACCTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgtaccg                                                          <  2:410486/140‑1 (MQ=255)
                                                                        aCTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgtctgtctc                                                       >  1:278279/1‑133 (MQ=255)
                                                                         cTCCTGCCACGGGGCCTGGAACCCCCCCTGCCCCTGAATCCCTTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTaa                                                      <  2:586865/140‑1 (MQ=255)
                                                                          tCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAg                                                     >  1:1553461/1‑140 (MQ=255)
                                                                             tcccACTGGGCCTGGAACCCCCCCTGCCCCTGAATCCTTTTTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGa                                                  <  2:15824/138‑1 (MQ=255)
                                                                             tcccACTGGGCCTGGAACCCCCCCTGCCCCTAAATCCCTTGTTTTGGTAAATCGGCCCCCTGCACCTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGa                                                  <  1:1041309/138‑1 (MQ=255)
                                                                             tGCCCCTGGGCCTGAAACCCCCCCTGCCCCTAAATCCTTTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGa                                                  <  2:546435/140‑1 (MQ=255)
                                                                             cgccccggggCCTGAAACCCCCCCTGCCCCTGAATCCTTTGTTTTCGTAAATCGCCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGa                                                  <  2:1553461/133‑1 (MQ=255)
                                                                              gCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATCCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGata                                                >  2:806008/1‑140 (MQ=255)
                                                                               ccccTGGGCCCGGAAATCTCCCCCGCCCCTGAATCCTTTTTTTTGGTAAATCGCCACCCTACAACTGCCAGTAAATTTGTTAAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGat                                                 <  2:1097097/137‑1 (MQ=255)
                                                                               ccACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGata                                                >  1:586865/1‑140 (MQ=255)
                                                                                    gggCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTg                                           >  2:407597/1‑140 (MQ=255)
                                                                                          ggAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATcg                                     >  1:226481/1‑140 (MQ=255)
                                                                                                cTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTg                                >  1:501403/1‑140 (MQ=255)
                                                                                                 tCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGa                               >  2:512915/1‑140 (MQ=255)
                                                                                                   ccctgcccctgAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGATc                             >  1:280228/1‑140 (MQ=255)
                                                                                                        cccctgAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGATCGCGAt                        >  2:440054/1‑140 (MQ=255)
                                                                                                            tgAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGATCGCGATTAAc                    >  2:1480352/1‑140 (MQ=255)
                                                                                                             gAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGATCGCGATTAAcc                   >  1:1050347/1‑140 (MQ=255)
                                                                                                                  ccTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGATCGCGATTAAcccgccc              >  1:384232/1‑140 (MQ=255)
                                                                                                                   cTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGATCGCGATTAAcccgcccg             >  1:1264603/1‑140 (MQ=255)
                                                                                                                   cTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGATCGCGATTAAcccgcccg             >  2:136202/1‑140 (MQ=255)
                                                                                                                   cTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGATCGCGATTAAcccgcccg             <  1:806008/140‑1 (MQ=255)
                                                                                                                     tGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGATCGCGATTAAcccgcccg             >  2:1420799/1‑138 (MQ=255)
                                                                                                                     tGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGATCGCGATTAAcccgcccg             <  1:1420799/138‑1 (MQ=255)
                                                                                                                        tttGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGATCGCGATTAACCCGCCCGGAGaa        >  2:1512314/1‑140 (MQ=255)
                                                                                                                            gTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCa    >  1:122187/1‑140 (MQ=255)
                                                                                                                              aaaTCGGCACCCTACAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCaaa  <  1:1209053/140‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                   
CCTGTGCTTGCCCGCTTAGCGTACTATCATTGGCGGTCAAACTAATTTGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATCGTGAGATCGCGCTGATCGCGATTAACCCGCCCGGAGAACTCAAA  >  NC_000913/2140793‑2141057

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