Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,956,003	G→A	D225N (GAT→AAT)	ilvA →	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,956,003	0	G	A	98.7%	245.5 / NA	77	D225N (GAT→AAT)	ilvA	threonine deaminase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (25/49);  total (27/50)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCG  >  NC_000913/3955764‑3956249                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         tGTCACAGCAGCAGGGGTTCACTTGGGGGCCGCCGGTCGACCAGCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAGCAACTGATGCCGCAAATCAAAGTGATCGCCGGATAAGCGTAAGACTCCGCGTGCCTGATAGCTGCGCTGTATGCGGGTCa                                                                                                                                                                                                                                                <  2:372889/250‑1 (MQ=255)    cacaGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCg                                                                                                                                                                                                                                            <  1:135993/251‑1 (MQ=255)    cacaGCAGCAGGGGTTCACCTGGGTGCCGCCGTGCGTCCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGGGGCGGTGCTGATCAAACAACTGATGCCGCAAAGCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATcc                                                                                                                                                                                                                                             <  2:56143/250‑1 (MQ=255)    cacaGCAGCAGGGGTTCACCTGGGTGCCGCCGGTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCg                                                                                                                                                                                                                                            <  1:248335/251‑1 (MQ=255)           gcagGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATc                                                                                                                                                                                                                                     <  1:363757/251‑1 (MQ=255)               gggTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGcc                                                                                                                                                                                                                                 <  2:230956/251‑1 (MQ=255)                 gTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCCGATCAAACAACTGATGCCGCAAATCCAAGTGATCGCCGGAGAAGCGGGAGACTCCGCCCGCCTGGCAGCCGCGCCGAATGCGGGTCCTCCGGGTGATCTGCcgc                                                                                                                                                                                                                               >  1:381840/1‑251 (MQ=255)                    cACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCACACACCTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGCCAACACCCCCTGCCTGAAAGCAGCGATGCATGCggtgcaaccgggtgctctgcccagcg                                                                                                                                                                                                                            >  1:388481/1‑227 (MQ=255)                                 ccgTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACCACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGGAGACTCCCCCCGCCTTAAAGCCGCGCTGGATGCGGGGCCTCCGGTTGGTTTGCCGCG‑CGTAGGGCTATTTgt                                                                                                                                                                                                              >  2:68730/1‑250 (MQ=255)                                     tCGACCAGCCGATGGTGATTGCCTGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGGCCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGTCTATTTGCTGa                                                                                                                                                                                                           <  2:68836/250‑1 (MQ=255)                                           aTCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTa                                                                                                                                                                                                    <  1:256164/251‑1 (MQ=255)                                             ccGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAg                                                                                                                                                                                                   <  2:78497/250‑1 (MQ=255)                                                   gTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTATGGCGATTGGCTGAAGGCGTAGCGGTaaa                                                                                                                                                                                            <  2:45959/251‑1 (MQ=255)                                                   gTGATTGCCGCGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGGATTTGTGCCAGCCGGCGGCGGCGGTCTGGCTGCTGGCGTAGCGGTACTGGTTAAAACACTGGAGGCGCAAACCAAAGTGATAGCCGTGGAAGCGGAAGAACACGCCTGCCTGACAAAAGCGGCGAATGCGGGGCCACCGGTTGAGTTGCCGCG‑GGTAGGGCTCTTTTGTGAAGGcagagcggcga                                                                                                                                                                                             >  2:248936/1‑240 (MQ=255)                                                         gCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCa                                                                                                                                                                                       <  1:11016/250‑1 (MQ=255)                                                         gCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATAAAACAACTGATGCCGCAAATCAAAGCGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTTCCGCGGCGT‑GGGCTATTTTGTGAAGGCGTAGCGGGAAAAAGCAt                                                                                                                                                                                      >  1:114208/1‑251 (MQ=255)                                                          ccGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGCCCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCAAACAACTGATGCCGCAAATCAAAAGGATCGCCGGAGAAGCGGAACAACACGCCTGCATGAAAGCAGCGGTGAATGCGGGTAATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGACGTAGCGGTACAACGCATc                                                                                                                                                                                     >  1:30996/1‑251 (MQ=255)                                                          ccGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATACAACAACTGATGCCGCCACTCAAAGTGATCCCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGAGGGGCATCCGGTTGCTCTGCCGCG‑AGTAGGGATATTTGCTGAAGGCGGAGCGGTAAAACGcctc                                                                                                                                                                                     >  1:231617/1‑248 (MQ=255)                                                             ggCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGt                                                                                                                                                                                 <  1:84479/251‑1 (MQ=255)                                                                 gggCACGCTGGCGCTGGAGCGGCTCCAGCAGGACGCCCATCTAGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGGCTGGCGGCTGGCGGGGCGGTGCGGATCAAACAACTGAGGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACGCCGCCTGTCTGTAAGCAGCGCTGAATGCGGGTCATCCGGTTGAGCTGGCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACg                                                                                                                                                                             <  2:388481/250‑1 (MQ=255)                                                                             gctggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAAGCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGtt                                                                                                                                                                 <  1:343067/251‑1 (MQ=255)                                                                                ggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTAt                                                                                                                                                              <  1:68730/251‑1 (MQ=255)                                                                                 gAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCCAAGTGATCGCCGCAGAAGCGGAAGACTCCGCCCGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGGTGATCTGCCGCG‑CGTAGGGCTATTTGCTTAAGGCGGAGCGGTAAAACGCATC‑GGTGGCGGAACCTTCCGTTTATg                                                                                                                                                             >  1:205957/1‑251 (MQ=255)                                                                                           agcagGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTAt                                                                                                                                                   <  1:401967/251‑1 (MQ=255)                                                                                             cagGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTATGGCTGCTGGCGTGGCGGCGCGGGTCACACAAATAATGCCGGAAAAAAAAATGATAACCGGAAAAGAGGGAGCCTCCGCACGCCTGTAAACCGGGCTGAATGCGGGGCATCCGGTGTAAATGCCGGG‑AGTAGGGGTATTTGGTGAAGGGGTAGCGGGAAAAAG‑ATAAGGGGACGAAAGCTTCCGATTATGCCAGGAGTGtct                                                                                                                                                 >  1:402198/1‑251 (MQ=255)                                                                                               ggACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTcg                                                                                                                                               <  1:268938/251‑1 (MQ=255)                                                                                               ggACGCCCATCTCGACCGCGTAGTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGAGCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCAGCCGGGTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTAtcgcg                                                                                                                                               <  1:268946/251‑4 (MQ=255)                                                                                                 aCGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTcgac                                                                                                                                             <  1:15159/251‑1 (MQ=255)                                                                                                  cGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTcgacg                                                                                                                                            <  1:308332/251‑1 (MQ=255)                                                                                                   gCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGGGATCGCCGTAGAAGCGGAAGACTCCCCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGGAAAAAGCATC‑GGTGACGAAACCTTACGTTTATGCCCGGAGAATCTcgacgc                                                                                                                                           >  1:124657/1‑250 (MQ=255)                                                                                                    cccATCTCGACCTCGTATTGGTGCCATTCGGCTGCTGCTGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACATCTGATGCCGCAAAACAAAGTGATCGCCGTAGAAGCTGAAGACACCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAATACGCATC‑GGTGACGAAACCTTCCGGTTCTGCCAGGAGTATCTCgacgac                                                                                                                                          <  2:276015/251‑1 (MQ=255)                                                                                                       aTCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGGGACGAAACCTTCCGTTTATGCCAGGAGGATCTCGACGAcatc                                                                                                                                       >  1:5851/1‑251 (MQ=255)                                                                                                           cGACCGCTTACTTGTGCGAGTCGGCGGCGGCGGTCTGGCTGCTGTCGTGGCGTTGCTGATCAAACAACTGATGCCGCAAATCAAAGTTATCGCCGTAGAAGTGGAAGACTCCGCCTTCCTGAAAGCATCTCGGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTTCTGAACGCGTATCGGTAGAACGCTTC‑GGTTACGAATCCTTCCGTTTATGGCTGTAGTATCTCGACGACAgcatca                                                                                                                                   <  2:18138/251‑1 (MQ=255)                                                                                                           cGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGTCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGTCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGAcatca                                                                                                                                      <  1:249348/248‑1 (MQ=255)                                                                                                           cGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGAcatca                                                                                                                                      <  1:52580/248‑1 (MQ=255)                                                                                                             aCCGCGTATTTGTGCCAGTCGGCGGCGGCGGTATGGCTGCTGGCGTGGCGGTGCTGATCTAACAAATGGTGCCGCAAATCAAATGGATCGCCGTAGAAACGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGAGGGTCATCCGGGTGATCTGCCGCG‑CGTAGGGCCATTTGCTGAAGGCGGAGCGGTAAAAAGCCTC‑GGTGACGAAACCTTCCGTTTATGCCAGGGGTATCTCGACGACATCATCAcc                                                                                                                                 >  1:3343/1‑251 (MQ=255)                                                                                                                      ttGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAg                                                                                                                         <  1:149840/250‑1 (MQ=255)                                                                                                                       tgtgCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCCAATCAAAGTGATCGCCGTTGAAGCGGAAGGCTCCGCCTGCCTGAAAGCAGAGCTGAATGCGGGTCCTCCGGGTTATCTGCCGCG‑CGTAGGGCTATTTGCCGAAGGCGTAGCGGTAAAACGCATC‑GGGGACGAAACCTTCCGTTTATGCCACGAGTATCTCGACGACATCATCACCGTCggtaacg                                                                                                                       >  1:142016/1‑245 (MQ=255)                                                                                                                          gCCAGTCGGCGGCGGCGGTCTGGCTGCGGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcgatg                                                                                                                    <  1:152670/251‑1 (MQ=255)                                                                                                                                 ggcggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCtg                                                                                                             <  1:314706/251‑1 (MQ=255)                                                                                                                                  gcggcggctgTCTGGCGGCGGGCGTGGCGGGGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCTCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACTAAACCTTCCGTTTATGCCAGGGGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCtgt                                                                                                            <  2:350068/251‑1 (MQ=255)                                                                                                                                   cggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCtgtg                                                                                                           <  1:267956/251‑1 (MQ=255)                                                                                                                                    ggcggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTgc                                                                                                          <  1:121085/251‑1 (MQ=255)                                                                                                                                       ggcggTCTGGCTGCTGGCGTGGCGGTGGTGATCCAACAACTGATGCCGCCAATCAAAGGGATCGCCGTAGAAGCGGAAAACTCCCCCCGCCTGAAAGCAGCGCTGAATGCGGGGCCTCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGGAGCGGTAAAACCCATC‑GGTGACGGAACCCTTCGTTTATGCCCGGAGGATCTCCGCGACATCATCACCGTCGATAGCGGTGCGATCTGTgcggg                                                                                                       >  2:85849/1‑250 (MQ=255)                                                                                                                                          ggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGAt                                                                                                    <  1:380083/251‑1 (MQ=255)                                                                                                                                                          gggcgGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCTTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTGTATGCCAGGAGTATCTCGACGACAGCATCACCGTCGATAGCGATGCGATCTGTGCTGCGATTAAGGATTTATTCGaa                                                                                    <  2:135015/250‑1 (MQ=255)                                                                                                                                                               gTGCTGATCAAACAACTGATGCCGCAAAACAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGAtgt                                                                               <  2:111348/251‑1 (MQ=255)                                                                                                                                                                 gCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGAtgtg                                                                              <  1:326297/250‑1 (MQ=255)                                                                                                                                                                    gATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTgcg                                                                            <  1:98093/249‑1 (MQ=255)                                                                                                                                                                       cAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGcgcgcg                                                                        <  2:372161/250‑1 (MQ=255)                                                                                                                                                                       cAAACAACTGAGGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGTTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTGCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGTAGGATTTATTCGAAGATGTGCGCGCgg                                                                       <  2:136384/251‑1 (MQ=255)                                                                                                                                                                           caacTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGc                                                                   >  1:380331/1‑251 (MQ=255)                                                                                                                                                                             acTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGTAGGATTTATTCGAAGATGTGCGCGCGGTGGCgg                                                                 >  2:129808/1‑251 (MQ=255)                                                                                                                                                                               tGATGCCGCAAATAAAAGTGATCGCCGTAGCAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGa                                                                >  1:68920/1‑250 (MQ=255)                                                                                                                                                                                gATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGGTCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGaa                                                               <  2:48846/250‑1 (MQ=255)                                                                                                                                                                                    ccGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCctc                                                          >  1:26156/1‑251 (MQ=255)                                                                                                                                                                                     cGCAAATCAAAGTGATCGCTGTAGAAGCGGAAGACGCCGACTGCCTGAAAGCAGCGCTGAATGCGGGTCATGCGTTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAGCGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGAGTTATTCGAAGATGTGCGCGAGTTGGCGGATCCctct                                                         <  1:277334/251‑1 (MQ=255)                                                                                                                                                                                        aaaTCAAAGTGCTCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAACGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGACGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGCTTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGgc                                                      >  2:167539/1‑251 (MQ=255)                                                                                                                                                                                                   aTCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggcgct                                           <  1:296124/251‑1 (MQ=255)                                                                                                                                                                                                     cGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGcgctggcgctgg                                         <  1:357298/251‑1 (MQ=255)                                                                                                                                                                                                     cGCCGGAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGcgctggcgctgg                                         <  2:26156/251‑1 (MQ=255)                                                                                                                                                                                                      gCCGTAGAAGCGGAAGACTCCGCCTGCCTGAATGCAGCGCTGAATGCGGGTCAGCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCACCACCGTCGATAGCGATGCGATCTGGGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCgctggcgctggc                                        <  2:142016/251‑1 (MQ=255)                                                                                                                                                                                                        cGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCgg                                      >  2:151220/1‑251 (MQ=255)                                                                                                                                                                                                        cGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCgg                                      >  1:64314/1‑251 (MQ=255)                                                                                                                                                                                                           aGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCATCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGaa                                   <  1:246029/251‑1 (MQ=255)                                                                                                                                                                                                             aaGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGCAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGAGGCGGACCCCTCTGGCGGGCTGGCGCTGGCGGGAATg                                 >  1:240687/1‑251 (MQ=255)                                                                                                                                                                                                              aGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGa                                <  1:45859/251‑1 (MQ=255)                                                                                                                                                                                                                cGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGCAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGaa                               >  1:383117/1‑250 (MQ=255)                                                                                                                                                                                                                   aaGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGaaaaa                            >  1:344611/1‑250 (MQ=255)                                                                                                                                                                                                                    aGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCAGC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAat                          <  1:145713/251‑1 (MQ=255)                                                                                                                                                                                                                         ccGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATc                     <  1:151220/251‑1 (MQ=255)                                                                                                                                                                                                                                   gagagtaacgcTGAATGCGGGTCATCCGGGTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTCGCGGTACAACGCATC‑GGTGACGAATCCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGTTAGCGCTGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGCTGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGGAGGAAAAACTACATGGCCCTGcac            >  2:358330/9‑250 (MQ=255)                                                                                                                                                                                                                                   gAAAGCAGCGCTGAATGCGGGTCATCCGGTTGGTCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGca             <  1:334845/249‑1 (MQ=255)                                                                                                                                                                                                                                         agcGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CTTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTc     <  2:356574/251‑1 (MQ=255)                                                                                                                                                                                                                                         agcGCTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTc     <  1:370607/251‑1 (MQ=255)                                                                                                                                                                                                                                          gcgcTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACAtt      >  1:22308/1‑249 (MQ=255)                                                                                                                                                                                                                                            gcTGAATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTcgcg  >  1:119710/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         TGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCG‑CGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATC‑GGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAACATTCGCG  >  NC_000913/3955764‑3956249

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