Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,598,833	G→A	12.2%	intergenic (‑72/+59)	dctB ← / ← bglX	C4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,598,833	0	G	A	12.2%	125.0 / 7.6	49	intergenic (‑72/+59)	dctB/bglX	C4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase
Reads supporting (aligned to +/- strand):  ref base G (22/21);  new base A (4/2);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-01

CGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGCAA  >  NC_002947/1598686‑1598981                                                                                                                                                                 |                                                                                                                                                     cGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTgc                                                                                                                                                     <  1:124818/149‑1 (MQ=255)       tGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCAcc                                                                                                                                               >  7:107041/1‑149 (MQ=255)                  cccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTTGGAGAAGGGGATGGGGAGAAGTCGGGCCCTGCCTCTTCGCGGGTAAACCCGCGCATACAGGAACTCTccagca                                                                                                                                                 >  1:271274/1‑149 (MQ=255)                     aGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGttt                                                                                                                                 >  2:168158/1‑149 (MQ=255)                              agcagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCcgcg                                                                                                                         <  3:146882/148‑1 (MQ=255)                              agcagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGCAGGCGCGGTTTTACCcgcg                                                                                                                         >  4:13038/1‑148 (MQ=255)                                cagGGTGATCAGCCCCAGACGCAGTGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTGCTGCAGAGTACCTGTAg                                                                                                                                          <  8:148282/129‑1 (MQ=37)                                cagGGTGATCAGCCCCAGACGCAGTGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTGCTGCAGAGTACCTGTAg                                                                                                                                          >  7:148282/1‑129 (MQ=37)                                       aTCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCg                                                                                                                <  8:107041/148‑1 (MQ=255)                                            ccccAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCAGatcgg                                                                                                          <  4:286715/149‑5 (MQ=255)                                            ccccAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTACGCGGGTAAACCCGCGCCTACAGGTAATCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGCCCGGAACAg                                                                                                          >  5:171561/1‑149 (MQ=255)                                              ccAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggc                                                                                                        >  7:6715/1‑149 (MQ=255)                                                  aCGCAGGGCGCGAAAGGAAATGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGt                                                                                                                                            <  8:307577/109‑1 (MQ=255)                                                  aCGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCAGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGt                                                                                                                                            >  7:307577/1‑109 (MQ=255)                                                       gggCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCttttt                                                                                               <  1:324549/149‑1 (MQ=255)                                                        ggCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCCCCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggc                                                                                                        >  5:8504/1‑139 (MQ=255)                                                        ggCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggc                                                                                                        <  6:8504/139‑1 (MQ=255)                                                         gcgcGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGCAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                             <  3:13038/149‑1 (MQ=255)                                                         gcgcGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGGCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                             <  7:109050/149‑1 (MQ=255)                                                             gAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGCCCTTTTTCAGTc                                                                                          >  7:39632/1‑148 (MQ=255)                                                              aaaGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGt                                                                                           >  3:147000/1‑146 (MQ=255)                                                              aaaGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGt                                                                                           <  4:147000/146‑1 (MQ=255)                                                                   aaaaGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCaa                                                                                   >  6:28893/1‑149 (MQ=255)                                                                       gggCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCa                                                                               <  3:96180/149‑1 (MQ=255)                                                                         gCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAATCCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCaaa                                                                             <  2:84425/149‑1 (MQ=255)                                                                              gggggCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGctc                                                                         <  8:6715/148‑1 (MQ=255)                                                                                 ggCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCCCCTGTAGGCGCGGTTTAACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                       >  4:217759/1‑131 (MQ=255)                                                                                 ggCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTAACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                       <  3:217759/131‑1 (MQ=255)                                                                                                       g‑tgg‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTg                                                        <  1:168158/149‑1 (MQ=255)                                                                                                           g‑‑‑gATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTggg                                                      <  5:28893/148‑1 (MQ=255)                                                                                                                     ggggCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGc                                              >  8:15855/1‑149 (MQ=255)                                                                                                                      gggCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGcc                                             <  7:15855/149‑1 (MQ=255)                                                                                                                       ggCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTGTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGccc                                            >  1:205139/1‑149 (MQ=255)                                                                                                                          ctgcctCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCGTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTAc                                         >  7:174771/1‑149 (MQ=255)                                                                                                                             cctCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTg                                      >  6:1577/1‑149 (MQ=255)                                                                                                                              ctctTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGc                                              >  4:299755/1‑140 (MQ=255)                                                                                                                              ctctTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGc                                              <  3:299755/140‑1 (MQ=255)                                                                                                                                    gctggTAAAACCGCGCCTACAGGTACTCTGAAGCAACTGTATGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTTCTGCTGCACCGCCTTGGAATCCAGCCCTACCTGGACATTg                               <  2:271274/146‑1 (MQ=255)                                                                                                                                      gggTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcgg                                                                                                      >  4:195506/1‑76 (MQ=255)                                                                                                                                      gggTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcgg                                                                                                      <  3:195506/76‑1 (MQ=255)                                                                                                                                          aaaaCCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGATCTGCTGCTGCACCGCCTgg                                                       >  3:52106/1‑119 (MQ=18)                                                                                                                                           aaCCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTccc                        >  5:326961/1‑149 (MQ=255)                                                                                                                                           aaCCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGCCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTccc                        <  8:39632/149‑1 (MQ=255)                                                                                                                                                gcgcCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGt                   >  4:316359/1‑149 (MQ=37)                                                                                                                                                   ccTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGa                              >  6:141317/1‑135 (MQ=25)                                                                                                                                                   ccTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGa                              <  5:141317/135‑1 (MQ=25)                                                                                                                                                     tACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCgg               >  4:128821/1‑148 (MQ=255)                                                                                                                                                              tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTg     >  2:34402/1‑149 (MQ=255)                                                                                                                                                                tGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                             <  6:36999/59‑1 (MQ=255)                                                                                                                                                                tGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                             >  5:36999/1‑59 (MQ=255)                                                                                                                                                                tGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGc    >  6:185268/1‑148 (MQ=255)                                                                                                                                                                 gcagcaCCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCGCCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGCaa  >  1:106303/1‑149 (MQ=255)                                                                                                                                                                 gcagcaCCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGAc                                   >  8:123611/1‑116 (MQ=255)                                                                                                                                                                 |                                                                                                                                                     CGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGT‑‑‑‑‑‑‑‑‑GG‑TGG‑‑‑GATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGCAA  >  NC_002947/1598686‑1598981

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