Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,609,858	C→A	8.9%	K125N (AAG→AAT) ‡	PP_3183 ←	SCO1/SenC family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,609,858	0	C	A	8.9%	161.2 / 6.5	56	K125N (AAG→AAT) ‡	PP_3183	SCO1/SenC family protein
Reads supporting (aligned to +/- strand):  ref base C (22/29);  new base A (3/2);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-01

GTTGCCGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTGGGCAGGAATCGGA  >  NC_002947/3609717‑3610004                                                                                                                                                |                                                                                                                                                   gTTGCCGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTc                                                                                                                                               <  1:93959/147‑1 (MQ=255)     ccGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCa                                                                                                                                           >  5:142281/1‑147 (MQ=255)     ccGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAAAAGTTCGGCAtc                                                                                                                                         <  1:119200/149‑2 (MQ=255)     ccGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGGTCGCCGGG‑CAGCAAGGGCCAAGCCTGGGCCGATGCCGAAAAGTTCGGCAt                                                                                                                                          <  1:189230/149‑1 (MQ=255)      cGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAg                                                                                                                                        >  8:227440/1‑149 (MQ=255)      cGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAg                                                                                                                                        >  2:179455/1‑149 (MQ=255)      cGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAg                                                                                                                                        >  1:303592/1‑149 (MQ=255)         tgatgaGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGt                                                                                                                                     <  3:193545/149‑1 (MQ=255)           atgaGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGttt                                                                                                                                   <  1:18200/149‑1 (MQ=255)              aGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGa                                                                                                                                 <  6:242142/148‑1 (MQ=255)                       ttGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTTGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCggg                                                                                                                       <  1:227161/149‑1 (MQ=255)                         gtgGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCgggg                                                                                                                      <  7:227440/148‑1 (MQ=255)                          gtgTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGTCCTGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCgggg                                                                                                                      <  1:256409/146‑1 (MQ=255)                                 tGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCg                                                                                                                <  6:238707/146‑1 (MQ=255)                                  ggAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTta                                                                                                             >  6:31140/1‑148 (MQ=255)                                  ggAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTta                                                                                                             <  6:165016/148‑1 (MQ=255)                                     gCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAgg                                                                                                         <  1:119480/149‑1 (MQ=255)                                          cGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATc                                                                                                     >  1:267231/1‑148 (MQ=255)                                               tCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGc                                                                                                <  6:57974/148‑1 (MQ=255)                                                 cAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGa                                                                                             >  6:273412/1‑149 (MQ=255)                                                     ccATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATggcctgga                                                                                          >  5:33731/1‑141 (MQ=255)                                                     ccATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGa                                                                                         <  4:314025/149‑1 (MQ=255)                                                      cATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAg                                                                                        >  2:279726/1‑149 (MQ=255)                                                        tCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAg                                                                                        >  3:227877/1‑147 (MQ=255)                                                          gATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTa                                                                                      >  6:114972/1‑147 (MQ=255)                                                           aTCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTagag                                                                                   <  7:96949/149‑1 (MQ=255)                                                            tCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAgg                                                                                  <  8:30266/149‑1 (MQ=255)                                                              cacaGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAgg                                                                                  >  8:118582/1‑147 (MQ=255)                                                                 aGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTATCGGCATAGCGTTTGAGGGGGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGAt                                                                             <  2:70052/149‑1 (MQ=255)                                                                      cAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTTGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCtt                                                                        <  2:267231/149‑1 (MQ=255)                                                                              ggcgTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCg                                                                                                                >  3:192689/1‑101 (MQ=255)                                                                              ggagTAGTTGCTCGATGTCGTCGGGTTGGGAGGT‑GAGCAAGGTCC‑AGCCTGGTGCTATGCCGATCTTTTCGTCATAGCTTTTGAGTGTGGCGGGGTTGTCg                                                                                                                <  4:192689/98‑1 (MQ=255)                                                                                cgTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCg                                                              <  2:303592/149‑1 (MQ=255)                                                                                         tCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCAGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGACGACCCGGTCACagag                                                     >  8:304775/1‑149 (MQ=255)                                                                                             tGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGAtt                                                 <  5:273412/149‑1 (MQ=255)                                                                                               tcgtcgGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGAtttt                                               <  5:114972/149‑1 (MQ=255)                                                                                                 gtcgGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGC‑GGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCt                                             <  6:142281/148‑1 (MQ=255)                                                                                                          gCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGc                                    >  4:302121/1‑149 (MQ=255)                                                                                                              gT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAg                                >  4:84478/1‑149 (MQ=255)                                                                                                                   gCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGgcgc                            >  8:182097/1‑149 (MQ=255)                                                                                                                       ggTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAg                        >  7:8975/1‑149 (MQ=255)                                                                                                                          cc‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGAGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTATGTACCTGGCGTAGGAGCGCAGtt                      >  8:279117/1‑148 (MQ=255)                                                                                                                           c‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTc                    >  3:159511/1‑149 (MQ=255)                                                                                                                               ccTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAg                                                                               <  2:111672/87‑1 (MQ=255)                                                                                                                               ccTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAg                                                                               >  1:111672/1‑87 (MQ=255)                                                                                                                                cTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCAc                 <  5:4868/148‑1 (MQ=255)                                                                                                                                   ggCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGTTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTg               <  1:279726/147‑1 (MQ=255)                                                                                                                                   ggCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTg               <  5:31140/147‑1 (MQ=255)                                                                                                                                     gctaTGCCGAAAAGTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGAtt                                                 >  6:107519/4‑111 (MQ=255)                                                                                                                                      cGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTGGGCAg          <  3:84478/149‑1 (MQ=255)                                                                                                                                       taTGCCGAAAAGTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTGGGCAg          >  2:327028/2‑148 (MQ=255)                                                                                                                                       taTGCCGAAAAGTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTGGGCAg          >  5:30451/2‑148 (MQ=255)                                                                                                                                           ccGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGGGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTGGGCAGGAATc     <  7:118582/149‑1 (MQ=255)                                                                                                                                             gAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGtt                      >  1:1177/1‑130 (MQ=255)                                                                                                                                             gAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGtt                      <  2:1177/130‑1 (MQ=255)                                                                                                                                              aaCTTTTCGGCATAGGGTTTGAGGTGGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTATAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTGGGCAGGAATctta  <  8:8975/149‑4 (MQ=255)                                                                                                                                                |                                                                                                                                                   GTTGCCGATGATGAGGCTCAGGTTGTGGTCTTTGGAGCGGCCGTTTTCCAGGCCATCGATCCACAGCCCCAGGCTGCGGCGTAGTTGCTCGATGTCGTCGGGTTCGCCGGT‑CAGCAAGGTCC‑AGCCTGGGCCGATGCCGAACTTTTCGGCATAGCGTTTGAGGGTGGCGGGGGTGTCGTTATAGGGATCGATGCTGATGGAGTAGAGGAAGATGTCTTTGCCGACCCGGTCACCGAGGATTTTCTGTACCTGGCGTAGGCGCGCAGTTTCCACTGGGCAGGAATCGGA  >  NC_002947/3609717‑3610004

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