Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,997,194	Δ1 bp	7.3%	intergenic (‑163/+37)	PP_5252 ← / ← PP_5253	hypothetical protein/arylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,997,194	0	C	.	7.3%	150.9 / 7.1	55	intergenic (‑163/+37)	PP_5252/PP_5253	hypothetical protein/arylesterase
Reads supporting (aligned to +/- strand):  ref base C (23/26);  new base . (2/2);  total (27/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 = 4.11e-01

GGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAG  >  NC_002947/5997047‑5997338                                                                                                                                                     |                                                                                                                                                 ggggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCt                                                                                                                                                 >  8:2673/1‑149 (MQ=255) ggggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCt                                                                                                                                                 >  3:78634/1‑149 (MQ=255)   ggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTg                                                                                                                                                >  2:21865/1‑148 (MQ=255)   ggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTg                                                                                                                                                >  6:227273/1‑148 (MQ=255)   ggCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTg                                                                                                                                                >  8:262206/1‑148 (MQ=255)    gCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTa                                                                                                                                              <  7:2673/149‑1 (MQ=255)     ccGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAg                                                                                                                                             <  7:262206/149‑1 (MQ=255)     ccGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAg                                                                                                                                             <  5:57358/149‑1 (MQ=255)     ccGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAg                                                                                                                                             <  6:251028/149‑1 (MQ=255)           tggggACTGCGAAGCATGCCGGTGGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGGTGAAGAAGGGTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCTCTCCTACA‑ACGGCTGTAGGAGCgg                                                                                                                                       <  1:50/148‑1 (MQ=39)           gggggACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGT‑TGTAGGAGCgg                                                                                                                                       <  4:234760/149‑1 (MQ=255)                 cTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTAc                                                                                                                                 <  5:290435/149‑1 (MQ=255)                  tGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTAcc                                                                                                                                >  7:157932/1‑149 (MQ=255)                  tGCGAAGCATACCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCGACA‑ACGTCTGTAGGAGCGGGTTTAcc                                                                                                                                >  1:257554/1‑149 (MQ=255)                    cGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGT‑TGTAGGAGc                                                                                                                                         >  8:146320/1‑138 (MQ=255)                    cGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGT‑TGTAGGAGc                                                                                                                                         <  7:146320/138‑1 (MQ=255)                         cATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAAt                                                                                                                         <  2:42476/149‑1 (MQ=255)                             ccGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGttt                                                                                                                     <  1:198707/149‑1 (MQ=255)                              cGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTc                                                                                                                    <  7:12554/149‑1 (MQ=255)                                        aaGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTAc                                                                                                          <  8:237870/149‑1 (MQ=255)                                           gCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGc                                                                                                       <  1:21865/149‑1 (MQ=255)                                                ggCTTGAATGTAGGTGCGGTGGGTGGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGa                                                                                                    <  4:254631/147‑1 (MQ=255)                                                     gAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACg                                                                                             >  1:181009/1‑149 (MQ=255)                                                     gAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGAGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACg                                                                                             >  1:212847/1‑149 (MQ=255)                                                      aaTGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCCCTCCTACC‑AAGTCTGTAGGGGCGGGTTTACCCGCGAATGTTTGCCACAATTACCGCTGAAGGGACGc                                                                                            >  4:262216/1‑149 (MQ=38)                                                               gCGGTCGGTTGAAGAAGGTTAGCCCTGCTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTTAACCCGCATCGACA‑ACGTAGTTAGGATCGGGTTTTCACGCGAATGTTTCCCACTATTACCGCTGAAGGAACGCCAGCAAGt                                                                                    >  6:95355/1‑148 (MQ=255)                                                                 ggTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCt                                                                                 <  4:78634/149‑1 (MQ=255)                                                                 ggTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCt                                                                                 <  5:120043/149‑1 (MQ=255)                                                                   tGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCa                                                                               <  5:227273/149‑1 (MQ=255)                                                                               ggTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCTGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCtgctgc                                                                   <  6:308025/149‑1 (MQ=255)                                                                                gtaaGCCCTGGACAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGTAAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGTTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCtgctgc                                                                   <  5:52346/145‑1 (MQ=255)                                                                                gTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTgcg                                                                  <  1:175218/149‑1 (MQ=255)                                                                                    gCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCAt                                                               <  2:212847/148‑1 (MQ=255)                                                                                      ccTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCCCTCCTACA‑ACGTCAGTAGGACCGGGGTTACCCGCGAATGTTTCCCACAATTACCGCAGAACGAACGCCAGCAAGACCTCATTAAGCTGCTGCGCATGAg                                                            >  6:63869/1‑149 (MQ=38)                                                                                       cTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGTGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGt                                                           <  4:112795/149‑1 (MQ=255)                                                                                       cTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTgcgc                                                                 >  1:171502/1‑143 (MQ=255)                                                                                       cTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTgcgc                                                                 <  2:171502/143‑1 (MQ=255)                                                                                            cAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACA‑ACGTGTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCCCAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTg                                                      >  2:95594/1‑149 (MQ=255)                                                                                                           cTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGTAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgcggcg                                       <  8:157932/149‑1 (MQ=255)                                                                                                            tGGGCCGGCCTCTTCGCGGGTAAACCCGCTACTAAC‑ACGTCTGTAAGAGCGGGTTTACCCCCGAATTTTTACACCCAATACCCCTGCACGAACCCCCGCAACTCCTCATTCACCTGCTTCGCATGATTCCATGCAAAACCGTGcgccgc                                      >  2:223415/1‑149 (MQ=16)                                                                                                               gccggccTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGcccc                                   >  6:139753/1‑149 (MQ=255)                                                                                                                ccggccTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCg                                  >  6:55388/1‑149 (MQ=255)                                                                                                                        ttCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATacac                          <  2:181009/149‑1 (MQ=255)                                                                                                                                 aaaCCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTc                 >  6:89065/1‑149 (MQ=255)                                                                                                                                 aaaCCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTc                 >  2:289128/1‑149 (MQ=255)                                                                                                                                  aaCCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTAGCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGTCCCGGCATACACCTTCAGCTc                 <  5:63869/148‑1 (MQ=255)                                                                                                                                    cccGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgc              >  1:275498/1‑149 (MQ=255)                                                                                                                                     ccGCTCCTACAGACGT‑TGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcg             >  6:303210/1‑149 (MQ=255)                                                                                                                                           cTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACg        >  4:104125/1‑148 (MQ=255)                                                                                                                                            tACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCtt                       <  6:88560/132‑1 (MQ=255)                                                                                                                                            tACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCtt                       >  5:88560/1‑132 (MQ=255)                                                                                                                                            tACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGaa      >  2:208719/1‑149 (MQ=255)                                                                                                                                            tACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGa       >  2:142573/1‑148 (MQ=255)                                                                                                                                              cA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAAt     >  3:171214/1‑148 (MQ=255)                                                                                                                                                 aCGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAg  <  4:171214/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GGGGCCGGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑ACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAG  >  NC_002947/5997047‑5997338

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