Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,297,782	G→A	P529S (CCT→TCT)	fdhF ←	formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,297,782	0	G	A	92.3%	131.7 / 0.4	56	P529S (CCT→TCT)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base A (20/28);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

AGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCG  >  NC_000913/4297548‑4298020                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                agaatgtctcactatccggcgcCGGCAGTGTCGCACAGTTCCCGGTCATCGAGCGGCAAGAGTAGTGACCACCTTCACGGACCGTTGACAGTATCATCGGGTAATCGTCGGTGAGTTTGTCGATTGGCGCTTCACAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAGGccc                                                                                                                                                                                                                                <  2:174985/233‑1 (MQ=255)     aGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGtt                                                                                                                                                                                                                            >  1:56017/1‑251 (MQ=255)      ggTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTc                                                                                                                                                                                                                           >  1:437003/1‑251 (MQ=255)       gTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCGTCGAACGGCAAGAGTATTGACCAACTTCCCGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGCAGAAATGCGCCAGACCGTTCGAGGTATAACACTTCTCTTTAAACAGATAAGAAGGCCCCTGATCGGCATCTGAAGTATCGAGGCCAGGCCCC‑TGAATGAAGCCCAGTTCg                                                                                                                                                                                                                          >  2:373729/1‑251 (MQ=255)        ttCATCAGCCAGCGCCGACAGTGCCGCACAGTTACCGGTCATCGAACGGCCAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAACCCAC‑TGAATGAAGCCCAGTTCGc                                                                                                                                                                                                                         >  2:234979/1‑251 (MQ=255)         tcatcaGCCAGCGCCGCCAGTGCCGCACCGTTACCGGTCATCGAACGGCAAGCGAAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGAGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCACGGCAAGACCAC‑TGAATGAAGCCCAGTTCGcc                                                                                                                                                                                                                        >  1:450512/1‑251 (MQ=255)         tcaccagccagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGc                                                                                                                                                                                                                         <  1:11729/246‑1 (MQ=255)          catcaGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGccc                                                                                                                                                                                                                       <  1:172590/251‑1 (MQ=255)                          ccgTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCg                                                                                                                                                                                                       <  1:327930/249‑1 (MQ=255)                                 gCACAGTTACCGGTCATCTAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTa                                                                                                                                                                                                <  1:107321/251‑1 (MQ=255)                                      gTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAAt                                                                                                                                                                                           <  1:342423/251‑1 (MQ=255)                                            ggTCATCGAACGGCAAGAGTAGGGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGc                                                                                                                                                                                     <  1:486112/251‑1 (MQ=255)                                             gTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATAGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAAATTCTCTTTAAACAGATAAGAAGTCCACTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGc                                                                                                                                                                                     >  1:233817/1‑250 (MQ=255)                                                        gCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGc                                                                                                                                                                          >  1:335004/1‑250 (MQ=255)                                                                agtGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATc                                                                                                                                                                  <  1:262408/250‑1 (MQ=255)                                                                      cAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGAtc                                                                                                                                                           <  1:440364/251‑1 (MQ=255)                                                                              cGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATACGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGGTCGCCCATTTTCTCGTAAGCCGCACCGTAGAAATCCGGGCACAGATGACGCACCTCATCCCAGCTCTCCTgggg                                                                                                                                                   >  2:166870/1‑250 (MQ=255)                                                                                        aCAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgttgt                                                                                                                                           >  1:130258/1‑249 (MQ=255)                                                                                         caggaCCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGt                                                                                                                                        <  2:3719/247‑1 (MQ=255)                                                                                          aGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACGTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTg                                                                                                                                       <  2:386689/251‑1 (MQ=255)                                                                                           gTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGTGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCAGCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGc                                                                                                                                      <  2:158283/251‑1 (MQ=255)                                                                                             aCCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCAt                                                                                                                                    <  2:94190/251‑1 (MQ=255)                                                                                             aCCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCAt                                                                                                                                    <  1:106408/251‑1 (MQ=255)                                                                                              ccATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATc                                                                                                                                   >  1:406561/1‑251 (MQ=255)                                                                                               cATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCg                                                                                                                                  >  1:260905/1‑251 (MQ=255)                                                                                                aTCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGAGAAGAACTGCGCCAGACCGTTCGGGGTATCCAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTGTCGCGGAAAACCCAC‑TGAATGAACCCCAGTTCGCCCATTTTCTCGTAAGTAGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCCCCTGGGTGTTGTTGTAGTGCATCgg                                                                                                                                 >  2:506769/1‑251 (MQ=255)                                                                                                  cGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGAt                                                                                                                               <  1:467200/251‑1 (MQ=255)                                                                                                     gtaaccggcggggAGTTTGGCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCTTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTCAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTTCATCGGATAAc                                                                                                                            <  2:53662/239‑1 (MQ=255)                                                                                                     gTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAAc                                                                                                                            <  1:337267/251‑1 (MQ=255)                                                                                                          cgtcgGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCCTTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACCGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGTTAACCCAt                                                                                                                        >  2:1389/1‑250 (MQ=255)                                                                                                                  aGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTTCCCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCCGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACTCAACCCATCCCAGATCTCCTGGGGGTTTTTGGAGTGCCTCGGATAACCCATACGGGTGGc                                                                                                               >  2:133956/1‑251 (MQ=255)                                                                                                                  aGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGTGGTATCAAACTTCTCTGTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTCGTTGTAGTGCATCGGATAACCCATACGGGTGGc                                                                                                               <  1:42998/251‑1 (MQ=255)                                                                                                                                            ggaGGTGAAGAACTGCGCCAGACCGCTCGGGGTATCAAACTTCGCTTTAAACAGATAAGAAGTCCCCTGTTCGGCATCTGACGTACCGCGGCAAGACCAC‑TGAAGGAAGCCCTGTGCGCCCATGTCCTCGTAAGTCGCACCGTAGAAAGCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTTTAGTGCATAGGAGAACCTATACGTGGGGCGATTTCACTGATGATTTGCCAGttcg                                                                                     <  2:297661/249‑4 (MQ=255)                                                                                                                                               ggTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGttt                                                                                  >  1:444951/1‑251 (MQ=255)                                                                                                                                                gTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGtttt                                                                                 <  1:471101/251‑1 (MQ=255)                                                                                                                                                     gaaCTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTc                                                                                >  1:55164/1‑247 (MQ=255)                                                                                                                                                      aaCTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTGTAAACAGATAAGAAGGCCCCTGATCGTCATCTGAATTATCGCGGCAAGACCAC‑TGGATGAAGCCCAGTTCGTCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCATCTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGGTTTCACTGATGATTTGCCCGACCGTCTTCAGATc                                                                           <  2:72273/251‑1 (MQ=255)                                                                                                                                                        ctgagcCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATccc                                                                         >  2:368989/5‑251 (MQ=255)                                                                                                                                                              catacgattCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCACGCAAGTCGAACCGTATAAATCCTGGCACATATGTCGCTACTCATCCCAGATCTCCTGGGTGTTGTTTTAGTTCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATCTCg                                                                   <  2:238811/244‑1 (MQ=255)                                                                                                                                                              cAGACCGTTCGGGGTATAAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGCTCGGCATCTGAAGTATCGCGGCAAGGc                                                                                                                                                                                                                                                >  2:376945/1‑79 (MQ=255)                                                                                                                                                              cAGACCGTTCGGGGTATAAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGc                                                                                                                                                                                                                                                <  1:376945/79‑1 (MQ=255)                                                                                                                                                                   cGTTCGGGGTATCAAACTTCTCTTTAATCAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTAACTGATGATTTGTCAGTCCGTTTTCAGATCCCATTTCGGTTc                                                               >  1:220845/1‑250 (MQ=255)                                                                                                                                                                        ggggTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAATCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGc                                                         >  1:229648/1‑251 (MQ=255)                                                                                                                                                                         gggTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGGGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGcc                                                        >  2:313139/1‑251 (MQ=255)                                                                                                                                                                             aTCAACCTTCACCTTAAACAGATAAGAAGACCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTGCGCCCATTTTCTCGTAAGTCGCACCTTAGAAATCCGGGCACAGATGACGCAACTCAGCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATCCGGGTGGCGATTTCACTGATGAGTTGCCAGTCCGTTTACAGATCCCATTTCGGTTCAACCGCCTTga                                                    <  2:460607/251‑1 (MQ=255)                                                                                                                                                                              tCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTga                                                    >  1:168429/1‑250 (MQ=255)                                                                                                                                                                                aaaCTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTga                                                    >  1:249214/1‑248 (MQ=255)                                                                                                                                                                                             aaCAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGTCAAGACCAC‑GGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTTTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGcc                                    <  1:267788/251‑1 (MQ=255)                                                                                                                                                                                             aaCAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGCCCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTGGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGcc                                    <  1:231638/251‑1 (MQ=255)                                                                                                                                                                                              aCAGATAAGAAGTCCCCTGATTGCCGTCGGAAGTACCGCGGCAATACCAT‑TGAATGAAGCCCAGTGCGCCCATTTTCTCGTAAGGCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCACCCCAGATCTCCTGGGTGTTGTTGTAGGGCTGCGTAGACCCCATACGTGTGGCGATTTCACTGATGATTTGCCAGTCCTTTTCCAGATCCCACTTCGTTTCATCCGCCTTGTAGAAACGCTGGAAGCCa                                   <  2:233817/251‑1 (MQ=255)                                                                                                                                                                                                 gATAAGAAGTCCCCTGAGCGGCATCTGAAGTATCGCGGCACGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTATCGTAAGTCGCACCGTATAAATCCGGGCACAGATGACGCAACTCATCCCAGAGCTCCTGGGGGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGCTTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGTAGAAACGCTGGACGCCCCgg                                >  2:369700/1‑251 (MQ=255)                                                                                                                                                                                                       aagtccctggaccGGCATCGGAAGTACCGCGGCAAG‑CCTCATGAAGGTAGCCCAGTACCCCAACTTTCCCGTAAGTTGCACCTTAGAAATCCGGGCACATATGACGCAACTCATCCCAGTTCGCCTGGTTGTTGTTTTAGGGCATCGGATAACCCATACGGGTGGCGATTTCCCTTATGATTTGCCGGTCCGTTTGCAGATCCCAGGTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACCGTCAGcc                          <  2:56017/239‑1 (MQ=255)                                                                                                                                                                                                              cTGATCGGCATCTGAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACCGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGCTAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGaa                                       >  1:241627/1‑231 (MQ=255)                                                                                                                                                                                                                     gCATCTGAAGTATCGCGGCAAGACCCT‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGCCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCAGCCCAGCTCTCCTGGGTGTGGTTGGAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGACCCCATTTCGGTTCACCCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGcc            <  2:437003/251‑1 (MQ=255)                                                                                                                                                                                                                           gAAGTATCGCGGCAAGACCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCAt       <  1:234979/250‑1 (MQ=255)                                                                                                                                                                                                                                aTCGCCGCAAGACCCC‑TGAATGAAGCCCAGTTCGCCCATTTGCGCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCTGATCTCCTGGGTGTTGTTGTAGGGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGACCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCg  <  2:147448/250‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                AGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCAC‑TGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCG  >  NC_000913/4297548‑4298020

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: ‑