Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,491,331	A→G	6.8%	T332A (ACC→GCC)	cysNC →	sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,491,331	0	A	G	6.8%	171.0 / 2.5	59	T332A (ACC→GCC)	cysNC	sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase
Reads supporting (aligned to +/- strand):  ref base A (29/26);  new base G (2/2);  total (31/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 = 9.15e-01

AGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCACCCACAAGGTCGATGTGAACACCCTGGAGCAGG  >  NC_002947/1491185‑1491478                                                                                                                                                   |                                                                                                                                                     agagCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCATGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCAcc                                                                                                                                                    <  6:46259/149‑1 (MQ=255)  gagCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCg                                                                                                                                                   >  8:83354/1‑149 (MQ=255)   agcagcCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGa                                                                                                                                                  >  1:287850/1‑149 (MQ=255)    gcagcCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGa                                                                                                                                                  <  3:133163/148‑1 (MQ=255)       gcCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAg                                                                                                                                              <  2:287850/149‑1 (MQ=255)        ccGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGt                                                                                                                                             >  1:1564/1‑149 (MQ=255)                 gTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCa                                                                                                                                    >  7:9028/1‑149 (MQ=255)                    cATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGc                                                                                                                                 >  6:247358/1‑149 (MQ=255)                      tCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTg                                                                                                                               <  8:9028/149‑1 (MQ=255)                         tCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGgtg                                                                                                                            >  4:25456/1‑149 (MQ=255)                              tACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTTGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCa                                                                                                                                    <  8:131412/136‑1 (MQ=255)                              tACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCa                                                                                                                                    >  7:131412/1‑136 (MQ=255)                                 gAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGc                                                                                                                    <  7:83354/149‑1 (MQ=255)                                    ggTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCga                                                                                                                 <  5:247358/149‑1 (MQ=255)                                           tcgAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTCGTACATCCCGACAACGTCCCGCCCGTCACCGACCAGTTc                                                                                                                                           >  8:92456/3‑116 (MQ=255)                                           tGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTc                                                                                                                                           <  7:92456/116‑1 (MQ=255)                                             gAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGAt                                                                                                        >  1:50303/1‑149 (MQ=255)                                                             gccaggccGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCGCCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGaaa                                                                                        >  3:89011/1‑149 (MQ=255)                                                                 ggccGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACg                                                                                    <  7:76921/149‑1 (MQ=255)                                                                  gccGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCGCCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGa                                                                                   >  3:59153/1‑149 (MQ=255)                                                                  gccGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACg                                                                                    >  7:201725/1‑148 (MQ=255)                                                                    cGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACa                                                                                 >  8:93255/1‑149 (MQ=255)                                                                      tGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACAt                                                                                >  5:115799/1‑148 (MQ=255)                                                                      tGACCCTGACCATGGAAGACGAGATCGACATCACCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATc                                                                               <  8:14960/149‑1 (MQ=255)                                                                         cccTGACCATTGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCaa                                                                             <  7:99207/148‑1 (MQ=255)                                                                              aCCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGa                                                                                   <  3:48834/137‑1 (MQ=255)                                                                              aCCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGa                                                                                   >  4:48834/1‑137 (MQ=255)                                                                                cATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGcca                                                                     >  1:305647/1‑149 (MQ=255)                                                                                       gACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCGCCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTa                                                              <  4:59153/149‑1 (MQ=255)                                                                                            gaTCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTCC‑GTGc                                                        >  5:76088/1‑149 (MQ=255)                                                                                              tCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCa                                                      <  8:201725/149‑1 (MQ=255)                                                                                              tCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCa                                                      <  2:305647/149‑1 (MQ=255)                                                                                                gACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCGCCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAg                                                     <  4:89011/148‑1 (MQ=255)                                                                                                      tcCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGAt                                              >  6:141890/1‑149 (MQ=255)                                                                                                        ccGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGAt                                              <  2:50303/147‑1 (MQ=255)                                                                                                            ggCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAg                                        >  2:170608/1‑149 (MQ=255)                                                                                                             gCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGc                                       <  2:1564/149‑1 (MQ=255)                                                                                                              cGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAg                                                                            <  1:280165/112‑1 (MQ=255)                                                                                                              cGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAg                                                                            >  2:280165/1‑112 (MQ=255)                                                                                                               gACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCAt                                     >  1:107611/1‑149 (MQ=255)                                                                                                                    gctgGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCAcc                                 >  1:126532/1‑148 (MQ=255)                                                                                                                    gctgGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCAcc                                 >  1:131187/1‑148 (MQ=255)                                                                                                                     ctgGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCAccc                                <  7:93255/148‑1 (MQ=255)                                                                                                                     ctgGTAAATGCCGTCAACGTCCCGCCGGTCACCTAACAGTTCGTCGCCATGGTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAATAAATACGACATCAAGCGTGCCACCAGCTACGGTGCCAGGTTCGATTGCCAGCATCAccc                                <  1:170608/149‑1 (MQ=255)                                                                                                                        gTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCACCcac                              >  4:123767/1‑147 (MQ=255)                                                                                                                        gTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCACCCACaa                            <  8:35563/149‑1 (MQ=255)                                                                                                                               ccGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCACCCACAAGGTCGAt                     <  5:141890/149‑1 (MQ=255)                                                                                                                               ccGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCACCCACAAGGTCGAt                     <  2:107611/149‑1 (MQ=255)                                                                                                                                   cAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCACCCACAAGGTCGATGTGa                 <  3:268313/149‑1 (MQ=255)                                                                                                                                       gTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAg                                        <  3:145146/122‑1 (MQ=255)                                                                                                                                       gTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAg                                        >  4:145146/1‑122 (MQ=255)                                                                                                                                         cccGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCACCCACAAGGTCGATCTGAACAccc           >  5:169051/1‑149 (MQ=255)                                                                                                                                           cgccgGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCACCCACGAGGTCGATGTGAACACCCTg         >  7:287070/1‑149 (MQ=255)                                                                                                                                             ccgGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCACCCACAAGGTCGATGTGAACACCCTg         >  8:256165/1‑147 (MQ=255)                                                                                                                                              cgGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCACCCACAAGGTCGATGTGAACACCCTGGAg      <  3:25456/149‑1 (MQ=255)                                                                                                                                                gTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCACCCACAAGGTCGATGTGAaca              <  5:193668/139‑1 (MQ=255)                                                                                                                                                gTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCACCCACAAGGTCGATGTGAaca              >  6:193668/1‑139 (MQ=255)                                                                                                                                                  cACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCACCCACAAGGTCGATGTGAACACCCTGGAGCAgg  >  4:119987/1‑149 (MQ=255)                                                                                                                                                  cACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCACCCACAAGGTCGATGTGAACACCCTGGAGCAgg  >  4:17929/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     AGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACATCTCCCGTGGCGACCTGCTGGTACATGCCGACAACGTCCCGCCGGTCACCGACCAGTTCGACGCCATGCTGGTGTGGATGGCCGAGGAGCCGATGCTCCCGGGCAAGAAATACGACATCAAGCGTGCCACCAGCTAC‑GTGCCAGGTTCGATTGCCAGCATCACCCACAAGGTCGATGTGAACACCCTGGAGCAGG  >  NC_002947/1491185‑1491478

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