Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,439,141	C→T	*123* (TAG→TAA)	yhdN ←	DUF1992 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,439,141	0	C	T	100.0%	231.0 / NA	73	*123* (TAG→TAA)	yhdN	DUF1992 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/36);  total (37/36)

TCGGTATATAGGCGAAACCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCACTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGG  >  NC_000913/3439006‑3439269                                                                                                                                        |                                                                                                                                 tCGGTATATAGGCGAAACCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAg                                                                                                                               >  2:1133335/1‑139 (MQ=255)             cGAAACCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGAtt                                                                                                                   <  1:1373161/139‑1 (MQ=255)               aaaCCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGAtttt                                                                                                                 <  1:429592/139‑1 (MQ=255)               aaaCCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGAtttt                                                                                                                 >  1:1535811/1‑139 (MQ=255)                 aCCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGt                                                                                                               >  1:902168/1‑139 (MQ=255)                            gTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTg                                                                                                    <  1:1228588/139‑1 (MQ=255)                                 cACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGAttt                                                                                                                  >  1:1461537/1‑120 (MQ=255)                                 cACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGAttt                                                                                                                  <  2:1461537/120‑1 (MQ=255)                                 cACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCTACAACTTGTCAGc                                                                                               >  2:1589714/1‑139 (MQ=255)                                    ttCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATa                                                                                            >  2:881624/1‑139 (MQ=255)                                    ttCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATa                                                                                            >  1:960667/1‑139 (MQ=255)                                      cATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACtt                                                                                                     <  2:562767/128‑1 (MQ=255)                                      cATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACtt                                                                                                     >  1:562767/1‑128 (MQ=255)                                         gCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATcgc                                                                                       >  2:830175/1‑139 (MQ=255)                                              atcatcTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGt                                                                                  >  2:758720/1‑139 (MQ=255)                                                catcTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTaa                                                                                <  2:1535811/139‑1 (MQ=255)                                                  tcTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTaaaa                                                                              <  1:1133335/139‑1 (MQ=255)                                                  tcTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTaaaa                                                                              >  2:1561826/1‑139 (MQ=255)                                                  tcTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTaaaa                                                                              >  2:407147/1‑139 (MQ=255)                                                     gctgTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTaa                                                                                >  2:202824/1‑134 (MQ=255)                                                     gctgTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTaa                                                                                <  1:202824/134‑1 (MQ=255)                                                       tgTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCg                                                                         >  2:253261/1‑139 (MQ=255)                                                          ttttCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTa                                                                      <  1:253261/139‑1 (MQ=255)                                                            ttCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACt                                                                    <  1:758720/139‑1 (MQ=255)                                                              cGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCcaac                                                                                                          <  1:189743/99‑1 (MQ=255)                                                              cGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCcaac                                                                                                          >  2:189743/1‑99 (MQ=255)                                                                   taaCGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCa                                                             >  2:390297/1‑139 (MQ=255)                                                                      cGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTcc                                                                                                             <  1:869741/88‑1 (MQ=255)                                                                      cGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTcc                                                                                                             >  2:869741/1‑88 (MQ=255)                                                                      cGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACtt                                                                                                     >  1:779243/1‑96 (MQ=255)                                                                      cGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACtt                                                                                                     <  2:779243/96‑1 (MQ=255)                                                                             tgtCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAAt                                                                           >  2:1602177/1‑115 (MQ=255)                                                                             tgtCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAAt                                                                           <  1:1602177/115‑1 (MQ=255)                                                                             tgtCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCg                                                                         >  2:1145777/1‑117 (MQ=255)                                                                             tgtCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCg                                                                         <  1:1145777/117‑1 (MQ=255)                                                                             tgtCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCa                                                   >  2:1265549/1‑139 (MQ=255)                                                                              gtCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTg                                                                                                    <  1:762304/89‑1 (MQ=255)                                                                              gtCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTg                                                                                                    >  2:762304/1‑89 (MQ=255)                                                                               tCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGc                                                 >  1:1237847/1‑139 (MQ=255)                                                                                    tgtTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACt                                                                    <  2:1448847/115‑1 (MQ=255)                                                                                    tgtTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACt                                                                    >  1:1448847/1‑115 (MQ=255)                                                                                     gtTACTTCCGCCATTTTTGCCAGCTCACCAATGCTATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAAt                                           <  2:1237847/139‑1 (MQ=255)                                                                                      ttACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAAtt                                          <  1:881624/139‑1 (MQ=255)                                                                                      ttACTTCCGCCATTTTTGCCAGCTCACCAATGCAATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAAtt                                          <  2:877804/139‑1 (MQ=255)                                                                                       tACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAAt                                           >  1:365362/1‑137 (MQ=255)                                                                                       tACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAAt                                           <  2:365362/137‑1 (MQ=255)                                                                                         cTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCa                                       <  1:830175/139‑1 (MQ=255)                                                                                             cGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGAtt                                                                                                                   <  2:127832/59‑1 (MQ=255)                                                                                             cGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGAtt                                                                                                                   >  1:127832/1‑59 (MQ=255)                                                                                              gCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAAt                                  >  1:910951/1‑139 (MQ=255)                                                                                               ccATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTaa                                                                                >  1:410083/1‑92 (MQ=255)                                                                                               ccATTTTTGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTaa                                                                                <  2:410083/92‑1 (MQ=255)                                                                                                     ttGCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATc                           <  1:1561826/139‑1 (MQ=255)                                                                                                       gCCAGCTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGa                         >  1:1404372/1‑139 (MQ=255)                                                                                                            cTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTc                                         <  1:1490359/118‑1 (MQ=255)                                                                                                            cTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTc                                         >  2:1490359/1‑118 (MQ=255)                                                                                                            cTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                    <  2:960667/139‑1 (MQ=255)                                                                                                            cTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                    <  2:513410/139‑1 (MQ=255)                                                                                                            cTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                    <  1:407147/139‑1 (MQ=255)                                                                                                            cTCACCAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                    >  2:1163589/1‑139 (MQ=255)                                                                                                                ccAATGCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAAcc                >  2:619708/1‑139 (MQ=255)                                                                                                                     gCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTc                                                                                                  <  1:83785/52‑1 (MQ=255)                                                                                                                     gCGATACATACATACTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTc                                                                                                  >  2:83785/1‑52 (MQ=255)                                                                                                                           catacatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTg     <  1:390297/139‑1 (MQ=255)                                                                                                                              acatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTc                                                                                                  <  2:907120/43‑1 (MQ=255)                                                                                                                              acatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTc                                                                                                  >  1:907120/1‑43 (MQ=255)                                                                                                                              acatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                    <  2:1162280/121‑1 (MQ=255)                                                                                                                              acatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCt                    >  1:1162280/1‑121 (MQ=255)                                                                                                                              acatacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCgg  <  2:902168/139‑1 (MQ=255)                                                                                                                               catacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTaa                                                                                >  2:443062/1‑60 (MQ=255)                                                                                                                               catacTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTaa                                                                                <  1:443062/60‑1 (MQ=255)                                                                                                                                   cTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCg                                                                         >  2:536268/1‑63 (MQ=255)                                                                                                                                   cTCCATTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCg                                                                         <  1:536268/63‑1 (MQ=255)                                                                                                                                        |                                                                                                                                 TCGGTATATAGGCGAAACCCACCTTCAGTACGCACTTCATGCTCCATCATCTGCTGTTTTTCGTAATAACGAATCGTGTCGGGTGTTACTTCCGCCATTTTTGCCAGCTCACCAATGCGATACATACATACTCCACTAGTTATCGTTGATTTTGTCCAACAACTTGTCAGCATAATCGCCGCGTAAAAAATCGGTACTCAATCCAGCTTGTCGCAGCTTCAATTCCAGTAATGACAATCGACGGCTAACCTCTTGATATTGCGG  >  NC_000913/3439006‑3439269

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