Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,608,544 A→G 7.2% L160S (TTG→TCG)  yneE ← bestrophin family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,608,5440AG 7.2% 168.3 / 3.2 59L160S (TTG→TCG) yneEbestrophin family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (26/27);  new base G (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

AGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAA  >  NC_000913/1608401‑1608686
                                                                                                                                               |                                                                                                                                              
aGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATa                                                                                                                                           >  6:44774/1‑149 (MQ=255)
    aaTTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGt                                                                                                                                       <  1:33717/149‑1 (MQ=255)
         gCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATCGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCt                                                                                                                                  <  3:116760/149‑1 (MQ=255)
           gTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAACATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTg                                                                                                                                 >  8:65569/1‑148 (MQ=255)
            tataGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCCACCATTCTCCCATTATTAACAAGATACGGTTAGCTg                                                                                                                                 >  1:79439/1‑147 (MQ=255)
            tataGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGga                                                                                                                               <  6:82190/149‑1 (MQ=255)
              taGGCAATGCGCTCGCATCCTGCAAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAACATACGGTTATCTGga                                                                                                                               >  4:149168/1‑147 (MQ=255)
                 gCAATGCGCTCGCATCCTGCCAGGCCCGCTGAAATATCATTAAGCAGATCGTTGAGGCTAATAAACAGGAGATCTGAAAGCTGTACATTGCGGCGATGAACCGCCAACCATGATCCCATTATTAACAAGATCCGGTTAGCTGGAGAGt                                                                                                                           >  7:174398/1‑148 (MQ=255)
                 gCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACTAGATTCGGTTAGCTGGAGAGtt                                                                                                                          >  1:207084/1‑149 (MQ=255)
                 gCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGtt                                                                                                                          >  8:156474/1‑149 (MQ=255)
                     tGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGaa                                                                                                                      <  7:156474/149‑1 (MQ=255)
                       cgcTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAg                                                                                                                     <  7:181707/148‑1 (MQ=255)
                       cgcTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGc                                                                                                                    <  5:44774/149‑1 (MQ=255)
                          tCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAg                                                                                                                 >  2:93657/1‑149 (MQ=255)
                                 cTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACgg                                                                                                                                        >  2:239534/1‑119 (MQ=255)
                                 cTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACgg                                                                                                                                        <  1:239534/119‑1 (MQ=255)
                                     cAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACGAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTg                                                                                                        >  5:154577/1‑147 (MQ=255)
                                     cAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACGAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTg                                                                                                        >  5:157954/1‑147 (MQ=255)
                                      cggACCGCTGAAATATAATTAAGACGATCGTTGATGCTAATAAACAGGATATATGAAAGCTGTCAATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGGTAGCTGGAGAGTTCGAAGCCAGTACACGCTTAAg                                                                                                     <  8:174398/148‑1 (MQ=255)
                                             cTGAAACATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAATCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCa                                                                                              <  7:65569/149‑1 (MQ=255)
                                                  atatCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACGAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCtt                                                                                         <  6:154577/149‑1 (MQ=255)
                                                  atatCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACGAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCtt                                                                                         <  6:157954/149‑1 (MQ=255)
                                                   tatCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCtt                                                                                         >  4:79463/1‑148 (MQ=255)
                                                     tCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTgag                                                                                      <  1:93657/149‑1 (MQ=255)
                                                             ccGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAg                                                                              >  2:50499/1‑149 (MQ=255)
                                                               gATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGc                                                                             >  8:220906/1‑148 (MQ=255)
                                                                    ttGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCAc                                                                       <  8:16378/149‑1 (MQ=255)
                                                                      gAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACtt                                                                     >  7:201322/1‑149 (MQ=255)
                                                                           taataaACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGcc                                                                >  5:185760/1‑149 (MQ=255)
                                                                            aataaACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCt                                                               >  7:153359/1‑149 (MQ=255)
                                                                            aagaaaCAGGATATCTGAAAGCTGTCCATTGCGGCGCTGCACCGCCAACCGTTCTCCCATTTTTAACCAGATACGGTGAGGTGGAGAGTTCGAAGCGAGGACCCGCGGATGATCTTCAGTCTGGAGTTAATGAGCCAGCACTTCCGACt                                                               >  6:49930/4‑149 (MQ=255)
                                                                             ataaACAGGATATCTGAAAGCTGTGGATGGCGGCGCTTAAGGGCCAACCATTCTCCCATTATTACCAAGATACGGTTAGCTGGAGCGTTCGAAGCCAGTACACGCTGTAGATCTTCAGTCTTGAGATTATGAGCCAGCACTTCCGCCt                                                               <  6:148779/148‑1 (MQ=255)
                                                                                aaCAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCa                                                                           >  6:156297/1‑133 (MQ=255)
                                                                                aaCAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCa                                                                           <  5:156297/133‑1 (MQ=255)
                                                                                aaCAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCAc                                                                       >  8:32726/1‑137 (MQ=255)
                                                                                aaCAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCAc                                                                       <  7:32726/137‑1 (MQ=255)
                                                                                  cAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCt                                                         >  3:16664/1‑149 (MQ=255)
                                                                                      atatCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGttt                                                     <  3:79463/149‑1 (MQ=255)
                                                                                         tCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGtt                                                      <  7:69920/145‑1 (MQ=255)
                                                                                         tCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGtt                                                      >  8:69920/1‑145 (MQ=255)
                                                                                         tCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCg                                                  >  6:120037/1‑149 (MQ=255)
                                                                                           tGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGc                                                                                                                                   <  6:192321/66‑1 (MQ=255)
                                                                                           tGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGt                                                 >  2:22071/1‑148 (MQ=255)
                                                                                            gAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTaa                                               >  1:85789/1‑149 (MQ=255)
                                                                                             aaaGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTaa                                               <  8:153359/148‑1 (MQ=255)
                                                                                                gCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGATGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTc                                           <  7:60734/149‑1 (MQ=255)
                                                                                                  tGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCAt                                         >  1:155178/1‑149 (MQ=255)
                                                                                                    tCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCAt                                         <  7:86891/147‑1 (MQ=255)
                                                                                                    tCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGc                                       <  2:85789/149‑1 (MQ=255)
                                                                                                       aTTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGt                                     >  4:61820/1‑148 (MQ=255)
                                                                                                               gcTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTGCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTgg                            <  2:203969/149‑1 (MQ=255)
                                                                                                                      cGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAg                     >  1:7373/1‑149 (MQ=255)
                                                                                                                           aCCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGAtt                <  8:201322/149‑1 (MQ=255)
                                                                                                                             cATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTg              <  7:220906/149‑1 (MQ=255)
                                                                                                                                    ccATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGt                              <  8:11078/126‑1 (MQ=255)
                                                                                                                                    ccATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGt                              >  7:11078/1‑126 (MQ=255)
                                                                                                                                    ccATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCgg       <  2:79439/149‑1 (MQ=255)
                                                                                                                                       ttattaACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTTATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGc     >  1:178750/1‑148 (MQ=255)
                                                                                                                                          taaaCAAGATACGGTTAGCTGGATAGTTCTAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCATCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGTCGAAAGCGATTTGCAGCCGGGCaaa  <  3:149168/146‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
AGGGAATTGGCGTATAGGCAATGCGCTCGCATCCTGCCAGGACCGCTGAAATATCATTAAGCCGATCGTTGAGGCTAATAAACAGGATATCTGAAAGCTGTCCATTGCGGCGCTGAACCGCCAACCATTCTCCCATTATTAACAAGATACGGTTAGCTGGAGAGTTCGAAGCCAGTACACGCTGAAGATCTTCAGTCTTGAGATAATGAGCCAGCACTTCCGCCTGTGGCTGTTTGCGTAATGTCATGCGTAAACAGTGGGCGAAAGCGATTTGCAGCCGGGCAAA  >  NC_000913/1608401‑1608686

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