Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,126,647 +C intergenic (‑101/+48) tdcG‑II ← / ← gcvP‑I L‑serine dehydratase/glycine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,126,6451.C88.2% 148.7 / 14.7 51intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (6/0);  new base C (17/28);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 5.61e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCG  >  NC_002947/1126497‑1126780
                                                                                                                                                            |                                                                                                                                         
cTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCCCGAATGCGGAGGTTTAGCCACCATCGCAATAGCGGGTaaa                                                                                                                                            >  2:89790/1‑149 (MQ=255)
   gAGGGGCCGATGCCGATCTTGAACAGGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcc                                                                                                                                         >  1:4764/1‑149 (MQ=255)
   gAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTCCCGCGAATGCGGAGGTATAGCCACCATAGCAATAGCGGGTAAAccc                                                                                                                                        >  4:73372/1‑148 (MQ=255)
   gAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGGTTTCCCCCCGATGGGGCGGGTGAGCCACCAGCGCGTTAGCGGGGGAAccc                                                                                                                                        >  4:12995/1‑148 (MQ=255)
     ggggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcccgc                                                                                                                                      >  1:1988/1‑146 (MQ=255)
      gggCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCGAATGCGGAGGGTTAGGCCCCAAGGCATTAGCGGGTAAAcccccg                                                                                                                                     >  6:288594/1‑145 (MQ=255)
         ccgatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCt                                                                                                                                   <  5:288594/149‑1 (MQ=255)
           gatgccgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTac                                                                                                                                 <  2:36274/149‑1 (MQ=255)
                cgatCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                             >  1:265347/1‑147 (MQ=255)
                       gatcAGGTCGTAGACGGCCAGTGACATGGGCACGTCCTGCCAGGTTGTGGGGGTGAGATTGACCCCTGTAGGCGCGGTTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTGGCGGGTAAACCCGCTCCTACACAGAGCCCg                                                                                                                     <  8:33726/144‑1 (MQ=25)
                        aaCAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTCCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCa                                                                                                                  >  6:397656/1‑149 (MQ=255)
                        aaCAGGGCGAAGACGCTCAGTGACATGGGGACCTCCGGCCCGGTTGTTGGGGGCACATTGAGCCCTGTAGGCGCGGGGTTACCCGCGAAGGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCa                                                                                                                  <  3:12995/149‑1 (MQ=37)
                          cAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTCCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTTAAACCGCTGCTACACATAGCCCGCCAcc                                                                                                                >  7:81597/1‑149 (MQ=255)
                              gTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGAGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAcc                                                                                                                                         <  1:138222/122‑1 (MQ=255)
                               tCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAAACGCCAATTCGGTCGGTTTTGCCACAAAGGAATAGCCGGTTAAAACGCACCCACACAAGACCCCcccacccag                                                                                                             >  7:319670/1‑148 (MQ=255)
                                gcaaGACGCTCACTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCGGCGAATGCGGAGGTATAGCCCCCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGc                                                                                                           <  5:40988/147‑1 (MQ=38)
                                cGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTac                                                                                                                                 <  1:216104/129‑1 (MQ=255)
                                       gCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCCCCATCGCATTAGCGGGTAAACCCGCTCCTACACAGAGCCCGCCCCCCAAGCCCGGTgg                                                                                                    >  6:292916/1‑149 (MQ=38)
                                         tCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCCCGAATGCGGAGCTATAACCACCAACGCCATAACCGGGAAACCCGCTACAACACACACCCCCCCCCCCCAGCCAACTGGCg                                                                                                  >  2:298136/1‑149 (MQ=255)
                                              gACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAg                                                                                             <  3:89303/149‑1 (MQ=255)
                                                 aTGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTCCCGCGAATGGGGGGGGATTGGCACCATCGCCTTGGCGGGGGACCCCGGTGCGAGACCGAGCCCGCCCCCCCCGCCCGGGGGGGGGGTgggg                                                                                          >  8:274130/1‑148 (MQ=255)
                                                  tgtgGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCt                                                                                                                                   <  7:184131/109‑1 (MQ=255)
                                                       aCCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg                                                                                    <  2:265347/149‑1 (MQ=255)
                                                       aCCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg                                                                                    <  2:1988/149‑1 (MQ=255)
                                                       aCCTCCTGCCAGGTTGGTGTGTGTACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAg                                                                                    <  8:277667/149‑1 (MQ=38)
                                                         ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTTACCCCCAAATCGGCGGGATAGCCACCATCGCCTTAACGGGTTAACCCCCTCCTAAACACAGCCCCCCACCCAACCAGGTTGCGGGGTGGGCTATAAGGCg                                                                                 >  7:125768/1‑149 (MQ=255)
                                                         ctcctGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGc                                                                                  <  3:363147/149‑1 (MQ=255)
                                                               ccAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTcc                                                                              <  3:73372/147‑1 (MQ=255)
                                                               ccAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTcc                                                                              <  5:397656/147‑1 (MQ=255)
                                                                  ggTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTAACCGCCAATGAGGAGGTATAGCCACCCTAGCATTACCCGGTAAAACCCCACCTACACAGAGCCCCCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGa                                                                           >  2:99163/1‑147 (MQ=255)
                                                                      gttgGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGct                                                                      <  5:292916/148‑1 (MQ=255)
                                                                       ttgGGTGGACATTGAGCCCTGTAGGCGCGGTTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctc                                                                     <  7:121120/148‑1 (MQ=255)
                                                                                cATTGAGCCCTTTAGGCGCGGGTTTACCCGCGAATGCGGAGGTAAAGCCACAATCGAATTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc                                                            <  8:319670/148‑1 (MQ=38)
                                                                                      gCCGTGTTGGCGCGGGTTTACCGGCGAATGCGGAGGTTTAGCCACCATGGCATTCGGGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGGAg                                                      <  3:129393/148‑1 (MQ=38)
                                                                                            tAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGca                                                 <  2:4764/147‑1 (MQ=255)
                                                                                               gagcgGGTTTACCCGCGAATGCGATGGTGGCTATACCTCCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCtctc                                                                   >  2:99086/3‑126 (MQ=18)
                                                                                               gagcgGGTTTACCCGCGAATGCGATGGTGGCTATACCTCCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCtctc                                                                   <  1:99086/124‑1 (MQ=18)
                                                                                                  ggggTTTACCCGCTAATGCGGAGGTATAGCCACCATCGAATTCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCAAAGCCTGGTGGCGGGTAGGACTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGt                                         <  1:298136/148‑1 (MQ=255)
                                                                                                  cgGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct                                                                    <  2:3581/122‑1 (MQ=255)
                                                                                                  cgGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctct                                                                    >  1:3581/1‑122 (MQ=255)
                                                                                                     gTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCCCCCCAGGTTGc                                      >  7:196692/1‑149 (MQ=255)
                                                                                                       tttCCCCTCGAATGCGGAGGTATAGCCACCATCGCATTCGCTGTTAAACCCGCTCCTACACAGATCCCGCCACCCAAGCCAGGTGTCGGGGAGGTCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCg                                     <  5:370744/149‑1 (MQ=255)
                                                                                                          aCCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGAc                                                            <  3:132186/123‑1 (MQ=255)
                                                                                                          aCCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGAGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCACGAACGAc                                                            >  4:132186/1‑123 (MQ=38)
                                                                                                            ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGTGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                >  2:291167/1‑149 (MQ=255)
                                                                                                             cgcgAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGACGcc                               >  2:260678/1‑149 (MQ=255)
                                                                                                                gAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGaa                            >  1:352782/1‑149 (MQ=255)
                                                                                                                     cGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgatcg                                                              >  7:272455/1‑110 (MQ=255)
                                                                                                                     cGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcgatcg                                                              <  8:272455/110‑1 (MQ=255)
                                                                                                                      ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCACCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGt                       >  1:97434/1‑149 (MQ=255)
                                                                                                                      ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGtt                      <  8:81597/149‑1 (MQ=255)
                                                                                                                      ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACAAAGAGCCCGCCAACCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACTGGCAGGCGCACACCAGGTTGCGGTCTCCGAACACGtt                      >  6:262816/1‑149 (MQ=255)
                                                                                                                           gTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCg                  <  5:175683/149‑1 (MQ=255)
                                                                                                                                 ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg            >  8:198217/1‑149 (MQ=255)
                                                                                                                                 ccaccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGAACGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAACCgg            >  6:19485/1‑149 (MQ=255)
                                                                                                                                  caccaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGc           <  7:274130/149‑1 (MQ=255)
                                                                                                                                       tCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGAcc      >  1:70375/1‑149 (MQ=255)
                                                                                                                                            aTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcg  >  1:201019/1‑149 (MQ=255)
                                                                                                                                                            |                                                                                                                                         
CTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCG  >  NC_002947/1126497‑1126780

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: