Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,123,923	+A	coding (534/1026 nt)	cytR ←	DNA‑binding transcriptional repressor CytR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,123,923	1	.	A	100.0%	122.3 / NA	34	Y178? (TAT→TAN)	cytR	DNA‑binding transcriptional repressor CytR
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (10/24);  total (10/24)

GCGCCGAGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAA‑TAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCGGCGG‑CAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGATATCCGT  >  NC_000913/4123678‑4124170                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                            gcgcCGAGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAGCCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCAATACAGCCAATCCGTTTATGCCCTTGCTCATataaaataa                                                                                                                                                                                                                                                          <  2:40006/250‑5 (MQ=255)  cgcCGAGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAAtt                                                                                                                                                                                                                                                        <  1:5168/251‑1 (MQ=255)    ccGAGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCGGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGTACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTa                                                                                                                                                                                                                                                      <  2:344962/251‑1 (MQ=255)        gcgcCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGc                                                                                                                                                                                                                                                  <  1:151658/251‑1 (MQ=255)          gcCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCAt                                                                                                                                                                                                                                                <  2:2513/251‑1 (MQ=255)          gcCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCAt                                                                                                                                                                                                                                                <  2:389581/251‑1 (MQ=255)                acaTCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATg                                                                                                                                                                                                                                          <  1:355864/251‑1 (MQ=255)                 caTCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCAGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATgc                                                                                                                                                                                                                                         <  1:36560/251‑1 (MQ=255)                              gaagaCAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTg                                                                                                                                                                                                                            <  1:198459/251‑1 (MQ=255)                                                 tgtgTAAGATCAAGCCGCTGCTGCAGCGCTTTGCGTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCAGAATGCCGCAGCGACGCAGCGCCTGGACATAGCCTGGCAGGCGGTTGTGACACAGCGGCCTCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAggg                                                                                                                                                                                                         <  2:216137/250‑2 (MQ=255)                                                            aaGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAAATGTAGG‑Gagctccagctcc                                                                                                                                                                                            >  1:107212/1‑251 (MQ=255)                                                                                    tCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCAc                                                                                                                                                                    <  1:398920/250‑1 (MQ=255)                                                                                             gAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCggcgg‑c                                                                                                                                                         >  1:22279/1‑251 (MQ=255)                                                                                                         gCCACGGGCGATGTCTTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCCGATTGT‑CGATATGAACTGt                                                                                                                                                                                                             >  1:103280/1‑190 (MQ=255)                                                                                                         gCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGt                                                                                                                                                                                                             <  2:103280/190‑1 (MQ=255)                                                                                                                  gATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAAATTTAGG‑CAGCTCCAGCTCCGGTGCC‑AATTCGTTCGCCATCCCCATCGGCCG‑CAGATTACGCTGTTCCTCAATg                                                                                                                                    >  1:100898/1‑251 (MQ=255)                                                                                                                         tGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AAtt                                                                                                                                                                                 <  2:83067/200‑1 (MQ=255)                                                                                                                                              cAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCGGCGG‑CAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCa                                                                                                        <  1:288798/251‑1 (MQ=255)                                                                                                                                               ggCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGGAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGCTTATGCCCTTGCTCATATAAAATAATTTACTGCATCCAATGCGGCGGGCCGAGTGTTCGTTA‑GAAATGTAGGGAAGCCCCA‑CTCCGGTGCCCAATT‑GTTCGCCCTCACCACTGGCGGGCGAATTAC‑CGCTTCCTCCATTCTGGCATCAAACGGCAGACCTGAACCCac                                                                                                       >  2:382111/2‑250 (MQ=255)                                                                                                                                                 cGACGCAGCGCCTGAACATAGCCTTGCAGTCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCAGCACCATCGGCGG‑CAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCGTGAACCCAGca                                                                                                     <  2:327799/251‑1 (MQ=255)                                                                                                                                                   aCGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCGGCGG‑CAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGcaac                                                                                                   <  1:237687/251‑1 (MQ=255)                                                                                                                                                      cAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCGGCGG‑CAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCaacaac                                                                                                <  1:7733/251‑1 (MQ=255)                                                                                                                                                              gAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCGGCGG‑CAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATc                                                                                        <  1:185548/251‑1 (MQ=255)                                                                                                                                                                         tGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCGGCGG‑CAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTg                                                                              <  2:323136/250‑1 (MQ=255)                                                                                                                                                                             ggCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCGGCGG‑CAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGtgat                                                                         <  2:113181/251‑1 (MQ=255)                                                                                                                                                                               cGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCGGCGG‑CAGAATACGCTGTTCCTCAATGCTGGCATCAAACGGCCGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGggatg                                                                        >  2:411640/1‑250 (MQ=255)                                                                                                                                                                                                       tCGGTACCGGCTATACAGCCAATCCATTTATGCCCTTGCTCAAATAAAATAATTTACGGCAACAAATGCGGCGGTCAGATTAT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCGGCGG‑CAGATTACGCTATTCCTCAATGCTGACATCAAACGACAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTTGTGATGATCAAATCGATAAAGGTTTTTTc                                                 <  1:109013/249‑1 (MQ=255)                                                                                                                                                                                                       tCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCGGCGG‑CAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTc                                                 <  1:109011/249‑1 (MQ=255)                                                                                                                                                                                                        cGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCGGCGG‑CAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCctg                                              <  1:382111/251‑1 (MQ=255)                                                                                                                                                                                                             ccGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCGGCGG‑CAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGAt                                         <  1:263206/251‑1 (MQ=255)                                                                                                                                                                                                                 cTATACATCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGAATGG‑CGATATGAACTGTAGG‑CAGCTCCAGATCCGGTGCA‑AATTCGTTTGCCATCACCATCGGCGG‑CAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCCGCAACAACATGCCATCCATTTGCTTGGTGATGATAAAATCGATAAAGGTTTTTTCCTGCTGTTTTTg                                     >  2:257723/1‑251 (MQ=255)                                                                                                                                                                                                                      cAGCCAATCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGCTTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCGGCGG‑CAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCCGCAACAACATGCCATCCATTTTCTTGGTGATGATCAAATCGATAAAAGGTTTTTCCTGCTGATTTTGATgcg                                >  2:93492/1‑251 (MQ=255)                                                                                                                                                                                                                           aaTCCGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATTAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCGGCGG‑CAGGTTTCGCTGTTCCTCAATTCTGGCCTCAAACCGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTTCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTTCTGATTTTTATgcgccca                            >  2:298223/1‑247 (MQ=255)                                                                                                                                                                                                                               cGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATc                                                                                                                                                                <  1:398733/115‑1 (MQ=255)                                                                                                                                                                                                                               cGTTTATGCCCTTGCTCATATAAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCCCCGGTGCC‑AATTCGTTCGCCATCACCATc                                                                                                                                                                >  2:398733/1‑115 (MQ=255)                                                                                                                                                                                                                                                   aaaa‑taaTTTACTGCATCAAATGCGGCGGTCAGCTTTT‑CGATATGAACTTTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCCTCGGCGG‑CCGATTCCGCTGTTCCTCCATTCTGGCATCCAACCGCAGCCTTGAACCCCGCAACAACATGCCATCAATTTTCTTGGTGATGATCACATCGATAAAGGTTTTTTCCTGCTGCTTTTGTTGCGCACCGTCGCCCATcctcccccgatatcctc  >  2:352760/2‑235 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                            GCGCCGAGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAA‑TAATTTACTGCATCAAATGCGGCGGTCAGATTGT‑CGATATGAACTGTAGG‑CAGCTCCAGCTCCGGTGCA‑AATTCGTTCGCCATCACCATCGGCGG‑CAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGATATCCGT  >  NC_000913/4123678‑4124170

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