Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,951,334	G→C	16.2%	intergenic (‑113/+48)	oguA ← / ← PP_2585	8‑oxoguanine deaminase/alpha‑ketoglutaric semialdehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,951,334	0	G	C	16.2%	124.6 / 18.0	55	intergenic (‑113/+48)	oguA/PP_2585	8‑oxoguanine deaminase/alpha‑ketoglutaric semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base C (5/4);  total (29/26)
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.44e-01

GAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTC  >  NC_002947/2951189‑2951481                                                                                                                                                  |                                                                                                                                                     gaggagCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTggcg                                                                                                                                                   >  3:301505/1‑149 (MQ=255)             tGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCcgc                                                                                                                                        <  7:267121/148‑1 (MQ=255)             tGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCcgc                                                                                                                                        <  6:115068/148‑1 (MQ=255)                  ggTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAagag                                                                                                                                  <  5:322651/149‑1 (MQ=255)                    tttGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGc                                                                                                                                >  8:228113/1‑149 (MQ=255)                       gTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCg                                                                                                                              <  1:226816/148‑1 (MQ=255)                                  ggTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTg                                                                                                                   >  5:22339/1‑147 (MQ=255)                                                cAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAATCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATca                                                                                                   <  3:112263/149‑1 (MQ=255)                                                   tACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCacac                                                                                                <  7:228113/149‑1 (MQ=255)                                                       cATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATc                                                                                                    >  2:151127/1‑141 (MQ=255)                                                       cATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATc                                                                                                    <  1:151127/141‑1 (MQ=255)                                                              gATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGgc                                                                                     <  4:301505/149‑1 (MQ=255)                                                                ttACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGgctg                                                                                   <  4:300100/149‑1 (MQ=255)                                                                    cGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATcaca                                                                                                 >  2:302851/1‑131 (MQ=255)                                                                    cGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATcaca                                                                                                 <  1:302851/131‑1 (MQ=255)                                                                          cATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGAATGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCt                                                                         >  1:66639/1‑149 (MQ=255)                                                                           aTAGGCTCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTc                                                                        <  8:99017/149‑1 (MQ=255)                                                                              ggCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCa                                                                      <  6:22339/148‑1 (MQ=255)                                                                                     ggTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCgg                                                              >  8:186881/1‑149 (MQ=255)                                                                                           gTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGAt                                                                                                     <  8:94204/104‑1 (MQ=255)                                                                                           gTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGAt                                                                                                     >  7:94204/1‑104 (MQ=255)                                                                                              agcaagcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGCATGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATccc                                                      <  1:296843/149‑1 (MQ=255)                                                                                                caagcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCAATTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGAAAACACCCCCTCCCGGCTGCTCACACCCTCCAGCCACCGGGCAAACCCCac                                                   >  6:245991/1‑148 (MQ=255)                                                                                                 aagcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCCCCCCCTCCAGCCACCGGGCAATCCCCAGc                                                  >  8:138744/1‑149 (MQ=255)                                                                                                   gcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAGCCGGGCAATCCCCAGCgg                                                <  2:52157/149‑1 (MQ=255)                                                                                                     tttttACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGt                                              >  2:64936/1‑149 (MQ=255)                                                                                                                 gTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCgg                                                >  8:254257/1‑135 (MQ=255)                                                                                                                  tataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGt                                              >  6:321924/1‑136 (MQ=255)                                                                                                                  tataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGt                                              <  5:321924/136‑1 (MQ=255)                                                                                                                    taCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGGAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGGGCACCCCCACCAGCAACCGGGCAATCCCCACCGGGTTAGCGGTTTTCAACg                               >  2:354551/1‑149 (MQ=255)                                                                                                                    taCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAAc                                <  7:186881/148‑1 (MQ=255)                                                                                                                     aCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACg                               >  4:186222/1‑148 (MQ=255)                                                                                                                     aCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCGGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGc                              >  8:298458/1‑149 (MQ=255)                                                                                                                               gTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGgcag                     >  7:229847/1‑148 (MQ=255)                                                                                                                               gTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGgcag                     >  3:217075/1‑148 (MQ=255)                                                                                                                                  cccGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAAc                                <  6:208101/134‑1 (MQ=255)                                                                                                                                  cccGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAAc                                >  5:208101/1‑134 (MQ=255)                                                                                                                                  cccGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGC‑GCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCaa                                 <  8:11373/132‑1 (MQ=255)                                                                                                                                  cccGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGC‑GCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCaa                                 >  7:11373/1‑132 (MQ=255)                                                                                                                                     gTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCg              >  4:152928/1‑149 (MQ=255)                                                                                                                                        ccGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCa           <  7:298458/149‑1 (MQ=255)                                                                                                                                        ccGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCa           <  1:64936/149‑1 (MQ=255)                                                                                                                                          ggCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAgg         <  5:245991/149‑1 (MQ=255)                                                                                                                                           gCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAggg        >  7:52601/1‑149 (MQ=255)                                                                                                                                            ccTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGt       <  3:152928/149‑1 (MQ=255)                                                                                                                                             ctctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAggg        >  4:354054/1‑147 (MQ=255)                                                                                                                                             ctctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCaa                                 <  6:153272/122‑1 (MQ=255)                                                                                                                                             ctctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCaa                                 >  5:153272/1‑122 (MQ=255)                                                                                                                                             ctctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAggg        >  4:354740/1‑147 (MQ=255)                                                                                                                                             ctctTCGCGGCTATAGCCGCCAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTa      >  6:327772/1‑149 (MQ=255)                                                                                                                                              tctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGtttt                                     <  5:24922/117‑1 (MQ=38)                                                                                                                                              tctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGtttt                                     >  6:24922/1‑117 (MQ=38)                                                                                                                                              tctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCa           >  5:68263/1‑143 (MQ=255)                                                                                                                                              tctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCa           <  6:68263/143‑1 (MQ=255)                                                                                                                                                ttGGCGGCTTGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGtt   >  5:149625/1‑149 (MQ=255)                                                                                                                                                 tGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTc  >  4:368238/1‑149 (MQ=255)                                                                                                                                                  ggcggcTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGtt   <  8:52601/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGG‑TTGCAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTC  >  NC_002947/2951189‑2951481

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