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

CCTTTAGAGATGAGTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGT  >  NC_000913/4139043‑4139531                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    ccTTTAGAGATGAGTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCg                                                                                                                                                                                                                                                 <  2:284756/250‑1 (MQ=255)        agaTGAGTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGc                                                                                                                                                                                                                                         <  1:394057/251‑1 (MQ=255)            gAGTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTATAACAGCGTTTGCAGCTCCTGAACCGTACGAAGGCCGTTTCCGCCCAGGGCATCt                                                                                                                                                                                                                                     >  2:12563/1‑251 (MQ=255)                 gtGCCAAATGCGGCCTTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTAcgcg                                                                                                                                                                                                                                <  2:165610/251‑1 (MQ=255)                   gCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGNATGCCGTTTCCGCCCTGGGCATCTACGCGGt                                                                                                                                                                                                                              >  2:242293/1‑251 (MQ=255)                    ccAAATGCGGCATTCCAGTCGTGTTCGAACTTCGCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGtt                                                                                                                                                                                                                             <  1:356355/251‑1 (MQ=255)                                  ccAGCCGTGTGCGAACTTCGCGATAGCTGACTCTACCGCAGGGGTGGTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGAGGGATTCCCATCCTGCCAGTAAACAGTCCAGCGCCTGACGCTGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTCCGCGGTTAACATACGGAGCa                                                                                                                                                                                                                <  2:421026/250‑1 (MQ=255)                                   cagtctggtTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGGGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAAc                                                                                                                                                                                                              <  1:102756/244‑1 (MQ=255)                                   cAGTCGTGTTCGAACTCCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAATCAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAAc                                                                                                                                                                                                              <  1:291047/251‑1 (MQ=255)                                          gtTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTAtttt                                                                                                                                                                                                       >  2:111892/1‑251 (MQ=255)                                          gtTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTAtttt                                                                                                                                                                                                       >  1:305895/1‑251 (MQ=255)                                           ttctaaCTTCTTTAGAGCGGACTCTCCCGCAGTGGTTTTTAGGATTTGTTGCGCGACATCAAATGTCGAGGTGTTGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCCCCCTTCGTTCTGGCTCGTTCATTTCTAACAGCGTGTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTAtttt                                                                                                                                                                                                       <  2:132440/246‑1 (MQ=255)                                                aCTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCt                                                                                                                                                                                                 <  1:170997/251‑1 (MQ=255)                                                         taGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCg                                                                                                                                                                                        >  1:85251/1‑251 (MQ=255)                                                           gCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGc                                                                                                                                                                                      <  1:46966/251‑1 (MQ=255)                                                            cTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCt                                                                                                                                                                                     >  1:155802/1‑251 (MQ=255)                                                              gACTCGACCACAGGGGTGTTGAGCATTTGTTGCGCTACATCGACGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGATTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGGGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATACCGTTTCCGCCCAGGGCATCTACGCGGTGAACATACGGAGCAACGTATTTTGCCACTGCCCGTGCGGCta                                                                                                                                                                                    >  2:313510/1‑250 (MQ=255)                                                               aCTCTCCCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCtaat                                                                                                                                                                                  <  1:418574/251‑1 (MQ=255)                                                                 tctACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCtaata                                                                                                                                                                                 >  1:417400/1‑250 (MQ=255)                                                                        cAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCAACATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGCACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGCATTTTGCCCCTGCCAGTGCGGCtaataaacattgt                                                                                                                                                                         >  2:423776/1‑244 (MQ=255)                                                                         aGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTtgtg                                                                                                                                                                        >  1:315778/1‑251 (MQ=255)                                                                         aGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTNGCCCCTGCCAGTGCGGCTAATAACCCTtgtg                                                                                                                                                                        >  2:257627/1‑251 (MQ=255)                                                                             gtgtTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTgcggc                                                                                                                                                                    <  1:290648/251‑1 (MQ=255)                                                                                gtGAGCATTTTGTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCTAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCTCCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCGCCTGAACCGTCCGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCAAGTGCGGCTATTAACCCGTGTGCGGCGCt                                                                                                                                                                 <  2:456587/250‑1 (MQ=255)                                                                                  gAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCtat                                                                                                                                                               >  2:448500/1‑251 (MQ=255)                                                                                          gttgCGCTACATCTAAGGGCAGGGTGATGGATTCACACCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATAcagc                                                                                                                                                        <  2:70883/250‑1 (MQ=255)                                                                                            agcgcTACATCTAATGAAAGGGCGATGGATTCACAGCCTGCCAGTAAACAGTCAAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGAGCATTTCTAACAGCGTTTGCAGCGCATGAACCGTACGAATACCGTTTCTGCCCTGGGCATATACGCGGTTAACAAACGGAGCAACGTATTTTGACACTGCCAGTGCGGCTAATCACCCTTGTGAGGCGCTATATACCGCAGt                                                                                                                                                     >  2:44395/2‑251 (MQ=255)                                                                                             gcgcTAAATCTAAGGGAAGGGTGATGGCTTCACATCCTGCCAGTAAACACTCCAGCGCCTGACGCGGCGTTTTAAAGATGGCTGCCAGCACCATGCTTTCTGGCGCGCGCATTTCTAACCGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCTTCTACGCGGTTAACAGACGGAGCAACGTACTTTGCCCCTGCCACCGCGCCTACTCACCCTTGGGCGGCGCTATCTAcatccgtg                                                                                                                                                    >  2:287489/1‑245 (MQ=255)                                                                                                  aCATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCa                                                                                                                                                 <  1:61886/249‑1 (MQ=255)                                                                                                  aCATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCGGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATAcagcacggccacg                                                                                                                                               >  2:286519/1‑243 (MQ=255)                                                                                                  aCATCTAAGGGCAGGGTGAGGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAg                                                                                                                                               <  1:75573/251‑1 (MQ=255)                                                                                                         aGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTTGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTgtagta                                                                                                                                        <  1:424241/251‑1 (MQ=255)                                                                                                          gggCAGGGTGATGGATTCACATCCTGCCAGTAAACCGTCCAGCGCCTTACGCGGCGTTTTAAAGCTGGCTGCCAGCCCCAGGCTTTCCGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGATCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTaa                                                                                                                                       <  2:421405/251‑1 (MQ=255)                                                                                                              aGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTACCATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATAcc                                                                                                                                   >  2:71006/1‑251 (MQ=255)                                                                                                                   gATGGATTCACAGCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCtt                                                                                                                              <  2:88693/251‑1 (MQ=255)                                                                                                                    aTGGTTTCACAGCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGt                                                                                                                                                     <  2:53321/227‑1 (MQ=255)                                                                                                                    aTGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGt                                                                                                                                                     >  1:53321/1‑227 (MQ=255)                                                                                                                    aTGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCttt                                                                                                                             <  1:481912/251‑1 (MQ=255)                                                                                                                      gtattccatcccttccAGTAAACAGTCCAGCGCCTGACGCGGCGGTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATGCCCTCttttt                                                                                                                           <  2:158717/237‑1 (MQ=255)                                                                                                                      ggATTCACACCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCtttt                                                                                                                            <  1:5812/250‑1 (MQ=255)                                                                                                                          cccaaTCCTGCCAGGAAACAGGCAAGCGCCTGACGCGGCGTATTAAAGTGGGCTGCAAGCCCCATGCTTTCTGGCGCGTGCATTTCTAACAGCGGTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGTCCTGGGCATCTACGCGGTTAACATACGGTGCAACGTATTTGGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTCTTTTCAg                                                                                                                       <  2:134905/247‑1 (MQ=255)                                                                                                                             caTCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAGCCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTAt                                                                                                                    <  1:257627/251‑1 (MQ=255)                                                                                                                              aTCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATTCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTAtt                                                                                                                   <  1:214369/251‑1 (MQ=255)                                                                                                                               tCCTGCCAGTAAACAGTCCAGTCCCTGACTCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCAGCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTAt                                                                                                                    <  1:438514/249‑1 (MQ=255)                                                                                                                                ccTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTAtttt                                                                                                                 <  1:266419/251‑1 (MQ=255)                                                                                                                                   gCCAGTAAACAGTCCAGCGCCCGACGCGGCGTTTTAAAGCCGGCTGCCAGCACCATGCTTGCTGGCGCGTGCATTTCTAACAGCGTTTGCAGTGCCTTAACCGTACGAATGCCGGTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGGGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTa                                                                                                                <  2:325572/249‑1 (MQ=255)                                                                                                                                             aGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCttttt                                                                                                                           <  2:80311/228‑1 (MQ=255)                                                                                                                                             aGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCttttt                                                                                                                           >  1:80311/1‑228 (MQ=255)                                                                                                                                              gTCCAGCGCCTGACGCGGCGTTTTAACGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCCGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGAc                                                                                                   >  2:433861/1‑251 (MQ=255)                                                                                                                                                    cgcCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAAt                                                                                                                                      >  2:23580/1‑210 (MQ=255)                                                                                                                                                    cgcCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAAt                                                                                                                                      <  1:23580/210‑1 (MQ=255)                                                                                                                                                         gACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGCAGTc                                                                                        >  1:253902/1‑251 (MQ=255)                                                                                                                                                             cggcgTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCTGAAGTCACCg                                                                                    >  1:169938/1‑251 (MQ=255)                                                                                                                                                                      aaaGCTGTCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGCTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTACCATACGTAGCAACGTATTTCGCCCCTGCCAGTGCAGCAACTAACCCTTGAGCGGCGCTATATCCAGCAGGGCCCAGTGTAGTCATTCCCTCTTTCTTCACAATTTTAATTACTTCCAGCCCTTCGTACGTCCCCCGGATTTTc                                                                            >  2:80492/1‑250 (MQ=255)                                                                                                                                                                                gCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATAcc                                                                 >  1:291039/1‑251 (MQ=255)                                                                                                                                                                                 ccAGCTCCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATAccc                                                                <  1:494904/251‑1 (MQ=255)                                                                                                                                                                                   aGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCgg                                                              >  1:251282/1‑251 (MQ=255)                                                                                                                                                                                              ttCTGGAGCGTGCATTTCTAACAGCGGTTGCAGCCCCTGACCCGTTCGAAGGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCATCGTATTTTGCCCCTGCCAGTGCGTCTAATAACCCTGGAGCGGCGCTATATACATCAGTGCTAAGTGTAGTAATACCCTCGGTTTTCAGTATTATAATTGCTGCCAGACCGTCGGAAGTCACCGGGATTTTcaccacgatgcccttgagagcgtct                                                    <  2:148225/250‑20 (MQ=255)                                                                                                                                                                                                  ggCGCGTGCTGTTCTAACGTCGTTTGCAGCTCCTTAACCGTACGAATGGCGTTTCCGCCCTTGGCGTCTACGCGGTTAACATACGGAGCATCGTATTTTGCCCCTGCCAGTGCGGCTAATACCCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCGTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCgcgc                                                 <  2:424105/249‑1 (MQ=255)                                                                                                                                                                                                  ggCGCGTGCATATCTACCAGCGGTTGCAGCTCCTGAACCGTACGAATGCCGTTTCAGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGGGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTGCGGATGTCACCGGGATTNTAACCACAATACCCGGAATAGCGTCgcgc                                                 >  2:256559/1‑249 (MQ=255)                                                                                                                                                                                                           aTTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCCGTGCGGCTCATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGTAAGTCACCGGGATTTTCACCACTATACCCGGAATAGCGTCTCGCAGGcgcttcgc                                      >  2:323981/1‑251 (MQ=255)                                                                                                                                                                                                             ttCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGcgcttcgct                                     <  2:483106/250‑1 (MQ=255)                                                                                                                                                                                                                  aCAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTcc                               <  1:230694/251‑1 (MQ=255)                                                                                                                                                                                                                       gTTTGCAGCTCCTGAACCGTACGAATGCCGTNTCCGCCCTGGGCCTCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCAt                          <  2:253902/251‑1 (MQ=255)                                                                                                                                                                                                                                    gAACAGTACGAAAGCCGGTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATGTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGGATTTTAAGTGCTGCCAGACCGTCGGAAGTCACCGGGATGTGCACCCCAATACCCGGAATAGCGTCGCGCAGGCGCGTCGCTTCTTCCACCATCCC‑TGCGCGTcgc            <  2:505208/251‑1 (MQ=255)                                                                                                                                                                                                                                        cGTACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCt        <  1:422192/251‑1 (MQ=255)                                                                                                                                                                                                                                           aCGAATGCCGGTTCCTCCCTGGGCATCTACGCGTTGAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCAt     <  2:504131/251‑1 (MQ=255)                                                                                                                                                                                                                                             gAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATgg   >  1:120191/1‑251 (MQ=255)                                                                                                                                                                                                                                              aaTGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGt  >  1:180380/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    CCTTTAGAGATGAGTAGTGCCAAATGCGGCATTCCAGTCGTGTTCGAACTTCTCTATAGCTGACTCTACCGCAGGGGTGTTGAGCATTTGTTGCGCTACATCTAAGGGCAGGGTGATGGATTCACATCCTGCCAGTAAACAGTCCAGCGCCTGACGCGGCGTTTTAAAGCTGGCTGCCAGCACCATGCTTTCTGGCGCGTGCATTTCTAACAGCGTTTGCAGCTCCTGAACCGTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGT  >  NC_000913/4139043‑4139531

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