Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	684,846:1	+A	7.0%	coding (1401/2400 nt)	cadA‑II ←	cadmium translocating P‑type ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	684,846	1	.	A	7.0%	222.7 / 10.7	57	coding (1401/2400 nt)	cadA‑II	cadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base . (25/28);  new base A (2/2);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

GCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGCCAGCCAC  >  NC_002947/684699‑684987                                                                                                                                                     |                                                                                                                                              gcCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑a                                                                                                                                              <  7:245768/149‑1 (MQ=255)          gCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGAGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGt                                                                                                                                     <  4:299500/149‑1 (MQ=255)           cGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTg                                                                                                                                    <  1:250946/149‑1 (MQ=255)           cGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTg                                                                                                                                    <  5:174645/149‑1 (MQ=255)                             ttGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCgg                                                                                                                   >  1:254155/1‑148 (MQ=255)                                     cAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCa                                                                                                          <  2:122096/149‑1 (MQ=255)                                     cAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCa                                                                                                          >  5:103220/1‑149 (MQ=255)                                       gcgcCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCCCGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCAt                                                                                                         >  2:266396/1‑148 (MQ=255)                                       gcgcCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCAt                                                                                                         >  3:12436/1‑148 (MQ=255)                                        cgcCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGa                                                                                                       <  6:103220/149‑1 (MQ=255)                                          ccTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATc                                                                                                     <  5:204932/149‑1 (MQ=255)                                              gCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCgcgg                                                                                                 >  8:214368/1‑149 (MQ=255)                                                  tgAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGcggcgg                                                                                              <  8:268541/148‑1 (MQ=255)                                                  tgAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCggcggg                                                                                             >  3:306363/1‑149 (MQ=255)                                                    aaCCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCggcgggc                                                                                            <  2:205327/148‑1 (MQ=255)                                                     accaccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGgtgt                                                                                          <  4:12436/149‑1 (MQ=255)                                                     accaccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGt                                                                                          >  4:289774/1‑149 (MQ=255)                                                       caccTGCGGGCTGCCGGAGGTGAGGGTGCAGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCCTCGGCATCCTTGATC‑AGGATGCCGTGCCGGGCGGCGACCCCGGAACCGGCCATCATCGCGGCGGCCGTgg                                                                                        >  7:97112/1‑149 (MQ=255)                                                         ccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGc                                                                                       >  1:145231/1‑148 (MQ=255)                                                          cTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCGGTGGCGGGCGGCGACCCCg                                                                                                                    >  2:311469/1‑118 (MQ=255)                                                          cTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCg                                                                                                                    <  1:311469/118‑1 (MQ=255)                                                             cGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAgg                                                                                   <  2:254155/148‑1 (MQ=255)                                                              gggCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGcc                                                                                 >  8:282798/1‑149 (MQ=255)                                                              gggCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGcc                                                                                 <  4:30208/149‑1 (MQ=255)                                                                gCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGa                                                                               <  5:60534/149‑1 (MQ=255)                                                                 cTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAg                                                                              <  5:37402/149‑1 (MQ=255)                                                                   gCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAggg                                                                            <  2:283063/149‑1 (MQ=255)                                                                       gAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCa                                                                        <  2:145231/149‑1 (MQ=255)                                                                          gTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAggg                                                                     <  7:236015/149‑1 (MQ=255)                                                                            gAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAgggg                                                                    >  3:71060/1‑148 (MQ=255)                                                                               ggTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAg                                                                 <  4:71060/148‑1 (MQ=255)                                                                                   ccGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGAt                                                            <  7:282798/149‑1 (MQ=255)                                                                                     ggTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGcggcgg                                                                                              >  4:195314/1‑113 (MQ=255)                                                                                     ggTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGcggcgg                                                                                              <  3:195314/113‑1 (MQ=255)                                                                                              gAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacgg                                                 >  8:1069/1‑149 (MQ=255)                                                                                                    tACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacgg                                           >  6:4161/1‑149 (MQ=255)                                                                                                      cgcgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacggcg                                         >  2:351537/1‑149 (MQ=255)                                                                                                      cgcgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacggcg                                         >  2:30917/1‑149 (MQ=255)                                                                                                              aCCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTggcgggcg                                                                                                                             >  1:269875/1‑57 (MQ=255)                                                                                                              aCCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTggcgggcg                                                                                                                             <  2:269875/57‑1 (MQ=255)                                                                                                              aCCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAg                                 >  4:159813/1‑149 (MQ=255)                                                                                                               ccGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAgc                                <  3:289774/149‑1 (MQ=255)                                                                                                                cGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAgcg                               >  8:43226/1‑149 (MQ=255)                                                                                                                  cATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGcgcg                             >  8:253863/1‑149 (MQ=255)                                                                                                                               cAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCgg                  <  1:351537/147‑1 (MQ=255)                                                                                                                               cAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTa                <  7:253863/149‑1 (MQ=255)                                                                                                                               cAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTa                <  7:214368/149‑1 (MQ=255)                                                                                                                                aGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAc               >  2:56374/1‑149 (MQ=255)                                                                                                                                aGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAc               <  7:1069/149‑1 (MQ=255)                                                                                                                                 gggCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGTGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAcc              <  1:266396/149‑1 (MQ=255)                                                                                                                                    cTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAccc             >  6:112773/1‑147 (MQ=255)                                                                                                                                       cGGCATCC‑TGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccagc         >  5:222133/1‑148 (MQ=255)                                                                                                                                        ggCATCC‑TGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacggcg                                         >  1:70295/1‑115 (MQ=255)                                                                                                                                        ggCATCC‑TGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacggcg                                         <  2:70295/115‑1 (MQ=255)                                                                                                                                        ggCATCC‑TGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccagcca       <  6:222133/149‑1 (MQ=255)                                                                                                                                         gCATGCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccagccag      >  6:132520/1‑149 (MQ=255)                                                                                                                                          cATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCcagccagc     <  2:62863/149‑1 (MQ=255)                                                                                                                                             ccTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGccagccac  >  4:316170/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                              GCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGCCAGCCAC  >  NC_002947/684699‑684987

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