Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,243,185	C→T	T156T (ACG→ACA)	folE ←	GTP cyclohydrolase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,243,185	0	C	T	100.0%	194.1 / NA	54	T156T (ACG→ACA)	folE	GTP cyclohydrolase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/29);  total (25/29)

GTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGCGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGG  >  NC_000913/2242951‑2243423                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               gTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGcc                                                                                                                                                                                                                                >  2:8808/1‑250 (MQ=255) gTGACGTTGCGCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGc                                                                                                                                                                                                                                 >  1:490842/1‑250 (MQ=255)   taCGTTGCTCTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGTACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTTAACAATCCACCAAGAGAAGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCATTAATTCTCCGCTTCGATCGAGACAGCCACGTTATTGGTTCCCATCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCAc                                                                                                                                                                                                                               <  1:568870/249‑1 (MQ=255)         gcgcTCCATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCgg                                                                                                                                                                                                                         <  2:452666/250‑1 (MQ=255)         gcgcTCCATTGTTCATGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCggc                                                                                                                                                                                                                        >  2:186637/1‑250 (MQ=255)              ccATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGc                                                                                                                                                                                                                     <  1:205595/249‑1 (MQ=255)              ccATTGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTg                                                                                                                                                                                                                   >  1:476967/1‑251 (MQ=255)                  tgTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGc                                                                                                                                                                                                                <  2:6806/249‑1 (MQ=255)                                   cAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAGCTGCACAATGc                                                                                                                                                                                               <  1:92122/250‑1 (MQ=255)                                           atgaCGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAAttt                                                                                                                                                                                      >  1:391826/1‑251 (MQ=255)                                             gaCGCGCAGCGCGCAGAAGCTCTTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGATGTCGTTGGCTTGGCTCTGGTTTCATCGCGGAGGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGAACGATACAGCCACGTTATTGGTGCCCGGCAGCGGTTGTAGCGCAATAAGGATTTGCTGTGTCAGACGTTGCTTCGCCGGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAAtttt                                                                                                                                                                                     <  1:51061/250‑1 (MQ=255)                                                   caGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCCCCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGAcgg                                                                                                                                                                               >  1:552248/1‑248 (MQ=255)                                                      cgcgcAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCg                                                                                                                                                                           <  2:413259/251‑1 (MQ=255)                                                        cgcAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGAt                                                                                                                                                                         <  1:549782/251‑1 (MQ=255)                                                         gcAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATc                                                                                                                                                                        <  1:108182/251‑1 (MQ=255)                                                         gcAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATc                                                                                                                                                                        >  1:416274/1‑251 (MQ=255)                                                                      gcgcGTATTCTGACTGGATTTGAACAATCCACCAAGAGCGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGCAATGCACCGCGTCGCTCGACACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCg                                                                                                                                                           >  1:399858/1‑251 (MQ=255)                                                                      gcgcGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGACTTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCg                                                                                                                                                           >  1:212512/1‑251 (MQ=255)                                                                        gcgTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACt                                                                                                                                                                                                        <  2:116988/204‑1 (MQ=255)                                                                        gcgTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACt                                                                                                                                                                                                        >  1:116988/1‑204 (MQ=255)                                                                          gTATTCTGACTGGATTTGAAAAATACACACAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGACTTCCCGCAGTACTCCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAACGAACTGCACAATGCGGTTAATTTTTGTCAGACCGATCACCGAATATTTCGGGat                                                                                                                                                       >  2:553714/1‑251 (MQ=255)                                                                                                aTCCACCAAGAGCGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGcc                                                                                                                                  >  1:131746/1‑250 (MQ=255)                                                                                                      caAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCCCCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATg                                                                                                                           >  2:224223/1‑251 (MQ=255)                                                                                                                    gTCGTGGCACGGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGgtg                                                                                                             <  2:312711/251‑1 (MQ=255)                                                                                                                         ggCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTcac                                                                                                        <  2:391826/251‑1 (MQ=255)                                                                                                                          gCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCACAGCAGCGTTTGTCGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTcac                                                                                                        >  1:240435/1‑250 (MQ=255)                                                                                                                                 ttGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGCCAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCACTAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGCACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTg                                                                                                >  2:417630/1‑251 (MQ=255)                                                                                                                                         cgGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAAGGCGTTTAATTTTTGACAGACCGAGCACCGAATCTTTCGGGATATAGGCCACCGTCGCGTGGCCATCGATGGGAACAAAATGGTGTTCACAGGTGCTGGTCAGAGt                                                                                        <  2:399858/251‑1 (MQ=255)                                                                                                                                          ggATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTg                                                                                       <  1:417630/251‑1 (MQ=255)                                                                                                                                            aTGCCCCGCACCTTCCCTCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTTCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTTACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGat                                                                                     <  2:466901/251‑1 (MQ=255)                                                                                                                                                     gcCTTCACGCAGTCATGCACCGCGACGACAGAGACAGAAAAGTTATCGGTGCCCAGCAGCGTCTGTAGCGCAATAAGAATTTGCTGTGACAGACGTACCTGCACCTGCGGACGCTGGGCACAGAACTGCACAATGCGGTTAATTTTTGACCGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTACCAAAATGGAGTTCACAGGTGCTGGCCCGAGGGATATCGCGCAc                                                                            >  2:533224/1‑251 (MQ=255)                                                                                                                                                       cTTCACGCAGTAATGCACCGCGTCGTTCGCGACAGCCACGTTATTGGTGCCAAGCAGCGTTTGTAGCGCCATAAGAACTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTCACCAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACgg                                                                          >  2:49894/1‑251 (MQ=255)                                                                                                                                                         tCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTg                                                                        >  1:21699/1‑251 (MQ=255)                                                                                                                                                           ccGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGTCGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGTATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTTTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGAc                                                                      <  1:81599/250‑1 (MQ=255)                                                                                                                                                                  aaTGCACCGCGTCGATCGCGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTc                                                                >  2:262826/1‑250 (MQ=255)                                                                                                                                                                   aTGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCAt                                                              <  1:189957/251‑1 (MQ=255)                                                                                                                                                                      cACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGCCCATTTCATCGa                                                           <  1:420876/251‑1 (MQ=255)                                                                                                                                                                       aCCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGAc                                                          <  2:2871/251‑1 (MQ=255)                                                                                                                                                                           cgTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTc                                                      <  1:118422/251‑1 (MQ=255)                                                                                                                                                                               gatcgaGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCAttt                                                  >  2:15274/1‑251 (MQ=255)                                                                                                                                                                                atcgaGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCAtttt                                                 <  1:453790/251‑1 (MQ=255)                                                                                                                                                                                        aGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCaa                                         <  1:498867/251‑1 (MQ=255)                                                                                                                                                                                         gCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGAGATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAAt                                        <  1:224223/251‑1 (MQ=255)                                                                                                                                                                                          ccACGTTATTGGTGCCCAGCAGCTTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATg                                       <  2:380887/251‑1 (MQ=255)                                                                                                                                                                                           cACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGa                                                                                                                                                                                  >  2:184752/1‑111 (MQ=255)                                                                                                                                                                                           cACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGa                                                                                                                                                                                  <  1:184752/111‑1 (MQ=255)                                                                                                                                                                                            aCGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTGTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAg                                     <  1:528979/251‑1 (MQ=255)                                                                                                                                                                                             cGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTGCGGGAGATAGGCCACCGTCGCTTTGCCATCGATGGGATCAAAATGGTGTTCACAGGTGCGGGTCAGAGTGAGATCTCGCGCGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAgg                                    <  1:352713/251‑1 (MQ=255)                                                                                                                                                                                                tATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTg                                 >  2:133936/1‑251 (MQ=255)                                                                                                                                                                                                          cagcagCGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAACTGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCggg                        >  2:456047/1‑250 (MQ=255)                                                                                                                                                                                                               gcGTTTGTAGCGCAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATGTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTg                  <  2:490842/251‑1 (MQ=255)                                                                                                                                                                                                                          gcAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATcc         >  1:99150/1‑249 (MQ=255)                                                                                                                                                                                                                          gcAATAAGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAg       <  1:549563/251‑1 (MQ=255)                                                                                                                                                                                                                               aaGAATTTGCTGTGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCgg  >  1:194966/1‑251 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               GTGACGTTGCGCTCCATGGTTCCTGCCTTTTAATCAGTTGTGATGACGCACAGCGCGCAGAAACTCGTGGCGCGTATTCTGACTGGATTTGAACAATCCACCAAGAGAGGTCGTTGTCGTGGCACTGGTTGCATCGCGGATGCCACGCGCCTTCACGCAGTAATGCACCGCGTCGATCGAGACAGCCACGTTATTGGTGCCCAGCAGCGTTTGTAGCGCAATAAGAATTTGCTGCGTCAGACGTTCCTGCACCTGCGGACGCTGGGCAAAGAACTGCACAATGCGGTTAATTTTTGACAGACCGATCACCGAATCTTTCGGGATATAGGCCACCGTCGCTTTGCCATCGATGGTAACAAAATGGTGTTCACAGGTGCTGGTCAGAGTGATATCGCGCACGGTGACCATTTCATCGACCTTCATTTTGTTTTCAATGAGGGTGATTTTCGGGAAATTGGCGTAATCCAGACCGG  >  NC_000913/2242951‑2243423

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