Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,628,6310TA35.3% 52.6 / 14.2 34intergenic (+84/+224)nadR/ettADNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base T (19/3);  new base A (0/12);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 8.30e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTT  >  NC_000913/4628402‑4628878
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        
tGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAACCAAGGTTTGAAAAAACCt                                                                                                                                                                                                                                    >  2:91999/1‑251 (MQ=255)
       cAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACggg                                                                                                                                                                                                                             >  2:99474/1‑251 (MQ=255)
        aaaGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGt                                                                                                                                                                                                                            >  1:76961/1‑251 (MQ=255)
            agTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGAcc                                                                                                                                                                                                                        <  2:366535/251‑1 (MQ=255)
             gTTCCAGAACTTGCTGGGGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGAGAccc                                                                                                                                                                                                                       >  2:50299/1‑251 (MQ=255)
                  aGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCAc                                                                                                                                                                                                                  >  1:282855/1‑251 (MQ=255)
                         gCTGGGGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACGGGTGCGGGAGATGAGGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGATTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTg                                                                                                                                                                                                           <  1:181810/251‑1 (MQ=255)
                             ggggAGATGCTAGAAGAGAACAATCTCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAgg                                                                                                                                                                                                        >  2:25959/3‑250 (MQ=255)
                                                        gattttgggggggTTAAAGAGGAAGATTACGAAAGTGGTTTCCGGCGCTGCGTGGAACTGGTGCGGGGGTTAATGGGGGAGAAAGAATAACGGTGATGAAACTGCTAAAAGGCGAGGTAAAAAATGAGTTTTTTTGATGAGTTGAAAC‑‑CCTCTGTGGAAGAGGCGGTCGAGATTAAACAAGGTTTAAAAAAACCTGCCGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACg                                                                                                                                                                            <  1:247582/241‑1 (MQ=255)
                                                        gAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACg                                                                                                                                                                            >  2:155217/1‑251 (MQ=255)
                                                          aTTTGTGCGGGTGGAAGAGGAAGATGACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTAAAAC‑‑CCTCTCTGGAAGAGGCTGTCGAGATAAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGgt                                                                                                                                                                          <  1:155217/251‑1 (MQ=255)
                                                           tttGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGgtg                                                                                                                                                                         >  2:117521/1‑251 (MQ=255)
                                                              gggggggTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACGGGTGGGGGAGATGAGGGGGGAGCAGAGATAACCGTGATGAAACTGCTAAAGGGCGAGGTAAAAAATGAGTTTTTTTGATGAGTTAAAAC‑‑CCTCTCGGGAAGAGGCGGTCGAAATAAAACAAGGTTTGAAAAACCCGGCAGGGGGGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGtt                                                                                                                                                                       <  2:130348/246‑1 (MQ=255)
                                                                        aaGAGGAAGATTACGACATTCGTTTCCTGCGCTGCGTGGAACTGGTGGGGGAGATGAGGGGGGAGCAGAGATAACCGTGATAAACCTGCTAAAGGGCGAGGTATAAAATGAGTTTTTTTGATGAGTGAAAA‑‑CCCTCTCTGGAAGAGGTTGTCGAGATAAAAAAAGGTTGAAAAAAACCTGCCGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATc                                                                                                                                                            <  1:316668/251‑1 (MQ=255)
                                                                          gagGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCg                                                                                                                                                          >  2:116678/1‑251 (MQ=255)
                                                                                    aCGACAGTCGTTTCCGGCGTGGCGGGGAACGGGTGGGGGAGATGAGGGGGGAGCAGAAATAACCGTGATAAAACTGCTCAAAGGCGAGATAAAAAATGAGTTTTTTTGATGAGTTAAAAA‑‑CCTCTCGGAAAGAGGCTGTCGAAATAAAACAAGGTTTAAAAAACCCTGCCGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTg                                                                                                                                                <  1:116678/251‑1 (MQ=255)
                                                                                                 ccTGCGCTGCGTGAAAGTGGTGGGGGAGATGTGGGGGGAGCAGAGTAACCCGGAATGAATCTGCCAAAAGGGGAGGAATAAAAGGATTTTTTGGGATGA‑‑TGAAAAACCCTCTCGGGAGGGGGCGGTCGAGATAAAAAAAGGTTGAAAAAACCCTGCGGGTTTGTCCCGCCACGAAAGTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCAGGCCGCATCCGTCATTCGGGGCTCAATGCCTGAt                                                                                                                                   <  2:249137/251‑1 (MQ=255)
                                                                                                      gctgcgTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGAc                                                                                                                              >  2:371050/1‑251 (MQ=255)
                                                                                                                   ggTGCGGGGGATGAGGGGGGAGGAGAGAAAACCGGAAAAAAACTGCTCAAAGGCGGGGAAAAAAATGATTTTTTGGAAGGGGTTAAAAC‑‑CCGCTCGGGAAGGGTCTGACGAAATAAAAAAAGGTTAAAAAAACCCGGAAGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCt                                                                                                                 <  2:280009/251‑1 (MQ=255)
                                                                                                                        gggAGATGATGGGGGATCAGATAAACCCGTGATAAATCGGATAAAAGGCGAGGTATAAAATGAGTTTTTTTGTTGAGTTGAAAA‑‑CCTCTTTGGAAGAGGTTGTCGAGTTAAAAAAAGGTTTGAAAAACCCTGCCTGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCa                                                                                                            <  1:226313/251‑1 (MQ=255)
                                                                                                                          gagaggaggggggAGAAGAGATAACCGGGATAAACCGGCTCAAGGGGGAGTTAAAAAATGATTTTTTTTGATGAGTGAAAA‑‑CCCTCTCGGGAAGGGGCTGTCGAAATAAAAAAAGGTTTAAAAAACCCGGCACGGGGGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGTAAACGGTGTTCACGCCGCAGCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAgg                                                                                                          <  1:50299/243‑1 (MQ=255)
                                                                                                                           agaTGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGc                                                                                                         >  1:126486/1‑251 (MQ=255)
                                                                                                                               gatgGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTAc                                                                                                     >  1:65142/1‑251 (MQ=255)
                                                                                                                                        gCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCt                                                                                                                 >  1:183127/1‑230 (MQ=255)
                                                                                                                                        gCAGAGATAACCGTGATAAAACGGATAAAGGGGGAGGTAAAAAATGATTTTTTTGGATGATTAAAAAA‑‑CCTCTCTGGAAGGGGTTGGCGAAATAAAAAAAGGTTAAAAAAACCCTGCAGGGGGGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACa                                                                                            <  2:65142/251‑1 (MQ=255)
                                                                                                                                                         atgaAACTGCTAAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGCTTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACACCGTATTGAAATTGc                                                                             >  2:19668/1‑251 (MQ=255)
                                                                                                                                                                             ggaataAAAGGGGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCGGGAAGAGGCGGGCGAAATAAAAAAGGGTTTAAAAAACCCTGCCGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATaa                                                         <  2:28700/248‑1 (MQ=255)
                                                                                                                                                                                                ggatgaGTTAAAAA‑‑CCTCTCGGGAAGAGGCGGTCGAGATAAAACAAGTTTTAAAAAACCCTGCCGGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATcc                                       <  1:109367/249‑1 (MQ=255)
                                                                                                                                                                                                 gatgaGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTACGGCTGTTGTAGACCGGAAACGGCGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCgg                                     >  1:307542/1‑251 (MQ=255)
                                                                                                                                                                                                           aaaa‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGCGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGc                                                                                                                                                        >  2:45798/1‑127 (MQ=255)
                                                                                                                                                                                                                    tctctGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTAc                      >  1:266841/1‑251 (MQ=255)
                                                                                                                                                                                                                    tctctGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTAc                      >  2:57858/1‑251 (MQ=255)
                                                                                                                                                                                                                       ctGGAAGGGGCTGGCGAAATAAAAAAAGGTTTAAAAAACCCGGCCGGGGGGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACt                     <  1:119585/249‑1 (MQ=255)
                                                                                                                                                                                                                                        gaTTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACtt  >  1:2019/1‑251 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        
TGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAA‑‑CCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTT  >  NC_000913/4628402‑4628878

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