Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,413,558	G→T	S61S (TCC→TCA)  noncoding (3/179 nt)	rcbA ← ralA →	double‑strand break reduction protein small regulatory RNA antitoxin RalA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,413,558	0	G	T	100.0%	134.9 / NA	37	S61S (TCC→TCA)  noncoding (3/179 nt)	rcbA ralA	double‑strand break reduction protein small regulatory RNA antitoxin RalA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/22);  total (15/22)

AAGGAGATACTCGTTTGAAGTGGAAGCCCTCCACCCAACAGTTCTGGCGGTATGCTTCAATTTGTCTGGCCCCCAGACCAGTGCGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTT  >  NC_000913/1413314‑1413794                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             aaGGAGATACTCGTTTGAAGTGGAAGCCCTCCACCCAACAGTTCTGGCGGTATGCTTCAATTTGTCTGGCCCCCAGACCAGTGCGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGa                                                                                                                                                                                                                                        <  1:436231/251‑1 (MQ=255)      gaTACTCGTTTGAAGTGGAAGCCCTCCACCCAACAGTTCTGGCGGTATGCTTCAATTTGTCTGGCCCCCAGACCAGTGCGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACCCCCGCTGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGttt                                                                                                                                                                                                                                   >  2:42209/1‑251 (MQ=255)         aCTCGTTTGAAGTGGAAGCCCTCCACCCAACAGTTCTGGCGGTATGCTTCAATTTGTCTGGCCCCCAGACCAGTGCGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTtctc                                                                                                                                                                                                                                >  2:296066/1‑251 (MQ=255)           cttccgatcttGGAAGCCCTCCACCCAACAGTTCTGGCGGTATGCTTCAATTTGTCTGGCCCCCAGACCAGTGCGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGc                                                                                                                                                                                                                              >  2:119833/11‑251 (MQ=255)            cgttttaaggcgaatcgcTCCCCCCAACANTTCGGGCTGTATTGTTCAATTTGGCGGACCCCCAGACCAGTGGGAAGCATCAGGCCGTATTCATCCATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCTCGATGTTGATTATTAATAATCAGCTATGAAGTTTTGATTTGAATACAATGCAATTCTCGATGACTGAAGTTtcacg                                                                                                                                                                                                                               <  2:230315/232‑4 (MQ=255)               ttGAAGTGGAAGCCCTCCACCCAACAGTGCTGGCTGTATGCTTCAATTTGGCTGGCCCCCAGACCAGTGCGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGTCTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCGCGATGACTGAAGTTTCTCGCAAt                                                                                                                                                                                                                           <  1:443438/250‑1 (MQ=255)                tGAAGTGGAAGCCCTCCACCCAACAGTTCTGGCGGTATGCTTCAATTTGTCTGGCCCCCAGACCAGTGCGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAAtt                                                                                                                                                                                                                          <  1:361588/250‑1 (MQ=255)                  aaGTGGAAGCCCTCCACCCAACAGTTCTGGCGGTATGCTTCAATTTGTCTGGCCCCCAGACCAGTGCGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTaaa                                                                                                                                                                                                                       <  1:373485/251‑1 (MQ=255)                      ggAAGCCCTCCACCCAACAGTTCTGGCGGTATGCTTCAATTTGTCTGGCCCCCAGACCAGTGCGAAGCATCAGGCCGTATTCAACCATCCACTCGTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAAtt                                                                                                                                                                                                                    <  2:137536/250‑1 (MQ=255)                                          ttCTGGCGGTATGCTTCAATTTGTCTGGCCCCCAGACCAGTGCGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGctc                                                                                                                                                                                               >  1:413127/1‑251 (MQ=255)                                              ggCGGTATGCTTCAATTTGTCTGGCCCCCAGACCAGTGCGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTg                                                                                                                                                                                           <  1:296066/251‑1 (MQ=255)                                                    aTGCTTCAATTTGTATGGACCCCAGACAAGTGCGAAGCATCAGGCCGTATTAAACCATCCACTCTTCATTCAAGATTACTTGTGCTATCGCATCACCTCTGGCAGTCGCCAATGTTAGACTGCAATTGACGCCCGATGTTGATTTTTAATAATGAGCTCTGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGATGTTTCTCGCAATTAAATTTTATCAGTTTTCCTTTCTGCTCTCTGGAAAct                                                                                                                                                                                     >  2:392102/1‑250 (MQ=255)                                                         tCAATTTGTCTGGCCCCCAGACCAGTGCGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGc                                                                                                                                                                                >  2:379944/1‑251 (MQ=255)                                                         tCAATTTGTCCAGCCCCCAGACCAGTACGCAGCATCAGGCCTTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATGCTCGATGACGGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCGCTCTGGAATCGCCTg                                                                                                                                                                                 <  2:307198/250‑1 (MQ=255)                                                              ttGTCTGGCCCCCAGACCAGTGCGAGGCATCAGGCCGTATTCATCCATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATGAATAAGCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACGGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGcttct                                                                                                                                                                            <  2:325857/250‑1 (MQ=255)                                                                           aGACCAGTGCGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGagag                                                                                                                                                              <  1:99054/251‑1 (MQ=255)                                                                                    cGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTc                                                                                                                                                     <  1:176084/251‑1 (MQ=255)                                                                                                   tATTCAACCATCAACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTANCTGAGAGCATTTTTTCGCATTCTGATTTCGt                                                                                                                                      >  2:228131/1‑251 (MQ=255)                                                                                                           cATCCACTCTTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTGCCTTAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAg                                                                                                                              <  2:426900/251‑1 (MQ=255)                                                                                                                   ctTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTg                                                                                                                                                                  >  1:78946/1‑207 (MQ=255)                                                                                                                   ctTCATTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTg                                                                                                                                                                  <  2:78946/207‑1 (MQ=255)                                                                                                                       aTTAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGACATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAACATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTCGTTTAGATTTTGAATATc                                                                                                                  >  1:475143/1‑251 (MQ=255)                                                                                                                        ttAAAGATTACTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCt                                                                                                                 >  1:239421/1‑251 (MQ=255)                                                                                                                                  cTTGTGCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGtt                                                                                                        >  2:406505/1‑250 (MQ=255)                                                                                                                                     gtgCTATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTtagt                                                                                                     >  1:472617/1‑250 (MQ=255)                                                                                                                                        ctaTCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAgg                                                                                                  >  1:142177/3‑250 (MQ=255)                                                                                                                                                        ggCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCttt                                                                                  <  1:401197/250‑1 (MQ=255)                                                                                                                                                          cAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCa                                                                               <  1:9397/251‑1 (MQ=255)                                                                                                                                                                      ttAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGNAGTGCCACCTTCCTTTTCAATAGTGGCgggg                                                                   >  2:245647/1‑249 (MQ=255)                                                                                                                                                                                               tGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCatt                                           <  1:119833/250‑1 (MQ=255)                                                                                                                                                                                                              aGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGtt                           <  1:478516/251‑1 (MQ=255)                                                                                                                                                                                                               gCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTATTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTc                          <  2:331329/251‑1 (MQ=255)                                                                                                                                                                                                                            aaTTTGAATACAATGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCt             <  1:278949/251‑1 (MQ=255)                                                                                                                                                                                                                                ttaaTACAATGCAATTCTCGATTACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTATCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCAGTTTTTCGCATTCTGATTTCGGTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTc          <  2:463535/248‑1 (MQ=255)                                                                                                                                                                                                                                       aaTGCAATTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGctt  <  1:66820/251‑1 (MQ=255)                                                                                                                                                                                                                                            aaTTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTAc                                                                                                                                                                     >  1:792/1‑83 (MQ=255)                                                                                                                                                                                                                                            aaTTCTCGATGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTAc                                                                                                                                                                     <  2:792/83‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                             AAGGAGATACTCGTTTGAAGTGGAAGCCCTCCACCCAACAGTTCTGGCGGTATGCTTCAATTTGTCTGGCCCCCAGACCAGTGCGAAGCATCAGGCCGTATTCAACCATCCACTCTTCATTAAAGATTACTTGTGCCATCGCATCACCTCTGGCAGGCGCCAATGTTAGACTGAAATTGACGCCCGATGTTGATTATTAATAATCAGCTATGAAGTTTTAATTTGAATACAATGCAATTCTCGAGGACTGAAGTTTCTCGCAATTAAAATTTATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTT  >  NC_000913/1413314‑1413794

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