Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,246,185	C→G	13.8%	M30I (ATG→ATC)	alr ←	alanine racemase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,246,185	0	C	G	13.8%	128.3 / 10.9	56	M30I (ATG→ATC)	alr	alanine racemase
Reads supporting (aligned to +/- strand):  ref base C (19/29);  new base G (2/6);  total (21/35)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

ATAGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGT  >  NC_002947/4246038‑4246325                                                                                                                                                    |                                                                                                                                               aTAGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCa                                                                                                                                               >  8:19512/1‑149 (MQ=255)   aGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGTTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATc                                                                                                                                             >  6:184911/1‑149 (MQ=255)     gCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGa                                                                                                                                           >  8:147847/1‑149 (MQ=255)     gCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGa                                                                                                                                           >  8:143115/1‑149 (MQ=255)           gCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGTCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCGATGGACAA‑Cgg                                                                                                                                    <  1:171240/149‑1 (MQ=255)             cTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCGATGGACAA‑Cggca                                                                                                                                  <  1:139705/149‑3 (MQ=255)             cTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTACCGTTGTCGATGGACAA‑Cggca                                                                                                                                  <  3:76557/149‑3 (MQ=255)              ttGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑Cggtgg                                                                                                                                 >  2:38363/1‑149 (MQ=255)                 aGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGc                                                                                                                              <  8:150757/149‑1 (MQ=255)                 aGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGc                                                                                                                              <  5:20480/149‑1 (MQ=255)                  gCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGcc                                                                                                                             <  7:41475/149‑1 (MQ=255)                   cACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTTCAGGGCGCTGGCGCTGACTTAGCCCCAAGCATTGCAGTTTTTGACCGGCAGGGTGTTGGTTCCCTTGTCCATAGACAA‑CGGTTGTGCCg                                                                                                                            >  1:98283/1‑149 (MQ=255)                    aCGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGc                                                                                                                           <  7:19512/149‑1 (MQ=255)                             aGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAgggg                                                                                                                  >  1:342601/1‑149 (MQ=255)                               cTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAgggg‑gg                                                                                                               <  3:92400/149‑1 (MQ=255)                                  ggACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGTTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGc                                                                                                              <  5:184911/147‑1 (MQ=255)                                        gCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCg                                                                                                       <  5:117094/148‑1 (MQ=255)                                        gCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGCGGTGCCGCATACAGGGG‑GGCCTGTcc                                                                                                        <  8:21059/147‑1 (MQ=255)                                         ccGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAACCGGTGGTGCCGCATACAGGGG‑GGCCTGTCCg                                                                                                       <  1:38363/148‑1 (MQ=255)                                           ggCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGa                                                                                                   >  8:259101/1‑149 (MQ=255)                                                 cTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGaa                                                                                             <  3:6371/149‑1 (MQ=255)                                                    ggCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAg                                                                                            <  2:98283/147‑1 (MQ=255)                                                          cAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCGATCGACAA‑CGGTGGTGcc                                                                                                                             <  4:230288/109‑1 (MQ=255)                                                          cAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGcc                                                                                                                             >  3:230288/1‑109 (MQ=255)                                                          cAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAg                                                                                      <  5:355803/147‑1 (MQ=255)                                                                 cGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGc                                                                             <  7:259101/149‑1 (MQ=255)                                                                         gtgCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAggg                                                                     >  3:165754/1‑149 (MQ=255)                                                                         gtgCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAgg                                                                      <  7:147847/148‑1 (MQ=255)                                                                         gtgCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAgg                                                                      <  7:271403/148‑1 (MQ=255)                                                                         gtgCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAgg                                                                      <  7:143115/148‑1 (MQ=255)                                                                            ctgcAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTAc                                                                  >  2:58145/1‑149 (MQ=255)                                                                             tgcAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGt                                                                                           >  2:39866/1‑123 (MQ=255)                                                                             tgcAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGt                                                                                           <  1:39866/123‑1 (MQ=255)                                                                                     gctggcgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGGCGGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAg                                                          <  1:98968/149‑1 (MQ=255)                                                                                       tggcgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCAGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCAAGGGATGCAGCCAGAAGGGTACGTCGAAAtgg                                                        >  4:272116/1‑148 (MQ=255)                                                                                         gcgctgACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGGCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCAt                                                     <  1:58145/149‑1 (MQ=255)                                                                                               aCTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGAt                                               <  2:342601/149‑1 (MQ=255)                                                                                                    gACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTTGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTc                                          <  7:199857/149‑1 (MQ=255)                                                                                                    gACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTc                                          <  4:165754/149‑1 (MQ=255)                                                                                                         aGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCgtgt                                      >  7:195160/1‑148 (MQ=255)                                                                                                                gCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTc                                          >  7:121358/1‑137 (MQ=255)                                                                                                                gCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTc                                          <  8:121358/137‑1 (MQ=255)                                                                                                                gCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCg                              >  1:183941/1‑149 (MQ=255)                                                                                                                  tGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTAc                                                                  <  7:227126/111‑1 (MQ=255)                                                                                                                  tGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTAc                                                                  >  8:227126/1‑111 (MQ=255)                                                                                                                         cACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTgcgc                      <  2:183941/148‑1 (MQ=255)                                                                                                                            ggTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCACAAGGGTACGGCGATAGGGCATGGTGATTACTCGTGTTCTTTTCGTAGTGCGCGCCt                  >  7:281243/1‑149 (MQ=255)                                                                                                                              tCAGGGTGTTGGTGCCGTTGTACATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGCTCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTa                <  8:283822/149‑1 (MQ=255)                                                                                                                                aGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTTCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAAc              <  3:272116/149‑1 (MQ=255)                                                                                                                                 gggTGGTGGTGCCGTTGTCGATCGTCAA‑CGGTGGTGCCGCATACAGGGG‑GGGCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGTGCATTGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAAc              <  8:281243/148‑1 (MQ=255)                                                                                                                                         gTGCCTGTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGGCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAg      <  8:195160/148‑1 (MQ=255)                                                                                                                                           gCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgt    >  5:280961/1‑148 (MQ=255)                                                                                                                                           gCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGgt    <  6:280961/148‑1 (MQ=255)                                                                                                                                           accGTTGTCGATGGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTc                               <  4:228904/120‑1 (MQ=255)                                                                                                                                           accGTTGTCGATGGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTc                               >  3:228904/2‑121 (MQ=255)                                                                                                                                            ccGTTGTCGATGGACAA‑CCGTGGTGCAGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGCATGCAGCAAGAAGGGTACGGCGAAAGGGCATGGTGATTAGTCGTGTTCTTGTCGAAGTGCGCTCATTAACAGCCGCATgtgt  >  4:190968/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                               ATAGGCATCGGCCTTGAGCACGGCGCACAGCTTGGACTTGCCGGCCAGCTCGGCCTGCAGCGTGCGGATGTTGTGCTGCAGGGCGCTGGCGCTGACTTCGACCCAGGCATTGCTGTTTTGCACGGTCAGGGTGTTGGTGCCGTTGTCCATCGACAA‑CGGTGGTGCCGCATACAGGGG‑GGCCTGTCCGGTGATCAGAAGTGCCAGGGATGCAGCCAGAAGGGTACGGCGAAAGGGCATGGTGATTACTCGTGTTCTTGTCGTAGTGCGCGCCTTAACAGACGCAGGTGT  >  NC_002947/4246038‑4246325

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