Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	167,798	C→A	Y105* (TAC→TAA)	fhuA →	ferrichrome outer membrane transporter/phage receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	167,798	0	C	A	100.0%	165.6 / NA	45	Y105* (TAC→TAA)	fhuA	ferrichrome outer membrane transporter/phage receptor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (20/25);  total (20/25)

TTAGCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGC‑GACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCA  >  NC_000913/167554‑168039                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  ttAGCGACAGGTCGGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCCCCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGC‑GACAGTCTGCTACCGGCTCTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGcc                                                                                                                                                                                                                                              <  1:313687/250‑1 (MQ=255)               gTTTATGCACAGGCAGCTGTTAAACCGAAAGAAGCCACTACCCCCTTTCCCGCTCCACCCGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGC‑GACAGTCGGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTg                                                                                                                                                                                                                                <  2:306317/250‑1 (MQ=255)                tttATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGC‑GACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGt                                                                                                                                                                                                                              >  2:379903/1‑251 (MQ=255)                       acaGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGC‑GACAGTCTGCTACCGGCACTAAAACCGATACGCCGATGCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTg                                                                                                                                                                                                                       <  1:384056/251‑1 (MQ=255)                                 gTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGC‑GACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTCCGCGTGGCGCATCCac                                                                                                                                                                                                             >  2:335077/1‑250 (MQ=255)                                    gAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCCTGGGGGCCTGCTGCAACTATTGCGGCGC‑GACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAacac                                                                                                                                                                                                          >  1:203935/1‑251 (MQ=255)                                    gAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGC‑GACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGCGTCTCTGTTGGTACGCGTGGCGCATCCAaca                                                                                                                                                                                                           <  1:176687/250‑1 (MQ=255)                                    gAACCGAAAGAAGACACCTTCCCCTGTACCGCCGCACCTGCGCCGGAAGTAGGCGCATGGGGGCCTGCTGCAACTATTGCGG‑GCAGACAGTCTGCGACCGCCCCTAAAATCTAGACGCCGATTCAAAAAGTGCCACCGTCTTTTTCTGTTGTGTCCGCCGCAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTTGTACGCGTGGCGCATCCAacat                                                                                                                                                                                                          <  2:93826/251‑2 (MQ=255)                                          aaaGAAGACGTTATCCCCGTTACCGCTGGACCTGCGGCGAAAGAAAGCGCTTGGGGGCCTGCTGCAACTATTGCGGGGC‑GACAGTCTGCTACCGGCCCTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGa                                                                                                                                                                                                    <  2:31894/251‑1 (MQ=255)                                              aagaCACTATCACCGTTCCCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGC‑GACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGaccac                                                                                                                                                                                                <  1:272862/251‑1 (MQ=255)                                                    caaTCAACGTTACCGCTGCACCTGCGCCGCAAGAACGAGCATGGGGCACTGCTACCACTATAGCGGAGC‑GACAGTCTGCTACCGGCACTAACACCGACACGCCGATTCAAAAAGAGCAACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTACAGGACGCGCTTAGCTAAACGCCGGGAGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGacagcctaat                                                                                                                                                                                           >  2:266148/3‑242 (MQ=255)                                                      aTCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGC‑GACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCa                                                                                                                                                                                         <  1:245792/250‑1 (MQ=255)                                                      aTCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGC‑GACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCAt                                                                                                                                                                                        >  2:366694/1‑251 (MQ=255)                                                              tACCGCTGCACCTGCGCCGCCAGAAAGCGCATGGGGGCCTGCTGCACCTATTGCGGCGC‑GCCAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGACGAGATGGCGCTGCATCAGCCGACGTCGGTAAAAGAAGCGCTTCGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTCTGCCCACCTGCTCATTCGCGGCt                                                                                                                                                                                >  2:43188/1‑251 (MQ=255)                                                                       aCCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGC‑GACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGc                                                                                                                                                                         >  2:79613/1‑249 (MQ=255)                                                                                  aaGAAAGCGCAGGGGGGCCTGCTGCAACTATTGCGGCGC‑GACAGTCTGCTACCGGCACTAAAACCGATACGGCGATTCAAAAAGTGCCACAGTCTAGTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCGTCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGc                                                                                                                                                            <  1:163817/251‑1 (MQ=255)                                                                                            aTGGGGGCCTGCTGCAACTATTGCGGCGC‑GACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACt                                                                                                                                                  >  2:247509/1‑251 (MQ=255)                                                                                                 ggCCTGCTGCAACTATTGCGGCGC‑GACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTg                                                                                                                                             >  2:202556/1‑251 (MQ=255)                                                                                                                       gc‑gACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGAGGAGCATCCAACACCTATGACCACCTGATCATTCGCGGCGTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGcc                                                                                                                       >  2:131715/1‑250 (MQ=255)                                                                                                                        c‑gACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCaa                                                                                                                      <  1:366694/251‑1 (MQ=255)                                                                                                                        c‑gACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCa                                                                                                                       >  1:360245/1‑250 (MQ=255)                                                                                                                            cAGTCTGCTCCCGTCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACtt                                                                                                                   <  1:233611/251‑1 (MQ=255)                                                                                                                            cAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACtt                                                                                                                   >  2:144780/1‑251 (MQ=255)                                                                                                                               tCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGGTGCAGGGCAACTTCt                                                                                                                 >  2:98944/1‑250 (MQ=255)                                                                                                                                  gCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCCCAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCtat                                                                                                               <  1:79613/249‑1 (MQ=255)                                                                                                                                        ggCACTAAACCCGATCCGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTATCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCATAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGc                                                                                                       <  2:135108/251‑1 (MQ=255)                                                                                                                                          cACTAAAACCGATCCGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGCATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCgg                                                                                                     >  2:125943/1‑251 (MQ=255)                                                                                                                                           aCTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGt                                                                                                    <  1:386002/251‑1 (MQ=255)                                                                                                                                              aaaaCCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGGCCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAAGGGCCTGAAGGTGCAGGGCAACTTCTATAACGATGCGGTc                                                                                                   <  2:394585/249‑1 (MQ=255)                                                                                                                                                  ccGATCCGCCGATTCAAAAAGGGCCACAGTCTATTTCTGTTGTGCCCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGAc                                                                                             <  2:87210/251‑1 (MQ=255)                                                                                                                                                     aTACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCg                                                                                          <  1:202556/251‑1 (MQ=255)                                                                                                                                                                     aGTGCCTCAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAAcgcg                                                                          >  1:114086/1‑251 (MQ=255)                                                                                                                                                                        gCCACAGTCTTTTGCGTTTGTGCCCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCt                                                                        <  2:261014/250‑1 (MQ=255)                                                                                                                                                                        gCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTg                                                                       <  1:102479/251‑1 (MQ=255)                                                                                                                                                                         ccACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTg                                                                       >  1:303355/1‑250 (MQ=255)                                                                                                                                                                             atcgtatttcggtgtagACCGCCTAAGAGATGGCGCTGCATCAGCCGAAGTCGGGAAAAGAAGCGCCTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACGATCTGAATGGCCTGTAGTTGCAGGGCAACTTCTATAACGAGGCTGTCATTGACCCGTATATGCTTGAACGCGCTGAAAtt                                                                  <  2:203935/235‑1 (MQ=255)                                                                                                                                                                              gTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTa                                                                 <  1:379903/251‑1 (MQ=255)                                                                                                                                                                                gtATTTCGGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGGCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGGCATTGACCCGTATATGCTGGAACGCGCTGAAATTAt                                                                <  2:164628/249‑1 (MQ=255)                                                                                                                                                                                   tttCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCCACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGc                                                              >  1:430854/1‑249 (MQ=255)                                                                                                                                                                                                                  gCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAg                               <  1:65725/249‑1 (MQ=255)                                                                                                                                                                                                                          cGAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTggc                     <  2:144098/251‑1 (MQ=255)                                                                                                                                                                                                                           gAAGTCGGTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGTCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCTGAAGGCCAAAGCCAGAATCACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTAGCACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTTCGGTAAAAGCAGTCCTggcg                    >  2:41718/1‑251 (MQ=255)                                                                                                                                                                                                                                 ggTAAAAGAAGCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGt              >  1:333469/1‑251 (MQ=255)                                                                                                                                                                                                                                          agCGCTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGAGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATgg     >  1:215989/1‑251 (MQ=255)                                                                                                                                                                                                                                             gcTTAGCTAAACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCAGTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCa  <  2:147939/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  TTAGCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGC‑GACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGCCACAGTCTATTTCTGTTGTGACCGCCGAAGAGATGGCGCTGCATCAGCCGAAGTCGGTAAAAGAAGCGCTTAGCTACACGCCGGGTGTCTCTGTTGGTACGCGTGGCGCATCCAACACCTATGACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCA  >  NC_000913/167554‑168039

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