Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,074,481 T→A 15.8% T172S (ACG→TCG)  tktA ← transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,074,4810TA15.8% 141.8 / 15.9 56T172S (ACG→TCG) tktAtransketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base T (25/22);  new base A (3/6);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01

CGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTG  >  CP009273/3074334‑3074627
                                                                                                                                                   |                                                                                                                                                  
cgTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTAc                                                                                                                                                  <  1:60006/150‑1 (MQ=255)
    aTGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCg                                                                                                                                               <  1:99261/149‑1 (MQ=255)
     tGACCGTCAATGGCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGcc                                                                                                                                             <  1:50132/150‑1 (MQ=255)
          gTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCagaga                                                                                                                                        <  1:66243/150‑1 (MQ=255)
          gTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCagaga                                                                                                                                        <  1:27907/150‑1 (MQ=255)
          gTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCaga                                                                                                                                          >  2:185967/1‑148 (MQ=255)
                gTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCGAGCTTCAGCGTACCCGCCAGAGAGCAAAc                                                                                                                                  >  2:42808/1‑150 (MQ=255)
                     gAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGt                                                                                                                             >  1:32942/1‑150 (MQ=255)
                            gTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGAt                                                                                                                      <  2:180956/150‑1 (MQ=255)
                               ccAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGcc                                                                                                                   >  2:618/1‑150 (MQ=255)
                                cAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCt                                                                                                                  <  2:182159/150‑1 (MQ=255)
                                  gCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTCCCCAGCTTCAGCGTCCCCGCCAGAGAGCCACCTTCGTGGGAGATGCCTTc                                                                                                                <  2:46248/150‑1 (MQ=255)
                                  gCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTc                                                                                                                >  1:190504/1‑150 (MQ=255)
                                    cGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCca                                                                                                              >  2:77725/1‑150 (MQ=255)
                                      tAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCcatc                                                                                                            <  1:124506/150‑1 (MQ=255)
                                      tAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCcat                                                                                                             >  1:81343/1‑149 (MQ=255)
                                              gAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGcc                                                                                                    >  2:83282/1‑150 (MQ=255)
                                                     aTTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGcc                                                                                              >  1:50836/1‑149 (MQ=255)
                                                     aTTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGcc                                                                                              >  1:176723/1‑149 (MQ=255)
                                                         cGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATg                                                                                          <  1:77725/149‑1 (MQ=255)
                                                                gtcgGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCg                                                                                    >  1:132846/1‑148 (MQ=255)
                                                                gtcgGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTa                                                                                  <  2:104852/150‑1 (MQ=255)
                                                                     tGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGgtgt                                                                             <  1:185967/150‑1 (MQ=255)
                                                                            gCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTc                                                                      >  1:152882/1‑150 (MQ=255)
                                                                             ccTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTc                                                                      >  2:49841/1‑149 (MQ=255)
                                                                                 cAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGAc                                                                   <  2:81343/148‑1 (MQ=255)
                                                                                 cAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACa                                                                  <  1:170380/149‑1 (MQ=255)
                                                                                  aaCGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTa                                                                            >  2:115508/1‑138 (MQ=255)
                                                                                   aCGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAg                                                                           <  1:115508/138‑1 (MQ=255)
                                                                                   aCGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATg                                                               <  1:83282/150‑1 (MQ=255)
                                                                                        aCCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGt                                                           <  2:32942/149‑1 (MQ=255)
                                                                                            tCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGccc                                                      <  2:98425/150‑1 (MQ=255)
                                                                                            tCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGccc                                                      >  1:69388/1‑150 (MQ=255)
                                                                                             cGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCg                                                     >  2:60793/1‑150 (MQ=255)
                                                                                                tAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGAc                                                  >  2:151025/1‑150 (MQ=255)
                                                                                                 aGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACg                                                 >  2:73896/1‑150 (MQ=255)
                                                                                                            tGTCGTCGTAGACGGCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTggc                                      >  1:75445/1‑148 (MQ=255)
                                                                                                             gtcgtcGTAGACGGCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTggc                                      <  1:194233/149‑3 (MQ=255)
                                                                                                             gTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTgcgc                                     >  2:157359/1‑150 (MQ=255)
                                                                                                              tcgtCGTAGACGGCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTggccc                                    <  2:165289/147‑5 (MQ=255)
                                                                                                                gtCGTAGACGGCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTggcccg                                   <  1:80840/148‑6 (MQ=255)
                                                                                                                   gTAGACGGCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAg                               >  2:176703/1‑150 (MQ=255)
                                                                                                                    tagacggCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGc                              <  1:139838/144‑1 (MQ=255)
                                                                                                                       acggCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCgcg                           <  2:94201/147‑3 (MQ=255)
                                                                                                                        aTGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGttt                          >  2:99941/1‑150 (MQ=255)
                                                                                                                         ggCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCcgcgcgcg                         >  1:109079/2‑146 (MQ=255)
                                                                                                                          gCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGttttt                        >  2:155328/1‑150 (MQ=255)
                                                                                                                            aaTCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAGCgcgcgttc                      <  2:39535/150‑8 (MQ=255)
                                                                                                                            aaTCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTc                      >  2:20401/1‑150 (MQ=255)
                                                                                                                                aGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGc                   <  2:190504/149‑1 (MQ=255)
                                                                                                                                  tttACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAAt                <  1:70251/150‑1 (MQ=255)
                                                                                                                                  tttACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAAt                >  2:11163/1‑150 (MQ=255)
                                                                                                                                    tACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATc               <  2:11936/149‑1 (MQ=255)
                                                                                                                                       ccAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCa           >  1:174197/1‑150 (MQ=255)
                                                                                                                                        cAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCa           >  2:117660/1‑149 (MQ=255)
                                                                                                                                                gCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAACCTGCGCCGCCAGCGTTTTTTCTGCCATCGCCATACCGACTg  <  2:158364/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTG  >  CP009273/3074334‑3074627

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