Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,486,387	T→C	F154S (TTT→TCT)	tsgA →	MFS transporter TsgA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,486,387	0	T	C	98.4%	236.3 / ‑4.9	64	F154S (TTT→TCT)	tsgA	MFS transporter TsgA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (25/38);  total (26/38)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTTTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCGGCTGTGAGTTCCCGGCGCTGGGCAAACATGCGCCAAAAACGGATGCT  >  NZ_CP009273/3486251‑3486517                                                                                                                                         |                                                                                                                                   gtCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGtctc                                                                                                                                  >  1:889786/1‑139 (MQ=255)      cGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTg                                                                                                                             >  1:505479/1‑139 (MQ=255)           tCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTg                                                                                                                             >  2:524679/1‑134 (MQ=255)           tCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTg                                                                                                                             <  1:524679/134‑1 (MQ=255)           tCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGgcgcg                                                                                                                        <  1:748139/139‑1 (MQ=255)           tCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGgcgcg                                                                                                                        <  2:33568/139‑1 (MQ=255)                  gTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCAtt                                                                                                                 >  2:520916/1‑139 (MQ=255)                     gATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAg                                                                                                              <  2:217916/139‑1 (MQ=255)                     gATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAg                                                                                                              <  1:541759/139‑1 (MQ=255)                     gATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAg                                                                                                              <  1:44194/139‑1 (MQ=255)                      aTTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCgcagca                                                                                                                   >  2:715424/1‑133 (MQ=255)                      aTTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCgcagca                                                                                                                   <  1:715424/133‑1 (MQ=255)                      aTTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGt                                                                                                             <  1:816867/139‑1 (MQ=255)                          gTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTa                                                                                                         <  1:1064530/139‑1 (MQ=255)                             cATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACt                                                                                                       <  2:382468/138‑1 (MQ=255)                             cATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACt                                                                                                       <  2:1002012/138‑1 (MQ=255)                                  cTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGtt                                                                                                  <  2:544386/138‑1 (MQ=255)                                    ggTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTTTCTACTGGCGCGCAGCCTTGAGTGGGACTGGGTTTAt                                                                                               >  2:872706/1‑139 (MQ=255)                                    ggTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGgcgcg                                                                                                                        <  1:1274872/114‑1 (MQ=255)                                    ggTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGgcgcg                                                                                                                        >  2:1274872/1‑114 (MQ=255)                                     gTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATg                                                                                              >  2:405703/1‑139 (MQ=255)                                     gTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATg                                                                                              >  2:583993/1‑139 (MQ=255)                                     gTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATg                                                                                              >  2:23461/1‑139 (MQ=255)                                     gTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATg                                                                                              >  2:625711/1‑139 (MQ=255)                                        acacAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGtt                                                                                                  >  2:998217/1‑132 (MQ=255)                                        acacAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGtt                                                                                                  <  1:998217/132‑1 (MQ=255)                                                 tATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCggg                                                                                   <  1:160596/138‑1 (MQ=255)                                                 tATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCggg                                                                                   <  1:1219934/138‑1 (MQ=255)                                                 tATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCggg                                                                                   <  2:1154533/138‑1 (MQ=255)                                                 tATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCggg                                                                                   <  2:1167397/138‑1 (MQ=255)                                                 tATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGc                                                                                  <  2:1274460/139‑1 (MQ=255)                                                     aaGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGgt                                                                              <  1:23461/139‑1 (MQ=255)                                                       gggCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCTTTGAGTGGTACTGTGTTTATGCCTGCATCGGGCTCgtgt                                                                            <  2:142253/139‑1 (MQ=255)                                                          cGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCggg                                                                                   >  1:1286562/1‑129 (MQ=255)                                                          cGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCggg                                                                                   <  2:1286562/129‑1 (MQ=255)                                                          cGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATg                                                                         >  1:95667/1‑139 (MQ=255)                                                           gTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGt                                                                        >  2:986111/1‑139 (MQ=255)                                                            tCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTTAGTATGGATGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCCTTGAGTGGTACTGGGTTTATGCCTGCCTCGGGCTGGTGTATGTc                                                                       >  2:20123/1‑139 (MQ=255)                                                              aGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGc                                                                     <  1:625711/139‑1 (MQ=255)                                                                cgcgGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTa                                                                   >  2:1103647/1‑139 (MQ=255)                                                                    gTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTAtttt                                                               >  2:820567/1‑139 (MQ=255)                                                                               ttATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCggg                                                                                   >  1:1159530/1‑108 (MQ=255)                                                                               ttATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCggg                                                                                   <  2:1159530/108‑1 (MQ=255)                                                                               ttATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCt                                                    >  2:119281/1‑139 (MQ=255)                                                                                  tttACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCg                                                 <  2:1206343/139‑1 (MQ=255)                                                                                        gACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCGGCtgtg                                           <  2:95667/139‑1 (MQ=255)                                                                                         aCTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAg                                                                                                              >  1:1017629/1‑71 (MQ=255)                                                                                         aCTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAg                                                                                                              <  2:1017629/71‑1 (MQ=255)                                                                                         aCTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCGGCTGTGa                                          <  1:384397/139‑1 (MQ=255)                                                                                                    aGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCGGCTGTGAGTTCCCGgcg                                <  2:505479/138‑1 (MQ=255)                                                                                                    aGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCGGCTGTGAGTTCCCGgcg                                <  2:843471/138‑1 (MQ=255)                                                                                                    aGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGACGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCGGCTGTGAGTTCCCGgcg                                <  1:1199167/138‑1 (MQ=255)                                                                                                       aTGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCGGCTGTGAGTTCCCGGCGCTgg                            >  2:835943/1‑139 (MQ=255)                                                                                                        tGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCGGCTGTGAGTTCCCGGCGCTggg                           >  1:538975/1‑139 (MQ=255)                                                                                                        tGGCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCGGCTGTGAGTTCCCGGCGCTggg                           >  1:656320/1‑139 (MQ=255)                                                                                                         ggCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAg                                                                                                              <  2:54741/55‑1 (MQ=255)                                                                                                         ggCTGGGATGATTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAg                                                                                                              >  1:54741/1‑55 (MQ=255)                                                                                                               gatgatTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCGGCTGTGAGTTCCCGGCGCTGGGCAAACAt                    >  2:294366/1‑139 (MQ=255)                                                                                                                atgatTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCa                                                                                        <  1:187415/70‑1 (MQ=255)                                                                                                                atgatTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCa                                                                                        >  2:187415/1‑70 (MQ=255)                                                                                                                   atTTTCCCAATGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCGGCTGTGAGTTCCCGGCGCTGGGCAAACATgcg                 <  2:683959/138‑1 (MQ=255)                                                                                                                             tGATCGCCGCGTCTCTACTGGCGCGCAGCATTGATTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCGTCGGCTGTGAGTTCCCGGCGCTGGGCAAACATGCGCCAAAAACGGa      <  1:595030/139‑1 (MQ=255)                                                                                                                             tGATCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCGGCTGTGAGTTCCCGGCGCTGGGCAAACATGCGCCAAAAACGGa      <  1:595018/139‑1 (MQ=255)                                                                                                                               aTCGCCGCGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCGGCTGTGAGTTCCCGGCGCTGGGCAAACATGCGCCAAAAACGGATg    <  1:820567/139‑1 (MQ=255)                                                                                                                                  gccgcGTCTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCGGCTGTGAGTTCCCGGCGCTGGGCAAACATGCGCCAAAAACGGATGCt  <  1:405703/138‑1 (MQ=255)                                                                                                                                         |                                                                                                                                   GTCAGCGGCATCACCATGTCGATTGGTACATTCCTGGTAACACAAATGTATGAAGGGCGTCAGCGCGGTTCCCGCCTGTTATTTACCGACTCCTTCTTCAGTATGGCTGGGATGATTTTCCCAATGATCGCCGCGTTTCTACTGGCGCGCAGCATTGAGTGGTACTGGGTTTATGCCTGCATCGGGCTGGTGTATGTCGCTATTTTTATTCTGACCTTCGGCTGTGAGTTCCCGGCGCTGGGCAAACATGCGCCAAAAACGGATGCT  >  NZ_CP009273/3486251‑3486517

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