Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,367,4000AC44.6% 104.4 / 14.9 92G361G (GGT→GGGmenEo‑succinylbenzoate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (28/23);  new base C (0/41);  total (28/64)
Fisher's exact test for biased strand distribution p-value = 8.69e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAG  >  NZ_CP009273/2367314‑2367485
                                                                                      |                                                                                     
aaCTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATAccc                                                                                    <  1:2965392/90‑1 (MQ=255)
     tttttCGGCAACGGGGATCACAAACCCCTGCAAACCCGCCGGATTTGCAGAAATTACGCGCTCGCCTTCTTCCGGCTGAATCCCCTCTcc                                                                               <  2:939181/86‑1 (MQ=255)
     cTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTcc                                                                               <  2:40900/90‑1 (MQ=255)
      ttGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTcc                                                                               <  1:790833/89‑1 (MQ=255)
      ttGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTcc                                                                               >  2:790833/1‑89 (MQ=255)
      ttGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCg                                                                              >  2:1525535/1‑90 (MQ=255)
        gTCGGCAAGGGGGACGTTAAACACCCGCAAAATCGCGGGATGTGCAGCAATTACCCGCTGGTCTTCTTCCGGCTGAATACCCTCTCCGcc                                                                            <  2:2360287/90‑1 (MQ=255)
        gTCGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTCCAGCAATTACGCGCTCGCCTTCTTCCGGCTGAATCCCCTCTCCGcc                                                                            <  2:169204/90‑1 (MQ=255)
          cGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCAc                                                                          <  2:1584149/90‑1 (MQ=255)
          cGGCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCAc                                                                          <  2:675518/90‑1 (MQ=255)
          cGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCAc                                                                          <  2:1232771/90‑1 (MQ=255)
          cGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCAc                                                                          >  2:3088905/1‑90 (MQ=255)
          cGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCAc                                                                          >  1:2229205/1‑90 (MQ=255)
            gCAACGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTg                                                                        <  2:1043976/90‑1 (MQ=255)
            gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTg                                                                        <  2:3381784/90‑1 (MQ=255)
            gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg                                                                        <  1:2050489/90‑1 (MQ=255)
            gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg                                                                        >  2:694075/1‑90 (MQ=255)
            gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg                                                                        <  2:2780535/90‑1 (MQ=255)
            gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg                                                                        >  2:859254/1‑90 (MQ=255)
            gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg                                                                        <  1:466299/90‑1 (MQ=255)
            gCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCCCTCCGCCACTg                                                                        <  2:3306362/90‑1 (MQ=255)
             aaaCGGGGACGATAAACACCTGCAAAACCGCAGTATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTga                                                                       <  1:2696908/89‑1 (MQ=255)
              aaCGGGGACGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTgaa                                                                      <  1:1700910/90‑1 (MQ=255)
              aaCGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTgaa                                                                      >  2:3104145/1‑90 (MQ=255)
                     aCGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGa                                                               <  1:3055015/90‑1 (MQ=255)
                      cGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGAt                                                              <  1:363798/90‑1 (MQ=255)
                      cGATAAACCCCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCCCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt                                                              <  1:1123884/90‑1 (MQ=255)
                      cGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt                                                              <  2:937763/90‑1 (MQ=255)
                      cGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGAt                                                              <  1:1004922/90‑1 (MQ=255)
                      cGATAAACACCTGCAAAACCGCAGGATGTACAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGAt                                                              <  1:1685511/90‑1 (MQ=255)
                        aTAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTg                                                            >  2:1340081/1‑90 (MQ=255)
                               ccTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAAc                                                     <  2:1492139/90‑1 (MQ=255)
                                ctccAAACCCGCAGGATGTGCAGCAATTCCGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACg                                                    <  2:872537/87‑1 (MQ=255)
                                cTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCAc                                                                          <  1:2477470/68‑1 (MQ=255)
                                cTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCAc                                                                          >  2:2477470/1‑68 (MQ=255)
                                cTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACg                                                    >  1:1081137/1‑90 (MQ=255)
                                cTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACg                                                    >  1:2720304/1‑90 (MQ=255)
                                 tGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGt                                                   >  1:1841568/1‑90 (MQ=255)
                                  gCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTc                                                  >  1:284591/1‑90 (MQ=255)
                                   cAAACCCGCAGGATGTGCAGCAATTACGCGCTCGCCTTCTTCCGGCTAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc                                                 <  2:2034395/90‑1 (MQ=255)
                                   cAAACCCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc                                                 <  1:2168631/90‑1 (MQ=255)
                                   cAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc                                                 <  2:2886922/90‑1 (MQ=255)
                                   cAAAACCGCAGGATGTGCAGCAATTACGCCCTCGCCTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTcc                                                 <  2:1755791/90‑1 (MQ=255)
                                   cAAAACCGAAGGATGTGCAGCAATTGCGCTCTGGTCTTTTTCCGGCTGAATCCCCTCCCCGCCACTGAAGAATAGATTGTCTAAACGTcc                                                 <  2:966025/90‑1 (MQ=255)
                                         cGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAAt                                           <  1:803662/90‑1 (MQ=255)
                                         cGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAAt                                           >  2:1459947/1‑90 (MQ=255)
                                             ggATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTc                                       <  2:2226756/90‑1 (MQ=255)
                                             ggATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTc                                       <  1:1402977/90‑1 (MQ=255)
                                                tgtgCAGCAATTCCGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGc                                    <  1:371913/90‑1 (MQ=255)
                                                 gtgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCt                                   <  2:1237192/90‑1 (MQ=255)
                                                 gtgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCt                                   >  2:2097379/1‑90 (MQ=255)
                                                 gtgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCt                                   >  1:1858809/1‑90 (MQ=255)
                                                 gtgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCt                                   >  1:1465211/1‑90 (MQ=255)
                                                  tgCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCtt                                  <  2:1623143/90‑1 (MQ=255)
                                                     agcaATTCCGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGcc                               <  1:99761/90‑1 (MQ=255)
                                                       caATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCAt                             <  2:864972/90‑1 (MQ=255)
                                                       caATTACGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCAt                             <  2:2686291/90‑1 (MQ=255)
                                                       caATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCAc                                                                          >  1:2969227/1‑45 (MQ=255)
                                                       caATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCAc                                                                          <  2:2969227/45‑1 (MQ=255)
                                                       caATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCAt                             <  2:1081137/90‑1 (MQ=255)
                                                        aaTTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTc                                       <  1:2508425/79‑1 (MQ=255)
                                                        aaTTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTc                                       >  2:2508425/1‑79 (MQ=255)
                                                         aTTACGCGCTCGATTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTa                           <  2:987411/90‑1 (MQ=255)
                                                         aTTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTa                           >  2:701139/1‑90 (MQ=255)
                                                            aCGCGCTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTg                                                                        <  1:1229444/42‑1 (MQ=255)
                                                            aCGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTg                                                                        >  2:1229444/1‑42 (MQ=255)
                                                            aCGCCCTCGCTTTTTTCCGGCTGAATCCCCTCCCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGc                        <  2:936795/90‑1 (MQ=255)
                                                               cgcTCGATTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtc                     <  2:347730/90‑1 (MQ=255)
                                                               cgcTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtc                     <  1:3278644/90‑1 (MQ=255)
                                                               cgcTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtc                     <  2:258174/90‑1 (MQ=255)
                                                               cgcTCGACTTCTTCCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtc                     <  1:2693403/90‑1 (MQ=255)
                                                               cgcTCGACTTCTTCCGGCTGAATCCCCTCCCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtc                     <  2:1149208/90‑1 (MQ=255)
                                                                gcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtct                    >  1:1513221/1‑90 (MQ=255)
                                                                gcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtct                    >  2:2272782/1‑90 (MQ=255)
                                                                gcTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCAtct                    >  2:488303/1‑90 (MQ=255)
                                                                   cGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCAc                 <  1:3088905/90‑1 (MQ=255)
                                                                   cGACTTCTTCCGGCTAAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCAc                 <  2:3001467/90‑1 (MQ=255)
                                                                      cttcttCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACcgc              <  1:2675129/90‑1 (MQ=255)
                                                                        tcttcCGGCTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTAACCGCGa            <  1:859254/90‑1 (MQ=255)
                                                                        tcttcCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTAACCGCGa            <  1:694075/90‑1 (MQ=255)
                                                                                cTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgt    <  2:3238020/90‑1 (MQ=255)
                                                                                cTGAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgt    <  1:1340081/90‑1 (MQ=255)
                                                                                cTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgt    <  1:1532886/90‑1 (MQ=255)
                                                                                  gAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag  <  2:52611/90‑1 (MQ=255)
                                                                                  gAATCCCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag  <  1:411469/90‑1 (MQ=255)
                                                                                  gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag  >  2:2047825/1‑90 (MQ=255)
                                                                                  gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag  >  2:2122355/1‑90 (MQ=255)
                                                                                  gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag  <  1:1459947/90‑1 (MQ=255)
                                                                                  gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag  <  2:3186087/90‑1 (MQ=255)
                                                                                  gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCgcgtag  >  1:3186087/1‑90 (MQ=255)
                                                                                  gAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCACGATCGCgcgtag  >  2:2622938/1‑90 (MQ=255)
                                                                                  aaaTACCCTCCCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCTCCGCGATCGCgcgtag  <  2:1932355/89‑1 (MQ=255)
                                                                                      |                                                                                     
AACTCCTTGTCGGCAACGGGGACGATAAACACCTGCAAAACCGCAGGATGTGCAGCAATTACGCGCTCGACTTCTTCCGGCTGAATACCCTCTCCGCCACTGAAGAATAGATTGTCTAAACGTCCGACAATGGTCAGCTTGCCATTATGCATCTCACCGCGATCGCGCGTAG  >  NZ_CP009273/2367314‑2367485

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