Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,894,839	T→C	L12P (CTC→CCC)	pabB →	aminodeoxychorismate synthase subunit 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,894,839	0	T	C	100.0%	110.3 / NA	35	L12P (CTC→CCC)	pabB	aminodeoxychorismate synthase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/11);  total (24/11)

CTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTC  >  NC_000913/1894612‑1895081                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    ctctTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAAt                                                                                                                                                                                                                              >  1:161342/1‑251 (MQ=255)    ttctGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAAtttt                                                                                                                                                                                                                           >  2:3931/1‑251 (MQ=255)                 gcgATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCCTGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCtt                                                                                                                                                                                                              >  2:70008/1‑251 (MQ=255)                   gagtGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTTCCCGGCACAGAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTTATTTCTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTaa                                                                                                                                                                                                            <  1:212812/248‑1 (MQ=255)                             ccGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGcc                                                                                                                                                                                                  >  1:207797/1‑251 (MQ=255)                                       tCCCTTGTGCCATCAGTTCCTGGATACGCTGGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATg                                                                                                                                                                                        >  2:284041/1‑251 (MQ=255)                                                           tGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGATTTTTATTTCTCCCGCTTACGCCACCTGCCGTGGGCGTTGCTTTTACACTCCGGCTATGc                                                                                                                                                                    >  2:86544/1‑251 (MQ=255)                                                           tGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGc                                                                                                                                                                    >  1:229127/1‑251 (MQ=255)                                                                        cAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGtat                                                                                                                                                       >  2:237618/1‑251 (MQ=255)                                                                                   tgctgTTCATGGGTGAGTGAAGGTAAACCTGCCAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTg                                                                                                                                             >  1:448295/1‑250 (MQ=255)                                                                                   tgctgTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTg                                                                                                                                             >  1:137441/1‑250 (MQ=255)                                                                                    gctgTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTg                                                                                                                                             <  1:237618/249‑1 (MQ=255)                                                                                        ttCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTCTAGCCGCTTTGATATtg                                                                                                                                       >  2:393039/1‑251 (MQ=255)                                                                                                 gagtgaAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCgccg                                                                                                                             >  2:350544/1‑251 (MQ=255)                                                                                                 gagtgaAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTAGCTCCCGCTGTGATTCCTTTACCCTGGCTTCAGGACGCCGCTGAATTTTATGTCTCCCGCTGAATCCACCTGCCGTGGGCGTTGCTTTTACACTCCGGCGAGGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTctcc                                                                                                                              <  2:459413/250‑4 (MQ=255)                                                                                                  agtgaAGGACAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCAGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTCATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGACACCTGCCGAGGGCGATGCTTTTACACTCCGGCTATGCCGATCAACAGTATAGCCGATTTGATATTTTG‑GTCgacg                                                                                                                             >  2:302989/1‑247 (MQ=255)                                                                                                     gaAGGTAATCCTGCAAACATTGTTAACTCCTGTGAATTTGTTGGGTCTAAGTGTTTCATGCTCCCCGGCACATCGCCAGGAGAGTCATGACTTATGAAGGCGTTATCTCCCGCTGGGAGTACTTTACCCGGGCGGCAGGACGCCGCTTAATTTTTTTTCTCCCGCTTAAGCCACCGGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATATCCGCGTTGATATTGTG‑GTCgccgagcc                                                                                                                         <  2:405407/251‑1 (MQ=255)                                                                                                           aaaCCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCGCCGAGCCGATTTg                                                                                                                   >  2:420370/1‑251 (MQ=255)                                                                                                              catGCAAACATTGTTAACTCCTGCCAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATCGCCCGTAGAGTCCGGGCTGATGATGACGTTATCTCCCGATGTGATTACTTTACCCTGGCGTCAGGAAGCCGCTGACTTTTATTTCTCCCGCTTAAGCCACCTGACGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATT‑TGTGTCGCCCACCCGAGTTGCAc                                                                                                                >  2:453483/3‑251 (MQ=255)                                                                                                                       aTTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCCGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCGCCGAGCCGATTTGCACTTTAACCAc                                                                                                       >  2:415322/1‑251 (MQ=255)                                                                                                                            tAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCGCCGATCCGATTTGCACTTTAACCACTTTCg                                                                                                  >  2:342507/1‑251 (MQ=255)                                                                                                                             aaCTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCGCCGAGCCGATTTGCACTTTAACCACTTTCg                                                                                                  <  2:366427/250‑1 (MQ=255)                                                                                                                                          ttgttgGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTc                                                                                                                                                                             >  1:242787/1‑163 (MQ=255)                                                                                                                                          ttgttgGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTc                                                                                                                                                                             <  2:242787/163‑1 (MQ=255)                                                                                                                                            gttgGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCgttgt                                                                                  >  1:137859/1‑251 (MQ=255)                                                                                                                                                    tAATTATTTCATGCTACCCGGCACATAGCCAGTCGAGTCAGGACTGATGAAGACGTTATCTCC‑GCTGTGATTACTTTACCCTGGCGTCAAGACACCGCTGCATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGaa                                                                           >  1:150261/1‑249 (MQ=255)                                                                                                                                                        tattGCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATGTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCGCGGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCg                                                                       <  1:284041/250‑1 (MQ=255)                                                                                                                                                           ttCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCCACCGCTGTGATTACTTTACCCTGGCGCAAGGACGCCGCTGAATTTTCTTTCTCCAGCTTAAGCCACCTGCAGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGaaa                                                                    >  1:372666/1‑250 (MQ=255)                                                                                                                                                                     ccGGCACATAGCCAGTAGAGTCAGGACTGATGAAGAAGTTCTCTCCCGCTGTGCTTACTTTACCCAGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTACGCCACCTGCCGAGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGACGCTTTGATATTGTG‑GTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACg                                                         >  1:262473/1‑251 (MQ=255)                                                                                                                                                                       ggCACATAGCCTGTAGAGTCAGGATTGATGAATACGTTATCTCCCTCTGTGATTACTTTGCCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACGCTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATTTTGTG‑GTCGCCGAGCCGATTTGCACTTGAACCACTTTCGGTAAAGAATCCTTTGTTAGGGAAAGCGAAAAACGCACAACGAc                                                       <  1:164889/251‑1 (MQ=255)                                                                                                                                                                                ccAGTAGAGTCAGGACTGATGAAGCCGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACtgatga                                              <  1:3931/251‑1 (MQ=255)                                                                                                                                                                                            ggACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCGGGCGTCAGGACGCCGCTGAATTTTATTCCGCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTTATATTGTG‑GTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGAGGACCCGCTACAGGt                                  <  2:372666/251‑1 (MQ=255)                                                                                                                                                                                                 gatgaAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCt                               <  2:112747/249‑1 (MQ=255)                                                                                                                                                                                                                    gCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCgcgc          >  2:351078/1‑251 (MQ=255)                                                                                                                                                                                                                         gATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACAGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGAc      >  1:379934/1‑250 (MQ=255)                                                                                                                                                                                                                             actttaccCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTc  <  1:70008/250‑1 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    CTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTG‑GTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTC  >  NC_000913/1894612‑1895081

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