Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,873,556	T→C	intergenic (+25/‑18)	dgcJ → / → yeaK	diguanylate cyclase DgcJ/mischarged aminoacyl‑tRNA deacylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,873,556	0	T	C	100.0%	126.2 / NA	37	intergenic (+25/‑18)	dgcJ/yeaK	diguanylate cyclase DgcJ/mischarged aminoacyl‑tRNA deacylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/16);  total (21/16)

CGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGTAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAAC  >  NC_000913/1873312‑1873799                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    cGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAAtctc                                                                                                                                                                                                                                               >  2:5332/1‑251 (MQ=255)             catcatTCTTGTCGATTCGAAGCCGCAAATTGCAGCACAACTGCCTGAACGGATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGGCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGaaa                                                                                                                                                                                                                                   <  2:330302/251‑1 (MQ=255)                   tCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACt                                                                                                                                                                                                                              <  2:114951/250‑1 (MQ=255)                           aTTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCTTGAAATGACTGAAATGGCt                                                                                                                                                                                                                     <  1:279044/251‑1 (MQ=255)                           aTTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCt                                                                                                                                                                                                                     <  1:368458/251‑1 (MQ=255)                               gACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTaa                                                                                                                                                                                                                   >  1:86149/1‑249 (MQ=255)                                               gcaCAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCa                                                                                                                                                                                                 <  1:306916/251‑1 (MQ=255)                                                    aCTGCCTGAACGTATAGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGt                                                                                                                                                                                            >  1:89105/1‑251 (MQ=255)                                                            aaCGTATCGAAAACCGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTATTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCCGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCt                                                                                                                                                                                    >  2:174193/1‑251 (MQ=255)                                                              cGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCttt                                                                                                                                                                                  >  2:313993/1‑251 (MQ=255)                                                                tATCGAAAAACGTCTGCCACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTttat                                                                                                                                                                                >  2:156779/1‑251 (MQ=255)                                                                    gAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATc                                                                                                                                                                            <  2:194051/251‑1 (MQ=255)                                                                               gCAACATATCGCGCCGCAGAAAGATATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGt                                                                                                                                                                  <  1:247012/250‑1 (MQ=255)                                                                                 aaCATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTCTTATCACAAGAAGGTGCt                                                                                                                                                               >  2:143279/1‑251 (MQ=255)                                                                                      atCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGGATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACtt                                                                                                                                                          <  1:143279/251‑1 (MQ=255)                                                                                                 gAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGtt                                                                                                                                                >  1:203914/1‑250 (MQ=255)                                                                                                  aaaGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTAc                                                                                                                                              >  1:85555/1‑251 (MQ=255)                                                                                                                   ttCCGGAATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCggtaggt                                                                                                                              >  2:294732/1‑250 (MQ=255)                                                                                                                                     gaaagaaaACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTc                                                                                                             >  1:172995/1‑249 (MQ=255)                                                                                                                                               gATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTgg                                                                                                 >  1:145227/1‑251 (MQ=255)                                                                                                                                                   cGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAACGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCAcc                                                                                             >  1:123911/1‑251 (MQ=255)                                                                                                                                                      tACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGct                                                                                          >  1:257787/1‑251 (MQ=255)                                                                                                                                                                   aaaaGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTTCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTg                                                                             <  2:178481/250‑1 (MQ=255)                                                                                                                                                                     aaGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCt                                                                           <  1:5332/251‑1 (MQ=255)                                                                                                                                                                          tCCGATGAGCGTTTATATGTCAATAAGGAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCTCTCTAGGCCAGGGTGCTAAgg                                                                       <  2:354565/250‑1 (MQ=255)                                                                                                                                                                           ccGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGc                                                                      <  1:313993/250‑1 (MQ=255)                                                                                                                                                                            cGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCAc                                                                    >  2:374973/1‑251 (MQ=255)                                                                                                                                                                                    gTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGt                                                            >  2:211885/1‑251 (MQ=255)                                                                                                                                                                                    gTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGATTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGCAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGTTCTAGTCCAGGGTGCTAAGGAACTGGTATGt                                                            >  2:211875/1‑251 (MQ=255)                                                                                                                                                                                      ttATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTa                                                           <  2:145227/250‑1 (MQ=255)                                                                                                                                                                                          atGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGt                                                       <  1:119831/250‑1 (MQ=255)                                                                                                                                                                                            gTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCa                                                     <  2:89105/250‑1 (MQ=255)                                                                                                                                                                                                           aaaaaCAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGc                                                <  1:16848/240‑1 (MQ=255)                                                                                                                                                                                                           aaaaaCAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGc                                                >  2:16848/1‑240 (MQ=255)                                                                                                                                                                                                            aaaaCAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCGTCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAAt                                    >  1:260838/1‑251 (MQ=255)                                                                                                                                                                                                                     gtttaTGATAACCTTCTGTGGATGTGTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATGAGCGGTTAATCGCTTTATTATCACAGGAAGGTGCTGACTTCCGCCTCGTTACTCATTAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCTCTCTAGGCCAGGGTGCTAAGGCGCTGGTATGTGAAGGCaaaggcaaaggctcaactgaacatgtt                           >  2:68896/5‑236 (MQ=255)                                                                                                                                                                                                                                      gtgGTTGTTTGCTTGCAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGc          >  1:291715/1‑251 (MQ=255)                                                                                                                                                                                                                                              ttgcttgcAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAAc  >  1:198475/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    CGATGAATTCTGCATCATTCTTGTCGATTCGACGCCGCAAATTGCAGCACAACTGCCTGAACGTATCGAAAAACGTCTGCAACATATCGCGCCGCAGAAAGAGATCGGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAAAAACAGCCGTTCATGATAACCTTCTGTGGTTGTTTGCTTGTAATCTCAGGAGCGTGAAATGACTGAAATGGCTAAAGGAAGCGTGACGCATCAGCGGTTAATCGCTTTATTATCACAAGAAGGTGCTGACTTCCGCGTCGTTACCCATGAAGCGGTAGGTAAATGTGAAGCGGTGTCAGAAATTCGTGGCACCGCTCTAGGCCAGGGTGCTAAGGCACTGGTATGTAAAGTCAAAGGCAATGGCGTAAATCAACATGTTCTGGCGATCCTCGCAGCCGATCAAC  >  NC_000913/1873312‑1873799

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