Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,645,313	G→C	8.8%	L222L (CTG→CTC) ‡	PP_4954 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,645,313	0	G	C	8.8%	173.8 / 7.1	57	L222L (CTG→CTC) ‡	PP_4954	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (26/26);  new base C (3/2);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01

GGTACAGCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGA  >  NC_002947/5645166‑5645455                                                                                                                                                    |                                                                                                                                                ggTACAGCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGa                                                                                                                                                <  4:288274/149‑1 (MQ=255)        cagAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCa                                                                                                                                           >  6:9241/1‑147 (MQ=255)        cagAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCa                                                                                                                                           <  5:9241/147‑1 (MQ=255)        cagAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGc                                                                                                                                         >  5:17126/1‑149 (MQ=255)           agCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCTGCGTTCAGCGCGGCC‑ac                                                                                                                                     <  5:106341/149‑2 (MQ=255)             ccTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGCGGTGGCCGCGCTGAACGCAGCC‑AAgg                                                                                                                                   <  1:8638/149‑1 (MQ=255)                 gCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGgaga                                                                                                                                <  6:172980/148‑1 (MQ=255)                 gCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGgaga                                                                                                                                <  4:283318/148‑1 (MQ=255)                     gACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACTGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGatcatc                                                                                                                           <  8:287798/149‑1 (MQ=255)                      aCTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGgaga                                                                                                                                >  3:263112/1‑143 (MQ=255)                      aCTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGgaga                                                                                                                                <  4:263112/143‑1 (MQ=255)                       cTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATcc                                                                                                                          <  7:276147/148‑1 (MQ=255)                             ttcGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCg                                                                                                                   >  5:90306/1‑149 (MQ=255)                                 aGCTGCTCAGCTACGTGAAGCTCAAGTAAATCCTGTTTTTTACGGTCCGCTCGCTCGACGTCTAGGCAAGGCGCCTGGAACTGGCCCAGCAGAAGGGGAAGGTGGTGGAGGGGCTGATCGGTGCC‑AAGGAGATCTTCCTGGAGCGGGTg                                                                                                               <  8:236219/149‑1 (MQ=255)                                                tGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGccaacc                                                                                                 >  6:60923/1‑148 (MQ=255)                                                 gAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGccaacca                                                                                                >  6:195988/1‑148 (MQ=255)                                                    gttCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACt                                                                                            >  6:90464/3‑149 (MQ=255)                                                      tCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTc                                                                                           >  8:252379/1‑148 (MQ=255)                                                      tCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTc                                                                                           >  2:43673/1‑148 (MQ=255)                                                        aaGGCCAGCCTGTTTGACCCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCCGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCg                                                                                          >  1:129732/1‑147 (MQ=255)                                                        aaGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGa                                                                                         <  5:195988/148‑1 (MQ=255)                                                        aaGGCCAGCCTGTCTGATCCCGTCCGGCCGCTCGACCACGACGCTCGGCGCCTGGAACTGGCCCAGCAGTACGGCAAGGTGGTGGACGCGCTGAACGCAGCC‑AAGCAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGaa                                                                                        <  6:164/149‑1 (MQ=255)                                                           gCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGt                                                                                     <  5:138587/149‑1 (MQ=255)                                                               gCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGc                                                                                 <  5:60923/149‑1 (MQ=255)                                                                       gAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCt                                                                         >  3:165729/1‑149 (MQ=255)                                                                         aCCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGCACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTc                                                                        >  4:172981/1‑148 (MQ=255)                                                                          ccGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTc                                                                        >  3:13846/1‑147 (MQ=255)                                                                           cGGTCCGGACGCTCGACGTCGAGGCACGGCGCCTGGAACTGTCCCAGGAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGTTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAg                                                                     <  6:90306/149‑1 (MQ=255)                                                                             gTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCGAGCCCACCTCGAAGGTTAGCCGCTGCCTCa                                                                       >  8:211064/1‑145 (MQ=255)                                                                             gTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCa                                                                       <  7:211064/145‑1 (MQ=255)                                                                             gTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCt                                                                   >  8:206521/1‑149 (MQ=255)                                                                                cGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGt                                                                 <  1:43673/148‑1 (MQ=255)                                                                                 gCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTAc                                                               <  3:294666/149‑1 (MQ=255)                                                                                     gCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGccaa                                                                                                   <  5:201418/109‑1 (MQ=255)                                                                                     gCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGccaa                                                                                                   >  6:201418/1‑109 (MQ=255)                                                                                       tCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTg                                                         >  1:312403/1‑149 (MQ=255)                                                                                              cGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTTGCCCAgg                                                  <  7:206521/149‑1 (MQ=255)                                                                                                 ggCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGt                                                                 >  5:204582/1‑131 (MQ=255)                                                                                                 ggCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGt                                                                 <  6:204582/131‑1 (MQ=255)                                                                                                    aCGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCCTAAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATc                                             >  1:249052/1‑149 (MQ=255)                                                                                                         gcCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAg                                       >  3:88431/1‑149 (MQ=255)                                                                                                             ggAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCgtg                                   >  5:141565/1‑149 (MQ=255)                                                                                                              gAACTGGCCCAGCAGAACGGCAAGGTTGTGGCCGCGCCGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCgt                                    <  2:129732/147‑1 (MQ=255)                                                                                                                      ccAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCt                          <  4:165729/149‑1 (MQ=255)                                                                                                                           agAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCAc                     <  5:90464/149‑1 (MQ=255)                                                                                                                           agAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCAc                     <  4:13846/149‑1 (MQ=255)                                                                                                                                     ggtggtGGCCGCGCCGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTAc            >  2:54891/1‑148 (MQ=255)                                                                                                                                      tccttggCTGCGTTCAGCGCGGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAc                  >  8:27213/5‑141 (MQ=255)                                                                                                                                      tccttggCTGCGTTCAGCGCGGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTAc                  <  7:27213/137‑1 (MQ=255)                                                                                                                                       tggtggCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAAc         >  7:9925/1‑149 (MQ=255)                                                                                                                                         ttggCTGCGTTCAGCGCGGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACg        >  2:97939/2‑148 (MQ=255)                                                                                                                                         ttggCTGCGTTCAGCGCGGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACg        >  4:136751/2‑148 (MQ=255)                                                                                                                                         gtggCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACg        >  2:310520/1‑148 (MQ=255)                                                                                                                                           ggCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCAc     >  5:71656/1‑149 (MQ=255)                                                                                                                                              cgcgCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCa                                >  4:11081/1‑119 (MQ=255)                                                                                                                                              cgcgCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCa                                <  3:11081/119‑1 (MQ=255)                                                                                                                                              cgcgCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGa  <  4:88431/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                GGTACAGCAGAGCCTGGCCAGACTGTTCTTCGAGCTGGGCAGCTACCTGAAGCTCAAGGCCAGCCTGTTTGAACCGGTCCGCACGCTCGACGTCGAGGCACGGCGCCTGGAACTGGCCCAGCAGAACGGCAAGGTGGTGGCCGCGCTGAACGCAGCC‑AAGGAGATCATCCTGCACCGGGTGGGCAACAGCCAACCAAACTCGAAGGTAAGCCGCTACCTCAAGCTGTACTTCCTGGCCCAGGATATCCACGAGCGTGTCAGTGCCTCGCACTACCCCTACAACGCACTGA  >  NC_002947/5645166‑5645455

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: ‑