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	93.6%	225.2 / ‑5.6	80	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/5);  new base A (36/37);  total (37/43)
Fisher's exact test for biased strand distribution p-value = 5.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCT‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTC  >  NC_000913/4558663‑4559151                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      aCCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTAc                                                                                                                                                                                                                                                >  1:363061/1‑251 (MQ=255)       gccggcggaggggAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGGACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGc                                                                                                                                                                                                                                          <  2:407378/239‑1 (MQ=255)              gAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCAtgtg                                                                                                                                                                                                                                   >  1:362110/1‑251 (MQ=255)               ggcggAGGGTGCCGCGGGCCAGCTGAAGACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCGGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACTGTACGTTGCGGTT‑AACGGGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTGTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTgg                                                                                                                                                                                                                                  <  1:275264/247‑1 (MQ=255)                       gggACGCGGGCCCGCGGAATACCCGCCTGAACGGCGCGGTCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGt                                                                                                                                                                                                                          <  2:176688/249‑1 (MQ=255)                       gTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCACTGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGt                                                                                                                                                                                                                          >  1:442230/1‑251 (MQ=255)                          aCGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGAc                                                                                                                                                                                                                       <  1:450789/251‑1 (MQ=255)                             cgTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTTAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                                                                                                                                                                                                    <  2:150147/251‑1 (MQ=255)                              gTGCCAGCGGCATAACCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGACTGCTCAAACAACGGTACGTTGCGGAT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCC‑TTTTACTGTCGCCCATGTGGAACACGGTGACGccc                                                                                                                                                                                                                   >  2:166996/1‑250 (MQ=255)                              gTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                                                                                                                                                                                                    >  1:48760/1‑250 (MQ=255)                                gCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCCTGTGGAACACGGTGACGCCAgg                                                                                                                                                                                                                 >  2:51053/1‑251 (MQ=255)                                gCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAgg                                                                                                                                                                                                                 >  1:288257/1‑251 (MQ=255)                                    gCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGttt                                                                                                                                                                                                              >  2:322777/1‑250 (MQ=255)                                     cGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGTACATCGCAGCTTTCCAGCAGGTCATAAATAGGCTGTAACGCAA‑TTTTACTGTCGCCCATGGGGAAAACGGTGGCGCTAGGTTTAc                                                                                                                                                                                                            >  2:251305/1‑251 (MQ=255)                                     cGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAc                                                                                                                                                                                                            >  1:338170/1‑251 (MQ=255)                                     cGGAATACCCGCCTGAACGGCGCGGGCAATACCTGCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTAc                                                                                                                                                                                                            <  2:146354/251‑1 (MQ=255)                                                     acgccgcggcgaaTACCTTCGGCAGGGCCGACCGGTTCGTCAAGGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCGGCTCAAACAACGGTGCGTTGCGGTT‑AACGGGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑GTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc                                                                                                                                                                                            <  2:145839/240‑1 (MQ=255)                                                     aCGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc                                                                                                                                                                                            >  1:99041/1‑251 (MQ=255)                                                      cGGCGCGGGCAATACCTTCGGCGGGGGCGACGGGTTCGTCAATGCGGCTGGTG‑TTATCGATGGTGCCGCCTTGGCGCGCGACCTCCCGCGCCTGCTCAAACAACGGTACGTGGCGGTT‑AACGGGGGTCGGCAGCAGCTTGCGGACCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTCCCGCCGAGCAAACCGCCa                                                                                                                                                                                           <  2:369776/251‑1 (MQ=255)                                                         cgcgGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAAcg                                                                                                                                                                                        <  1:288795/251‑1 (MQ=255)                                                                  tACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCt                                                                                                                                                                               >  2:100631/1‑251 (MQ=255)                                                                   cccTTCGGCAGGGGCGACCGGTTCGTCAAGGCTGCTGGTG‑ATATCGATGGTGCCGCCTTGGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCt                                                                                                                                                                               <  2:357167/249‑1 (MQ=255)                                                                    ccTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCGTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTgc                                                                                                                                                                             <  2:143288/251‑1 (MQ=255)                                                                       tCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTgcc                                                                                                                                                                            <  2:150660/249‑1 (MQ=255)                                                                            aGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATAGCGAGGGTGCCGCCTTTGCGCGCGAACACCAGCGCCTGCGCAAACAACGGTACGTTGCGGTT‑AACGTGTGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACGGTCGCCCATGGGGGACACGGGGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCatat                                                                                                                                                                     <  1:452847/251‑1 (MQ=255)                                                                             ggggggACCGGTCCGGCAATGCTGCTGGTG‑AAATGGATGGGGCCCCCTTGGCGCGCGAACTCCAGCGCCGGCTAAACCAACGGTACGTGGCGGTT‑AACGTGGGTCGGCAGCAGTTTGCTGATCGCAACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGAAACCCGGTGCCGCCAGTTTTACCGCCGAACAAACCGCCAACGCGGGATTCTGCCGCCATAtt                                                                                                                                                                    <  2:132317/251‑1 (MQ=255)                                                                             ggggCGCCCGGTCCGTCAATGCTGC‑GGTGGATAGCGATGGGGCCGCCTTTGCGCGCGACCCCCAGCCCCGGCTCAAACAACGGTACGTGGCGGTT‑ACCGTGGGCCGGCAGCAGTTTGCTGATCGGCACATCGCAGTTTGCCAGCGGGTCATAAATAGGCTGTAACGCCT‑TTTTACTGTCGCCCATGTGGAACGCGTGGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATAtt                                                                                                                                                                    <  2:53127/251‑1 (MQ=255)                                                                                  gACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGcc                                                                                                                                                                >  1:403936/1‑250 (MQ=255)                                                                                           gTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATCGGCCAGGTGATaa                                                                                                                                                       >  2:339331/1‑250 (MQ=255)                                                                                            tCAATGCTGCTGGTG‑ATATCGATGGAGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCACATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATAATGGCCAGGTGATAAAc                                                                                                                                                     >  1:224109/1‑251 (MQ=255)                                                                                                 gctgctGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCTCCGGCACAAACAACGGTACGTGGCGGTT‑AACGTGGGGCGGCAGCAGCTTGCGGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCAGGTGGAACGGGGTGACGCCAGGTTTACCGCCGGGCAAGCCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTcc                                                                                                                                                 <  2:276716/250‑1 (MQ=255)                                                                                                 gctgctGGTG‑ATATAGATGGTGCAGACCTTGCGCGAGAACTCCAGCGTCTGCTCATACAACGGTACGATGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGAACATCGCAGTTTTCCAGCAGGTACTAAATCGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCACCGCGGGATTATGCCGCCATATTGGCCAGGTGATAAACGACCg                                                                                                                                                >  2:458406/1‑251 (MQ=255)                                                                                                              atCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAAc                                                                                                                                    >  2:360534/1‑251 (MQ=255)                                                                                                                cGATGGTGCCGCCTTTGCGCGCGACCTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACg                                                                                                                                   <  1:51053/250‑1 (MQ=255)                                                                                                                     gTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACTTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCATGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATc                                                                                                                             <  1:419223/251‑1 (MQ=255)                                                                                                                     gTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATc                                                                                                                             >  1:202721/1‑251 (MQ=255)                                                                                                                               gccccccgAACCCCAGCGCCTGCTCAAACAGCGGTACGTTGCGGTT‑ACCTTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCAAAAATAGGCTTTACCCCCT‑TTTTACTGTCGCCCATGTGGAACACGGTGCCGCCAGGTTTACCGCCGAGCAAACCGCCACCGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGgc                                                                                                                    <  2:322196/244‑1 (MQ=255)                                                                                                                                        aCTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGGGGGTCGGCAGCAGCTTGCTGATCGACACATCGCAGTTTTCCAGCAGGTCATAAATAGGTGGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCATCGCGGGAGTCGGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACg                                                                                                          <  1:430982/251‑1 (MQ=255)                                                                                                                                            cggcgcCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGAAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCaa                                                                                                      <  1:20095/249‑1 (MQ=255)                                                                                                                                              gcgcCTGCTCAAACAACG‑TAGCCTGGGGTTGAACGGGGGTGGCCGGCAGCTGGCTGATCGGCCAATCGCAGTTTTCCAGCAGGTAATAAAAAGGCTGTAACGCCA‑TTTTGCTGTCCCCCAGGGGGACCACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATc                                                                                                    <  2:55461/251‑1 (MQ=255)                                                                                                                                               cgcCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCGGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAAttc                                                                                                   >  2:347635/1‑249 (MQ=255)                                                                                                                                                gcCTGCTCAACCAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGAGCGGCCCACCGCATTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATcac                                                                                                  <  1:360789/251‑1 (MQ=255)                                                                                                                                                    gCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGCCAGAAGCTTGCGGATCGGCACATCGCAGTTTTCCAGCAGGTAATAAATAGGCTGTAACGCCA‑TTTTACTGGCCCCCATGGGGAACCCGGTGACGCCAGGTTTCCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcg                                                                                               <  2:311586/250‑1 (MQ=255)                                                                                                                                                      ccAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcgatc                                                                                            <  1:111918/250‑1 (MQ=255)                                                                                                                                                       cAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcgatc                                                                                            <  1:383048/250‑1 (MQ=255)                                                                                                                                                        aaaCAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCCGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCg                                                                                                                                                                                                                                           >  1:68281/1‑106 (MQ=255)                                                                                                                                                        aaaCAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCg                                                                                                                                                                                                                                           <  2:68281/106‑1 (MQ=255)                                                                                                                                                              cGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCg                                                                                    >  2:437230/1‑251 (MQ=255)                                                                                                                                                                   cGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACg                                                                               >  2:127830/1‑251 (MQ=255)                                                                                                                                                                         ggTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCtttt                                                                         >  1:261487/1‑251 (MQ=255)                                                                                                                                                                         ggTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCtttt                                                                         >  1:391445/1‑251 (MQ=255)                                                                                                                                                                         ggTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCtttt                                                                         >  2:137304/1‑251 (MQ=255)                                                                                                                                                                             gaacgtgggccggcggcatctgccggATCGGCACATCGCAGTTTTCCAGCAGGTCATAAAAGGGCTGTAACGCAT‑TTTTCCGTTCGCCCATGTGGAACAGGGTGACGCCAGTTTTCCCGCCGAGCAATCCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGGGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTcc                                                                      <  2:71864/226‑1 (MQ=255)                                                                                                                                                                              aaCGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGGGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTccc                                                                     >  2:122219/1‑250 (MQ=255)                                                                                                                                                                                    gggcGGCCGCAGCTTGCTGATCGGCACAGCGCAGTTTCCCAGCAGGTCAAAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCTGGTGGAACACTGAGAGGCCGGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGTCCAGGTGTTAATCGTCCGGTGCTGCAGAACGGTGATCATAGATGGCGCATTTCACGCCAATCTCACGATCGATAATCGCCACGTCTTTTTCCACGGAAc                                                               <  1:458406/248‑1 (MQ=255)                                                                                                                                                                                    ggTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACCAGATCGATAATCGCCACGCCTTTTTCCACGGAAc                                                               >  2:355124/1‑251 (MQ=255)                                                                                                                                                                                       cGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATGGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAAccc                                                             <  1:137961/250‑1 (MQ=255)                                                                                                                                                                                        ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCCGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGt                                                           >  1:401914/1‑251 (MQ=255)                                                                                                                                                                                        ggCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGt                                                           >  1:327893/1‑251 (MQ=255)                                                                                                                                                                                        gcgggctgcttcctgtgcgccccaTCGCAGTTTTCCAGGAGGTAATAAAGGGGC‑GTAAACCCTATTTTACGGTCGCCCGTGGGAAACAGGGTGCCGCCGTGTTTCCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATGTCACTCCAATCACACGATCGATAATCGCCACGTCTTTGTCCACGGAACCCGt                                                           <  2:288257/229‑1 (MQ=255)                                                                                                                                                                                           agcagcTGGCTGCTCGGCCCATCGCAGTTTTCCGGCAGGTCATAAAAGGGCTGTAACGCCA‑TTTTACTGTCCCCCCGGTGGAACACGGGGACGCCAGGTTTACCGCCGAGCAACCCGCCACCGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAAt                                                        <  2:333981/251‑1 (MQ=255)                                                                                                                                                                                              cgctagctgacccgCACACCGCAGTTTTCCAGCAGGTCATAATTAGTCTGTAACGCAC‑TTTTACTGTCCCCCATGGGGAACCAGCTGACGCCAGGTTCTCCGCCGAAGAAACCGCCATGGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGTTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGACCCCGTAATGGt                                                     <  2:1074/238‑1 (MQ=255)                                                                                                                                                                                                cTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGttc                                                   >  2:8107/1‑249 (MQ=255)                                                                                                                                                                                                             ccaccGCAGTTTTCCAGCAGGTCTTAAATAGGCTGTAACGCCA‑TTTTACGGCCGCCCATGTGGAACACTGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATg                                      <  1:251305/247‑1 (MQ=255)                                                                                                                                                                                                                  gCAGTTTTCAGGAAGGTCATAAATATGCCTTAACCCCT‑TTTTCCTGCCCCCCTTGTGGAACACGGTGACGCCGGGTTTACCGCCGAGCAAACCGCCACCGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAg                                 <  2:160610/251‑1 (MQ=255)                                                                                                                                                                                                                   catttttCCAGCAGTTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACCCGGTGACGCCAGGTTTCCCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAgc                                <  1:347635/248‑1 (MQ=255)                                                                                                                                                                                                                    aGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAgcg                               >  2:120670/1‑251 (MQ=255)                                                                                                                                                                                                                     gTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAgcgc                              >  2:42022/1‑251 (MQ=255)                                                                                                                                                                                                                     gTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCA‑TTTTACTCTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGCTCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGACGGGACATGATAAgcgc                              >  1:186076/1‑251 (MQ=255)                                                                                                                                                                                                                         cccAGCAGGTCATAAAAAGGCTTTACCGCCA‑TTTTCCTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCGGCCGCCATATTGGCCAGGTGATGAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTGTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCgg                           <  1:447615/249‑1 (MQ=255)                                                                                                                                                                                                                           cagcggGTCAGAAAGAGGCTGTAACGCCA‑TTTTGCTGTCGCCCATGTGGAACAGGGGGACGCCGGGTTTACCGCCGAACAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCa                        <  1:392749/251‑1 (MQ=255)                                                                                                                                                                                                                                   cATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGa                                     >  2:88024/1‑230 (MQ=255)                                                                                                                                                                                                                                   cATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGa                                     <  1:88024/230‑1 (MQ=255)                                                                                                                                                                                                                                   cATAAATAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAg                <  1:127830/251‑1 (MQ=255)                                                                                                                                                                                                                                     aaaaTAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGAAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGGCATGATAAGCGCCGGTCAGCATCCAGgc              <  1:120670/250‑1 (MQ=255)                                                                                                                                                                                                                                     aaaaTAGGCTGTAACGCCA‑TTTTACTGCCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGgc              <  1:437230/250‑1 (MQ=255)                                                                                                                                                                                                                                      aaaTAGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGgcg             >  1:426016/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaaGGCTGTAACGCCA‑TTTTACTGTCGCCAATGTGGACCACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCt           <  1:4548/249‑1 (MQ=255)                                                                                                                                                                                                                                          aGGCTGTAACGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTg          >  1:401506/1‑250 (MQ=255)                                                                                                                                                                                                                                                 aaCGCCA‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCcttc  >  1:189450/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      ACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTG‑ATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTT‑AACGTGGGTCGGCAGCAGCTTGCTGATCGGCACATCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCT‑TTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTC  >  NC_000913/4558663‑4559151

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