Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,773,173	Δ1 bp	coding (893/1914 nt)	mtlA →	mannitol‑specific PTS enzyme IICBA component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,773,173	0	G	.	100.0%	202.0 / NA	46	coding (893/1914 nt)	mtlA	mannitol‑specific PTS enzyme IICBA component
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (27/19);  total (27/19)

CAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTAAAATCAT  >  NC_000913/3772942‑3773422                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          caggtccaggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTCCGGGCGATGCGGCAATCATCCACTTCCTGGGGGGTATCCACGCAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCCCGCTGACTATCCTGGGCGGAGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTCCTgg                                                                                                                                                                                                                                       >  1:155145/1‑251 (MQ=255)                    gctgctGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTTTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGc                                                                                                                                                                                                                    <  1:116306/251‑1 (MQ=255)                            cGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTcc                                                                                                                                                                                                            >  1:174651/1‑250 (MQ=255)                                       tttGGTCGTGGAAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATcgcg                                                                                                                                                                                                 >  2:53564/1‑251 (MQ=255)                                          ggTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGTGGTATCCACGAAATCTACTTCCCGGATTTGCTGATGAATCCGCGTCTGAGCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGGTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGgt                                                                                                                                                                                              <  1:60752/251‑1 (MQ=255)                                                 gTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTAGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCATCGGCGGTATGACTGGAGTGTTCACGCCGCCTACCCTGGGCGGTGGTCTGGTTTCTCCGGCCTCTCCGG‑TTTTATCCTTGCTGTACTTGCGATGACACCACAAGGAGCTTACTTCGCTAACATCGCGGGtgtgtgtg                                                                                                                                                                                       >  1:48937/1‑251 (MQ=255)                                                    gcgcTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTgcgg                                                                                                                                                                                    <  1:279073/251‑1 (MQ=255)                                                         aaaCAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGAATGAGCTGATGCATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACAGGCCTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGCTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGGGCTTACTTCGCTAACATCGCGGGAGCGAGTGCGGCGa                                                                                                                                                                                 >  1:132656/1‑249 (MQ=255)                                                           aCAGTCTGCGGGCGGTGCGGCAATCAGCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCGCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGc                                                                                                                                                                             <  1:325622/251‑1 (MQ=255)                                                                   cgggcggTGCGGCCATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGCTCCTCGCCGTCATCCTCGCCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGGTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGGGCTTACTTCGCTAACATCGCTGGTGTGTGTGCGGCGATGGCGTTCTCCt                                                                                                                                                                     >  2:189618/1‑251 (MQ=255)                                                                    gggcggTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCATACTTCGCTAACATCGCGGGTGGGTGTGCGGCGCTGGCTGGCTCCtt                                                                                                                                                                    >  2:235096/1‑251 (MQ=255)                                                                        ggTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTAGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATAGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGtt                                                                                                                                                                >  1:31203/1‑251 (MQ=255)                                                                         gTGCGGCACTCATCCACTTCCTGGGGGCTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGCGTTTCCGCTGACTATCTTGGGCGGGGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTCTCCTTGCTGTACTGGCGATGACACCAACAGGTGCTTCCTTCGCTAACATCTCGGGCGAGTGTggggcgcgggttgtcctcttcgttt                                                                                                                                                               >  2:177858/1‑232 (MQ=255)                                                                                 atcatcCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCt                                                                                                                                                       >  1:225676/1‑251 (MQ=255)                                                                                              tGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTAtt                                                                                                                                                    >  2:92738/1‑241 (MQ=255)                                                                                              tGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTAtt                                                                                                                                                    <  1:92738/241‑1 (MQ=255)                                                                                                    gTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGGTGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTCTTTTGCTGAAAAcccgcc                                                                                                                                    >  2:130486/1‑247 (MQ=255)                                                                                                      aTCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTgcg                                                                                                                                                                                     <  1:307018/200‑1 (MQ=255)                                                                                                      aTCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTgcg                                                                                                                                                                                     >  2:307018/1‑200 (MQ=255)                                                                                                          aCGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGAGAGCTCACGCTGACTATACTGGACGGTGGACTGGTTTCTCCGGCAACTCCGG‑TTTTATCCTTGCTGTACTGGCGATGAAACCAAAAGATGCTTACTTCGCTAACATCGAGGGTGTGTGTGCGGCGATGGCTGACTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTg                                                                                                                               >  1:197421/1‑250 (MQ=255)                                                                                                             aaaacTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGTTATGACTGGCGTGTTCGCGCTTTCTACCCTGGGCGGTGGTCCGGTTTCTCCGGCATCTCCTG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGAGGCTATGGCTGTCTCCATCGTTGTCACTGCTATTTTGCTGAATACCAGCAAAGTGAaag                                                                                                                           <  2:115595/247‑1 (MQ=255)                                                                                                                   aCTTCCCGTCTCCGCTGCTGAATCCGCGTCTGATCCTCGCAGTCCTCCTCGGCGGTACGCCTGGCGTGTTCCCGCTGACTCTCCTGGGCGGTGGTCTGGTTTCTCCGGCCTCTCCGG‑TTTTATCCTTGCTGTCCTGGCGATGACACCCAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTCTTTTTCTGACAACCCGCAACGTGAaagaagc                                                                                                                       >  2:317755/1‑248 (MQ=255)                                                                                                                           tATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAACGTGAAAGAAGAAGATGATAtt                                                                                                             >  1:46436/1‑251 (MQ=255)                                                                                                                               tgCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCCCGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGCTGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGa                                                                                                           >  1:52869/1‑249 (MQ=255)                                                                                                                                 cTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCAGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCCAAGGGAAAGAAGAAGATTATATTGAagc                                                                                                        >  2:103043/1‑250 (MQ=255)                                                                                                                                      gaATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCaa                                                                                                   <  1:162112/250‑1 (MQ=255)                                                                                                                                        aTCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACtc                                                                                                <  1:147181/251‑1 (MQ=255)                                                                                                                                        aTCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGGCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACtc                                                                                                >  1:283434/1‑251 (MQ=255)                                                                                                                                             cgTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTcgtcgt                                                                                           <  2:174651/251‑1 (MQ=255)                                                                                                                                                   aTCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGc                                                                                       >  1:83284/1‑249 (MQ=255)                                                                                                                                                        cGCAGTCATCCGCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGAc                                                                                  <  2:103227/249‑1 (MQ=255)                                                                                                                                                               aTCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCCAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCt                                                                         >  2:270315/1‑251 (MQ=255)                                                                                                                                                                     ggcggTATGACTGGAGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAACAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGGGTCt                                                                   >  2:141827/1‑251 (MQ=255)                                                                                                                                                                             gACTGGCGTGTTCACGCTGTCTATCCTGTGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGgcg                                                            <  2:100031/250‑1 (MQ=255)                                                                                                                                                                                tGGCGCGTTCACGCTGACTGTAATGGAAGGAGGCCTGGTTTCTACAGAATCTCCGG‑TATTATCCTTGCTGTACTGGCGCTGTCACCCAACGGTGCTTACTTCGCTAACATAGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTAGTTGTCTCTGCTATTTTGCTGAATAGCAGCAAAGTGAAAGAAAAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACTTGAAAGCTGCGTCTAAAGgcgaatt                                                        >  2:296084/1‑247 (MQ=255)                                                                                                                                                                                   cGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTcc                                                     <  2:323113/251‑1 (MQ=255)                                                                                                                                                                                    gagtTCACGCTGACTATCCTGGGCTGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGGGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAGGTGAAAGAAGAAGATGATATTGAATAAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTAGCAAGGCGCATCTCCg                                                    >  2:279290/3‑251 (MQ=255)                                                                                                                                                                                      gtTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGc                                                   <  1:270315/250‑1 (MQ=255)                                                                                                                                                                                        tCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGCGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGt                                                >  1:121837/1‑251 (MQ=255)                                                                                                                                                                                          aCGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTct                                              <  1:109368/251‑1 (MQ=255)                                                                                                                                                                                                         gtcGGTGTTCTGGTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACt                               <  1:288812/249‑1 (MQ=255)                                                                                                                                                                                                               ggtCTGGTTTCTCCGGCATCTCCGG‑TGTTATCCCTTCTGTACTGGCGATGACACCAAAAGGTGCTTACGTCGCTATCATCGCGGGGGTGTGGGCGGCGATGGCGGTCTCCTTAGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGAATATATTGAAGCATCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCACTGTCTGCTGGCGATGTTACTATCGAc                         <  2:288977/251‑1 (MQ=255)                                                                                                                                                                                                                    ggTTTCTCCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGa                     <  1:35906/250‑1 (MQ=255)                                                                                                                                                                                                                      tttttcCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCCGGAAATGAAAGCTGAgatgccggagaaa                                                          >  1:215196/5‑199 (MQ=255)                                                                                                                                                                                                                      tttttcCGGCATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACGCGTCGTATGCAGGACATGAAAGCttagatgccggagaaa                                                          <  2:215196/207‑16 (MQ=255)                                                                                                                                                                                                                               cATCTCCGG‑TTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCacgtacgt          >  1:155982/1‑250 (MQ=255)                                                                                                                                                                                                                                   ctccggTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTGACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAACCCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGTCATGAGCCACGTAAGTaaa       >  2:184401/5‑250 (MQ=255)                                                                                                                                                                                                                                     ccggTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGTCCTGAGCCACGTACGTaaa       <  2:72174/246‑1 (MQ=255)                                                                                                                                                                                                                                       ggTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGCAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTCAAatcat  >  1:266845/1‑251 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          CAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTAAAATCAT  >  NC_000913/3772942‑3773422

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