Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,139,288	C→T	D141N (GAC→AAC)	fsaB ←	fructose‑6‑phosphate aldolase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,139,288	0	C	T	100.0%	200.7 / NA	55	D141N (GAC→AAC)	fsaB	fructose‑6‑phosphate aldolase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/26);  total (29/26)

GGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTCT  >  NC_000913/4139070‑4139533                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       ggCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACg                                                                                                                                                                                                                        >  1:190910/1‑251 (MQ=255) ggCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACg                                                                                                                                                                                                                        >  1:326862/1‑251 (MQ=255)  gCATTCCAGTCGTGTTCGAACTACTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTGGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACgg                                                                                                                                                                                                                       <  2:41492/251‑1 (MQ=255)               gtTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTAtttt                                                                                                                                                                                                          >  2:497504/1‑251 (MQ=255)                 ccgacgttcgctATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGc                                                                                                                                                                                                        <  2:270521/241‑1 (MQ=255)                  cGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGc                                                                                                                                                                                                        >  1:270521/1‑250 (MQ=255)                    agcTTCTCTATAGTGGACTTTACCGCAGGGTGGTTGAGCATTGGTTGCGCTACATCTAAGGGCAGGGTGATGGATTAACAGCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGTATCAACGTATTTTGcccc                                                                                                                                                                                                     <  2:448488/249‑1 (MQ=255)                     ccTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCt                                                                                                                                                                                                    <  1:163176/250‑1 (MQ=255)                         ctctATAGTTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGTTTTCTGGCGCGTGCATTTGTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTTCGCCCTGGGCATCTCAGCGGATAACATACGGAGCAAAGTATTTTGCCACTGcc                                                                                                                                                                                                 >  2:200554/1‑250 (MQ=255)                           ctATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGt                                                                                                                                                                                              >  2:361415/1‑251 (MQ=255)                                       ctACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCtaata                                                                                                                                                                                    <  1:462214/249‑1 (MQ=255)                                         aCCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAAc                                                                                                                                                                                  <  1:218798/249‑1 (MQ=255)                                           cGCAGGTGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGTATTCACATCCTGCCAGTAAACAGGCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGTCGTTTCCGCCCTGGGCATCTACTCGGTTAACATACTGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCtt                                                                                                                                                                              <  2:376536/251‑1 (MQ=255)                                           cGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACACCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCtt                                                                                                                                                                              <  2:476972/251‑1 (MQ=255)                                              aGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACAGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTtgtg                                                                                                                                                                           >  2:328340/1‑251 (MQ=255)                                                      tGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCtc                                                                                                                                                                   >  1:43144/1‑250 (MQ=255)                                                      tGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCt                                                                                                                                                                    <  1:95932/250‑1 (MQ=255)                                                        tgCATTGGTTGCGCTACAACAAAGTGCCGGGTGATGGATTCCCATCCTGCCAGTAAACAGTCCTGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCCGCCCCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTAAACATACGGAGCATCGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCtat                                                                                                                                                                  <  2:316667/249‑1 (MQ=255)                                                        cgCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCtat                                                                                                                                                                  >  2:408079/2‑250 (MQ=255)                                                           aTTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATAc                                                                                                                                                              <  1:5249/251‑1 (MQ=255)                                                              tgttgCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCACCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACACACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATAcagc                                                                                                                                                           >  2:98592/1‑251 (MQ=255)                                                               gttgCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATAcagca                                                                                                                                                          <  1:81496/251‑1 (MQ=255)                                                                     cTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTATCATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCCGTGcc                                                                                                                                                     >  2:480753/1‑250 (MQ=255)                                                                       aCATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGACCCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCa                                                                                                                                                    >  2:153456/1‑249 (MQ=255)                                                                                 gaggggTGAGGGATTCAAATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCCCCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATTCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATa                                                                                                                                        <  2:476718/248‑1 (MQ=255)                                                                                     ggtgtctgATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCCGGCTGCCAGCACCATGCTGTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCAGCTACGCGGTTATCATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCCGTGCCAAGTGTAGTAATACCct                                                                                                                                    <  2:315330/244‑1 (MQ=255)                                                                                     ggTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATACCACTTGTGCGGCGCTATATACAGCAGTGCCAAGGGTAGTAATACCct                                                                                                                                    >  2:304888/1‑251 (MQ=255)                                                                                               tCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTCCCAGCACCATGCTTTCTGGCGCGTGCCTTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGCTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATCTACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTc                                                                                                                            >  1:382341/1‑249 (MQ=255)                                                                                                cacaTCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAg                                                                                                                          <  1:304311/250‑1 (MQ=255)                                                                                                   aTCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTAtt                                                                                                                      >  1:427626/1‑251 (MQ=255)                                                                                                           aGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGAACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCGCCTGAACCGTACGACTGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCACGTGTAGTAATACCCCCTTTTTTCAGTATTTTAATtgc                                                                                                              >  2:314229/1‑251 (MQ=255)                                                                                                               aaCAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGcc                                                                                                          <  1:209978/251‑1 (MQ=255)                                                                                                                 cAGTCCAGCCCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAg                                                                                                        <  2:51718/251‑1 (MQ=255)                                                                                                                   gTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGAc                                                                                                      <  1:304888/251‑1 (MQ=255)                                                                                                                        gcggcTGACGAGACGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGCGCATTTCTAACAGCGTTTGCAGCTCCTTAACCGTACGAAAGCCCATTCCGCCCTGAGCATCTACGCGGTTAATATACGGTGCAAAGGATTTTGCCACTACCAGTGAGGCGAATAACCCTTGTGCGGCGCTATATCCAGCCGTGCAAAATGAAGTCATACACTCTTTTTTCAGTATTTAAATTGCTGCCAGACATTCg                                                                                                 >  2:354710/5‑251 (MQ=255)                                                                                                                                          tAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGAtttt                                                                               <  1:480753/250‑1 (MQ=255)                                                                                                                                            aaGCTAGCTGTCAGCCCAATGCTGTCTGGCGCGTCAAACTTTAACAGCGTTTGCAGCTCCTGAACCGTACGCATGCCGTTTCAGCCCTGGGCATCTACGTGGTTAACATACGGAGAAACGTATTTTGCCCCTGCCAGTACGGCTCCTAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGCAGCAATACACTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGT‑ACCAGGGATTTTcac                                                                            >  2:376929/1‑251 (MQ=255)                                                                                                                                              gctggcttCCAGCACCATTCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGACCCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTACTAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCCCC‑GGGATTTTcacc                                                                           >  2:131579/1‑250 (MQ=255)                                                                                                                                               ctggctgCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATAcagc                                                                                                                                                           <  1:273200/170‑1 (MQ=255)                                                                                                                                               ctggctgCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATAcagc                                                                                                                                                           >  2:273200/1‑170 (MQ=255)                                                                                                                                                     gCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCACAATAcc                                                                   >  1:493683/1‑251 (MQ=255)                                                                                                                                                      ccAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAATATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCACAATAccc                                                                  >  1:489106/1‑251 (MQ=255)                                                                                                                                                              aTGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCACAATACCCGGAATAGc                                                          >  1:500452/1‑251 (MQ=255)                                                                                                                                                                     cTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCACAATACCCGGAATAGCGTcgc                                                     >  2:209478/1‑249 (MQ=255)                                                                                                                                                                     cTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCACAATACCCGGAATAGCGTcgc                                                     <  1:209478/249‑1 (MQ=255)                                                                                                                                                                      tGGCGCGTGCATTTCTAACCGTGTTTGCAGCTCTTGAACCGTCCGACTGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACCGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATATCCTCTTTTTTCAGTATTTTACTTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCACAATACCCGGAATAGCGTCgcgc                                                   >  2:129514/1‑250 (MQ=255)                                                                                                                                                                                aTTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGcgcttcgc                                        <  1:498848/251‑1 (MQ=255)                                                                                                                                                                                     tAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATCCGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCACAATACCCGGAATa                                                            <  2:395965/226‑1 (MQ=255)                                                                                                                                                                                     tAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCACAATACCCGGAATa                                                            >  1:395965/1‑226 (MQ=255)                                                                                                                                                                                     tAACACCGATTGCAGCTCCAGAACCGTACGAATGCCGTTTCCGCACTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATCCAGCAGTGCCAAGTGTAGTACTCCCCTCATTTTTCAGTATTTTACTTGCTGCCAGACCTTCGGGAGTCACC‑GGGATTTTAACCACAATACCCGGAATAGCGTCTCGCAGGCGCTTCGcttctt                                   >  2:157235/1‑251 (MQ=255)                                                                                                                                                                                      aaccgCGTGTGCAGCTCCTGAACCGTTCGATTGCCGTTTCCGCCCTGGGCATCTACGCGGGTAACATACGGATCATCGTAGTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACATCAGTGCCAAGTGTAGTAATTCCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCGCAATACCCGGAATAGCGTCGCGCAGGCGCTTCGcttctt                                   <  2:28845/246‑1 (MQ=255)                                                                                                                                                                                          gCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTccacc                              <  1:497504/251‑1 (MQ=255)                                                                                                                                                                                                  aGCTCCTGACCCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTGACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATccctgc                      >  1:80506/1‑248 (MQ=255)                                                                                                                                                                                                          aaCCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTcgcg             <  1:323360/251‑1 (MQ=255)                                                                                                                                                                                                                     tGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCt  >  1:469757/1‑251 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       GGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACC‑GGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTCT  >  NC_000913/4139070‑4139533

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