Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,558,907	T→A	K207M (AAG→ATG)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,558,907	0	T	A	92.0%	245.3 / ‑4.6	89	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/7);  new base A (45/36);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 5.08e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCT  >  NC_000913/4558669‑4559149                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       gCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATCCCCGCCTGAACGGCGCGGGCAATACATTCGGCAGGGGCGAACGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGc                                                                                                                                                                                                                                            >  1:498916/1‑251 (MQ=255) gCCGTCAGAGCTGAGGGTGACGCGAGCCAGCGGACTACCCGCCTGAACGGCGCGGGCAATACCTTCGGACGGGGCGACCGGTTCGTCA‑ACGCTGCTGGTGATATCGCTGGTGCCGCCTTTGCGCGAGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGAAACGCCATTTTACTGTCGc                                                                                                                                                                                                                                            >  2:252180/1‑251 (MQ=255)  ccGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGcc                                                                                                                                                                                                                                           <  1:93556/251‑1 (MQ=255)   cGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCAGCCTGCACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑CTGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGCTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGccc                                                                                                                                                                                                                                          >  1:481769/1‑251 (MQ=255)                   gACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGa                                                                                                                                                                                                                         <  1:425967/251‑1 (MQ=255)                       cgTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGcc                                                                                                                                                                                                                     <  1:516296/251‑1 (MQ=255)                        gTGCCAGCGGCATACCCGCCTGAACGGCACGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGCTATCGATGGTGCCGCCTTTGCGCGCGACCTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCCGATCGGCACATCGCAGTTTACCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCCTGTGGacctcggttaggcc                                                                                                                                                                                                                      >  2:537248/1‑237 (MQ=255)                        gTGCCAGCGGAATACCCGCCTGACCGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AAAGGTACGTTGCGGTTAACGTGGGTCGGCAGCCGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGcc                                                                                                                                                                                                                     >  2:206989/1‑250 (MQ=255)                        gTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCa                                                                                                                                                                                                                    <  1:104245/251‑1 (MQ=255)                            cAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATATATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGtt                                                                                                                                                                                                                <  1:79221/251‑1 (MQ=255)                             ggCGGAATCCCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGttt                                                                                                                                                                                                               <  2:113778/250‑1 (MQ=255)                             aGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGttt                                                                                                                                                                                                               <  1:157419/251‑1 (MQ=255)                                ggAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑CTGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCACAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCC‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCCTGTGGAA‑CACGGTGACGCCAGGTTT‑Acc                                                                                                                                                                                                           >  2:233750/1‑251 (MQ=255)                                   aTACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCCGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCCGGTTT‑Accgcc                                                                                                                                                                                                        >  2:438088/1‑251 (MQ=255)                                        cGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCa                                                                                                                                                                                                 <  1:142866/251‑1 (MQ=255)                                         gCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCACGTGGAA‑CACGGTGACGCCAGGTTT‑CCCGCC‑G‑AGCa                                                                                                                                                                                                 >  1:346916/1‑250 (MQ=255)                                                            aaCCTTCGGCAGGGGCGTCCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGTTACGTGTCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACAGCGCATGTTTCCA‑GCAGGTCATAAATAGGCGGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGTTCACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCt                                                                                                                                                                             <  1:512302/250‑1 (MQ=255)                                                             aCCTTCGGCAGGGGCGACCGGTTAGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTg                                                                                                                                                                            >  1:453348/1‑251 (MQ=255)                                                             aCCTTCGGCAAGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCt                                                                                                                                                                             <  1:524325/250‑1 (MQ=255)                                                                 tCGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTgcc                                                                                                                                                                          <  1:363960/249‑1 (MQ=255)                                                                 tCGGCAGGGGCGACCGGCTCGACA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCCAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCCTAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTgccgc                                                                                                                                                                        >  2:319945/1‑251 (MQ=255)                                                                  cGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTgccgcc                                                                                                                                                                       >  1:360463/1‑251 (MQ=255)                                                                  cGGCAGGGGCCACCGGTTC‑TGACATGCGGCTGGTGATACCGTTGGTGCCGCCTTGTAGTGGGAACTCCAGCGCCTGCGCAAAC‑AACGGTACGTTGCGGGTAACGTGGGTCGGCAGAAGCTTGCTGATCGGTACATCCCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCTTGTTACTGTCGCCCATGTGGAA‑CGCGGTGACGCCGGGTTT‑CACGCC‑G‑AGCAAACCGCCAACGCGGGATTCTgccgcc                                                                                                                                                                       <  2:330790/251‑1 (MQ=255)                                                                      aGGGGCGCACGGTTCGTCA‑TTGTTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AGCGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTTATCGGAACATCGCAGTTTTCCA‑GCAGTTCATAAATAGGCGGTACCGCCATTTTACTGTCGCCCATTTGGAA‑CACGTTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCatat                                                                                                                                                                   <  2:69046/251‑1 (MQ=255)                                                                          gCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGc                                                                                                                                                               >  1:292213/1‑251 (MQ=255)                                                                          gCGACCGGTTCGTCA‑ATGCTGCTGGCGATATCGATGGAACCGCCTTAACGAGCAAACTCCAGCGCCTGCTAAAAC‑CACGGTACGTTGCGGTTAACGTGGGTAGGCAGCAGCTTGCTGATCGGCACATCGCAGTCTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATCTTACTGTAGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCACACCGCCAACGCGGGATTCTGCCGCCATATTgcc                                                                                                                                                               >  1:486709/1‑249 (MQ=255)                                                                               cGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAggg                                                                                                                                                          >  1:292214/1‑250 (MQ=255)                                                                                  ttCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGAt                                                                                                                                                       >  1:12142/1‑251 (MQ=255)                                                                                     gTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTACCGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATaaa                                                                                                                                                    <  1:393108/251‑1 (MQ=255)                                                                                     gTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATaa                                                                                                                                                     >  1:220033/1‑250 (MQ=255)                                                                                     gTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACGCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGGGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGAGTCTGCCGCCATATTGGCCAGGTGATaaa                                                                                                                                                    <  1:438088/251‑1 (MQ=255)                                                                                      tCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAAc                                                                                                                                                   >  1:530699/1‑251 (MQ=255)                                                                                        caaggctgctGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTGGCGGTTAACGGGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAAc                                                                                                                                                   <  1:233750/246‑1 (MQ=255)                                                                                        a‑aTGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGt                                                                                                                                                 >  1:484922/1‑251 (MQ=255)                                                                                          atgcggctgggtataccGATGGTGCCGCCTTGGCGCGCGAACCCCAGCGCCTGCCCAAAC‑AACGGTACGTTGCGGTTAACGGGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAAAAGGCTGTAACGCCATTTTCCTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑CCCGCC‑G‑AGCAAACCGCCACCGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTc                                                                                                                                                <  2:128292/235‑1 (MQ=255)                                                                                               gctgGTGATATCGATGGTGCCGCTTTGGCGCGCGAACTCCGGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTAATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCGGTAACGCCATTTTACTGTCGCCCTTGTGGAA‑CACGTTGACGCCAGGTTT‑GCCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCgg                                                                                                                                             <  1:497393/249‑1 (MQ=255)                                                                                                 gggTGATATCGATGGTGCCGCCTGTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑ACCGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CCCGGGGACGCCGGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCg                                                                                                                                         <  2:187583/250‑1 (MQ=255)                                                                                                   gTGATATCGATGGCGCCGCCTTTACGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGCGATAAACGTCCGGTGCGGc                                                                                                                                       >  2:281592/1‑251 (MQ=255)                                                                                                     gATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGc                                                                                                                                       >  1:407445/1‑249 (MQ=255)                                                                                                       tatCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGCCGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGa                                                                                                                                    >  2:275072/1‑250 (MQ=255)                                                                                                       tatCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGa                                                                                                                                    >  1:165686/1‑250 (MQ=255)                                                                                                        cgcGATGGCGCCGCCGTTGCGCGCCAACCCCAGCGCCTGC‑CAAACCACCGGTCCGTTGCGGTTACCGTGGGTCGGCAGCAGCTTCCGGTTCGGCCCATCGCAGTTTTCCA‑ACAGGTCATAAATAGGCGGTAACGCCATTTTACTTTCCCCCATGTGGAC‑CTCGGGGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTTCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAAc                                                                                                                                  <  2:467317/249‑1 (MQ=255)                                                                                                               gTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATc                                                                                                                           >  1:56654/1‑251 (MQ=255)                                                                                                                      ctgtGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGTCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCGGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGAt                                                                                                                                                                                <  2:304555/188‑1 (MQ=255)                                                                                                                      cTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGAt                                                                                                                                                                                >  1:304555/1‑191 (MQ=255)                                                                                                                         tGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGgc                                                                                                                  <  1:275072/250‑1 (MQ=255)                                                                                                                         tGCGCGCGAACCCCAGCGCCGGCTAAAAC‑AACGGGCGGTGGCGGTTAGCGTGGGTCGGCAGGAGCTTGCTGATCGGCACAGCGCATTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGTTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGgc                                                                                                                  <  2:167842/250‑1 (MQ=255)                                                                                                                            gcgcgAACCCCAGCGCCGGCTCAAAC‑AACGGTACGTTGCGTTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCATTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGGGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCAt                                                                                                              <  2:165686/251‑1 (MQ=255)                                                                                                                                 aaCTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCAc                                                                                                         >  2:210683/1‑251 (MQ=255)                                                                                                                                  aCTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACg                                                                                                        >  1:7455/1‑251 (MQ=255)                                                                                                                                              gCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcg                                                                                             >  1:363772/1‑250 (MQ=255)                                                                                                                                                 aaaaC‑AACGGGACGTTGCGGTTAACGGGGGCGGGCGGCGGTTTGCGTATCGGCAAATCGAGGTTTCCCG‑GCAGGTCATAAATGGGCTGAAACGCCTTTTTAGTGCCGCCCATGGGGAC‑CACGGTGACGC‑AGGTTTTACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcgatc                                                                                          <  2:265516/249‑1 (MQ=255)                                                                                                                                                   aac‑aacTGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCCCATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcgatcgat                                                                                       >  1:70944/1‑251 (MQ=255)                                                                                                                                                     c‑aacGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATaa                                                                                     <  1:206442/251‑1 (MQ=255)                                                                                                                                                       aacGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAAt                                                                                    >  1:67621/1‑251 (MQ=255)                                                                                                                                                           gTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACCTCGCAGTTTTCCA‑GCAAGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑CGCAAACCGCCAACGCGGGATTCTGACGCCATCTTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGcc                                                                                >  1:212362/1‑251 (MQ=255)                                                                                                                                                                 cttcggacctgTGGGTCGGCAGCAGCTTGCTGATCGGCACACCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGGGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGt                                                                            <  1:412410/239‑1 (MQ=255)                                                                                                                                                                  gCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGGGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCtt                                                                         <  1:220136/251‑1 (MQ=255)                                                                                                                                                                           gTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTACCGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCA‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATACACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACgg                                                                >  1:244901/1‑251 (MQ=255)                                                                                                                                                                           gTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGt                                                                            >  2:412410/1‑239 (MQ=255)                                                                                                                                                                           gTGGGCCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACgg                                                                <  1:418522/251‑1 (MQ=255)                                                                                                                                                                            tGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACgg                                                                >  2:529801/1‑250 (MQ=255)                                                                                                                                                                             gggTCGGCAGCAGCTTGCTGTTCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCCCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGaa                                                              <  1:395499/251‑1 (MQ=255)                                                                                                                                                                                 cGGCAGCAGCTTGCGGATCGCCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTCCTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAAccc                                                           <  1:168929/250‑1 (MQ=255)                                                                                                                                                                                   gcagcaGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCCCGTAATAAATAGGGTGTAACGCCCTTTTACTGTCGCCCATGCGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCACACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCCGAACGGTGATCAGACATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGt                                                         >  1:71981/1‑250 (MQ=255)                                                                                                                                                                                   gcagaagtttGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTa                                                        <  1:298257/243‑1 (MQ=255)                                                                                                                                                                                       cagcTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATgg                                                    >  1:32011/1‑251 (MQ=255)                                                                                                                                                                                         gcTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTg                                                  >  1:485614/1‑251 (MQ=255)                                                                                                                                                                                            tGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑CGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCgg                                               >  1:32043/1‑251 (MQ=255)                                                                                                                                                                                            tGCTGATCGGCACATCGCAGTTTCCCA‑GCAGGTCATAAATAGGCGGTAACGCCATTTTACGGTCGCCCATGTGGAA‑CACGGTAACGCCGGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGc                                                 <  1:327546/249‑1 (MQ=255)                                                                                                                                                                                               tGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGa                                             >  2:150977/1‑250 (MQ=255)                                                                                                                                                                                                gttgggCACATC‑CCGTTTTCCAAGCAGTTCATAAATAGGCTGTACCGCTTTTTTACTGAGGCCGATGTGGAA‑CACGGTGACTCCAGTTTT‑GCCGCC‑G‑AGCAACCCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGaa                                            <  2:463066/246‑1 (MQ=255)                                                                                                                                                                                                 aTCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAgg                                          >  2:441523/1‑251 (MQ=255)                                                                                                                                                                                                      cacaTCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACAt                                     >  2:347510/1‑251 (MQ=255)                                                                                                                                                                                                        caTCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATg                                    >  1:436111/1‑250 (MQ=255)                                                                                                                                                                                                            ccAGTGTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCATTGTACTGCCGCCCATGTGGAA‑CACGGTGACGCCGGGTTT‑ACCGTC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGAGAAACGTCCGGGGCGGCAGAACGGTGATCAGAGAGGGCGCATTGCACGCCAATCACACGGTCGATAAGCGCCACGTCTTTTTCCTCGGAACCCGTAATGGTGCGGGAAGGGACATGATAAg                               <  1:252180/250‑1 (MQ=255)                                                                                                                                                                                                              aGTTTTCCA‑GCAGGTCATAAATAAGCTGTAACGCCATATTACTGTCGACAATGTGGAA‑CACGGAGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGCTTCTACCGCCATATTGGCCAAGTGATAAACGTCCGGTGCGGCAGGACGGGGATCAGAGATGGCGCATTTCACGCCAATCAAACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGGGCGGGAAGGGACATGATAAgcg                             >  1:331286/1‑251 (MQ=255)                                                                                                                                                                                                                  tcccA‑GCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGCCGCCAGGTTT‑CCCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCTCGTCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGATCCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCg                          <  1:516688/248‑1 (MQ=255)                                                                                                                                                                                                                     ca‑gcagGTAAAAAATAGGTTGTAACCCATTTTTACTGCCCCCCAGGTGAAA‑CACGGTGACGCCAGGTTT‑CCCGCC‑G‑AGCAAACCGCCACCGCGGGATTCTGCCGCCATATTGGCGAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGAAATGGTGCGGTAAGGGACATGATAAGCGCCGGTCa                      <  2:67621/251‑1 (MQ=255)                                                                                                                                                                                                                           ggccATAAATAGGCTGTAACGCGATTTTACGTTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTC‑GCCGCC‑G‑AGCAAACCGGCAACGCGGGAt                                                                                                                                                                                <  2:518109/89‑1 (MQ=255)                                                                                                                                                                                                                           ggTCATAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGAt                                                                                                                                                                                >  1:518109/1‑92 (MQ=255)                                                                                                                                                                                                                              cATAAATAGGCTGAAAGCCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCTCC‑G‑AGCAAACCGCAAACGGGGGATTCTGCCTCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCa               <  1:474038/250‑1 (MQ=255)                                                                                                                                                                                                                                tAAATAGGCTGTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACAGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGCACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCCCGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGCTAAGCGCCGGTCAGCATCCAGgc            >  2:489270/1‑251 (MQ=255)                                                                                                                                                                                                                                aaaaTAGGCTGAAACGCCTTTTTACTTTCGCCCATGGGGAA‑CTGGTTGGCGCC‑GGTTT‑TACGCC‑GGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGgc            <  2:407445/250‑1 (MQ=255)                                                                                                                                                                                                                                   atgggCTGTAACGCCTTTTTCCGGTCGCC‑GTGGGGAAACCCGGTGACGCC‑GGTTT‑ACCCCCCG‑AGCAAACCGCCACCGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCt         <  2:220033/248‑1 (MQ=255)                                                                                                                                                                                                                                      ggcggTAACGCCATTTTGCTGGCGCCCGTGTGGAA‑CACGGTGCCGCCAGGTTT‑ACCGCC‑G‑CGAAAACCGCCAACGCGGGATTCTGCCGCCATATGGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGAGCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGAt      <  2:56654/247‑1 (MQ=255)                                                                                                                                                                                                                                          gTAACGCCTTTTTACGGGCGCCCATGTGGAA‑CTCGGTGCCGCCGGGTGT‑CCCGCC‑G‑AGCAACCCGCCACCGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGcc   <  2:386876/250‑1 (MQ=255)                                                                                                                                                                                                                                          gTAACGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCct  >  1:166010/1‑251 (MQ=255)                                                                                                                                                                                                                                             aCGCCATTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGAt      >  1:321931/1‑244 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       GCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCA‑ATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAC‑AACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCA‑GCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAA‑CACGGTGACGCCAGGTTT‑ACCGCC‑G‑AGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCT  >  NC_000913/4558669‑4559149

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