Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,226,699	G→A	G958R (GGG→AGG)	metH →	cobalamin‑dependent methionine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,226,699	0	G	A	97.0%	210.7 / ‑2.6	67	G958R (GGG→AGG)	metH	cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (29/35);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

ATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTA  >  NC_000913/4226457‑4226945                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       atgatTTTGTCGCCCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAATCCGCGCACACCACCGGGCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTGGACGGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTAt                                                                                                                                                                                                                                                <  2:334173/251‑1 (MQ=255)           cGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACAGCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCAtt                                                                                                                                                                                                                                       <  1:288691/250‑1 (MQ=255)              ccgtccccggtaGGATTACGAAACCGTACGTATTCTGCCCGGGCTCAAGATACCGCGCACACCGCCGGTCTCGCTGGAAAGGGCTCGCGAGAACGTTTTCGCTGTTGACTGGCAGGCTTACACGACGCCGGTGGCGCACCGTCGCGGTGTGCAGGAAGTCGAAGCCAGCATCGAAACGCCGCGTAATTACATCGTCTGGTCACCGTTCGTTATGGCCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTgg                                                                                                                                                                                                                                   <  2:43715/240‑1 (MQ=255)                 tACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAAATACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGa                                                                                                                                                                                                                                  >  1:286007/1‑249 (MQ=255)                     cGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGa                                                                                                                                                                                                                            <  1:5362/251‑1 (MQ=255)                     cGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGa                                                                                                                                                                                                                            <  1:147394/251‑1 (MQ=255)                       cAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCGCCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGa                                                                                                                                                                                                                            >  1:232219/1‑249 (MQ=255)                        aaGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAgg                                                                                                                                                                                                                         >  1:123774/1‑250 (MQ=255)                            aGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTAGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAgtggtg                                                                                                                                                                                                                     <  2:388087/251‑1 (MQ=255)                                 gAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCCGCATCGAAACGCTGCGTAATTACATCGCCTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGACGTATCCGCGCATTCTGGCAGATGAAGGGGTGGGCGt                                                                                                                                                                                                                >  2:189506/1‑251 (MQ=255)                                  aaaCCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCggt                                                                                                                                                                                                               >  1:363408/1‑249 (MQ=255)                                  aaaCCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGtt                                                                                                                                                                                                               >  1:443203/1‑251 (MQ=255)                                   acccGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCACATCCAGCATCGAAACGCTGCGTAATTACATCGCCTGGCCACCGTTCTTTATGACCTGGTCGCTGGCCCGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGGGGGCGTTg                                                                                                                                                                                                              >  2:467449/3‑251 (MQ=255)                                                cAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCTATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCgtct                                                                                                                                                                                                 <  2:102581/251‑4 (MQ=255)                                                   cACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGtt                                                                                                                                                                                              <  1:205396/251‑1 (MQ=255)                                                        gcgcAAGAAACCGCGCACACCACCGGTCATGCTGGAAGCGGCGCGCGATAACGATTGCGCTTGTGACTGGCAGGCTTACACGCCGCCGGTGGCGCGCCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACGGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAg                                                                                                                                                                                         <  1:37588/251‑1 (MQ=255)                                                          gcAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGAc                                                                                                                                                                                       >  1:34333/1‑251 (MQ=255)                                                               aaaCCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGGGGGCGTTGAGGCGCAGCGGCTGATTCAAGACGCCa                                                                                                                                                                                   >  2:113271/1‑250 (MQ=255)                                                                        cccccaccGGTCCCGCTGGAAGCGGCGCGCGATAACTATTTCGCTGTTGAATGGCAGGCTTTCACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGTATCGAAACGCTGCGTAATTACATCGACTGGACCCCGTTCTTTATGACCTGTTAGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCTCAGCGGCTGTTTAAAGACGCCAACGACATGCt                                                                                                                                                                         <  2:316665/248‑1 (MQ=255)                                                                        acacCCCCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCt                                                                                                                                                                         <  1:254369/251‑1 (MQ=255)                                                                          accaccGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCAACGTCTCGGCGTGCAGGAAGTCGACACCAGCCTCGAAACGCTGAGTAATTAAATCGCCTGGACACCGTTATTTATGACCCGGTCGCTGGCAAGGAAGTATACGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGAGGCTGTTTAAAGAAGCCAACGACCTGccgg                                                                                                                                                                       >  2:370271/1‑248 (MQ=255)                                                                               cGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATaa                                                                                                                                                                   >  1:9265/1‑250 (MQ=255)                                                                                ggTCACGCTGGAGGCGGCGCGCGACAACGATTTCTCTTTTGACTGGCAGGCTTACACGCCGCCGGGGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCTGCTGTTTATAGACGCCAACGACATGCTGGATAAAt                                                                                                                                                                 <  1:363523/251‑1 (MQ=255)                                                                                 gTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGGGGCGCAGCGGCTGGTTAAAGACGCCAACGACATGGTGGAAAAAtt                                                                                                                                                                >  2:141498/1‑251 (MQ=255)                                                                                      gCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTTGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTAGAAGCCAGCATCGAAACGCAGCGTAATTACATCGGCTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCAGCGCATTCTGGAAGATGAAGTGGTGGGCTTTGAGGCGCAGCGGCTGGTTAAAAGCGCCACCGGAATGATGGATCAATTAAgcg                                                                                                                                                           >  2:310232/1‑251 (MQ=255)                                                                                              cggcgCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTTAAGACGCCAACGACATGGTGGATAAATTAAACGCCGAGaa                                                                                                                                                    >  2:409462/1‑250 (MQ=255)                                                                                                   cgcgTTAACGATTTCGCTTTTGGCTGGCAGGTTTACACGCCGACGGTGGCGCACCGTCGCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCGTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCt                                                                                                                                              <  2:39131/251‑1 (MQ=255)                                                                                                         aaCGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTTGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCa                                                                                                                                                                                                       <  1:425739/188‑1 (MQ=255)                                                                                                         aaCGAATTCGCTTTTAACTGGCAGGCTTACACGCCGCCCGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGCAGCCAGCATCGAAACGCTGCGTAATTAAATCGCATGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCa                                                                                                                                                                                                       >  2:425739/1‑188 (MQ=255)                                                                                                             aTTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCcgcgg                                                                                                                                    >  1:127245/1‑250 (MQ=255)                                                                                                                       tGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGAc                                                                                                                                                                                       <  1:129166/190‑1 (MQ=255)                                                                                                                       tGACTGGCAGGCTTACACGACGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCCACGACATGCTGGATAAATTAAGCGCCGAGAAAACCCTGAATCCCcgtggggtggggg                                                                                                                          >  2:41134/1‑248 (MQ=255)                                                                                                                              cAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGtt                                                                                                                   <  1:199129/251‑1 (MQ=255)                                                                                                                               aGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTc                                                                                                                  <  1:101586/251‑1 (MQ=255)                                                                                                                                ggctgacccgccgcCGGGGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTcc                                                                                                                 <  2:443203/243‑1 (MQ=255)                                                                                                                                 gCTTACACTCCGACGGTTGCGCACCGGCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTccc                                                                                                                <  1:206436/251‑1 (MQ=255)                                                                                                                                           cgccgGTGGCGCACCGTCGCGGCGTGCAGGAAGTCGACGCCTGCACCGAAACGGGGCGTAAGTACATCGACGGGACACCGTTGTTTATGACCTGGTCGCTGGACAAGAAGTATCCGCGCATTATGGAAGATGAAGTGGTGGGCTGTGAGGCGCAGCGGCTTTTTAAAGACGCCAACGACATGCTGGATAAATTAGGCGCCGAGAAAACGCTGAACCCGAGTTGTGTGGTGGGCGTGTTCCCGCCAAACCgt                                                                                                      <  2:21219/251‑1 (MQ=255)                                                                                                                                             ccgGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCgtgt                                                                                                    <  1:162076/251‑1 (MQ=255)                                                                                                                                                gTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTTGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAGGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTggg                                                                                                 <  2:407320/251‑1 (MQ=255)                                                                                                                                                atGGCGCACCGTCTCGGCGTGAAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGAGATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGGGGTGGGCCTGTTCCCGGCAAACCGTGgggg                                                                                                 >  2:108177/2‑251 (MQ=255)                                                                                                                                                      cACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATg                                                                                            >  1:157620/1‑250 (MQ=255)                                                                                                                                                              cGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGa                                                                                     <  1:386205/249‑1 (MQ=255)                                                                                                                                                                          aGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGCCATTGAAATCTACCGTGACg                                                                       >  1:234310/1‑251 (MQ=255)                                                                                                                                                                           gTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCTACTGGCCACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCTTCCTGGAAGGTGAAGTGGTGGGCGTTGAGGCAAAGCGACTGTTTAAAGACGCCCCCGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAGTCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCCTGTGGGCGGTGGCATTGTAATCTACCGGGACgg                                                                      >  1:467751/1‑250 (MQ=255)                                                                                                                                                                             cGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGACGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGGGGTGGGCATGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGaa                                                                     >  1:456330/1‑250 (MQ=255)                                                                                                                                                                             cGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGaa                                                                     <  2:135496/250‑1 (MQ=255)                                                                                                                                                                               aaGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGGGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAAcg                                                                  >  1:713/1‑251 (MQ=255)                                                                                                                                                                                  ccAGCATCGAAACGCTGAGTAATTACATCGACTGGCCACCGTTCTTTCTGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAGATACGCCAACGACATGCTGGGTCAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGAGGGCCTGTTCCCGTCAACCCGTGTGGGCGATGGCATTGACATCTACCGTGCCGAAAcacg                                                                >  2:55365/1‑247 (MQ=255)                                                                                                                                                                                         tgAAACGCTGCGTAATTACATCTACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTTCCGTGACGAAACGCGTACCCAt                                                         <  2:134548/249‑1 (MQ=255)                                                                                                                                                                                          gAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCAtgt                                                       >  1:130321/1‑251 (MQ=255)                                                                                                                                                                                          gAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCAtgt                                                       >  2:116083/1‑251 (MQ=255)                                                                                                                                                                                           aaaCGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGCCGAAACGCGTACCAAtgtg                                                      >  1:25071/1‑251 (MQ=255)                                                                                                                                                                                               gCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCa                                                  <  1:130404/251‑1 (MQ=255)                                                                                                                                                                                                   cGTAATTACATCGACTTGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAAGTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGt                                              <  1:303061/251‑1 (MQ=255)                                                                                                                                                                                                   cGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTATGGAAGATGACGTGGTGGGCGTTGAGGCGTAGCGGCTGTTTAAAGACGCCAACGACAGGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGCTACCGGCAAACCGTGTGGGCAGTGACATTGAAATCTACCGTGCCGAAACGCGAACCCATATGATCAACGt                                              >  2:214189/1‑251 (MQ=255)                                                                                                                                                                                                       aTTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGGAATCTACCggtacga                                                                      >  2:410678/1‑217 (MQ=255)                                                                                                                                                                                                       aTTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGa                                                                      <  1:410678/223‑1 (MQ=255)                                                                                                                                                                                                         tACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGcc                                         <  1:23680/250‑1 (MQ=255)                                                                                                                                                                                                            aTCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGTCCTGTTCCCGGCAAACCGTGTGGGCGATGACCTTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGccacc                                      >  1:287492/1‑250 (MQ=255)                                                                                                                                                                                                                        ccGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAAc                          >  1:386932/1‑250 (MQ=255)                                                                                                                                                                                                                         cGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAAGCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAg                        <  1:55365/251‑1 (MQ=255)                                                                                                                                                                                                                               ttATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCg                    <  2:109528/249‑1 (MQ=255)                                                                                                                                                                                                                                 aTGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGaaaa                <  1:95194/251‑1 (MQ=255)                                                                                                                                                                                                                                           cggtGGCCAGGAAGTATCCGCGCATTGTGGAGGATGAAGTGGTGGGCGTGGAGGCGCAGCGGCTGTGTAAAGACGCCAACGACATGCTGGATGAAGTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTCCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACTCGTACCCATGTGATCACCGTCAGCCACCATCTGCGTCAACAGTCCGAAAAAACAggatgc      <  2:246185/248‑5 (MQ=255)                                                                                                                                                                                                                                            gCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCg     <  1:41134/251‑1 (MQ=255)                                                                                                                                                                                                                                               ggCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCt   >  1:263741/1‑250 (MQ=255)                                                                                                                                                                                                                                               ggCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTGGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTa  <  2:34333/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       ATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTA  >  NC_000913/4226457‑4226945

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