Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,074,485 C→T 16.1% A170A (GCG→GCA tktA ← transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,074,4850CT16.1% 136.4 / 18.4 56A170A (GCG→GCAtktAtransketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (26/21);  new base T (3/6);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

ATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTT  >  CP009273/3074338‑3074631
                                                                                                                                                   |                                                                                                                                                  
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                    >  2:11163/1‑150 (MQ=255)
                                                                                                                              tttACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAAt                    <  1:70251/150‑1 (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)
                                                                                                                                                aCCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGtt  >  2:37575/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTTTTTCTGCAATCGCCATACCGACTGCGTT  >  CP009273/3074338‑3074631

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: