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%	276.5 / NA	73	D141N (GAC→AAC)	fsaB	fructose‑6‑phosphate aldolase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/33);  total (40/33)

TAGAGATGAGTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTC  >  NC_000913/4139047‑4139532                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     tAGAGATGAGTCGTGCCAAAGGCGACATGCCAGACGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTCACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCt                                                                                                                                                                                                                                              <  2:46799/250‑1 (MQ=255)                ccAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGtt                                                                                                                                                                                                                              <  1:485262/251‑1 (MQ=255)                   aaTGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAAc                                                                                                                                                                                                                           <  1:102726/251‑1 (MQ=255)                      gCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCATTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACAt                                                                                                                                                                                                                         <  1:329020/250‑1 (MQ=255)                       cGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATAc                                                                                                                                                                                                                       >  2:506646/1‑251 (MQ=255)                              ccAGTCGTGTTCGAACTTCTCTATAGCTTACTCTACCGCAGGGGTGTTGATCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGCTTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCa                                                                                                                                                                                                                 <  1:56740/250‑1 (MQ=255)                                aGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACg                                                                                                                                                                                                              >  1:346135/1‑251 (MQ=255)                                 gTCGTGTTCGAACTTCGCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGt                                                                                                                                                                                                             <  2:341424/251‑1 (MQ=255)                                  tCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGCATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGt                                                                                                                                                                                                             >  2:406003/1‑250 (MQ=255)                                    gtgtTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTAtt                                                                                                                                                                                                          <  1:53503/251‑1 (MQ=255)                                           aaCTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGACCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGcccc                                                                                                                                                                                                   >  1:225424/1‑251 (MQ=255)                                                  ctATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGt                                                                                                                                                                                            >  1:1049/1‑251 (MQ=255)                                                    ataGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGAGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGc                                                                                                                                                                                          >  1:34583/1‑251 (MQ=255)                                                      aGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCgg                                                                                                                                                                                        >  2:313615/1‑251 (MQ=255)                                                         tGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCt                                                                                                                                                                                      >  1:319885/1‑250 (MQ=255)                                                         tGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGAGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCt                                                                                                                                                                                      >  2:371626/1‑250 (MQ=255)                                                          gtctctACCGCTGGGGTGGTGAGCATTTTTTGCGCTACATCTAAGGGCAGGTTTATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCt                                                                                                                                                                                      <  1:173529/247‑1 (MQ=255)                                                          gACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCta                                                                                                                                                                                     >  1:68664/1‑250 (MQ=255)                                                                cccGCAGGGGAGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAAccc                                                                                                                                                                              <  2:98814/250‑1 (MQ=255)                                                                  cGCAGGGGGGATGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCCGCCAGCACCATGCTTTATGGCGCGTGCTTTTCTAACAGGGTTTGCAGCTCCTGAACCGTACGAAGGCCGTTTCCGCCCTGGGCATCTACGAGATTAACAGACGGAGCAACGTAATTTGCCCCTGCCAGTGTGGCTAATAAccatt                                                                                                                                                                            >  2:487757/1‑248 (MQ=255)                                                                         gtgtTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTgcggc                                                                                                                                                                     >  2:422734/1‑251 (MQ=255)                                                                             ttagcattgtttgCGCAACATCTAAGGGCGGGGTTATGAATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTGAAAGCGGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGTGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCt                                                                                                                                                                  <  2:44023/240‑1 (MQ=255)                                                                                 cATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCtatat                                                                                                                                                              >  1:389937/1‑250 (MQ=255)                                                                                      gttgCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATAcagca                                                                                                                                                        <  1:125536/251‑1 (MQ=255)                                                                                                  cTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCACgtgt                                                                                                                                             >  2:427611/1‑250 (MQ=255)                                                                                                    aaGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCtat                                                                                                                                                                <  2:403044/229‑1 (MQ=255)                                                                                                    aaGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCtat                                                                                                                                                                >  1:403044/1‑229 (MQ=255)                                                                                                    aaGGGCAGGGTGAGGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTgtagt                                                                                                                                          <  2:300842/251‑1 (MQ=255)                                                                                                      gggCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTgtagta                                                                                                                                         >  1:401718/1‑250 (MQ=255)                                                                                                       ggCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAAt                                                                                                                                       <  1:56729/251‑1 (MQ=255)                                                                                                       ggCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAAt                                                                                                                                       <  1:96598/251‑1 (MQ=255)                                                                                                          aGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCACGTGCAGTAATAcc                                                                                                                                    >  2:518390/1‑251 (MQ=255)                                                                                                          aGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATAcc                                                                                                                                    >  1:223922/1‑251 (MQ=255)                                                                                                            ggTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCct                                                                                                                                  <  1:472096/251‑1 (MQ=255)                                                                                                             gTGATTGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCctc                                                                                                                                 <  2:270360/251‑1 (MQ=255)                                                                                                             gTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGCATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCctc                                                                                                                                 >  2:414074/1‑251 (MQ=255)                                                                                                              tGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCctct                                                                                                                                >  2:63551/1‑251 (MQ=255)                                                                                                                 tGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGAGCCAAGTGTAGTAATACCCTCtttt                                                                                                                             >  2:68813/1‑251 (MQ=255)                                                                                                                  ggATCCACACCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCttttt                                                                                                                            <  2:346135/251‑1 (MQ=255)                                                                                                                      tCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTc                                                                                                                          >  1:501311/1‑249 (MQ=255)                                                                                                                          aTCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTAtt                                                                                                                    <  1:422734/251‑1 (MQ=255)                                                                                                                                  aGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATtgc                                                                                                            >  2:222275/1‑251 (MQ=255)                                                                                                                                    tAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATtgctg                                                                                                          >  2:218770/1‑251 (MQ=255)                                                                                                                                     aaaCAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATtgctgc                                                                                                         <  1:166528/251‑1 (MQ=255)                                                                                                                                      aaCAGTCCAGCGCCTTGCGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCGGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATGCCCTCTTTTTTCAGTATTTGAATTGCTGcc                                                                                                        <  2:195908/251‑1 (MQ=255)                                                                                                                                              aGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGGTTGAAGCTCCTGAACCGTACGAACGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATCCGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAAGAACCCTTGTGCGGCGCTATATTCAGCAGTGCCAAGTGTAGTAATCCTCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTc                                                                                                <  2:464/251‑1 (MQ=255)                                                                                                                                                 gcCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGGGCATTTCTAACAGCGGTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCTTCTACGCGGTTAACATACGGATCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGGAGTAATACCCTCTTTTTGCAGTATTTTAATTGCTGCCAGACCTTCGGa                                                                                             <  2:79568/251‑1 (MQ=255)                                                                                                                                                     gACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTc                                                                                         <  1:455900/251‑1 (MQ=255)                                                                                                                                                     gACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTAGTTTAATTGCTGCCAGCCCTTCGGAAGt                                                                                          <  1:362904/250‑1 (MQ=255)                                                                                                                                                          ggcgTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCgg                                                                                    <  1:99194/251‑1 (MQ=255)                                                                                                                                                            cgTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCgggt                                                                                  >  2:104657/1‑250 (MQ=255)                                                                                                                                                               tttAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGAttt                                                                               >  2:293852/1‑251 (MQ=255)                                                                                                                                                               tttAAAGCGGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGAtt                                                                                <  2:14411/250‑1 (MQ=255)                                                                                                                                                                 tAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTg                                                                                                                                                     >  1:181583/1‑179 (MQ=255)                                                                                                                                                                 tAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTg                                                                                                                                                     <  2:181583/179‑1 (MQ=255)                                                                                                                                                                 tAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTc                                                                             >  1:277746/1‑251 (MQ=255)                                                                                                                                                                    aGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTcacc                                                                          >  2:36508/1‑251 (MQ=255)                                                                                                                                                                               aGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCgg                                                               >  1:507206/1‑251 (MQ=255)                                                                                                                                                                                 cACCCTGCTTTCTGGCGCATGCATTTCTAACAGCGTTTGCAGCTCATGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCCGTGCCAAGTGTAGCAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGCCACCGGGATTTTACCCACAATACCCGGa                                                              >  2:105958/1‑250 (MQ=255)                                                                                                                                                                                  aCCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCggcat                                                            >  1:522959/1‑248 (MQ=255)                                                                                                                                                                                  aCCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAAt                                                            >  1:116753/1‑251 (MQ=255)                                                                                                                                                                                   ccATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAAt                                                            >  2:283572/1‑250 (MQ=255)                                                                                                                                                                                                      cATTTCTAACAGCGTTTGCAGCTCCTGCACCGTACGAATGCCGTTTCCTCCCTGGGCATCTACGCGGTTAACAGACGGATCAACGTATTGTGACCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCCGTGCCAAGTGTAGTAATACCCTCTCTGTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGcgcttcg                                        >  2:439806/1‑251 (MQ=255)                                                                                                                                                                                                         ttCTAACAGCGTTTGCAGCACCTGAACCGTCCGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGAATTTTGCCCCTGCCAGTGCGGCTAATAACCATTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGcgcttcgctt                                     >  2:7252/1‑251 (MQ=255)                                                                                                                                                                                                           cTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATCCAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATCGCGTCGCGCAGGCGCCTCGcttct                                   >  2:321458/1‑251 (MQ=255)                                                                                                                                                                                                                            tCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGCTTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGcgcg                 >  2:267297/1‑251 (MQ=255)                                                                                                                                                                                                                                gAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCCGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCTCCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTc               <  2:239235/249‑1 (MQ=255)                                                                                                                                                                                                                                gAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTc               <  2:401718/249‑1 (MQ=255)                                                                                                                                                                                                                                 aaCCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTcgcg            >  1:475134/1‑251 (MQ=255)                                                                                                                                                                                                                                    cGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCt         <  2:319885/251‑1 (MQ=255)                                                                                                                                                                                                                                     gTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTc        <  2:116753/251‑1 (MQ=255)                                                                                                                                                                                                                                         gAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATgg    >  2:40066/1‑251 (MQ=255)                                                                                                                                                                                                                                           aTGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTc  <  1:406003/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     TAGAGATGAGTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTC  >  NC_000913/4139047‑4139532

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