Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	284,586	T→A	8.5%	intergenic (+33/+30)	betT‑II → / ← atsK	choline/carnitine/betaine transporter/alpha‑ketoglutarate‑dependent sulfonate dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	284,586	0	T	A	8.5%	151.6 / 3.6	56	intergenic (+33/+30)	betT‑II/atsK	choline/carnitine/betaine transporter/alpha‑ketoglutarate‑dependent sulfonate dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (25/26);  new base A (2/3);  total (27/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 = 9.89e-01

ACCTGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATT  >  NC_002947/284438‑284729                                                                                                                                                      |                                                                                                                                                aCCTGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCt                                                                                                                                                 <  3:77143/149‑1 (MQ=255)   cTGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTg                                                                                                                                               <  7:85857/149‑1 (MQ=255)        cGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCcgc                                                                                                                                          >  7:132600/1‑149 (MQ=255)           aGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                       <  6:137134/149‑1 (MQ=255)             gCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCCGCGAaga                                                                                                                                     <  8:132600/149‑1 (MQ=255)               aGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCCGCGAAGAgg                                                                                                                                   <  5:293868/149‑1 (MQ=255)                gCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                  <  2:229228/149‑1 (MQ=255)                gCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                  <  4:177937/149‑1 (MQ=255)                  cAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGc                                                                                                                                  >  4:285153/1‑147 (MQ=255)                   aGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCGCTCCTGTTCTGGCCTCTTCGCGGGCAAGCCCGCGAAGAGGCCAg                                                                                                                               <  8:86259/149‑1 (MQ=255)                        tATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAgg                                                                                                                          >  5:45738/1‑149 (MQ=255)                             ccTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAAc                                                                                                                      >  1:86944/1‑149 (MQ=255)                               tggtggGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTAGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGc                                                                                                                   >  2:255164/1‑149 (MQ=255)                               tggtggGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTATTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGc                                                                                                                   >  5:113907/1‑149 (MQ=255)                               tggtggGCTACAGCAAGGAGCACATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGCTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGCGCTTGCCCGCGAAGCGGCCAGACAGGCAACCg                                                                                                                    >  6:136775/1‑148 (MQ=255)                                       tACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCaa                                                                                                                       >  3:307214/1‑137 (MQ=255)                                         cAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTg                                                                                                         >  5:169718/1‑149 (MQ=255)                                            cAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCAt                                                                                                      >  4:250319/1‑149 (MQ=255)                                              agtcgCATATCATCAACCACATCGTCTACCAGTACTT‑GGGGAACTTGCTGTACGTCCACCTGGTCGGCTTTCAAAGGTCTCTTGTCGGGGCCTCTTGGGGGGCTGCCCGGCTATTTGGCCAGACAGGCAACCGCTAGAAAGGTGCATcc                                                                                                    <  4:29606/145‑1 (MQ=255)                                                  gCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCGCTTAGCGGGCATGGCCCCGCACAGGCCAGCCAAGCACCCCCGAGAAAAGTGCaacaccc                                                                                                 >  5:223920/1‑142 (MQ=255)                                                      atcatcAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGAt                                                                                            >  5:74831/1‑149 (MQ=255)                                                         atcaACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGt                                                                                         <  1:255164/149‑1 (MQ=255)                                                          tcaACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTc                                                                                        >  5:336920/1‑149 (MQ=255)                                                           caACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAgg                                                                                              <  6:110574/142‑1 (MQ=255)                                                           caACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAgg                                                                                              >  5:110574/1‑142 (MQ=255)                                                                 cATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGt                                                                                 <  8:223534/149‑1 (MQ=255)                                                                    ccTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGa                                                                                                                                        >  2:254950/1‑91 (MQ=255)                                                                    ccTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGa                                                                                                                                        <  1:254950/91‑1 (MQ=255)                                                                          ccAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGc                                                                                                        >  5:273214/1‑117 (MQ=255)                                                                          ccAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGc                                                                                                        <  6:273214/117‑1 (MQ=255)                                                                           cAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAg                                                                                               >  1:120344/1‑125 (MQ=255)                                                                           cAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACccg                                                                       >  7:84219/1‑149 (MQ=255)                                                                             gTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACccgcc                                                                     >  8:85476/1‑149 (MQ=255)                                                                                 gA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGGTGGTCGCCCGGTGCATCACCCCCCgggt                                                                 >  6:27354/1‑147 (MQ=255)                                                                                     cGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcg                                                               <  6:113907/148‑1 (MQ=255)                                                                                        cacaTGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcggtcg                                                           <  6:169718/149‑1 (MQ=255)                                                                                         acaTGCAGTACCTGCACCGGGTCCGCTGACCAAGCTCGCCTGTCCTGGCCGCTGGGGGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcggtcgg                                                          <  5:27354/149‑1 (MQ=255)                                                                                            tccAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGCCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTaa                                                       <  3:161286/147‑1 (MQ=255)                                                                                                    gagagcttGGTCAGCGGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTaa                                                       <  1:16653/134‑1 (MQ=255)                                                                                                    gagagcttGGTCAGCGGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTaa                                                       >  2:16653/8‑141 (MQ=255)                                                                                                        aCCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGCGGCCAGACAGGCCGCCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGc                                            >  1:288949/1‑148 (MQ=255)                                                                                                                cGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTg                                   >  3:83947/1‑149 (MQ=255)                                                                                                                 gCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTgg                                  <  6:336920/149‑1 (MQ=255)                                                                                                                     aCCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCa                              >  3:115505/1‑149 (MQ=255)                                                                                                                      ccAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCAc                             <  3:285153/149‑1 (MQ=255)                                                                                                                       cAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCAcg                            >  8:302088/1‑149 (MQ=255)                                                                                                                        aaGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCAcgc                           <  5:171755/149‑1 (MQ=255)                                                                                                                          gCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCgg                         <  4:83947/149‑1 (MQ=255)                                                                                                                          gCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCgg                         <  2:86944/149‑1 (MQ=255)                                                                                                                           ctgtcgtgtccggGCCTCTTCGGGGGGTTCCCGGCTATGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCgg                         <  6:45738/142‑1 (MQ=255)                                                                                                                           ctctCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATc                                                                            <  4:239236/97‑1 (MQ=255)                                                                                                                           ctctCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATc                                                                            >  3:239236/1‑97 (MQ=255)                                                                                                                                   tcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCa                 <  3:300397/148‑1 (MQ=255)                                                                                                                                   tcctgGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCa                 <  3:250319/148‑1 (MQ=255)                                                                                                                                        gCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCTCCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGc           <  5:136775/149‑1 (MQ=255)                                                                                                                                             ttCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGt      >  5:299392/1‑149 (MQ=255)                                                                                                                                                 cgGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCAtt  <  7:85476/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                ACCTGAGCGAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGA‑GCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATT  >  NC_002947/284438‑284729

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