Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,312,219	A→G	L258P (CTA→CCA)	alsR ←	DNA‑binding transcriptional repressor AlsR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,312,219	0	A	G	100.0%	170.9 / NA	47	L258P (CTA→CCA)	alsR	DNA‑binding transcriptional repressor AlsR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/15);  total (32/15)

GCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCAT  >  NC_000913/4311975‑4312457                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               atgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTa                                                                                                                                                                                                                                          >  1:141842‑M2/83‑250 (MQ=255)    tggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAatat                                                                                                                                                                                                                                      >  1:96955‑M2/80‑251 (MQ=255)                    tattacgggcactatttattcaaaactctgacgaaaaacaggatgtcgtcagtcttgacgtgaGGTAAGGGGACA‑GTAACTACGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCCTATTAATATTAGCTTGATCAATGTTGAGCTGGGCGACAGCGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGAAtt                                                                                                                                                                                                                       >  2:417339‑M2/64‑250 (MQ=255)                        acgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGAc                                                                                                                                                                                                                  >  1:309764‑M2/60‑251 (MQ=255)                                              tctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGc                                                                                                                                                                                            >  1:304933‑M2/38‑251 (MQ=255)                                                         acaggctgtcgtctgttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGGATTACCTTTTTCTTATTACAGCGCGCCTGGTGAGAAGAAATAAACAATTGCGCCGGTTTTTTCCATATTATTATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTATTGCCGAGGGATTAAGACCTAATAACGGCGTTTCCGGGGctggtgggaaa                                                                                                                                                                                  >  2:129‑M2/27‑245 (MQ=255)                                                             gctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAaataat                                                                                                                                                                             <  1:499150‑M2/229‑1 (MQ=255)                                                                      agttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGc                                                                                                                                                                    <  1:125738‑M2/238‑1 (MQ=255)                                                                         tttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGcc                                                                                                                                                                   >  2:442872‑M2/11‑249 (MQ=255)                                                                           tgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGtt                                                                                                                                                               >  1:20637‑M2/9‑251 (MQ=255)                                                                             acgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCGTGGTGATAATAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCTTCCAGCAAGGTTGGTTGTAATATTCTTGGCGTTGCATTGCGTCCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTc                                                                                                                                                             <  1:32064‑M2/245‑1 (MQ=255)                                                                               gtgaggTAAGGGGACAG‑TAACTCGGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGGGCGAAGAAATCAACAATTGCGACGGTTTTTTCCGTATTAATATTAGATTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCAAGCAAGGTTGTTTGTAATATTCCTGCCGAGGCAATACGACCTAATAACTGCGTTTGCGGGGCTGGTGAGCAAATAAAATGATCGGGCCGTTTCGc                                                                                                                                                           >  2:571413/5‑251 (MQ=255)                                                                                            gacagtacctcagcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGAGAGAAGAAATCAACACTTGCGCAGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGATGGGCGCCAGAGACAAAAAAAGCGTACAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGTCTGGGGAGCAAATAATATGACAGGCCAGCTTCTCTATCGGTGagtg                                                                                                                                               >  2:305131/13‑248 (MQ=255)                                                                                             acagTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATgg                                                                                                                                              >  1:470519/4‑251 (MQ=255)                                                                                             acagTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATgg                                                                                                                                              >  1:355594/4‑251 (MQ=255)                                                                                               agTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGt                                                                                                                                             >  1:493651/2‑250 (MQ=255)                                                                                               agTAACGCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGACTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTTGGCGACAGAGGCAAAAAAAGCGGCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGGTTCCGGGGCTGGTGAGCAAATAATATGATCGTCCAGTTTCGCTATCGGTGAAGGGTa                                                                                                                                            <  2:166123/250‑1 (MQ=255)                                                                                                               cGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACACTTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGACTGGTAGCTATGGGTTATACa                                                                                                                             >  2:355488/1‑250 (MQ=255)                                                                                                                  aaTTACCTTTTTCTTATTTCAGCGCGACTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGGTGGTTGGAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAaat                                                                                                                          >  1:127829/1‑250 (MQ=255)                                                                                                                             tCTTATTTCAGCGCGCCTGGTGAGAAGAAAACAACAATTGCGCCGGTTTTTTCCCTCTTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGGTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCCAATAATATGATCGGGCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGcccc                                                                                                              >  1:547695/1‑251 (MQ=255)                                                                                                                                   ttCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAA‑GCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCtt                                                                                                        <  1:234586/250‑1 (MQ=255)                                                                                                                                    tCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCttt                                                                                                       >  1:65257/1‑251 (MQ=255)                                                                                                                                     cAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGCGGCATTACGACATAATAACGGCGATTCCGGGGCTGGCGAGCAAATAATATGATCGGCCAGTTTCGCCATCGGTGAATGGTAGATATGGGTTATACAAATAACATTTGCCCCGTTCtttt                                                                                                      >  2:1565/1‑251 (MQ=255)                                                                                                                                        cgcgcCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTCTACAAATAATCTTTGCCCCGTTCTTTTTTg                                                                                                   >  1:542304/1‑251 (MQ=255)                                                                                                                                                   agaagaAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACg                                                                                        >  1:491661/1‑251 (MQ=255)                                                                                                                                                     aagaaATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGc                                                                                      >  1:333443/1‑251 (MQ=255)                                                                                                                                                         aaTCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCt                                                                                  >  1:231216/1‑251 (MQ=255)                                                                                                                                                            caacaaTTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCttt                                                                                <  1:129390/250‑1 (MQ=255)                                                                                                                                                                 aTTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGGTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTGTACGGCCGCTTTTACATc                                                                          <  2:493651/251‑1 (MQ=255)                                                                                                                                                                    gcgcCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGGGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACTGCCGCTTTTACATccct                                                                       >  2:197090/1‑248 (MQ=255)                                                                                                                                                                                           tttaGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGAATTACGACCTAATAACGGCGTTTCAGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGACTGGGAGCTATGGGTTATACAAATAAACTTTGCCCCGTTCTTTTTTGCCAATTCTACGGCCGCTTTTACATCAATGGTCCGCCCGGAATGGGTCACt                                                 >  2:127920/2‑250 (MQ=255)                                                                                                                                                                                           attaGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGAATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACt                                                 >  2:68817/1‑250 (MQ=255)                                                                                                                                                                                             taGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCTGAATGGGTCACTAc                                               <  2:217725/250‑1 (MQ=255)                                                                                                                                                                                                cTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGc                                           >  1:580051/1‑251 (MQ=255)                                                                                                                                                                                                     taaaTGTTGAGCTGTGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATACTTGCCGAGGCATTACGACCTAATAACGTCGTTTCCTGGGCTGTTGAGCAAATAATATGATCGGCCAGTTTCGCTATCTGTGAATGGTAGCTATGGGTTATACAAAGAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTGTTACATCACTGGCTCGCCCGGAATAGGTCACGACCAGCacaa                                       <  2:518795/248‑1 (MQ=255)                                                                                                                                                                                                      cAATGTTGAGCTGTGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCTCCCGTAATGGGTCACTACCAGCacaac                                      <  1:129/250‑1 (MQ=255)                                                                                                                                                                                                             gAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCtcct                               <  1:569049/250‑1 (MQ=255)                                                                                                                                                                                                                  gggCGACAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCtccttcct                           <  1:342226/249‑1 (MQ=255)                                                                                                                                                                                                                       aCAGAGACAAAAAAAGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTaaca                     >  1:341370/1‑250 (MQ=255)                                                                                                                                                                                                                           agaCAAAAAAAGCGTCCATCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGAGTTTCCGGGGCTGGTGAGCAAATAATATGATCGACCAGTTTCGCTATCGGAGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGa                 >  2:495174/1‑250 (MQ=255)                                                                                                                                                                                                                                   aaaGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGAACTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGAcat        <  2:388653/251‑1 (MQ=255)                                                                                                                                                                                                                                    aggCGTCCAGCAAGGTTGGTTGTAATAGTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCGGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCCACAGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGAcatc       <  2:90389/249‑1 (MQ=255)                                                                                                                                                                                                                                    aaGCGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGACACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGAcatc       >  1:289934/1‑251 (MQ=255)                                                                                                                                                                                                                                       cGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCCGCACAACAtct                                  >  2:15440/1‑221 (MQ=255)                                                                                                                                                                                                                                       cGTCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACAtct                                  <  1:15440/221‑1 (MQ=255)                                                                                                                                                                                                                                         tCCAGCAAGGTTGGTTGTCATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACACATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGAcatcatca   >  2:130560/1‑250 (MQ=255)                                                                                                                                                                                                                                         tCCAGCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGAcatcatcat  >  2:73763/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               GCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCAT  >  NC_000913/4311975‑4312457

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