Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,396,772	T→C	F236F (TTT→TTC)	amiB →	N‑acetylmuramoyl‑L‑alanine amidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,396,772	0	T	C	100.0%	236.1 / NA	62	F236F (TTT→TTC)	amiB	N‑acetylmuramoyl‑L‑alanine amidase B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/34);  total (28/34)

CGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTG  >  NC_000913/4396526‑4397012                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 cGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAg                                                                                                                                                                                                                                              <  2:322434/251‑1 (MQ=255) cGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGGTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCAGCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCGCCATCGCCAGCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGAGGTTCAAAg                                                                                                                                                                                                                                              <  2:27527/251‑1 (MQ=255)    aCCGGCGCGAAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTAGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCg                                                                                                                                                                                                                                           >  1:292293/1‑251 (MQ=255)             cAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAAccc                                                                                                                                                                                                                                  <  1:259940/251‑1 (MQ=255)               aTCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCgg                                                                                                                                                                                                                                >  1:386648/1‑250 (MQ=255)               aTCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGATCAATGACGATCCGATGGTCAAAAGCGTTTTAACCCGt                                                                                                                                                                                                                                >  2:330603/1‑251 (MQ=255)                         aaCGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGTGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGAtt                                                                                                                                                                                                                      <  2:82410/251‑1 (MQ=255)                         aaCGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGAtt                                                                                                                                                                                                                      <  2:192010/251‑1 (MQ=255)                         aaCGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGAtt                                                                                                                                                                                                                      <  1:7855/251‑1 (MQ=255)                         aaCGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGAtt                                                                                                                                                                                                                      <  1:324021/251‑1 (MQ=255)                          aCGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTa                                                                                                                                                                                                                     <  1:306128/251‑1 (MQ=255)                               aaGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTa                                                                                                                                                                                                                <  2:154459/251‑1 (MQ=255)                               aaGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTa                                                                                                                                                                                                                <  1:49820/251‑1 (MQ=255)                                aGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGGTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTAt                                                                                                                                                                                                               <  2:388693/251‑1 (MQ=255)                                   aaCCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTAtct                                                                                                                                                                                                             <  1:25393/250‑1 (MQ=255)                                       gCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTg                                                                                                                                                                                                        <  2:156917/251‑1 (MQ=255)                                            aCGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCGCGGTGATg                                                                                                                                                                                                     <  2:149388/249‑1 (MQ=255)                                              gggTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCCTCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGgc                                                                                                                                                                                                 >  1:377116/1‑251 (MQ=255)                                              gggTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGAAGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGAGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTAATTTGCTCAATGACCATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATggggg                                                                                                                                                                                                 >  1:40510/1‑250 (MQ=255)                                                  gtTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTCACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGGACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGgcgc                                                                                                                                                                                               >  1:242625/1‑249 (MQ=255)                                                  gtTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGGAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGGTTACTTTATCTCGGTGATGGGgcgc                                                                                                                                                                                               >  2:412227/1‑249 (MQ=255)                                                          cagTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTgg                                                                                                                                                                                     <  2:292293/251‑1 (MQ=255)                                                            gTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCa                                                                                                                                                                                   <  1:138091/251‑1 (MQ=255)                                                                  aggTACCGCGTCCGGCGGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAg                                                                                                                                                                             <  2:177252/250‑1 (MQ=255)                                                                   ggTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGGTTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGc                                                                                                                                                                            >  2:323312/1‑251 (MQ=255)                                                                   ggTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCAGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCCCATAAGc                                                                                                                                                                            >  2:208512/1‑251 (MQ=255)                                                                        cgcgTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACg                                                                                                                                                                                                                            <  2:116702/198‑1 (MQ=255)                                                                        cgcgTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACg                                                                                                                                                                                                                            >  1:116702/1‑198 (MQ=255)                                                                           gTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCCTCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGGGATGTGGCACGTAAGCAAAACGcc                                                                                                                                                                    >  1:347134/1‑251 (MQ=255)                                                                                             cGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTTGCGCTATAGGCCCCGGGGGTACGCGGGGGAAAAATATCACCATCGCCCTCGCGCGGAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGGGACGGGGATTACTTTATCTCGGTGATTGGGCGCAGTGATGTGGCACGTAAACAAAACGCCAATTTCCTCGTGTCGAtt                                                                                                                                                  >  2:59525/1‑251 (MQ=255)                                                                                              ggCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTc                                                                                                                                                 <  1:1623/251‑1 (MQ=255)                                                                                              ggCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCAGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGGAAATTTCGTACTTTGCTCAATGACGATCCGATGTTCCAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGGGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTTTCGATTc                                                                                                                                                 >  2:199341/1‑251 (MQ=255)                                                                                                   acacTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCt                                                                                                                                            <  1:59525/251‑1 (MQ=255)                                                                                                                           cTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGGACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCg                                                                                                                    >  2:105167/1‑251 (MQ=255)                                                                                                                             aTTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCg                                                                                                                    <  1:177156/249‑1 (MQ=255)                                                                                                                               tGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTg                                                                                                                >  1:346823/1‑251 (MQ=255)                                                                                                                                   gCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGGTGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGgcgc                                                                                                            >  2:154646/1‑251 (MQ=255)                                                                                                                                      ggACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCCTCGCCCTCGCGCGTAAATTTCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTCAACCCGTGACGGGGATTACTTTATCTCGGTGGTGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCCCGCTGATGCCGCCCCGAACCGCAGGGCGACTGGCGCTTc                                                                                                         >  2:279403/1‑251 (MQ=255)                                                                                                                                           cggcggTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACATAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTAt                                                                                                    <  1:208512/251‑1 (MQ=255)                                                                                                                                               ggTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGCCTGGCGCTTCCGTATgggg                                                                                                >  2:68283/1‑250 (MQ=255)                                                                                                                                                  cAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGcc                                                                                                                                      <  1:306594/210‑1 (MQ=255)                                                                                                                                                  cAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGcc                                                                                                                                      >  2:306594/1‑210 (MQ=255)                                                                                                                                                  cAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGGTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGct                                                                                             <  1:199341/251‑1 (MQ=255)                                                                                                                                                   aGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTGCTTTATCTCGGTg                                                                                                                                                                                                        <  1:333528/143‑1 (MQ=255)                                                                                                                                                   aGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTg                                                                                                                                                                                                        >  2:333528/1‑143 (MQ=255)                                                                                                                                                     gACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGctctc                                                                                          >  1:248561/1‑251 (MQ=255)                                                                                                                                                      aCCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGctctct                                                                                         <  2:92006/251‑1 (MQ=255)                                                                                                                                                                     gCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTACCCGTCGCGCCAAc                                                                          >  2:92244/1‑251 (MQ=255)                                                                                                                                                                      ccccGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACa                                                                         <  2:386648/251‑1 (MQ=255)                                                                                                                                                                       cccGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAg                                                                        <  2:226300/251‑1 (MQ=255)                                                                                                                                                                        ccGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCTCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGt                                                                       <  1:371248/251‑1 (MQ=255)                                                                                                                                                                          ggtggtACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTg                                                                      >  1:83172/1‑250 (MQ=255)                                                                                                                                                                                      gagaAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAgctg                                                         <  1:323312/251‑1 (MQ=255)                                                                                                                                                                                         aaaaaTGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAgctggct                                                      <  1:412227/251‑1 (MQ=255)                                                                                                                                                                                                        gccatcgcGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACgaga                                        >  2:93371/1‑250 (MQ=255)                                                                                                                                                                                                         ccatcgcGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGaa                                       >  2:346112/1‑250 (MQ=255)                                                                                                                                                                                                                 cgTAAATTGCGTACTTTGCGCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGa                              <  2:318090/251‑1 (MQ=255)                                                                                                                                                                                                                    aaaTTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAgct                           >  1:206311/1‑251 (MQ=255)                                                                                                                                                                                                                         gCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCtgg                      <  1:323458/251‑1 (MQ=255)                                                                                                                                                                                                                         gCGTACTTTGCTCAATGACGATCAGATGTTCAAAGGCGGTTTAACCCGGGACGGGGATTACTGTATCGCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAGAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACGGGCGCTTCCGTATGGGTGCGATCGAACCGTCGCGCGAACAGTGAAATGGCAAGCTGGCAGGAGCAGCACGAGAAACAGTCGGAGCTGCtgg                      >  1:96914/1‑251 (MQ=255)                                                                                                                                                                                                                                     cAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTg          <  1:154646/251‑1 (MQ=255)                                                                                                                                                                                                                                             aTCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTg  >  1:46971/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 CGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCAGGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTG  >  NC_000913/4396526‑4397012

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