Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,226,790	G→A	S988N (AGC→AAC)	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,790	0	G	A	100.0%	197.2 / NA	58	S988N (AGC→AAC)	metH	cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/26);  total (32/26)

GCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAAT  >  NC_000913/4226506‑4227084                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                       gCACGGGCGCAAGAAACCGCGCACACCAACGGTCACGCTGGAAGCGGCGCGTGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCAGGTGGCGCACCGTCTAGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAGTTACATCGACTGGACACCGCTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCt                                                                                                                                                                                                                                                                                         >  2:146396/1‑300 (MQ=255)      ggCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTNCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATc                                                                                                                                                                                                                                                                                    <  1:513828/300‑1 (MQ=255)      ggCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTAGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATc                                                                                                                                                                                                                                                                                    >  1:641492/1‑300 (MQ=255)              aaaCCGCGCACACCACCGGTCACGCTGGAAGNGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGNAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCgcgtggc                                                                                                                                                                                                                                                                            >  2:551731/1‑300 (MQ=255)               aaCCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCgcgtggcg                                                                                                                                                                                                                                                                           <  1:728076/300‑1 (MQ=255)                   gcgcACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCgtggt                                                                                                                                                                                                                                                                       >  1:577912/1‑300 (MQ=255)                    cgcACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCgtggtg                                                                                                                                                                                                                                                                      <  1:200872/300‑1 (MQ=255)                     gcACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGtggtgg                                                                                                                                                                                                                                                                     <  1:413934/300‑1 (MQ=255)                        caccacCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGcc                                                                                                                                                                                                                                                                  >  1:412353/1‑300 (MQ=255)                        caccacCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTAGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGcc                                                                                                                                                                                                                                                                  >  1:762098/1‑300 (MQ=255)                         accaccGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCt                                                                                                                                                                                                                                                                 <  1:106546/300‑1 (MQ=255)                          ccaccGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTg                                                                                                                                                                                                                                                                <  1:268733/300‑1 (MQ=255)                                  cACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGTAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGc                                                                                                                                                                                                                                                        <  1:170903/300‑1 (MQ=255)                                  cACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAATTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGc                                                                                                                                                                                                                                                        <  2:247355/300‑1 (MQ=255)                                    cGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGTTTACACGCCGCCGGTGGCGCACCGTCTAGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAGTTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATNCGCGCATTCTGAAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCAAACGAAATGCTGGATAGATTAAAGGCCAAGAAGACGCTGAATCAGTGTGACGTGGTGGGcacgttcacggtaa                                                                                                                                                                                                                                                      >  2:169325/1‑287 (MQ=255)                                                                           ggCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGaaa                                                                                                                                                                                                               <  1:614499/300‑1 (MQ=255)                                                                               ggCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATNCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAANGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCAGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAAcgcg                                                                                                                                                                                                           >  2:601474/1‑300 (MQ=255)                                                                               ggCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGCGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAAcgcg                                                                                                                                                                                                           >  2:495693/1‑300 (MQ=255)                                                                               ggCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAAcgcg                                                                                                                                                                                                           >  1:56502/1‑300 (MQ=255)                                                                               ggCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAAcgag                                                                                                                                                                                                           >  1:420365/1‑298 (MQ=255)                                                                               ggCTTACACGCCGCCGGTGGCGCACCGTCTAGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAAcgcg                                                                                                                                                                                                           >  2:309143/1‑300 (MQ=255)                                                                               ggCTTACACGCCGCCGGTGGCGCACCGTCTAGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAAcgcg                                                                                                                                                                                                           >  1:237191/1‑300 (MQ=255)                                                                                    acacGCCGCCGGTGGCGCACCGCCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCGGGGAAGTATCCGCACATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCCTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTAccc                                                                                                                                                                                                      <  2:21537/300‑1 (MQ=255)                                                                                          cgccgGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGa                                                                                                                                                                                                <  1:143811/300‑1 (MQ=255)                                                                                               gTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAANCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAAc                                                                                                                                                                                           >  1:546933/1‑300 (MQ=255)                                                                                                 ggCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGNGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGt                                                                                                                                                                                         <  1:563856/300‑1 (MQ=255)                                                                                                      aCCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGcc                                                                                                                                                                                    <  1:617988/300‑1 (MQ=255)                                                                                                         gTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGGCACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCACATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGATGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGGCATTGAAATCTACCGTGACGAAACGCGTACCCATGTGGTCAACGTCAGtcacc                                                                                                                                                                                 >  2:211478/1‑300 (MQ=255)                                                                                                         gTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCCTGTGATCAACGTCAGccacc                                                                                                                                                                                 >  1:747353/1‑300 (MQ=255)                                                                                                         gTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCAAACGACATGCTGGATAAATTACACGCCGCGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCTATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGccacc                                                                                                                                                                                 >  1:52209/1‑300 (MQ=255)                                                                                                                 gTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGt                                                                                                                                                                         >  1:18559/1‑300 (MQ=255)                                                                                                                     aGGAAGTCGAAGCCAGCATCTAAACGCTGCATAACTACATCGACTGGACGCCGTTCTTTATGACCTGGTCGCTGGCTGGGAAGTATCCGAGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGTAGCGGCTGTT‑TAAAGACGCCAACGACATGTTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGGCAGCCACCATCGGCGTCAAc                                                                                                                                                                     <  2:2620/300‑1 (MQ=255)                                                                                                                     aGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAAc                                                                                                                                                                     <  1:737280/300‑1 (MQ=255)                                                                                                                                     cATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAgg                                                                                                                                                     >  2:436992/1‑300 (MQ=255)                                                                                                                                                  cGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGGTGAGGCGCAGCGGC‑GTTCTAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAGATCTACCGTGACGAAACGCGTACCTATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAactac                                                                                                                                        <  2:256726/300‑1 (MQ=255)                                                                                                                                                    tAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTg                                                                                                                                      >  1:699112/1‑300 (MQ=255)                                                                                                                                                    tAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTg                                                                                                                                      >  2:109843/1‑300 (MQ=255)                                                                                                                                                    tAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAACCACCATCTGCGTCAACAGACCGAAACAACAGGCTTCGCTAACCACTg                                                                                                                                      >  2:110141/1‑300 (MQ=255)                                                                                                                                                    tAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCAGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTAACGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTg                                                                                                                                      >  2:286062/1‑300 (MQ=255)                                                                                                                                                       ttACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTTGCCGGTAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGtct                                                                                                                                   <  2:205521/300‑1 (MQ=255)                                                                                                                                                       ttACATCGACTAGACACCGTTCTCTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTTGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACTAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGtct                                                                                                                                   <  2:88355/300‑1 (MQ=255)                                                                                                                                                               aCTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACt                                                                                                                           >  1:27354/1‑300 (MQ=255)                                                                                                                                                                    acacCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGtt                                                                                                                      >  1:649220/1‑300 (MQ=255)                                                                                                                                                                            cTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGaa                                                                                                              >  2:136100/1‑300 (MQ=255)                                                                                                                                                                                               gCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGAt                                                                                           >  2:250001/1‑300 (MQ=255)                                                                                                                                                                                                       gTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCgg                                                                                   >  1:85118/1‑300 (MQ=255)                                                                                                                                                                                                              cgcATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCAttt                                                                            <  1:704741/300‑1 (MQ=255)                                                                                                                                                                                                                  ttCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTAGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCg                                                                        >  2:211968/1‑300 (MQ=255)                                                                                                                                                                                                                        aaGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTg                                                                  >  1:679199/1‑300 (MQ=255)                                                                                                                                                                                                                            tGAAGTGGTGGGCGTTGAGGCGCAGCTGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCTTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCgg                                                              <  1:246446/300‑1 (MQ=255)                                                                                                                                                                                                                            tGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCNCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCgg                                                              <  1:531975/300‑1 (MQ=255)                                                                                                                                                                                                                             gAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTAGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCggg                                                             >  2:740284/1‑300 (MQ=255)                                                                                                                                                                                                                                    tggGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAag                                                      <  1:570743/300‑1 (MQ=255)                                                                                                                                                                                                                                      ggCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAagag                                                    <  1:204057/300‑1 (MQ=255)                                                                                                                                                                                                                                                       gCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGc                                   <  1:717438/300‑1 (MQ=255)                                                                                                                                                                                                                                                                 aGACGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAgc                          <  1:361927/300‑1 (MQ=255)                                                                                                                                                                                                                                                                   aCGCCAACGACATGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAgcgc                        <  1:299960/300‑1 (MQ=255)                                                                                                                                                                                                                                                                              aTGCTGGATAAATTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAGCGCAGCACgatgat             >  1:324010/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                         aTTAAACGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAAt  >  2:324143/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                       GCACGGGCGCAAGAAACCGCGCACACCACCGGTCACGCTGGAAGCGGCGCGCGATAACGATTTCGCTTTTGACTGGCAGGCTTACACGCCGCCGGTGGCGCACCGTCTCGGCGTGCAGGAAGTCGAAGCCAGCATCGAAACGCTGCGTAATTACATCGACTGGACACCGTTCTTTATGACCTGGTCGCTGGCCGGGAAGTATCCGCGCATTCTGGAAGATGAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTT‑TAAAGACGCCAACGACATGCTGGATAAATTAAGCGCCGAGAAAACGCTGAATCCGCGTGGCGTGGTGGGCCTGTTCCCGGCAAACCGTGTGGGCGATGACATTGAAATCTACCGTGACGAAACGCGTACCCATGTGATCAACGTCAGCCACCATCTGCGTCAACAGACCGAAAAAACAGGCTTCGCTAACTACTGTCTCGCTGACTTCGTTGCGCCGAAGCTTTCTGGTAAAGCAGATTACATCGGCGCATTTGCCGTGACTGGCGGGCTGGAAGAGGACGCACTGGCTGATGCCTTTGAAGCGCAGCACGATGATTACAACAAAAT  >  NC_000913/4226506‑4227084

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