Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,007,693	G→A	*610* (TAG→TAA)	recQ →	ATP‑dependent DNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,007,693	0	G	A	100.0%	231.7 / NA	72	*610* (TAG→TAA)	recQ	ATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/37);  total (35/37)

CTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCC  >  NC_000913/4007563‑4007818                                                                                                                                   |                                                                                                                              cTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCagcagc                                                                                                                       >  2:919333/1‑139 (MQ=255)     aCAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGACGAGTAATcactgtctctt                                                                                                                   >  2:417649/1‑130 (MQ=255)      cAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATaaa                                                                                                                  <  1:890358/139‑1 (MQ=255)         aTGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATaaaaaa                                                                                                               >  2:554362/1‑139 (MQ=255)              gATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          <  1:925610/139‑1 (MQ=255)                    cGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATg                                                                                                    >  1:1247264/1‑139 (MQ=255)                    cGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATg                                                                                                    >  2:223152/1‑139 (MQ=255)                       cAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                           <  2:939663/113‑1 (MQ=255)                       cAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                           >  1:939663/1‑113 (MQ=255)                       cAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAg                                                                                                 <  2:1247264/139‑1 (MQ=255)                        aGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGa                                                                                                <  1:349448/139‑1 (MQ=255)                         gCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          >  2:1000283/1‑128 (MQ=255)                         gCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          <  1:1000283/128‑1 (MQ=255)                                               cGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGa                                                                         <  1:9091/139‑1 (MQ=255)                                                 ttGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                       <  1:1581923/139‑1 (MQ=255)                                                      aTGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAat                                                                  <  1:223152/139‑1 (MQ=255)                                                        gcgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                           >  2:1473107/1‑80 (MQ=255)                                                        gcgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                           <  1:1473107/80‑1 (MQ=255)                                                         cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          <  2:1117897/96‑1 (MQ=255)                                                         cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          >  2:1040377/1‑96 (MQ=255)                                                         cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          <  1:1040377/96‑1 (MQ=255)                                                         cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          >  1:1117897/1‑96 (MQ=255)                                                         cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGc                                                                                                           <  2:669444/95‑1 (MQ=255)                                                         cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGc                                                                                                           >  1:669444/1‑95 (MQ=255)                                                         cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                            <  2:181542/126‑1 (MQ=255)                                                         cgcAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                            >  1:181542/1‑126 (MQ=255)                                                                tGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGaa                                                                                                  <  2:694996/97‑1 (MQ=255)                                                                tGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGaa                                                                                                  >  1:694996/1‑97 (MQ=255)                                                                   aaCGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGtt                                                     >  1:577339/1‑139 (MQ=255)                                                                    aCGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAg                                                                        >  2:409584/1‑119 (MQ=255)                                                                    aCGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAg                                                                        <  1:409584/119‑1 (MQ=255)                                                                        tttGGCAAACCGTTTATGTCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGGCTCCGTAAACGTTTCCCCCGCGAGTCAAATGTAGGTTGAGGTTATTtc                                                >  1:758095/1‑139 (MQ=255)                                                                        tttGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                           >  2:1572375/1‑64 (MQ=255)                                                                        tttGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATcag                                                                                                                           <  1:1572375/64‑1 (MQ=255)                                                                        tttGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          >  2:485829/1‑81 (MQ=255)                                                                        tttGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGcc                                                                                                          <  1:485829/81‑1 (MQ=255)                                                                            gCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCAc                                            >  2:572331/1‑139 (MQ=255)                                                                                 ccGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCg                                       <  1:919333/139‑1 (MQ=255)                                                                                 ccGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCg                                       <  1:554362/139‑1 (MQ=255)                                                                                   gTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGcc                                     <  2:577339/139‑1 (MQ=255)                                                                                    tttATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                       <  1:1068954/104‑1 (MQ=255)                                                                                    tttATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                       >  2:1068954/1‑104 (MQ=255)                                                                                           cgcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGca                             <  1:572331/139‑1 (MQ=255)                                                                                            gcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATGAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGcac                            >  2:458655/1‑139 (MQ=255)                                                                                            gcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATg                                                       >  1:1693812/1‑112 (MQ=255)                                                                                            gcTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATg                                                       <  2:1693812/112‑1 (MQ=255)                                                                                              tGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                            >  2:1131589/1‑89 (MQ=255)                                                                                              tGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                            <  1:1131589/89‑1 (MQ=255)                                                                                              tGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                       >  2:882982/1‑94 (MQ=255)                                                                                              tGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                       <  1:882982/94‑1 (MQ=255)                                                                                                  tCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGaa                                                                                                  >  1:22836/1‑63 (MQ=255)                                                                                                  tCGTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGaa                                                                                                  <  2:22836/63‑1 (MQ=255)                                                                                                   tgTGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATAAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTgg                     <  2:758095/138‑1 (MQ=255)                                                                                                     tGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATtg                        >  1:1485583/1‑134 (MQ=255)                                                                                                     tGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATtg                        <  2:1485583/134‑1 (MQ=255)                                                                                                     tGCGCATGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGgcg                   >  1:3801/1‑139 (MQ=255)                                                                                                          aTGTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGT‑GCGCTTATg             >  2:1693395/1‑139 (MQ=255)                                                                                                            gTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtct                                               >  2:288688/1‑104 (MQ=255)                                                                                                            gTTGATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTtct                                               <  1:288688/104‑1 (MQ=255)                                                                                                               gATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATtg                        >  2:1795401/1‑124 (MQ=255)                                                                                                               gATGGCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATtg                        <  1:1795401/124‑1 (MQ=255)                                                                                                                  ggCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCg      <  1:527854/139‑1 (MQ=255)                                                                                                                   gCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGtt                                                           <  2:376368/85‑1 (MQ=255)                                                                                                                   gCGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGtt                                                           >  1:376368/1‑85 (MQ=255)                                                                                                                    cGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                       <  1:1343693/72‑1 (MQ=255)                                                                                                                    cGATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGt                                                                       >  2:1343693/1‑72 (MQ=255)                                                                                                                     gATGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGT‑GCGCTTATGAGCCCCGGTcc  <  1:1693395/139‑1 (MQ=255)                                                                                                                      aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                            <  1:1313369/65‑1 (MQ=255)                                                                                                                      aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCcg                                                                            >  2:1313369/1‑65 (MQ=255)                                                                                                                      aTGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTcc  >  1:1773618/1‑139 (MQ=255)                                                                                                                       tGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGgcg                   <  1:901342/121‑1 (MQ=255)                                                                                                                       tGACGAAGAGTAATCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGgcg                   >  2:901342/1‑121 (MQ=255)                                                                                                                                   |                                                                                                                              CTGAACAGATGCCGATCACCGCCAGCGAAATGCTCAGCGTTAACGGCGTTGGGATGCGCAAGCTGGAACGCTTTGGCAAACCGTTTATGGCGCTGATTCGTGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCCGGTCC  >  NC_000913/4007563‑4007818

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑