Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,637,814 A→G G373G (GGA→GGG creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,637,8140AG94.3% 94.9 / ‑2.0 35G373G (GGA→GGGcreCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (10/23);  total (10/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.36e-01

CGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCG  >  NC_000913/4637569‑4638035
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                             
agTTTTACCGACAACATTCAGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGTAAGACCGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGTGAAGCGCTCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACCGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAAttt                                                                                                                                                                                                                          <  1:371520/250‑1 (MQ=255)
               aTTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGTAGCAGGCGCTGGGGAATTTACTGGATAACGCCAt                                                                                                                                                                                                           <  1:120509/251‑1 (MQ=255)
                 tCTGCCGCAAAATGCGCGGATGCAGGTATTGGCTGAAACGTTCCTACGCTATGCAAGACTGGAGACTCTTCATGAAGTTGTTCTGACTGCTGCTGATGTGGAGGCATTATCCCGCCGCGTCTGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAACCACTCTGCAGGTTTCGCCCACCGAGGTTAACGCTGCTGCTGAACCTGCGTTACTGGCTCAGGCGCTGGGAAATTTAATGGATAACGTCATc                                                                                                                                                                                                          <  2:99140/250‑1 (MQ=255)
                  cTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGCCTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTTCCCGCGTCTGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTCCATGTTATGCCAACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTCCTGGTTAACCCCATCg                                                                                                                                                                                                         >  1:22941/1‑250 (MQ=255)
                          aaaTGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTAccc                                                                                                                                                                                                >  1:47727/1‑251 (MQ=255)
                            aTGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTAccccc                                                                                                                                                                                              <  1:3141/251‑1 (MQ=255)
                                       cAGGCATTGGGAGCAACGTTCCTACTCCAGGCAAGCCTGTAGAATCGTCCGGAAGTCTTTCTGCCTGTTGTTGATGTGGCGGCATTAGTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTGCTTGATAACGCCATCGATTTTACCCCCGAGATCGGtc                                                                                                                                                                                    <  2:315553/250‑2 (MQ=255)
                                          gCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCAt                                                                                                                                                                                <  1:209287/251‑1 (MQ=255)
                                                gTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCt                                                                                                                                                                          <  1:197893/251‑1 (MQ=255)
                                                               cGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGa                                                                                                                                                           <  1:41988/251‑1 (MQ=255)
                                                                 ccAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTgggtc                                                                                                                                                         >  1:104248/1‑248 (MQ=255)
                                                                                                    tGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGt                                                                                                                      >  1:202727/1‑251 (MQ=255)
                                                                                                                       ggcATTATTTCGCCGCGGCAGCGAATCGATCACCGTGCAGTTGGCAGAAAAAAATATCCCTTTACATGTTATGCCCTCTGAGGTTAATTTTGATGCTGAACCGGCTTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCTCCCTTAAGGTGCTGTAGACCGGTAGTGGTATTCCTGACTAcg                                                                                                   <  2:63590/251‑1 (MQ=255)
                                                                                                                             atttcgccgcGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAac                                                                                                                                                   <  1:365738/193‑1 (MQ=255)
                                                                                                                             atttcgccgcGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCAGAAGTGGATCAGGAac                                                                                                                                                   >  2:365738/5‑197 (MQ=255)
                                                                                                                                 cgccgcGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACg                                                                                         <  1:248712/251‑1 (MQ=255)
                                                                                                                                   ctgcgTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCTACCGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCTCTGGGGAATTTACTGGATAACGCCAGCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTTACTACGCGCTGTCACg                                                                                         <  2:37153/247‑1 (MQ=255)
                                                                                                                                           gCGAAGCGCGCACCGTGCAGTTTGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGGTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGATGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGaa                                                                               <  2:104248/251‑1 (MQ=255)
                                                                                                                                             gAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATGCCGGTAGAGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACg                                                                             >  1:263088/1‑251 (MQ=255)
                                                                                                                                             gAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACg                                                                             >  1:371013/1‑251 (MQ=255)
                                                                                                                                                 cgcgcACCGTGCAGTTGGCAGAAAAAAAAATCAATCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCttt                                                                         <  1:59037/251‑1 (MQ=255)
                                                                                                                                                  gcgcTCCGGGCAGTTGGCTGAAAAAAGCATCACGCGACATGTTATGCCCGCTGAGGTTAATGTGGCTGCTGACCCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCtttt                                                                        <  2:114885/251‑1 (MQ=255)
                                                                                                                                                                      aaaaaaacaTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCTGTCGCATAACTCTAAGCGCCGAAGTGGATCATGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAAt                                                    <  1:240021/251‑1 (MQ=255)
                                                                                                                                                                         aaaaaaTCACTCTACATGTTATGCCCACTGAGGTTAATTTTGCTGCTGAACCGTCGTTACTTGAGCAGGCTCTGTGGAATTTACTTGATAACGCCATCGATTTTACCCCCGAGAGCGTTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAATGTGCTGGATACCGGTAGTGGTATTCCTGACTACTCGCTGTCACATATTTTTGAACGCTTTTACTCTTTGCCGCATGCAAATggg                                                 <  2:173717/251‑1 (MQ=255)
                                                                                                                                                                          aaacaTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCACATGGGc                                                >  1:143313/5‑251 (MQ=255)
                                                                                                                                                                                   tctacATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAagc                                        >  1:4952/5‑250 (MQ=255)
                                                                                                                                                                                      acATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCg                                    >  1:151729/2‑251 (MQ=255)
                                                                                                                                                                                      acATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCg                                    <  1:159609/250‑1 (MQ=255)
                                                                                                                                                                                       cATGTTATTCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGGTACTGGAGCAGGCGCTGGGGAATTTACTGGAGAACGCCATCGATTTTACCCCCGAGAGCGGTCGCAGAACGCTAAGCGCCGAAGTGTATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTGTGAACGCTTTTACGCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCg                                    <  2:236118/250‑1 (MQ=255)
                                                                                                                                                                                             atgCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTgg                              <  1:93393/248‑1 (MQ=255)
                                                                                                                                                                                                                  gctgctGACCCGGCGTTATTGTCGCAGGCGCTGGGGAATTTACTGGCTAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCGAAGCGCCGAAGTGGATCCGGAACACGTCGCCCTTAAGGTGTTGGATTCCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAATCTGCCGTCTGGGGTTGGCGTTCGGCAGTGAGGt        <  2:276643/251‑1 (MQ=255)
                                                                                                                                                                                                                    gtcggAACCGTCGTTACTGGAGCAGGCGCTTGGGAATTTTCTGTATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATAACGCTAAGTGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCTCTGTCACGTATTTTTGAATGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTc       <  2:186163/246‑1 (MQ=255)
                                                                                                                                                                                                                     gctgAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGc     <  1:195276/251‑1 (MQ=255)
                                                                                                                                                                                                                      cagAACCGGCGTTACTGGAGCAGGCGCGGGGGAATTTACTGTAGAACTCCTTCCATTTTCCCCCCTAGAGCGGGCGCAGAACGCTAAGCGACGAAGTGGACGAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCGGTCACGTATTTTTGAACGCTTGTACTCTTTGCCGCGTGCATGTGGTCAAAAAAGCATCGTTCTGGGGGTGGCGTTCGTCAGTGTGGTCgac    <  2:249851/249‑3 (MQ=255)
                                                                                                                                                                                                                        gAACCGGCGTTACGGGAGCGGGCGCTTGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCGGTCACGGATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATTGGCAAAAAAGCAGCTGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCg  <  2:179516/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                             
CGTTTTACTGACAACATTCTGACGCAAAATGCGCGTATGCAGGCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCG  >  NC_000913/4637569‑4638035

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