Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,008,3920GT31.4% 52.6 / 10.6 35W39L (TGG→TTG) yedEputative selenium transporter YedE
Reads supporting (aligned to +/- strand):  ref base G (12/12);  new base T (10/1);  total (22/13)
Fisher's exact test for biased strand distribution p-value = 2.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGC  >  NC_000913/2008165‑2008641
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                            
gagaAATTTGATCTAACTATTCTTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTTTCTCTTCTTTCTTTTTTGGCCTTTCTTGGACCCTTTTGGCTGTCACCGGTGTAttttacc                                                                                                                                                                                                                                      >  1:114430/1‑247 (MQ=255)
         gATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTTTTTTCTTCACTTTCCATTTTTGCATTACTGGCCCCTTTTTGGCTGTCCACGGGGAATTTACCCGTTGGggcg                                                                                                                                                                                                                            >  2:109515/1‑251 (MQ=255)
              aaCTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGcc                                                                                                                                                                                                                         <  1:305302/249‑1 (MQ=255)
                cTATTATTTTTGTGAATTATTACCCCGAAAATCAATTTGCCTGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTTGCCTTTATGGCCCCCTTTTGGGTTTCACGGGTGCATTTAccctcttggggggccacgc                                                                                                                                                                                                                     >  2:173207/1‑235 (MQ=255)
                             gAATTATTACCCCGAAACTAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGATCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTTCCTCCTCTTTTTGCCTTACTTGCCCCTTTTTGGGTGGCACGGGGGAATTTTTCCGTTGGGGGGGcccggcccctccactgttt                                                                                                                                                                                                        >  2:479057/1‑234 (MQ=255)
                             gAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTc                                                                                                                                                                                                       <  2:131848/251‑1 (MQ=255)
                                         cGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAAACATCTCAGTTCAGCGGCTACATTCATGCCATGGCCGCACTTCAAACACGCCCGGTTGATTACATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGCCTTACTCTCTACATACCATTTTTGCAATTCTGGCACCTTTTTGGCAGCTACGGGTGGATTTACCCGTTGGGGGGGCCAGCCCCTGCACCTGTTCGGCGtacctgct                                                                                                                                                                                           >  1:503469/1‑244 (MQ=255)
                                                aaatTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTTTCTCTACTTACTATTTTGGGATTACTGGCACCTTTTTGGCTGTCCCGGGGGAATTTACCCGTTTGGGCGGCCCGCTCCTGCCACTGTTCGGCGTCCCTGCTGCagagg                                                                                                                                                                                    >  1:232491/1‑250 (MQ=255)
                                                          ggTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTTTTTTTGCCTTTCTGGCACCTTTTTGGGTGTCCCGGGTGGATTTTCCCGTTCGGGCGGCCAGCTCCTGCCACTGTTGGGCGTCCATGCTGGAGAGTggggggtctt                                                                                                                                                                          >  2:449686/1‑245 (MQ=255)
                                                           gtgtATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTCCTTACTATTTTTGGATTACTGGCACCTTTTGGGGTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACCGTTCGGGGGCCATGCTGAAGAGGGGGGTTACttt                                                                                                                                                                         >  1:181032/1‑251 (MQ=255)
                                                             gtATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCTTCCATGCTGAAGAGTGGGGTTACTTTa                                                                                                                                                                        <  1:203962/250‑1 (MQ=255)
                                                                    tCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTAt                                                                                                                                                                 <  1:110163/250‑1 (MQ=255)
                                                                        gACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGGCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATc                                                                                                                                                            <  1:192350/251‑1 (MQ=255)
                                                                                                  cAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTACATTCTGGGCCCCCATCACTGCGGTCATCGCGGCGGGGATTTTCTCCACTTTCTATTTTTGGATTACTTGCACCTTTTTGGCTGTCACGGGGGGATTTTCCCGTTTGGGCGGCACGCTCCTACAACTGTTCGGCGGCCATGCTTCAGAGTGGGGGTACTTTTAAATTATCCCTCTGGGAGGATTGCCATTAAACCGCacc                                                                                                                                  >  2:452294/1‑249 (MQ=255)
                                                                                                  cAGAGGCTACATTCATGTAATGGCCGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCAATCACTGCGGTCATCGCGGCGGGGATTTTCTCTTATTACCATTTTTGCATTACTTGCCCCTTTTTGGCTTTTACGGGGGAATTTTCCCGTTGGGGGGGCCAGCTCCCGGAACTGTTTGGGGGTCCTGCTGCAGAGTGGGGTTTATTTTAAATTAACCATCTGGAAgggtacgcatatacccgcatc                                                                                                                                  >  2:354705/1‑232 (MQ=255)
                                                                                                               catgtcatgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGgatgatga                                                                                                                     <  2:232491/251‑1 (MQ=255)
                                                                                                                    catgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTAATTTCTCTACTTACTATTTTTGCATTTCTGGCACCTTTTTGGGTTTCCCGGGGGGATTTTCCCGTTGGGGGGGCCAGCTCCTGCAACTGTTCGGCGGCCCTGCTGAAGAGTGGGGTTACTTTATAATTATCCATCTGGAAAGATCGCCCTTAACCACCCTCGACGGGATGATGATcccc                                                                                                                >  2:120963/1‑249 (MQ=255)
                                                                                                                                tCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTgg                                                                                                     <  1:449686/250‑1 (MQ=255)
                                                                                                                                  aaaCACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCCCGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGGTCGGCGGCCATGCTGAAGAGTGGGGGTACTTTTAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGCCGGGGTGATGATCCTCGGGATGTTTggcgg                                                                                                  >  1:10320/1‑251 (MQ=255)
                                                                                                                                            ggTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGCTTCTCCCTACTTCCTATTTTGGCATTACTGGCCCCCTTTGGGCCGTCAAGGGGGAATTTAACCGTTGGGGCGGCCAGCCCCTGCCACTGTTCGGCGTCCATGCTGAAGAGTGGGGGTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCCACGGGATGGTGATCCTCGGTATGTTTGGCGGCTGCTTTGcc                                                                                        >  1:267958/1‑251 (MQ=255)
                                                                                                                                              ttGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCATTAATTAATATTTTGGCAGTACTGGCACCCTTTTGGCTGTCACCAGGGAATTTACCACGTTGGGCGGGCAGCTCCTGCCACTGTTTGGGGTCCAAGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGGTCGCCACTAAACCCCCTCCACGGGTAGATGATCCTCGGTATGTTTGGCGGCTTCTTTGcagc                                                                                      >  1:43199/1‑248 (MQ=255)
                                                                                                                                                            ggCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGGGAATTTACCCGTTGGGGCGGCCCGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGGTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGccc                                                                        >  1:387138/1‑250 (MQ=255)
                                                                                                                                                                cccATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTTCTATTTTTGCATTACTGGCAACCTTTTGGCTGTCACGGGTGAATTTAACCGTTGGGGCGGCCAGCTCCTGCCACTGTTCGGCGTCCATGCTGAAGAGTGGAGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAAACCGCATCGAAGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCCGTGGGCcaaca                                                                     >  1:28102/1‑250 (MQ=255)
                                                                                                                                                                              cATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGgct                                                                                                <  1:138950/209‑1 (MQ=255)
                                                                                                                                                                              cATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCAACTGTAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTGGTCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGa                                                       <  1:354705/250‑1 (MQ=255)
                                                                                                                                                                                           aTTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCg                                         >  2:140512/1‑251 (MQ=255)
                                                                                                                                                                                             tctctcTACTTACGTTTTTGGCGTTACTGGCACCTTTGTGCCTGTAACGGGTGAATTTTCCCGTTGGGGCGGCCAGCTCCGGCAACTGTTCGTCTTCCATGCGGAGGAGTGGTGTTACTATAAAATTATCCATCTGGAAGGATCGCCATTAGCCCGCATCGACGGGATGATGATCCTCTGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGGGGGCTAACAATGTGAAGCTGCGAATGGCGGGCAGCCGt                                        <  2:452430/250‑1 (MQ=255)
                                                                                                                                                                                                ctctACTTACTATTTTGGCATTACTGGCACATTTTGGGCTGTCACGGGAGAATTTACCCGTTGGGGCGGCCAGCTCCTGCCACTGTTCGGCGTCCCTGCTGAAGAGTGGGGTTACTTTAACATTCTCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGACTGCCGCGCAGCCGTATcc                                    >  2:348004/1‑251 (MQ=255)
                                                                                                                                                                                                       tactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCAtt                               <  1:120963/249‑1 (MQ=255)
                                                                                                                                                                                                        actATTTTGGCCTTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGACGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCACCAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCTCATTATg                            >  2:358431/1‑251 (MQ=255)
                                                                                                                                                                                                                ggCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCcat                    <  1:56666/251‑1 (MQ=255)
                                                                                                                                                                                                                   aTTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCcatcat                 >  2:436730/1‑251 (MQ=255)
                                                                                                                                                                                                                     tACTGGCACCTTTTGGGCTGTCACGGGCGAATTTACCCATTGGGGCGCCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGACGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGTTCGCCATTAACCCGCATCGACGGGCTGCTGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCCACAATGTCAAACTGCGAATGCCGCGCAGACGTATCCGCATTATGCAGGCCATCATTg               >  1:237832/1‑251 (MQ=255)
                                                                                                                                                                                                                       cTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTAt                                                                                                                                                                 >  2:387879/1‑103 (MQ=255)
                                                                                                                                                                                                                       cTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTAt                                                                                                                                                                 <  1:387879/103‑1 (MQ=255)
                                                                                                                                                                                                                                     gggCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGc  <  1:407201/250‑1 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                            
GAGAAATTTGATCTAACTATTATTTTTGTGAATTATTACCCCGAAAATAAATTTGCATGGTGTATGATTCGCGACTCTTGTTTCAAACCTTCTGAGTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGC  >  NC_000913/2008165‑2008641

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