Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	372,646	C→G	6.7%	intergenic (+66/‑113)	dgcA → / → dgcB	dimethylglycine dehydrogenase subunit/dimethylglycine dehydrogenase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	372,646	0	C	G	6.7%	188.7 / 2.8	59	intergenic (+66/‑113)	dgcA/dgcB	dimethylglycine dehydrogenase subunit/dimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base C (31/24);  new base G (2/2);  total (33/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 = 9.61e-01

TACCGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTC  >  NC_002947/372500‑372791                                                                                                                                                     |                                                                                                                                                  tACCGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCAc                                                                                                                                                 >  4:123836/1‑149 (MQ=255)   ccGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGA‑TCAACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACt                                                                                                                                                >  8:76201/1‑148 (MQ=255)       aTCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGacaa                                                                                                                                           <  4:74662/149‑1 (MQ=255)           gCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACt                                                                                                                                       >  6:72097/1‑149 (MQ=255)             gACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATTCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTg                                                                                                                                      >  1:53231/1‑149 (MQ=255)              aCTGCGTGGCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGtt                                                                                                                                    <  7:314314/149‑1 (MQ=255)                      gCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAcc                                                                                                                            >  4:20998/1‑149 (MQ=255)                      gCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAc                                                                                                                             <  7:76201/148‑1 (MQ=255)                      gCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAc                                                                                                                             <  3:123836/148‑1 (MQ=255)                        ccAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCgg                                                                                                                          <  5:72097/149‑1 (MQ=255)                               aaCGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCAtt                                                                                                                   >  3:299486/1‑149 (MQ=25)                                  gTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTcgc                                                                                                                <  4:299486/149‑1 (MQ=21)                                     cATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg                                                                                                              <  1:326093/148‑1 (MQ=255)                                      aTGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTcgc                                                                                                                <  4:339315/145‑1 (MQ=255)                                      aTGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTcgc                                                                                                                >  3:339315/1‑145 (MQ=255)                                      aTGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGt                                                                                                            <  2:53231/149‑1 (MQ=255)                                          ggcgA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGt                                                                                                            >  4:18382/1‑145 (MQ=255)                                          ggcgA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGt                                                                                                            <  3:18382/145‑1 (MQ=255)                                              a‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGc                                                                                                    >  1:293771/1‑149 (MQ=255)                                                    cgacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTa                                                                                               <  4:84422/149‑1 (MQ=255)                                                       cgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACAt                                                                                            <  3:20998/149‑1 (MQ=255)                                                          cTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTc                                                                                                                  >  7:264508/1‑124 (MQ=18)                                                          cTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTc                                                                                                                  <  8:264508/124‑1 (MQ=18)                                                                 cTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGtc                                                                                  >  6:79031/1‑149 (MQ=255)                                                                    tgCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTAACCCCCAATTCCGAAACCAATCAAAACAATGTTGCATGACCGGACCAATACCCGGTTGAACCCGCTGTTACATGGGCCCataaac                                                                                >  1:269908/1‑144 (MQ=255)                                                                      cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc                                                                                                                        <  1:153639/106‑1 (MQ=255)                                                                      cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc                                                                                                                        >  2:153639/1‑106 (MQ=255)                                                                      cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGaa                                                                                                         >  3:168827/1‑121 (MQ=255)                                                                      cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGaa                                                                                                         <  4:168827/121‑1 (MQ=255)                                                                             ttCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGt                                                                      >  5:339444/1‑149 (MQ=255)                                                                                tGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCgg                                                                   >  8:162435/1‑149 (MQ=255)                                                                                tGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCgg                                                                   >  6:165445/1‑149 (MQ=255)                                                                                      ttAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATggg                                                                                         >  2:305663/1‑121 (MQ=255)                                                                                        aaCAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGaa                                                           >  2:206526/1‑149 (MQ=255)                                                                                             tCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAg                                                      >  5:314226/1‑149 (MQ=255)                                                                                                   ccATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGttt                                                <  1:206526/149‑1 (MQ=255)                                                                                                      tCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATct                                                        >  2:253000/1‑137 (MQ=255)                                                                                                      tCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATct                                                        <  1:253000/138‑2 (MQ=255)                                                                                                      tCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgtg                                             <  2:293771/149‑1 (MQ=255)                                                                                                         ggCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCaa                                                              <  2:7863/129‑1 (MQ=255)                                                                                                         ggCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCaa                                                              >  1:7863/1‑129 (MQ=255)                                                                                                                tCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATg                                                                                           <  1:265497/93‑1 (MQ=255)                                                                                                                tCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATg                                                                                           >  2:265497/1‑93 (MQ=255)                                                                                                                   cTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATg                                <  7:162435/149‑1 (MQ=255)                                                                                                                      tAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTg                             <  2:269908/149‑1 (MQ=255)                                                                                                                         gcgCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGACCCCGCACTTACATGGGGCCAGTCTCCCCCACTGGTCGGCGCAAGAAGCCAGCTGGTAGTCGGAGCCTCCCATGTAGAac                          >  3:346063/3‑149 (MQ=255)                                                                                                                         gAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGc                                                                                                    >  7:94749/1‑75 (MQ=255)                                                                                                                         gAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGc                                                                                                    <  8:94749/75‑1 (MQ=255)                                                                                                                         gAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAac                          >  3:53333/1‑149 (MQ=255)                                                                                                                           gCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAacac                        >  1:206488/1‑149 (MQ=255)                                                                                                                           gCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAac                          >  7:208354/1‑147 (MQ=255)                                                                                                                           gCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAac                          <  8:208354/147‑1 (MQ=255)                                                                                                                               gTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCtt                    >  8:254805/1‑149 (MQ=255)                                                                                                                                      ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgt                                              <  6:150935/117‑1 (MQ=255)                                                                                                                                      ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgt                                              >  5:150935/1‑117 (MQ=255)                                                                                                                                      ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTAcc               >  8:10445/1‑148 (MQ=255)                                                                                                                                                gTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGt    >  8:289488/1‑149 (MQ=255)                                                                                                                                                 tCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGtt   >  7:102913/1‑149 (MQ=255)                                                                                                                                                   aGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTc  <  6:314226/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TACCGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGA‑TCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTT‑ACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTC  >  NC_002947/372500‑372791

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