Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	684,846:1	+A	7.1%	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.1%	217.8 / 10.2	56	coding (1401/2400 nt)	cadA‑II	cadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base . (26/26);  new base A (2/2);  total (28/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 = 8.39e-01

GCCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGCCAGCCAC  >  NC_002947/684699‑684987                                                                                                                                                     |                                                                                                                                              gcCGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑a                                                                                                                                              <  6:277914/149‑1 (MQ=255)   cGGCCAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑Agg                                                                                                                                            <  4:118830/149‑1 (MQ=255)      ccAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATg                                                                                                                                         <  3:220186/149‑1 (MQ=255)       cAGGCGGTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGc                                                                                                                                        >  5:177905/1‑149 (MQ=255)            ggTAAAGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCC‑TGATCAAGGATGCCgaag                                                                                                                                   <  1:131782/149‑4 (MQ=255)                 aGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTg                                                                                                                                    >  4:185284/1‑143 (MQ=255)                 aGGTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTg                                                                                                                                    <  3:185284/143‑1 (MQ=255)                   gTCGGCGCTGTTGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGgcgggcgg                                                                                                                            <  2:152861/149‑1 (MQ=255)                             ttGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCgg                                                                                                                   >  7:173110/1‑148 (MQ=255)                             ttGCCGTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCgg                                                                                                                   >  4:166195/1‑148 (MQ=255)                                  gTCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCgg                                                                                                             >  6:264395/1‑149 (MQ=255)                                   tCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGc                                                                                                            >  5:64079/1‑149 (MQ=255)                                   tCCAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGc                                                                                                            >  7:150154/1‑149 (MQ=255)                                     cAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGgcgggcgg                                                                                                                            <  5:3192/131‑1 (MQ=255)                                     cAGCGCCTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGgcgggcgg                                                                                                                            >  6:3192/1‑131 (MQ=255)                                          ccTGGCTGTGAACCACCTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATc                                                                                                     >  3:63301/1‑149 (MQ=255)                                                      ccaccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTg                                                                                         >  5:197995/1‑149 (MQ=255)                                                       caccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGc                                                                                                            <  6:184594/129‑1 (MQ=255)                                                       caccTGCGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGc                                                                                                            >  5:184594/1‑129 (MQ=255)                                                             cGGGCTGCCGGAGGTGAGGGTGCCGGTTTTGTCTAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGc                                                                                  <  6:177905/149‑1 (MQ=255)                                                                         ggTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAgg                                                                      <  4:164745/149‑1 (MQ=255)                                                                          gTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAggg                                                                     >  5:39496/1‑149 (MQ=255)                                                                          gTGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAggg                                                                     <  6:197995/149‑1 (MQ=255)                                                                           tGAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAgggc                                                                    >  1:262556/1‑148 (MQ=255)                                                                            gAGGGTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAgggg                                                                    >  8:130525/1‑148 (MQ=255)                                                                               ggTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAgg                                                                <  7:130525/149‑1 (MQ=255)                                                                                gTGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGc                                                               <  4:63301/149‑1 (MQ=255)                                                                                 tGCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCg                                                              <  4:23671/149‑1 (MQ=255)                                                                                  gCCGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGa                                                             >  6:136270/1‑149 (MQ=255)                                                                                   ccGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGAt                                                            <  5:264395/149‑1 (MQ=255)                                                                                    cGGTTTTGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATg                                                           <  3:166195/149‑1 (MQ=255)                                                                                          tGTCGAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAg                                                      <  6:39496/148‑1 (MQ=255)                                                                                              gAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGGCGAGGGCGCAGGGGCAGGCGATGACCAGCacgg                                                 <  8:173110/149‑1 (MQ=255)                                                                                              gAACACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATCATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacgg                                                 >  5:224633/1‑149 (MQ=255)                                                                                               aaCACTACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacgg                                                 <  1:120302/148‑1 (MQ=255)                                                                                                    tACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacgg                                           >  4:239119/1‑149 (MQ=255)                                                                                                    tACGCGGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacg                                            >  7:71144/1‑148 (MQ=255)                                                                                                      cgcgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCacggcgacg                                            <  8:71144/146‑1 (MQ=255)                                                                                                       gcgGTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGt                                        >  4:73599/1‑149 (MQ=255)                                                                                                         ggTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGAcc                                                        <  2:98045/131‑1 (MQ=255)                                                                                                         ggTTGACCGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGAcc                                                        >  1:98045/1‑131 (MQ=255)                                                                                                               ccGCATGGGCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAgc                                >  5:38240/1‑149 (MQ=255)                                                                                                                     gggCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGa                                                                               <  2:25809/96‑1 (MQ=255)                                                                                                                     gggCACGCTCCAGGGCTTCGGCATCCTTGAAC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGa                                                                               >  1:25809/1‑96 (MQ=255)                                                                                                                      ggCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGttt                         >  2:25349/1‑149 (MQ=255)                                                                                                                       gCACGCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGTCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGTGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTc                        <  3:239119/149‑1 (MQ=255)                                                                                                                           gCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGc                    <  6:64079/149‑1 (MQ=255)                                                                                                                           gCTCCAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGc                    <  8:150154/149‑1 (MQ=255)                                                                                                                               cAGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTa                >  4:281766/1‑149 (MQ=255)                                                                                                                                aGGGCTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGCTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTAc               >  5:162498/1‑149 (MQ=255)                                                                                                                                    cTTCGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCcca           <  5:136270/149‑1 (MQ=255)                                                                                                                                     cacGGCATCC‑TGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccag          >  7:257593/3‑149 (MQ=255)                                                                                                                                      acGGCATCC‑TGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGCGGCAGGCGATGACCAGCACGGCGACGGCGATGATCAGCGCGGTTTCCAGCGCTACAccaga         >  8:196393/2‑148 (MQ=255)                                                                                                                                       cGGCATCCTTGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccagc         <  7:200776/148‑1 (MQ=255)                                                                                                                                       cGGCATCCTTGATC‑AGGATGCAGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGAAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccagc         <  7:197225/148‑1 (MQ=255)                                                                                                                                         gCATCC‑TGATCAAGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCccagccag      <  8:257593/149‑1 (MQ=255)                                                                                                                                               ttGATC‑AGGATGCCGTGGCGGGCGGCGACCCCGGTACCGGCCATGATCGCGGCGGGCGTGGCCAGGCCGAGGGCGCAGGGGCAGGCGATGACCAGCACGGCGACGGCGTTGATCAGCGCGGTTTCCAGCGGTACCCCAGccagccac  >  7:63634/1‑147 (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: ‑