Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,853,215	T→C	7.4%	intergenic (+60/‑64)	relA → / → mazG	ATP:GTP 3'‑pyrophosphotransferase/modified nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,853,215	0	T	C	7.4%	150.1 / 2.8	54	intergenic (+60/‑64)	relA/mazG	ATP:GTP 3'‑pyrophosphotransferase/modified nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (23/27);  new base C (2/2);  total (25/29)
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.99e-01

CCAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGACCTAAAGCAG  >  NC_002947/1853069‑1853362                                                                                                                                                    |                                                                                                                                                        ccAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTGCCTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAAt                                                                                                                                                         <  1:80108/148‑1 (MQ=255) ccAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCt                                                                                                                                                       >  7:301621/1‑149 (MQ=255)     gCCTGGACGCGCTGGGGCGCCTGCAGGGGCGGATCTCGCAGTTGCCGAACATCATCCAGACGCGGCGTAATCGTACACCTTGAGCCCGGGGCGTATG‑CTTAGCGGGTAAACCCGCTGCTACAGGATGGAAGCAAGTTTGAAATCTGCAc                                                                                                                                                   >  4:18242/1‑149 (MQ=255)         ggACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTa                                                                                                                                               >  1:143222/1‑149 (MQ=255)              cgcTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCtgt                                                                                                                                          <  6:89701/149‑1 (MQ=255)                cTGGGGCGCCTGCTGGGGCGGATCTCGGAGTTGCGGAACATCATCGAGACGCGGCGTAATCGTAGCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAAGCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTgg                                                                                                                                        <  7:4104/149‑1 (MQ=255)                   gggCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAg                                                                                                                                     >  5:179223/1‑149 (MQ=255)                    ggCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGc                                                                                                                                    <  5:264345/149‑1 (MQ=255)                       gcCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCggg                                                                                                                                 <  6:337965/149‑1 (MQ=255)                               gggCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCAcc                                                                                                                                                  >  5:24226/1‑124 (MQ=255)                               gggCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCAcc                                                                                                                                                  <  6:24226/124‑1 (MQ=255)                                   ggATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCcgcg                                                                                                                       >  6:104314/1‑147 (MQ=255)                                   ggATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTA‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCcgcg                                                                                                                       >  8:202862/1‑147 (MQ=255)                                        tcGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGc                                                                                                                                                     <  1:329032/112‑1 (MQ=255)                                          gCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGcca                                                                                                              >  4:110157/1‑149 (MQ=255)                                          gCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTCTGGGAGCGGGCAAGCCCGCGAAGAAGcc                                                                                                               <  1:200636/148‑1 (MQ=255)                                                  cGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGTAAGCCCGCGAAGAAGCCACCA‑CCg                                                                                                       >  2:257564/1‑147 (MQ=255)                                                                  gcggcgTAATCGTAACCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGa                                                                                   <  5:104314/149‑1 (MQ=255)                                                                      cgTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGTAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTa                                                                               <  1:257564/149‑1 (MQ=255)                                                                      cgTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTa                                                                               <  6:179223/149‑1 (MQ=255)                                                                          aTCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACAcc                                                                           <  4:304184/149‑1 (MQ=255)                                                                                ccccTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCTCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTgac                                                                     <  1:262985/149‑1 (MQ=255)                                                                                   cTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTgacga‑c                                                                 >  3:152518/1‑149 (MQ=255)                                                                                      gagaCCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑cc                                                                <  8:301621/147‑1 (MQ=255)                                                                                       agaCCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑Cctgc                                                             >  2:290557/1‑149 (MQ=255)                                                                                         aCCGGGGCGTCTG‑CTTAGCGGGTTAACCCGCACCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑Cctgctg                                                           >  2:351394/1‑149 (MQ=255)                                                                                           cGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGGATGACCTACACCCTTGACGA‑CCTGCTGCa                                                         >  4:144585/1‑149 (MQ=255)                                                                                             gggCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGAcc‑tgc                                                                                            >  3:312051/1‑114 (MQ=255)                                                                                             gggCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGAcc‑tgc                                                                                            <  4:312051/114‑1 (MQ=255)                                                                                                     g‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGcc                                               >  1:124441/1‑149 (MQ=255)                                                                                                       cTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGccc                                              <  3:110157/149‑1 (MQ=255)                                                                                                             gggTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTgc                                        <  7:202862/149‑1 (MQ=255)                                                                                                                  aaCCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGAc                                   >  2:268616/1‑149 (MQ=255)                                                                                                                   aCCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGAcc                                  >  1:175994/1‑149 (MQ=255)                                                                                                                     ccGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑Cctgc                                                             <  6:363537/120‑1 (MQ=255)                                                                                                                     ccGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑Cctgc                                                             >  5:363537/1‑120 (MQ=255)                                                                                                                      cGCTCCCACAGGTACGGTGCAGATTTCAAACCTGCATCCATCCTGTAGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCg                                <  6:7749/148‑1 (MQ=21)                                                                                                                      cGCTCCCACAGGTACGGTGCAGATTTCAAACCTGCATCCATCCTGTAGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCg                                >  5:7749/1‑148 (MQ=21)                                                                                                                              cAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACgg                        <  1:198402/148‑1 (MQ=255)                                                                                                                              cAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACgg                        <  2:124441/148‑1 (MQ=255)                                                                                                                              cAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACg                         <  1:140188/147‑1 (MQ=255)                                                                                                                                  atggatgCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTgcgc                                      <  2:322032/130‑1 (MQ=255)                                                                                                                                  atggatgCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTgcgc                                      >  1:322032/1‑130 (MQ=255)                                                                                                                                     ggtgCAGATTTCAAACCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGAc                                   <  1:269937/128‑1 (MQ=255)                                                                                                                                     ggtgCAGATTTCAAACCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGAc                                   >  2:269937/3‑130 (MQ=255)                                                                                                                                      atgCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGATCCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGt                <  1:268616/149‑1 (MQ=255)                                                                                                                                      atgCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTg               >  1:48873/1‑149 (MQ=255)                                                                                                                                        gCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTgg              >  1:1020/1‑148 (MQ=255)                                                                                                                                         cAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGa            <  1:351394/149‑1 (MQ=255)                                                                                                                                         acggTTTTAAATATGCACCGTACCTGTGGTTGGGGGCAAGCCCGCGAATAATCCACCATCCCAATCCCCGGGCCACCTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGt                <  1:290557/146‑1 (MQ=255)                                                                                                                                          aGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACgg                        <  5:94663/136‑1 (MQ=255)                                                                                                                                          aGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACgg                        >  6:94663/1‑136 (MQ=255)                                                                                                                                           ggTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGAcc          <  2:175994/149‑1 (MQ=255)                                                                                                                                             tttGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGACCt         >  2:183815/1‑148 (MQ=255)                                                                                                                                                  aaTCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGACCTAAAGCa   <  2:186507/149‑1 (MQ=255)                                                                                                                                                    tCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGACCTAAAGCAg  >  4:82853/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                        CCAGGCCTGGACGCGCTGGGGCGCCTGCTGGGGCGGATCTCGCAGTTGCCGAACATCATCGAGACGCGGCGTAATCGTACCCCTTGAGACCGGGGCGTCTG‑CTTCGCGGGTAAACCCGCTCCTACAGGATGGATGCAGGTTTGAAATCTGCACCGTACCTGTGGGAGCGGGCAAGCCCGCGAAGAAGCCACCA‑CCGAAGGACC‑TGCC‑TGTATGACCTACACCCTTGACGA‑CCTGCTGCACCTCATGGCCCGCCTGCGCGACCCGCAGTACGGCTGCCCGTGGGACCTAAAGCAG  >  NC_002947/1853069‑1853362

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