Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 86,245 A→C 11.7% M158R (ATG→AGG)  PP_0075 ← choline sulfate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294786,2450AC11.7% 148.0 / 11.4 59M158R (ATG→AGG) PP_0075choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base C (4/3);  total (32/27)
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.82e-01

GCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGC  >  NC_002947/86097‑86383
                                                                                                                                                    |                                                                                                                                          
gCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCAt                                                                                                                                           >  1:843122/1‑150 (MQ=255)
gCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCAt                                                                                                                                           >  2:518305/1‑150 (MQ=255)
  tCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCATGGCCCTGc                                                                                                                                         >  2:982970/1‑150 (MQ=255)
          cGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCATGGCCCTGCCcagccag                                                                                                                                 <  1:982970/150‑4 (MQ=255)
            tGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGtt                                                                                                                               >  2:142439/1‑150 (MQ=255)
                    agagTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCt                                                                                                                       >  1:1147048/1‑150 (MQ=255)
                              ggCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTg                                                                                                                              <  2:1192833/133‑1 (MQ=255)
                              ggCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTg                                                                                                                              >  1:1192833/1‑133 (MQ=255)
                               gCCCCACCGATCACGACGGCCGCCAACCCGACCAGCAACGTCGGTACCTCGAGCTCGCCTGCGTTACTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGAGTAGCTGGCCCAAGGCa                                                                                                            >  2:919992/1‑150 (MQ=255)
                                ccGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCcagcag                                                                                                                                  <  2:783816/127‑1 (MQ=255)
                                ccGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCcagcag                                                                                                                                  >  1:783816/1‑127 (MQ=255)
                                 cGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAAt                                                                                                          <  2:1147048/150‑1 (MQ=255)
                                   cACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGa                                                                                                        <  1:496203/150‑1 (MQ=255)
                                    aCCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGAc                                                                                                       <  2:688096/150‑1 (MQ=255)
                                    aCCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGAc                                                                                                       <  1:518305/150‑1 (MQ=255)
                                    aCCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGAc                                                                                                       <  2:843122/150‑1 (MQ=255)
                                           acgacgGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAg                                                                                                               >  2:827084/1‑135 (MQ=255)
                                           acgacgGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAg                                                                                                               <  1:827084/135‑1 (MQ=255)
                                           acgacgGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATggcg                                                                                                <  2:1105669/150‑1 (MQ=255)
                                                         gCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGaa                                                                                  >  1:1072029/1‑150 (MQ=255)
                                                           ccAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGTTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAg                                                                                                               <  2:1073423/119‑1 (MQ=255)
                                                           ccAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAg                                                                                                               >  1:1073423/1‑119 (MQ=255)
                                                            cagcagCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGcc                                                                               >  2:917400/1‑150 (MQ=255)
                                                              gcagcaAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGa                                                                             <  2:1072029/150‑1 (MQ=255)
                                                                  caAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGcac                                                                         <  1:115700/150‑1 (MQ=255)
                                                                    aggaaggTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCcagca                                                                                                                                   <  2:178283/90‑1 (MQ=255)
                                                                    aggaaggTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCcagca                                                                                                                                   >  1:178283/1‑90 (MQ=255)
                                                                      gaaggTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACg                                                                     >  1:364501/1‑150 (MQ=255)
                                                                          gtaggtCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAAt                                                                 >  1:754333/1‑150 (MQ=255)
                                                                              gtCGACCTCGCCTGCGTGACTGGCCAGGTCCTCCACCGTCATCAGCGCGGTGGCCTCGCGGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGga                                                             <  1:919992/150‑1 (MQ=255)
                                                                               tCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAg                                                                                 <  1:436882/129‑1 (MQ=255)
                                                                               tCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAg                                                                                 >  2:436882/1‑129 (MQ=255)
                                                                               tCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGACATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGAACCACGGAATGGgag                                                            >  2:29386/1‑150 (MQ=255)
                                                                                cGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGgaga                                                           >  2:49992/1‑150 (MQ=255)
                                                                                    cTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGaa                                                       >  2:156363/1‑150 (MQ=255)
                                                                                    cTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCT‑GCTGGGCATGGCCCTGCCcagcca                                                                                                                                  <  2:110372/74‑3 (MQ=25)
                                                                                    cTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCT‑GCTGGGCATGGCCCTGCCcagcca                                                                                                                                  >  1:110372/1‑72 (MQ=25)
                                                                                        ccTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATggcgg                                                                                               >  1:889148/1‑106 (MQ=255)
                                                                                        ccTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATggcgg                                                                                               <  2:889148/106‑1 (MQ=255)
                                                                                          tGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGATGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGAc                                                 <  1:142439/150‑1 (MQ=255)
                                                                                                     ccAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGa                                      >  1:1066953/1‑150 (MQ=255)
                                                                                                         gTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATg                                                                >  1:939451/1‑120 (MQ=255)
                                                                                                         gTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATg                                                                <  2:939451/120‑1 (MQ=255)
                                                                                                         gTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAgcg                                  <  1:1139767/150‑1 (MQ=255)
                                                                                                              gCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCCCCACGGAATGGga                                                             >  1:974036/1‑118 (MQ=255)
                                                                                                              gCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGga                                                             <  2:974036/118‑1 (MQ=255)
                                                                                                                  ggTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGGTGGGCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGccc                         >  2:1007841/1‑150 (MQ=255)
                                                                                                                    tcatcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGa                       >  1:732384/1‑150 (MQ=255)
                                                                                                                    tcatcaGCGCGGTGGCCTGGCTGGGCAGAGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGa                       >  2:150025/1‑150 (MQ=255)
                                                                                                                     catcaGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGag                      <  2:364501/150‑1 (MQ=255)
                                                                                                                           cgcgGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGagcagcag                <  2:275526/150‑1 (MQ=255)
                                                                                                                                 caacctgctggGCATGGCCCTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGagcag                   <  2:727141/135‑1 (MQ=255)
                                                                                                                                 caacctgctggGCATGGCCCTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGagcag                   >  1:727141/7‑141 (MQ=255)
                                                                                                                                    ctggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATa                             >  1:1007476/1‑128 (MQ=255)
                                                                                                                                    ctggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATa                             <  2:1007476/128‑1 (MQ=255)
                                                                                                                                     tggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGAc                                                 <  1:731705/107‑1 (MQ=255)
                                                                                                                                     tggctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGAc                                                 >  2:731705/1‑107 (MQ=255)
                                                                                                                                      tgctggGCATGGCCCTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATgc     >  1:843300/2‑150 (MQ=255)
                                                                                                                                         tggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATgccgc  >  1:853839/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
GCTCGGCCAGCGTGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGC  >  NC_002947/86097‑86383

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