Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,373,506 A→T 30.0% N90K (AAT→AAA pbuG ← Guanine/hypoxanthine permease PbuG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,373,5060AT30.0% 82.0 / 48.7 63N90K (AAT→AAApbuGGuanine/hypoxanthine permease PbuG
Reads supporting (aligned to +/- strand):  ref base A (20/24);  new base T (11/8);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.96e-01

GCATTAATAATAACTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTA  >  USA300TCH1516_ALE/2373375‑2373644
                                                                                                                                   |                                                                                                                                           
gCATTAATAATAACTTCCCTAAAAAGTCTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATccc                                                                                                                                    >  2:760554/1‑141 (MQ=255)
       aaataaCTGCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTg                                                                                                                             <  1:387546/140‑1 (MQ=255)
         ataaCTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATATATCCACTTGTTAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGt                                                                                                                           <  1:573626/141‑1 (MQ=255)
         ataaCTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGt                                                                                                                           <  2:604687/141‑1 (MQ=255)
         ataaCCTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGt                                                                                                                           <  1:313598/141‑1 (MQ=255)
          taaCTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTg                                                                                                                          >  2:216651/1‑141 (MQ=255)
           aaCTTCCCTAAACCATGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGc                                                                                                                         >  2:536484/1‑141 (MQ=255)
                  cTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGACTGCTTTTAATCCCATACCTGGTGCTAACGCa                                                                                                                 >  1:264708/1‑141 (MQ=255)
                    aaaCCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAAt                                                                                                               >  2:749874/1‑141 (MQ=255)
                     aaCCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCCACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAAtt                                                                                                              >  1:139207/1‑141 (MQ=255)
                      aCCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTg                                                                                                             >  2:15338/1‑141 (MQ=255)
                      aCCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTg                                                                                                             <  1:312737/141‑1 (MQ=255)
                       cccTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTCCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTgc                                                                                                            <  2:694480/141‑2 (MQ=255)
                        ccTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGa                                                                                                           <  1:9846/141‑1 (MQ=255)
                        ccTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATACTTTTAATCCCATACCTGGTGCTAACGCAATTGGa                                                                                                           <  2:217618/141‑1 (MQ=255)
                             gACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGagcttt                                                                                                      >  1:487912/1‑136 (MQ=255)
                              aCGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTa                                                                                                     >  2:670072/1‑141 (MQ=255)
                              aCGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAAATTGCCAAGGAATACCCATGGTTAACAAAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTggcgacgga                                                                                                     >  2:118795/1‑134 (MQ=255)
                               cGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAg                                                                                                    >  2:349437/1‑141 (MQ=255)
                                 gTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCt                                                                                                  <  1:544438/141‑1 (MQ=255)
                                    aaTATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCGAACGAAATTGGATATTTAGCTAtt                                                                                               >  2:750108/1‑141 (MQ=255)
                                      tatCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAg                                                                                             <  1:36666/141‑1 (MQ=255)
                                      tatCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAg                                                                                             <  2:785624/141‑1 (MQ=255)
                                             aaGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATg                                                                                      <  1:642526/141‑1 (MQ=255)
                                               gAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGaa                                                                                    >  1:277289/1‑141 (MQ=255)
                                                    taCCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTGGCCCCATGAATAGCg                                                                               >  1:123969/1‑141 (MQ=255)
                                                    taCCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCg                                                                               >  1:123915/1‑141 (MQ=255)
                                                        tGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGcc                                                                           <  1:503784/141‑1 (MQ=255)
                                                        tGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGcc                                                                           <  1:443516/141‑1 (MQ=255)
                                                        tGACAATACAACGGCTTGTAAACCAACTTGCCAAGGAATACCTCTTGTTACCACAACAGTAAATGCAAAGAATGCTTTTAATCCCATTCCTGTTGCTAACTCAATTGGATATTTAGCTATTATCCCCATGAATAGCGAGcc                                                                           <  2:141350/141‑1 (MQ=255)
                                                         gagaATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCt                                                                          <  1:270242/141‑1 (MQ=255)
                                                          agaATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTa                                                                         <  2:282114/141‑1 (MQ=255)
                                                           gaATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTAc                                                                        >  2:188293/1‑141 (MQ=255)
                                                             aTAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACaa                                                                      <  2:85842/141‑1 (MQ=255)
                                                                      cTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTa                                                             <  1:135593/141‑1 (MQ=255)
                                                                       tGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTaa                                                            <  1:475659/141‑1 (MQ=255)
                                                                        gTTATACCAACTTGCAAAGGAATACCCATGGTTAACACCACAGTAACTGCAAAGAATGCATTCAATCCCATACCTGATGCCAACGCAATTGGATGTTTAGCTATTAGCCCCATGAATAGCGAGCCTCCAAATGCTGCCAAt                                                           >  2:220159/1‑141 (MQ=255)
                                                                                  cTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGACTACAAGTTCTGCTAATGCAGTCGCTa                                                 >  1:212707/1‑141 (MQ=255)
                                                                                    tGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCAATGAATAGCGAGCCTACAAATGGTGCTAATGCAGTCGCGACa                                               >  2:551581/1‑141 (MQ=255)
                                                                                     gCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACaa                                              >  2:711031/1‑141 (MQ=255)
                                                                                     gCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACaa                                              >  1:285384/1‑141 (MQ=255)
                                                                                             attcccATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCa                                      <  1:743743/138‑1 (MQ=255)
                                                                                             aTACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCa                                      <  1:412127/141‑1 (MQ=255)
                                                                                                   aTGGTTAACACAACAGTAAATGTAAAGAAGGCATTTAATCCCTACTTTGTACCAAGCTCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTgg                                <  1:771188/141‑1 (MQ=255)
                                                                                                   aTGGTTAACACAACAGTAAACGCAAACAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTgg                                <  1:457185/141‑1 (MQ=255)
                                                                                                    tGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGt                               <  1:760554/141‑1 (MQ=255)
                                                                                                     ggTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATACCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTc                              <  1:223000/141‑1 (MQ=255)
                                                                                                        tAACACAACAGTAAATGCAATGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAt                           <  1:15338/141‑1 (MQ=255)
                                                                                                        tAACACAACAGTAAATGCAAAGAATGCTTTTTATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAt                           >  1:219931/1‑141 (MQ=255)
                                                                                                        tAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAt                           >  2:129882/1‑141 (MQ=255)
                                                                                                          acacAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAttt                         >  2:11382/1‑141 (MQ=255)
                                                                                                          acacAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCCACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCAttt                         >  2:671824/1‑141 (MQ=255)
                                                                                                              aacaGTAAATGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCAGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCat                     >  2:404544/1‑141 (MQ=255)
                                                                                                                   tAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAtctt                <  1:418058/141‑1 (MQ=255)
                                                                                                                   aaaaTGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAtctt                <  2:568526/140‑1 (MQ=255)
                                                                                                                    aaaTGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCACCATGAATAGCGAGCCTACAAACGCTGCGAATGCAGTCGCTACAAAAATGGCACCTTGGTACATTCTCATAtcttc               >  2:278243/1‑141 (MQ=255)
                                                                                                                    aaaTGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAtcttc               >  1:677173/1‑141 (MQ=255)
                                                                                                                     caTGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATGTTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTACTATTTCAGTCGCTACAAACATGGCACCTTGGTCCATTTTCATAtgttct              <  2:479737/140‑1 (MQ=255)
                                                                                                                     aaTGCAAAGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATAtcttct              >  1:629449/1‑141 (MQ=255)
                                                                                                                       tGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTCGCTATTAGCCACATGAATAGCGAGCCTACCAATGCTGCTAATGCAGTCGCTGCAGAAATGGCACCTTGGTCCATTTTCATATATTCTGa            >  1:194294/1‑141 (MQ=255)
                                                                                                                            aGAATGCTTTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGc       >  1:517771/1‑141 (MQ=255)
                                                                                                                            aGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGc       <  1:87968/141‑1 (MQ=255)
                                                                                                                                 gCATCTAATCCAAGACCTCGTGCAAACGCAATTAGATATTTAGCTATTAGCCCCAGGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCATCTTGGTCCATTTTCATATCTTCTGATACGCCTTTa  <  1:430065/141‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                           
GCATTAATAATAACTTCCCTAAACCCTGTGACGGTTAATATCGCAAAGAATATACCTGAGAATAAAACGCCTGTTAAACCAACTTGCCAAGGAATACCCATGGTTAACACAACAGTAAATGCAAAGAATGCATTTAATCCCATACCTGGTGCTAACGCAATTGGATATTTAGCTATTAGCCCCATGAATAGCGAGCCTACAAATGCTGCTAATGCAGTCGCTACAAAAATGGCACCTTGGTCCATTTTCATATCTTCTGATACGCCTTTA  >  USA300TCH1516_ALE/2373375‑2373644

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