Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,705,915 G→A 7.4% R193C (CGC→TGC)  PP_4162 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,705,9150GA 7.4% 155.5 / 2.1 54R193C (CGC→TGC) PP_4162membrane protein
Reads supporting (aligned to +/- strand):  ref base G (27/23);  new base A (2/2);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

CCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGCACGCCGCTGCCGATCTGCGCGCTGCGCGCCAGCAACCAGTTGCC  >  NC_002947/4705771‑4706062
                                                                                                                                                |                                                                                                                                                   
ccAACCGCAGCACGATCGCCCCGGGCACGCGGGCTATCAGGAAGCCGATCAGTGCGAGCAGGGCCTGCGCGGCCGCCGTGCCGATCTCCCCGTTGTCCACCCGCTGCTCGGTAGCGGCTACCAGTTCCAGGTAATACTCGGCGCGCTc                                                                                                                                                  <  4:180561/148‑1 (MQ=255)
ccAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCg                                                                                                                                                 <  4:79367/149‑1 (MQ=255)
    ccccAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGc                                                                                                                                                <  8:241917/146‑1 (MQ=255)
    ccccAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGc                                                                                                                                                >  7:241917/1‑146 (MQ=255)
    ccccAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCcca                                                                                                                                             <  6:217055/149‑1 (MQ=255)
       cAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCccac                                                                                                                                            >  6:366885/1‑147 (MQ=255)
         gCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGc                                                                                                                                        >  1:64773/1‑149 (MQ=255)
           aCGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCg                                                                                                                                      <  8:282488/149‑1 (MQ=255)
               tCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTgc                                                                                                                                  <  1:338989/149‑1 (MQ=255)
                     cGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAg                                                                                                                            >  7:159907/1‑149 (MQ=255)
                             ccGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCAc                                                                                                                    >  2:295559/1‑149 (MQ=255)
                                  gATCAGGAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATg                                                                                                               >  6:205072/1‑149 (MQ=255)
                                  gATCAGGAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATg                                                                                                               >  5:339847/1‑149 (MQ=255)
                                    tCAGGAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATgcc                                                                                                             >  1:311746/1‑149 (MQ=255)
                                       ggAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATgccgc                                                                                                           >  8:290115/1‑148 (MQ=255)
                                         aaGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCaa                                                                                                        <  6:339847/149‑1 (MQ=255)
                                             cGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGt                                                                                                    <  3:230141/149‑1 (MQ=255)
                                                tCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGt                                                                                                 <  8:159907/149‑1 (MQ=255)
                                                      ccAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGTCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTcg                                                                                           <  7:265814/149‑1 (MQ=255)
                                                      ccAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCACTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTcg                                                                                           >  2:130705/1‑149 (MQ=255)
                                                      ccAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCACTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTcg                                                                                           >  4:245099/1‑149 (MQ=255)
                                                       cagcagGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTcgc                                                                                          <  7:95939/149‑1 (MQ=255)
                                                            gggCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTag                                                                                     <  2:311746/149‑1 (MQ=255)
                                                                   cgcgGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTagaagaag                                                                               <  3:320048/148‑1 (MQ=255)
                                                                      ggCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTc                                                                                                <  2:128134/128‑1 (MQ=255)
                                                                      ggCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTc                                                                                                >  1:128134/1‑128 (MQ=255)
                                                                        ccgccgTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCACTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAgaacacg                                                                         <  3:245099/149‑1 (MQ=255)
                                                                        ccgccgTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCACTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAgaacacg                                                                         <  1:130705/149‑1 (MQ=255)
                                                                          gccgTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAgaacacgaa                                                                       <  1:311053/149‑1 (MQ=255)
                                                                              tGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAgaacacgaacgc                                                                    >  5:146326/1‑146 (MQ=255)
                                                                                       ccccGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAgg                                                          >  5:3165/1‑149 (MQ=255)
                                                                                        cccGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAggg                                                         >  2:68348/1‑149 (MQ=255)
                                                                                              gACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCTTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGt                                                    <  1:85732/148‑1 (MQ=255)
                                                                                              gACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCa                                                            <  4:11149/140‑1 (MQ=255)
                                                                                              gACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCa                                                            >  3:11149/1‑140 (MQ=255)
                                                                                               accaccCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAgg                                                                <  2:321347/135‑1 (MQ=255)
                                                                                               accaccCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAgg                                                                >  1:321347/1‑135 (MQ=255)
                                                                                                  accCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAg                                               <  5:366885/149‑1 (MQ=255)
                                                                                                      gCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAgaac                                                                            >  8:340770/1‑116 (MQ=255)
                                                                                                      gCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAgaac                                                                            <  7:340770/116‑1 (MQ=255)
                                                                                                             ggTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGCACGCCGCTGc                                    >  5:66500/1‑149 (MQ=255)
                                                                                                              gTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGCACGCCGCTGcc                                   >  5:135847/1‑149 (MQ=255)
                                                                                                                   ggCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGCACGCCGCTGCCGATct                              <  3:178177/149‑1 (MQ=255)
                                                                                                                     cTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGCACGCCGCTGCCGATctgc                            >  6:244494/1‑149 (MQ=255)
                                                                                                                           ggtccaggtAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGCACGCCGCTGCCGATctgcgcgctg                      >  2:75658/1‑149 (MQ=255)
                                                                                                                                 ggtAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGCACGCCGCTGCCGATCtgcgcgctgcgcgc                 >  1:57236/1‑148 (MQ=255)
                                                                                                                                 ggtAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGCACGCCGCTGCCGATCTGCGCGCTGCGCGcc                <  6:66500/149‑1 (MQ=255)
                                                                                                                                  gtAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGc                                              <  7:267363/118‑1 (MQ=255)
                                                                                                                                  gtAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGc                                              >  8:267363/1‑118 (MQ=255)
                                                                                                                                       aCTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGCACGCCGCTGCCGATCTGCGCGCTGCGCGCCAGCAAc          >  3:102958/1‑149 (MQ=255)
                                                                                                                                       aCTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGCACGCCGCTGCCGATCTGCGCGCTGCGCGCCAGCAAc          >  2:190679/1‑149 (MQ=255)
                                                                                                                                         tCGGCGCGCTCGCCCACCCGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGCACGCCGCTGCCGATCTGCGCGCTGCGCGCCAGCAAc          >  5:155062/1‑147 (MQ=255)
                                                                                                                                         tCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCGCGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGCACGCCGCTGCCGATCTGCGCGCTGCGCGCCAGCAAc          >  3:165907/1‑147 (MQ=255)
                                                                                                                                         tCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGCACACGAACACCAGGCTCAGGGTCAGTTCCAGCACGCCGCTGCCGATCTGCGCGCTGCGCGCCAGCAAc          >  7:219594/1‑147 (MQ=255)
                                                                                                                                            gcgcgcTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTcc                                                 >  5:238906/1‑105 (MQ=255)
                                                                                                                                            gcgcgcTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTcc                                                 <  6:238906/105‑1 (MQ=255)
                                                                                                                                               cgcTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGCACGCCGCTGCCGATCTGCGCGCTGCGCGCCAGCAACCAGTTGcc  <  5:244494/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CCAACCCCAGCACGATCGCCCCGGGCACCCCGGCGATCAGGAAGCCGATCAGTGCCAGCAGGGCCTGCGCGGCCGCCGTGCCGATCACCCCGTTGACCACCCGCTGCACGGTAGCGGCTACCAGGTCCAGGTAATACTCGGCGCGCTCGCCCACCAGCCGGTGCAGCAAGCGCAGCACGAATGCCGCCAAGCGTGGTCCGTCGCGGTAGAAGAAGAACACGAACACCAGGCTCAGGGTCAGTTCCAGCACGCCGCTGCCGATCTGCGCGCTGCGCGCCAGCAACCAGTTGCC  >  NC_002947/4705771‑4706062

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