Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,635,889	T→C	D123D (GAT→GAC)	creA →	PF05981 family protein CreA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,635,889	0	T	C	97.1%	86.6 / ‑4.3	34	D123D (GAT→GAC)	creA	PF05981 family protein CreA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (17/16);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.01e-01

TTGATGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGATGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATA  >  NC_000913/4635654‑4636130                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   ttGACGATCCCGATGTGAATAATGTCAACAGTTAAGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTTTGGCGGAATATACCTCCGATGCTGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGGCGTGCGCTTTTATGACGCCAaacgcaacgc                                                                                                                                                                                                                                      <  2:41117/250‑1 (MQ=255)  tgacgatCCCGATGTGAAAAATGTCCCCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGTGCCGATTGAACTGTCGGATCGTATTAAAAACGTCAAAGCTCAGGTCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGTTCGTGCGCTTTTATGACGCCAAACGCAAcgcg                                                                                                                                                                                                                                     <  1:155492/246‑1 (MQ=255)                                       aGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAgg                                                                                                                                                                                              >  1:33071/1‑251 (MQ=255)                                         ccGGGCGAAACCAGGTGGTATTAAAGGTGGATTGGGTCTGGCGGAAGATACCTCAGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGCGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTAATCTGGCTTAATCCGACAAAGT‑TGtgtaagggt                                                                                                                                                                                            >  2:148652/1‑243 (MQ=255)                                             gCGAAAACAGGTGGTATTAAAGGGGGATTTGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGcc                                                                                                                                                                                        <  1:204809/251‑1 (MQ=255)                                                   acagagggtATTAAGTGGGTATTGGGTCGGGCGGAAGATTCCTACGATGAGTCGATATCATGTCAGGAAGTCGGGCCTATGTAATTGTCGGAGCGTATTAAAAACGGCAAAGCTCAGGACGAGGGAGTAGTCAAAAAACGCACGTCACTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTGTATGACGTCAAACGCAGCGCGCTCGATAATCTGGCTTACTCCGACAAAGT‑TGTATAAGTTTCGCCGaaaaa                                                                                                                                                                                  <  2:117753/245‑1 (MQ=255)                                                                               tGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTATCGCGGTTCCTGTCATGCCg                                                                                                                                                      >  2:274046/1‑251 (MQ=255)                                                                                           ccTCCGATGCGGCCGTTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAGACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGAACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGATGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGTAGGTTCACCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCaaaa                                                                                                                                           >  2:143261/1‑248 (MQ=255)                                                                                                gATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACagag                                                                                                                                     <  1:285272/251‑1 (MQ=255)                                                                                                    cggctatttcaggTCAGCAAGTCGAGCCGATTGAACGGTCGGATCGTATTAAAAACGGCAAAGCGCAGGGCGAGGTAGTATGCAAAAAACGCACGTCCCTGGTCTGTAAGTCGTTACAGTTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCGTAGCTGGCTTACGCCGACAAAGT‑TGTAGAAGGTTCTCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGa                                                                                                                                 <  2:353354/239‑1 (MQ=255)                                                                                                           ttCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAACCGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAACCGCGATTAGCGCGGTTCCTGTCTTGCCGTGGCGGCAATAACAGAGGCGTTTTATGc                                                                                                                          >  2:128164/1‑251 (MQ=255)                                                                                                                       aGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGACATAACAGAGGAGATTTATgcaacgggcaacgg                                                                                                              >  2:165746/1‑251 (MQ=255)                                                                                                                            ggCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGGGGCGGCCATAACAGAGGCGATTTATGCAACGGGAACCGGTGTgg                                                                                                         >  2:167327/1‑251 (MQ=255)                                                                                                                                  ttGAACTGTCGGATCGTATTAAAAACGGCAAAGATCAGGGCGAGGTAGTATTCACAAAATGCAAGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAATGCGCTCGCTGATCTGGCTTACTCCGACAA‑GTGTGTAGAAGGATCGCCGAAAAACGAGATTAGCGCGGTTCCTGTTAGGCCGTGGCTGGAATAACAGACGCGCTTTCAGCAACGGGAACCgggatgggtcgtg                                                                                                   >  2:279421/1‑240 (MQ=255)                                                                                                                                      aCTGTCGGATCGTAGGAAAGCAGGCAACGCTCAGGGCGCGGTAGAATTCAAAAAAAGCACGTCCCTGGTATTTAAGTCGTTACAGGACGTGCGCTTTTATGCCGCCAAACGCAACGCGCACGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCAGAAAAACGCGAGTCTCGCGGTTCATGTCATGCCGTAGCGGCACTTAAAGAAGCGATTTATGCCAGAGGAAAAGGTCTGGtaggttgat                                                                                                >  2:218150/1‑242 (MQ=255)                                                                                                                                        tGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCgg                                                                                                                                                                    <  1:288422/180‑1 (MQ=255)                                                                                                                                             gctCGTATTAAAAACGGAAAAGCTCAGGGCGAGGTAGTATTCACAAAACGCACCTCCCTGGTCTATAAACCGTTACAGGACGTGCGCTTTTATGACGCCAAACGCAACGCGATCGCTTATATGACTTACTCCGACAAAGT‑TGTAGAAGATTCGCCGAAAAACGAGCTTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAAAAAAGGCGCTTTATACAACTGGAAACGGGCTGGTGAGTggtaggtggga                                                                                         >  2:373323/3‑241 (MQ=255)                                                                                                                                                       aaaaCGGCAAATCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGcc                                                                              <  1:335084/251‑1 (MQ=255)                                                                                                                                                         aaCGGGAAAGCTCAGGGCGAGGTAGGATGCAAAAAACGCCCGGCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGTCAAACGCAACGCGCTCGCTTACCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGa                                                                            <  2:78063/251‑1 (MQ=255)                                                                                                                                                            ggCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGGTAGTGGAAGATGAGCAAGGGATAGCCGacac                                                                         >  1:332323/1‑251 (MQ=255)                                                                                                                                                               aaaGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCt                                                                      <  1:268075/251‑1 (MQ=255)                                                                                                                                                                   cccAGGGCGAGTTAGTATGCAAAAAGCGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCGTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAATGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATGTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTc                                                                  <  2:137477/249‑1 (MQ=255)                                                                                                                                                                            aGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTTTGCAACGGGGAACGGTCTGGTTATTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCttactgttg                                                         >  2:67820/1‑243 (MQ=255)                                                                                                                                                                              gtagtaTTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTgc                                                        >  1:157901/1‑250 (MQ=255)                                                                                                                                                                                agtaTTCAAAAAACGCACGTCCCTGGACTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTgaagc                                                     >  1:112833/1‑251 (MQ=255)                                                                                                                                                                                           aaCGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCACAGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACAAGATGGTCTACATGTTGAAGCAGGAAGGgttt                                          >  2:309708/1‑251 (MQ=255)                                                                                                                                                                                            aCGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTg                                         <  1:167327/251‑1 (MQ=255)                                                                                                                                                                                              gCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGcc                                       >  1:150012/1‑251 (MQ=255)                                                                                                                                                                                                    cccTGGGCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGa                                  <  1:201809/250‑1 (MQ=255)                                                                                                                                                                                                       tGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTc                              >  1:203734/1‑251 (MQ=255)                                                                                                                                                                                                                 gTCGTTACTGGCCGTGGCCTTTTATGACGACAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCTATTAGCGCGGTTCCGGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTGAGTGGAAGAAGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCATGAAGGCTTTGCCGTCGAGGTCTTTGAGCGAg                    <  2:157901/251‑1 (MQ=255)                                                                                                                                                                                                                   cGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGGTAGTGGAAGATGAGCAAGGGATAGCCGACCCGCTGGTCTAATTGTTGCAGCAGGAAGGTTTTGCCGTCGGGGTCTTTGAGCGAGGc                  >  2:20254/1‑251 (MQ=255)                                                                                                                                                                                                                        cagttcattcgcTTTTATGACGCGAAACGCAACGCGCCCGCTTATCTGGCTTACTCCGACAAAGG‑TGTAGAAGGTTCGCCGAAAAACTCGATTTGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGTAAACGGTCTGGTTAGTGGAAGATGATCAAGGGATAGCCGGCACGCTGGTTTACATTTTGCAGCAGGAAGTTTTTGCCGTCGAGGTCTGTGAGCGAGGCTTGcc             <  2:227202/242‑1 (MQ=255)                                                                                                                                                                                                                                     tttATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATa  <  1:218150/249‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                   TTGATGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGATGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGT‑TGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAGGCGATTTATGCAACGGGAAACGGTCTGGTTAGTGGAAGATGAGCAAGGGATAGCCGACACGCTGGTCTACATGTTGCAGCAGGAAGGTTTTGCCGTCGAGGTCTTTGAGCGAGGCTTGCCGGTGCTGGATA  >  NC_000913/4635654‑4636130

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