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	96.6%	173.2 / ‑2.5	59	G958R (GGG→AGG)	metH	cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (37/20);  total (38/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

GTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAAC  >  NC_000913/4226454‑4226934                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            gTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAg                                                                                                                                                                                                                                        <  1:4971/251‑1 (MQ=255)      gagtttGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCCCCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATcc                                                                                                                                                                                                                                   <  1:157439/248‑1 (MQ=255)                   gTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGA‑AACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGa                                                                                                                                                                                                                       >  1:421491/1‑249 (MQ=255)                                  aCGAAACCGTACGTATTCAGCACGGTCGCAAGAAACCGCGCACACCACCGGTCACGCTGAAAGCGGCGCGCGATACCGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCCTCGAAACGCTGCGTAATTATATCGACTGGACACCGTTCTTTATGCCCTGGGCGCTGGACAGGACGTATCCGCGCATTTTGGGCGATGAAGTGGTgggg                                                                                                                                                                                                       >  2:470916/1‑250 (MQ=255)                                      aaCCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGACGTATCCGCGCATTCTGGAAGATGACGTGGTGGGCGTTg                                                                                                                                                                                                   >  1:155090/1‑251 (MQ=255)                                          gtacgtaTTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGGAGATGAAGTGGGGGGCGTTGGGgc                                                                                                                                                                                               >  1:127543/1‑251 (MQ=255)                                                   cAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCTGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCt                                                                                                                                                                                      <  1:79187/251‑1 (MQ=255)                                                   cAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACCTCGAATGGCCACCGGTCTTTATGACCTGGTCGCTGGCCAGGAAGATTCCGCGAATTCTGTAAGa                                                                                                                                                                                                                    >  2:268501/1‑221 (MQ=255)                                                   cAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGa                                                                                                                                                                                                                    <  1:268501/221‑1 (MQ=255)                                                     gCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGt                                                                                                                                                                                    <  1:125111/251‑1 (MQ=255)                                                     gCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGt                                                                                                                                                                                    >  1:458596/1‑251 (MQ=255)                                                         gggCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGAAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGGGGCGCCGCGGCTGGTTa                                                                                                                                                                                 >  2:1155/1‑250 (MQ=255)                                                            ggcAAGAAACCGCGCACTCCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCGTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGa                                                                                                                                                                             <  1:56924/250‑1 (MQ=255)                                                            cgcAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCCTACACGCCGCCGGGGGCGCACCGTCTCGGCGGGCAGGAAGTCGAAGCCAGCATCGAAACGATGCGTAATTACCTCAACTGTACACCGTTCTTTATAGCCTGGTCGCTGGCCAGGAAAGATCCGCGCAGTCTGGGAGATGGAGGGGGAGGCGTTGAGGAGCAGCGGCTGtttca                                                                                                                                                                                >  2:210947/1‑246 (MQ=255)                                                            cgaaaGAAAACGCGGGCACCTCCGGTCACGCTGTAAGCGGCGAGCGATATCGATTCCGCTTTTGACTGGCGGGCTTCAACGGCGCCGGTGGCGTAACTTCTTGGCGTGCAGGAAGTAGAAGCACGCATCGAAACGCCGCGTAATTACATCGACTGGACACCGTGCTCTATGACCTGGTCGCTGGCGAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGa                                                                                                                                                                             <  1:89533/248‑1 (MQ=255)                                                                 gAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGccc                                                                                                                                                                        >  1:250902/1‑250 (MQ=255)                                                                   aaCCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAAc                                                                                                                                                                      >  1:360137/1‑251 (MQ=255)                                                                    aCCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGGCTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGCAGCGCAGCGGCTGGTTTAAGACGCCAAcc                                                                                                                                                                     >  2:112912/1‑250 (MQ=255)                                                                    aCCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACg                                                                                                                                                                     >  1:50236/1‑251 (MQ=255)                                                                    aCCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACg                                                                                                                                                                     >  1:110951/1‑251 (MQ=255)                                                                         gcACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCCGGCTTACACGCCGCCGGGGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAAATACATAGACAGGCCACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGAGGATTCTGGAAGATAAAGTGGGGGGCGTTGGGGGGCAGCGGCTGtttgaaagcgcccacgagatg                                                                                                                                                                >  2:122107/1‑233 (MQ=255)                                                                             ccccccGGTCACGCTGGAAGCGGCGCGCGA‑AACGATTTCTCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGa                                                                                                                                                           <  2:421491/247‑1 (MQ=255)                                                                              ccaccGGTCACGCTGGAATCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTGACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGa                                                                                                                                                           <  2:175486/251‑1 (MQ=255)                                                                              ccaccGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTgg                                                                                                                                                            >  1:445635/1‑250 (MQ=255)                                                                                        cGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGAAGCCAACGACATGCTGGATAAATTAAgc                                                                                                                                                 >  1:16124/1‑251 (MQ=255)                                                                                              aaGCGGCGCGCGATAACGAGTTCGCTTTTGACGGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCgat                                                                                                                                           <  2:315262/251‑2 (MQ=255)                                                                                                 cggcgCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTACCCCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCAACTGGACACCGTTCTTTATGACCTCGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGCTGAAGTGGTGGGCGTTGGGGCGCAGCGGCTGTTTAAAGACGACAACGGCATGCTGGATAAATTCAGCGACGAGaa                                                                                                                                         >  1:219819/1‑250 (MQ=255)                                                                                                                aTTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGCCTGGGCACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGGGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGGGCCGCGaacaagctggatacgcgg                                                                                                                         >  2:319130/1‑235 (MQ=255)                                                                                                                  ttCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCgcgtg                                                                                                                        <  1:162/250‑1 (MQ=255)                                                                                                                  tgcgctgtggacgggCAGGCTTGCACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGTCGAGAAAACGCTGAAGCCgcgtg                                                                                                                        <  2:105580/237‑1 (MQ=255)                                                                                                                                                  ggTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTgg                                                                                       >  1:131136/1‑251 (MQ=255)                                                                                                                                                   gTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGGGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGGGGCGTGGTGGGCCTGTTCCCGGCAAACCGGGTggg                                                                                      >  1:100051/1‑251 (MQ=255)                                                                                                                                                    tGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTAGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGCGGCGCAGCGGCTGTTTAAAGACGCAAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGGGTGGTGGGCCTGTTCCCGGAGACCCGTGTgggg                                                                                     >  2:241359/1‑250 (MQ=255)                                                                                                                                                       cgcACCGTCTCGGCGTGCAGGAAGTCGA‑GCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGGATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATg                                                                                 >  1:58767/1‑251 (MQ=255)                                                                                                                                                          aCCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGAc                                                                               <  1:319130/251‑1 (MQ=255)                                                                                                                                                               ctcGGCGTGCAGGAAGTCGAAGCCAGCATCGCAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTg                                                                           >  1:412573/1‑250 (MQ=255)                                                                                                                                                                     gTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCGTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTa                                                                    <  2:50236/251‑1 (MQ=255)                                                                                                                                                                      tGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTAc                                                                   >  1:375512/1‑251 (MQ=255)                                                                                                                                                                              gTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAgg                                                                                                                                                                                                <  1:128221/118‑1 (MQ=255)                                                                                                                                                                              gTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAgg                                                                                                                                                                                                >  2:128221/1‑118 (MQ=255)                                                                                                                                                                              gTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGGGACg                                                            >  1:257340/1‑250 (MQ=255)                                                                                                                                                                              gTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGGGTGGGCGATGACATTGAAATCAACCGTGACg                                                            >  2:74505/1‑250 (MQ=255)                                                                                                                                                                                cGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGCATACGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGaa                                                          >  2:227635/1‑250 (MQ=255)                                                                                                                                                                                 gacgCCAGCATCGAAACGCTGCGTAATTGCATCGACGGGACACCGTTCTTTATGACCTGTTCGCTGGTCAGTAAGTAACCGAGCATTCTGGAAGATGAAGTTGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGAAATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAAGCTATCGTGACGAAAc                                                        <  2:436951/248‑1 (MQ=255)                                                                                                                                                                                 gAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAAc                                                        <  1:406200/251‑1 (MQ=255)                                                                                                                                                                                   aGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGCCGCCAACGACATGCTGGATAAATTAAGCGCCGGGAAAACGCTGAATCCGCGTGACGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGGGCCGAAAcgc                                                      >  2:480080/1‑251 (MQ=255)                                                                                                                                                                                         cATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTAccc                                                >  1:296462/1‑251 (MQ=255)                                                                                                                                                                                                       gTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGGCATGCTGGATAAATTAAGCGCCGAGCAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGCCCTTGAAATCTACCGGGACGAAACGCGTACCCATGTGATCAACGTc                                  >  2:180739/1‑251 (MQ=255)                                                                                                                                                                                                        tAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGCGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTc                                  >  2:163356/1‑250 (MQ=255)                                                                                                                                                                                                             aCATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGccac                            >  1:363982/1‑251 (MQ=255)                                                                                                                                                                                                             aCATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGccac                            >  1:123922/1‑251 (MQ=255)                                                                                                                                                                                                                  gACTGGACACCTTTCTGTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGGGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCt                       <  2:313021/251‑1 (MQ=255)                                                                                                                                                                                                                  gACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGAGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCAACATCt                       >  1:73014/1‑251 (MQ=255)                                                                                                                                                                                                                     tGGACACCGTTCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCg                    >  1:96109/1‑251 (MQ=255)                                                                                                                                                                                                                      ggACACCGTTCTTTATGACCTGGGCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGt                   <  2:445635/251‑1 (MQ=255)                                                                                                                                                                                                                               tCTTTATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGGGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGAcc          >  2:471472/1‑251 (MQ=255)                                                                                                                                                                                                                                  ttATGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGGGGGCGGTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGTGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGa        <  2:378604/250‑1 (MQ=255)                                                                                                                                                                                                                                    aTGACCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGaaa      >  1:56156/1‑250 (MQ=255)                                                                                                                                                                                                                                       aCCTGGTCGCTGGCCAGGAAGTATCCGCGCATTCGGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACAGGCTGGATAAATTAAGCGCCGAGAAAACGCGGAAGCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAAc  <  2:250902/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            GTGATGATTTTGTCGCTCGTACCCGCAAGGAGTACGAAACCGTACGTATTCAGCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAAC  >  NC_000913/4226454‑4226934

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