Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	557,635	A→T	L36Q (CTG→CAG)	folD ←	bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	557,635	0	A	T	100.0%	303.2 / NA	82	L36Q (CTG→CAG)	folD	bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (28/54);  total (28/54)

CCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTC  >  NC_000913/557393‑557873                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                ccGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGt                                                                                                                                                                                                                                         >  1:128460/1‑251 (MQ=255)                 cccttgtcgtcgttattcatcccccGCCGGTAACGGCAGTTGAACCAGAATGCCCTCGATGGTGGTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTg                                                                                                                                                                                                                         <  2:113489/229‑1 (MQ=255)                   actgtgtgcttttaccactccccGCCGGTAACGGCAGTGGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGAGAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCa                                                                                                                                                                                                                       <  2:566032/231‑1 (MQ=255)                        tgtaccgttcaagcaCCCGCCGGTAACGGCCGTTGAACCATAAGGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCCCCGCTTCGCTGGTGGTTTCCGGGATGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGGCATAAATTTGCGATGCAGGGTTACTACCCTCCAGCACAACGGACTGTCCTGGTGCCCGCAGTCCTGCTGCAATAc                                                                                                                                                                                                                   <  2:165239/236‑1 (MQ=255)                        tcaaatttctcactccccGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATAc                                                                                                                                                                                                                   <  1:156255/235‑1 (MQ=255)                                    aTACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACttt                                                                                                                                                                                                      >  1:463799/1‑251 (MQ=255)                                     cccccGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACtttt                                                                                                                                                                                                     <  1:39191/249‑1 (MQ=255)                                              gTAACGGCTGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATACGCTCCAGCAGCTCCGCGTCGCTGGTGGTTTCCGTGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTTCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCCCAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGATCGTTTTGAGCAACt                                                                                                                                                                                            <  2:560968/251‑1 (MQ=255)                                              gTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACt                                                                                                                                                                                            <  1:261311/251‑1 (MQ=255)                                                aaCGGCAGTTGATCCAGAATGCCAGCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTCCCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACtt                                                                                                                                                                                           <  2:413619/250‑1 (MQ=255)                                                 aCGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTc                                                                                                                                                                                          >  1:332051/1‑250 (MQ=255)                                                       tttGAACCAGAATGCCATCGATGGGGTTGTCGGCATTCAGCGTATCGCTAAGCTACAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGCCTTCTTCACAAGCCGTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTAATACACACCAGAACAACAGCCTGACCGGGGGCCCGCAGTCCGACTGCAATACGCGCCTGCACTTTTTGAGCAACTTCACAgcg‑c                                                                                                                                                                                  >  2:305074/2‑251 (MQ=255)                                                          gAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CAcc                                                                                                                                                                               <  1:256142/251‑1 (MQ=255)                                                           aaCCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACct                                                                                                                                                                              <  1:254624/251‑1 (MQ=255)                                                              cAGAATGCCATCGATGTTGTTGTCGGCTTTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACctgct                                                                                                                                                                           <  2:40897/251‑1 (MQ=255)                                                                gAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCtgctgc                                                                                                                                                                         <  1:47640/251‑1 (MQ=255)                                                                gAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCtgctgc                                                                                                                                                                         <  1:5660/251‑1 (MQ=255)                                                                    gCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCaa                                                                                                                                                                     <  1:136342/251‑1 (MQ=255)                                                                    gCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCaa                                                                                                                                                                     <  1:581345/251‑1 (MQ=255)                                                                    gCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCa                                                                                                                                                                      <  1:198586/250‑1 (MQ=255)                                                                          gATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGttt                                                                                                                                                               <  1:416118/251‑1 (MQ=255)                                                                           agggggttgtCGGCTTTCCGCGTGTCTACAAGCGACAGCCGCTCCGCTTGGTGGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCCGCTGCGCAATCGtttt                                                                                                                                                              <  2:309242/246‑1 (MQ=255)                                                                           aTGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGtttt                                                                                                                                                              >  1:90149/1‑251 (MQ=255)                                                                           aTGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGtttt                                                                                                                                                              <  2:90149/251‑1 (MQ=255)                                                                                  tgtCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCCATCGTTTTACCGTc                                                                                                                                                        >  2:525129/1‑250 (MQ=255)                                                                                     cGGCATTCTGCGTATCGATAAGCCCCAGCAGCTCCTCTTCGCTGGTGGTTGCCGGGAGGTAATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCaata                                                                                                                                                    <  2:332029/251‑1 (MQ=255)                                                                                     cGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCaat                                                                                                                                                     <  1:543360/250‑1 (MQ=255)                                                                                     cGGCATTCAGCGTATCGATAAGCGCCAGCAGCTGCGCTTCGCGGGTGGTTTCCGGGCGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGGGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCaata                                                                                                                                                    <  2:400537/251‑1 (MQ=255)                                                                                     cGGCATGCAGCGTCTCGATAAGCCTCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGTGAGACGATCCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCCCCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCaat                                                                                                                                                     <  2:513091/250‑1 (MQ=255)                                                                                      ggCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCACTCGTTTTACCGTCaata                                                                                                                                                    >  2:358646/1‑250 (MQ=255)                                                                                      ggCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCaata                                                                                                                                                    >  1:237640/1‑250 (MQ=255)                                                                                         aTTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCt                                                                                                                                                >  1:519849/1‑251 (MQ=255)                                                                                         aTTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCt                                                                                                                                                >  1:399518/1‑251 (MQ=255)                                                                                            cAGCGTATCGACAAGCACCAGCAGCCCCGCTTCGCTGGTGGTTTCCGGGAGGTCATCAGAGCGGGAGACGACCCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAACTTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTACCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAAATTTTTGAGCAACTTCAGAGCGACACCCNC‑GCGCAATCGTTTTACAGTCactatactttg                                                                                                                                             >  2:291964/1‑241 (MQ=255)                                                                                             aGCGTATCGATAAGATCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCCGCAATACGCGCCTGAACTTTTTGAGCACCTTCCGAGCG‑CACCTGCTGCGCACTCGTTTTACCGTCAACAATATTtgc                                                                                                                                            >  2:440352/1‑251 (MQ=255)                                                                                              gCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTtgct                                                                                                                                           <  1:463909/251‑1 (MQ=255)                                                                                              gCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTtgcc                                                                                                                                           >  2:392298/1‑250 (MQ=255)                                                                                                 tATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATCAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGACCTGGTACCCGCAGTCCGGCTGCAATACCCGCCTGAACTTTTTGAGCCAATTCAGAGCG‑CAGCAGCTGCGCAATCGTTTTACCGTCCATAATCTTTGCTGcc                                                                                                                                        >  2:481881/1‑251 (MQ=255)                                                                                                 tATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGcc                                                                                                                                        >  1:483497/1‑251 (MQ=255)                                                                                                  aTCGATAATCGCCAGCCGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATATGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGTCATAAATGTGCGATGCAGCGTTACTACCCACCAGCACAACGTCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGACCTTTTTGAGCATCTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGCCAATAATCTCTGCTGcc                                                                                                                                        <  2:540529/250‑1 (MQ=255)                                                                                                    cGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCAt                                                                                                                                      <  1:3203/250‑1 (MQ=255)                                                                                                        aaGCTCCAGCAGCTCCGCTTCGCTGGTGGTNTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCaga                                                                                                                                  <  2:282574/250‑1 (MQ=255)                                                                                                         aGCTCCAGCAGCGCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGGGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAAGTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CTCCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCagaga                                                                                                                                <  2:178554/251‑1 (MQ=255)                                                                                                               agcagcTCCGCTGCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGGCAATAATCTTTGCTGCCATCAGAGAGAGGAt                                                                                                                          <  1:425027/251‑1 (MQ=255)                                                                                                                   gcTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGCGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTcc                                                                                                                       >  2:565380/1‑250 (MQ=255)                                                                                                                       cGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTg                                                                                                                  <  1:448036/251‑1 (MQ=255)                                                                                                                        gCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGt                                                                                                                 >  1:98640/1‑251 (MQ=255)                                                                                                                        gCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAAACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGGGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGt                                                                                                                 <  1:473931/251‑1 (MQ=255)                                                                                                                           tCGCTTGTGGTTTCCGGTAGGTCATAAGAGCTGGAGTCGAACCCGAGTTCTTCGCAAGCTTTGCGTTGTCTTGCGACATAAATTTGCGATGCAGGGTTACTCCCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCTGGCTGCAATACGCTCCTGAACTTTTTGATCAACGTCATAGCG‑CACCTGCTGCGCAATCGTTTTAGCGTCTATAATCTGTTCTGCCATCAGAGAGAGGATTGCATCTGtt                                                                                                                <  2:329509/249‑1 (MQ=255)                                                                                                                             gCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGTCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGt                                                                                                            <  1:293614/251‑1 (MQ=255)                                                                                                                              cTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGt                                                                                                            >  2:39812/1‑250 (MQ=255)                                                                                                                                   ggtTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCCGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGt                                                                                                            >  2:559274/1‑245 (MQ=255)                                                                                                                                   ggtTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATc                                                                                                       <  1:389351/250‑1 (MQ=255)                                                                                                                                        ccGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCCCAAGCCTTGCGTTTTCATGCAACATAAATTTGCGATGCAGGGTTGCTACCCACCAGCACACCGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGTACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAgg                                                                                                  <  1:436390/250‑1 (MQ=255)                                                                                                                                        ccGGGAGGTCATAAGAGCGGGAGACGAACCCGACTGCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAggg                                                                                                 <  2:519849/251‑1 (MQ=255)                                                                                                                                          gggAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGa                                                                                               <  1:358646/251‑1 (MQ=255)                                                                                                                                             aGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATgc                                                                                            <  1:572139/251‑1 (MQ=255)                                                                                                                                                     agagCGGGAGACGAACCCGACTTCTTCACAAGCCGTGCGTTTGCGTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGTCTGTCCTGGTGCCCGCAGTCCGGCTGCAATTCGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAttt                                                                                    <  2:58088/251‑1 (MQ=255)                                                                                                                                                            gagaCGAACACGACTTCTTCACAAGCCTTGCGTTTGCTTGAGACATAAATTAGCGATGCAGGGTTACTACCCACCAGCTCAACGGCCTGTCCTGGTGCCCGCAGTCCGCCTGCAATACGCGCCTGACCTCTTTGTGCAACTTCAGAGCG‑CACCAGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGCGCGAGGATTCCATCTGTTACGTAGCTCGCAGGGGATGCGCCTAATTTGTc                                                                                >  2:309538/1‑248 (MQ=255)                                                                                                                                                               aCGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGc                                                                          <  1:55739/251‑1 (MQ=255)                                                                                                                                                                 gAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCgg                                                                        <  1:123196/251‑1 (MQ=255)                                                                                                                                                                               cacaAGCCTTGCGTTTGCCTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCTGGGCGCGCTGTCAgg                                                          <  2:332051/251‑1 (MQ=255)                                                                                                                                                                                acaAGCCTTGCGTTTGCGTGCGACATAAATTTGCGATGCAGGGATACTACCCACCAGCACAACGGCCTGTCCTGGTCCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGAG‑CACCTGCTGCGCAATCGTTTTACCGTCACTAATCTTTGCTGCCATCAGAGAGCAGATTCAATCTGTTACGGAAATCGCAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGGTGTGAggg                                                         >  2:415633/1‑250 (MQ=255)                                                                                                                                                                                  aaGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgttt                                                       >  1:199804/1‑251 (MQ=255)                                                                                                                                                                                   agacgtGCGGTTGCTTGCGACATAAATTTGGGATGCAGGGTTTCTACCCACCAGCACACCGGCCTGGCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttc                                                      <  2:39068/246‑1 (MQ=255)                                                                                                                                                                                        tgcggtggcttgcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttcgttt                                                  <  1:481831/243‑1 (MQ=255)                                                                                                                                                                                        tGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttcgtttc                                                 >  1:1798/1‑251 (MQ=255)                                                                                                                                                                                         gcggttgcttgcGACATAAATTTGCGATGCAGGGTTCCTACCCCCCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCATAAGCGGGGCGCGCTGTCAGgtttcgttgca                                                <  2:363328/247‑4 (MQ=255)                                                                                                                                                                                            gttgattgcGGCATAACTTTGCGAGGCAGGGTTCCTCCACACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGTATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGAt                                             <  2:1798/246‑1 (MQ=255)                                                                                                                                                                                              tgcttgcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttcgtt                                                   >  2:463800/1‑243 (MQ=255)                                                                                                                                                                                              tgcttgcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttcgtt                                                   <  1:463800/243‑1 (MQ=255)                                                                                                                                                                                               gcttgcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGAttt                                           <  1:434972/250‑1 (MQ=255)                                                                                                                                                                                                   gcGACATAAATTTGCGATGAAGGGGTACTACCCACCAGCCCACCGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATTTGTTACGTAGAGCGATGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCTATTCAGATTTATcg                                       <  2:555659/250‑1 (MQ=255)                                                                                                                                                                                                        aTAAATTTTCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAg                                                                           <  2:448170/209‑1 (MQ=255)                                                                                                                                                                                                        aTAAATTTTCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAg                                                                           >  1:448170/1‑209 (MQ=255)                                                                                                                                                                                                          aaaTTTGCGAGGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CCCCTGCTGCGCAATCGTTTTACCGGCAATAATCTTTGCTGCCATCAGAGAGAGGCTTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGc                               <  2:455043/251‑1 (MQ=255)                                                                                                                                                                                                                  gATGCAGGGTTACTACACACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGCAGCGACCTCtt                       >  2:330133/1‑251 (MQ=255)                                                                                                                                                                                                                         ggTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAgg                >  1:249582/1‑251 (MQ=255)                                                                                                                                                                                                                             actacCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAg            <  1:465183/251‑1 (MQ=255)                                                                                                                                                                                                                                tcccACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGAGTCCATCTGTTACGTAGATCAAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAgg           <  2:406728/248‑1 (MQ=255)                                                                                                                                                                                                                                acCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCATGAACTTTTTGAGCACCTTCAGAGCG‑CACCTGCTGCGCTATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGCTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTATTGGGAAGGTgtagcg         >  2:345225/1‑246 (MQ=255)                                                                                                                                                                                                                                       aGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTc  <  1:574534/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                CCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCG‑CACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACCGTC  >  NC_000913/557393‑557873

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