Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	881,531	A→T	I269F (ATC→TTC)	dacC →	D‑alanyl‑D‑alanine carboxypeptidase DacC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	881,531	0	A	T	100.0%	126.4 / NA	39	I269F (ATC→TTC)	dacC	D‑alanyl‑D‑alanine carboxypeptidase DacC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (21/18);  total (21/18)

GGCATTGCTGGGTAAAGCATTGATCCACGATGTGCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTATCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTTATACGTTAACGGAACCGCAGCTTACCGCACCGCTGAAAAAAGGTCA  >  NC_000913/881287‑881772                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  ggCATTGCTGGGTAAAGCATTGATCCACGATGTGCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGtt                                                                                                                                                                                                                                             <  2:281355/251‑1 (MQ=255)      tGCTGGGTAAAGCATTGATCCACGATGTGCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAAt                                                                                                                                                                                                                                        >  1:447928/1‑251 (MQ=255)                      gATCCACGATGTGCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATTGCATGAAGACATGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTGTAAGCTCCGTAGTGGTTGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGACTGACaaagta                                                                                                                                                                                                                         <  2:80271/250‑4 (MQ=255)                      gATCCACGATGTGCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAGCCGACCGTTTCCGTTTTAATGAGTCTGAGAAatta                                                                                                                                                                                                                         <  1:129430/250‑1 (MQ=255)                         ccACGATGTGCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGACATTATTg                                                                                                                                                                                                                      >  1:440054/1‑250 (MQ=255)                              aTGTGCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTgg                                                                                                                                                                                                                <  1:504234/251‑1 (MQ=255)                                   ccGGAAGAGTACACCATTCATAAAGAGAAAGAGTTCACCATCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACATGGGGttt                                                                                                                                                                                                           >  1:169533/1‑251 (MQ=255)                                     ggAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTcc                                                                                                                                                                                                         >  2:6960/1‑251 (MQ=255)                                        agagTACGCCATTCATATCGAGAAAGAGTTCAACTTCAACAAACTTCGTCAGCATAACCGTAACCGTCTGCTCTGGAGCAGCAATCTGAATGTTGACGGCATGAAGACAGGAACCGCGGCAGGCGCGGGATCTAATCTGGTTGCATCGGCTACCCGGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAAACGACCGTTTCCGTTTTATTGAGTCTTAGAAATAATTGACATGGGGTTTCAGct                                                                                                                                                                                                      >  2:45438/1‑251 (MQ=255)                                           gTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGcttct                                                                                                                                                                                                   <  2:122107/251‑1 (MQ=255)                                               gCCATTCATAATGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCGACCCCGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGGTTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTg                                                                                                                                                                                                >  2:336275/1‑250 (MQ=255)                                                                      cACCTTCAACAACATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGCTGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTg                                                                                                                                                                        >  2:294382/1‑251 (MQ=255)                                                                        ccTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGAt                                                                                                                                                                      >  1:129869/1‑251 (MQ=255)                                                                              acaaAATTCTTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAAGCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTTATTCCAcc                                                                                                                                                                <  1:408663/251‑1 (MQ=255)                                                                                                                   gAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGGTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACAGTGACGCCAATTAAATCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGGGATAAGAGCg                                                                                                                           >  1:173149/1‑251 (MQ=255)                                                                                                                      cagcaATCTGAATGTTGATGGCATGAAGACAGGCACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGa                                                                                                                          >  1:204317/1‑249 (MQ=255)                                                                                                                                       aTGGCATGAAGACAGGACCCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGCCCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGAGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGt                                                                                                       >  2:535122/1‑251 (MQ=255)                                                                                                                                         ggCATGAAGACAGGAACCAAGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAAGTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACGCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGa                                                                                                     <  2:46002/251‑1 (MQ=255)                                                                                                                                            aTGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCggg                                                                                                             <  2:164755/240‑1 (MQ=255)                                                                                                                                            aTGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCggg                                                                                                             >  1:164755/1‑240 (MQ=255)                                                                                                                                             tGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGCCCGTTTCCGTTTTAATGAGTCTGCGCAATTATTGAACTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCg                                                                                                 >  1:582935/1‑251 (MQ=255)                                                                                                                                                    aGGAACGACGGCCGGCGCGGGATATAATCTGGTTGCATCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGCGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGACTGGTTTGGTGATAAGAGCGAAGTGAATATCGGGGCAGGTGACGCgggatc                                                                                            >  2:425720/1‑246 (MQ=255)                                                                                                                                                    aGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTc                                                                                            >  1:460509/1‑249 (MQ=255)                                                                                                                                                              gCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATAcc                                                                                <  1:446226/251‑1 (MQ=255)                                                                                                                                                                    gcgGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTgg                                                                          <  1:187840/251‑1 (MQ=255)                                                                                                                                                                           ataATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCt                                                                    <  1:352485/250‑1 (MQ=255)                                                                                                                                                                                  ggTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGTCCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGaaa                                                                <  2:125811/247‑1 (MQ=255)                                                                                                                                                                                       cTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGaa                                                        >  2:434569/1‑250 (MQ=255)                                                                                                                                                                                           ggCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGTGCGAAGTGAATCTCGGGGCAGGTGAAGCGTGCTCAGTGACAATACCGCGTGGGCGGCTGAAAAACCTGAAAGCgc                                                   >  2:227742/1‑250 (MQ=255)                                                                                                                                                                                                       cGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACAGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACATTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGAGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTtatcagtta                                        >  2:433172/1‑244 (MQ=255)                                                                                                                                                                                                        gATATGCGTTTAATCGCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTGAATGAGTCTGAGAAAGTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGAGGCCACCTTTGTGACGCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGGTCAGTGACCATTCCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTTATACGTTAAc                                      <  2:544648/251‑1 (MQ=255)                                                                                                                                                                                                               gTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGATTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAACCCTGAAAGCGAGTTATACGTTAACGGAACCg                               >  1:316301/1‑251 (MQ=255)                                                                                                                                                                                                                        ccGTAGTGTTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTTATACGTTAACGTAACCGCAGCTTacc                      <  1:150222/251‑1 (MQ=255)                                                                                                                                                                                                                        ccGTAGTGGTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTTAGAAATTAGTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCTCCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTTATACGTTAACGGACCCGCAGCTTacc                      <  2:440054/251‑1 (MQ=255)                                                                                                                                                                                                                               gtTGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGCTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTTATACGTTAACGGAACCGCAGCTTaccgcaccg                <  2:19079/250‑1 (MQ=255)                                                                                                                                                                                                                                 tGGGGGCGAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTTATACGTTAACGGAACCGCAGCTTACCGCCCCGCTg             >  1:521343/1‑251 (MQ=255)                                                                                                                                                                                                                                        gAAAACCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTTATACGTTAACGGATCCGCAGCTTACCGCACCGCTGAAAAAAg      <  2:303058/251‑1 (MQ=255)                                                                                                                                                                                                                                         aaaaCCGACCGTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTTATACGTTAACGGAACCGCAGCTTACCGCACCGCTGAAAAAAgg     >  1:587770/1‑251 (MQ=255)                                                                                                                                                                                                                                            accgaccgTTTCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTTATACGTTAACGGAACCGCAGCTTACCGCACCGCTGAAAAAAGGTCa  <  2:447928/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  GGCATTGCTGGGTAAAGCATTGATCCACGATGTGCCGGAAGAGTACGCCATTCATAAAGAGAAAGAGTTCACCTTCAACAAAATTCGTCAGCCTAACCGTAACCGTCTGCTGTGGAGCAGCAATCTGAATGTTGATGGCATGAAGACAGGAACCACGGCAGGCGCGGGATATAATCTGGTTGCTTCGGCTACCCAGGGCGATATGCGTTTAATCTCCGTAGTGTTGGGGGCGAAAACCGACCGTATCCGTTTTAATGAGTCTGAGAAATTATTGACCTGGGGTTTCCGCTTCTTTGAAACCGTGACGCCAATTAAACCTGATGCCACCTTTGTGACTCAGCGCGTCTGGTTTGGTGATAAGAGCGAAGTGAATCTCGGGGCAGGTGAAGCGGGCTCAGTGACCATACCGCGTGGGCAGCTGAAAAACCTGAAAGCGAGTTATACGTTAACGGAACCGCAGCTTACCGCACCGCTGAAAAAAGGTCA  >  NC_000913/881287‑881772

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