Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	389,601	C→T	D43N (GAC→AAC)	hemB ←	porphobilinogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	389,601	0	C	T	98.7%	300.1 / NA	79	D43N (GAC→AAC)	hemB	porphobilinogen synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (32/46);  minor base G (0/1);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTATTCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TCGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACAT  >  NC_000913/389354‑389843                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     gtgtATTCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAg‑ttgtc                                                                                                                                                                                                                                                   >  2:212070/1‑251 (MQ=255)     aTTCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAAt                                                                                                                                                                                                                                                <  1:416671/250‑1 (MQ=255)                cAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTcaac                                                                                                                                                                                                                                   >  1:406449/1‑251 (MQ=255)                  ggtgccggACATAACGATCATTTCTGGCACGGTCGTCGTGCAGATGCGCGACATACGCGCCGCCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCAGCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCGCC‑TTcaacaa                                                                                                                                                                                                                                 <  2:114135/244‑1 (MQ=255)                     gtCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTGG‑TTGTCAATTTCTTC‑TTcaacaa                                                                                                                                                                                                                                 >  2:256327/1‑248 (MQ=255)                          aCATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCCCGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGc                                                                                                                                                                                                                         >  2:242403/1‑251 (MQ=255)                          aCATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGc                                                                                                                                                                                                                         <  1:392142/251‑1 (MQ=255)                           cATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTACATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGGATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCCTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGc                                                                                                                                                                                                                         >  2:194903/1‑250 (MQ=255)                            aTAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCa                                                                                                                                                                                                                        >  2:393776/1‑250 (MQ=255)                               aCGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAacac                                                                                                                                                                                                                    <  1:81823/251‑1 (MQ=255)                                 gctCATTTCTGGCACGGTCTCCTTGCAGATCCGCTACATCCGCGCCCACAGTCCATCTTCCCGCAAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCAAGGAACGAATACCGGCGTTGGCGATGCGTTCCATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGACATGGCTTAAACGTCTTTGTAG‑ATGTCAATTTCTTC‑TTCAACAAAGCTCGTCAACAcc                                                                                                                                                                                                                   >  2:390026/3‑250 (MQ=255)                                    cATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGt                                                                                                                                                                                                               >  2:58287/1‑251 (MQ=255)                                    cATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGt                                                                                                                                                                                                               >  1:99880/1‑251 (MQ=255)                                     aTTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTc                                                                                                                                                                                                              <  2:301133/251‑1 (MQ=255)                                         cTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑tt                                                                                                                                                                                                          >  2:12791/1‑250 (MQ=255)                                         cTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑tt                                                                                                                                                                                                          >  2:304050/1‑250 (MQ=255)                                          tGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTa                                                                                                                                                                                                         >  1:166598/1‑250 (MQ=255)                                           ggCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTaa                                                                                                                                                                                                        <  2:154695/250‑1 (MQ=255)                                               cGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTa                                                                                                                                                                                                   <  1:291670/251‑1 (MQ=255)                                                 gTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCAGCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGTTCG‑TTAAGGCTAAg                                                                                                                                                                                                 <  2:174237/251‑1 (MQ=255)                                                    tgcgtGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAgtgt                                                                                                                                                                                              >  2:230693/5‑251 (MQ=255)                                                     gCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAgt                                                                                                                                                                                                <  2:307711/248‑1 (MQ=255)                                                        tGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGtgttgt                                                                                                                                                                                           <  1:304050/250‑1 (MQ=255)                                                         ccAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGTCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGtct                                                                                                                                                                                         <  1:152589/250‑1 (MQ=255)                                                           aGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTctct                                                                                                                                                                                       >  1:404076/1‑251 (MQ=255)                                                             aTGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTc                                                                                                                                                                                     <  1:104166/251‑1 (MQ=255)                                                              tGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTc                                                                                                                                                                                     >  2:21810/1‑250 (MQ=255)                                                               gcgcGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCa                                                                                                                                                                                    <  1:389704/250‑1 (MQ=255)                                                                cgcgACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCaaa                                                                                                                                                                                  <  2:281014/251‑1 (MQ=255)                                                                   gACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACAt                                                                                                                                                                               >  2:77340/1‑251 (MQ=255)                                                                   gACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTA‑Gtgttgttttcttcaacaa                                                                                                                                                                               >  1:249866/1‑238 (MQ=255)                                                                      aTACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAgc                                                                                                                                                                            >  1:389038/1‑251 (MQ=255)                                                                      aTACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAgc                                                                                                                                                                            <  1:88553/251‑1 (MQ=255)                                                                      aTACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAgc                                                                                                                                                                            <  1:254614/251‑1 (MQ=255)                                                                                  aGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAgg                                                                                                                                                                <  1:12791/251‑1 (MQ=255)                                                                                           tCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCTCAGCGCAGGAGATTTgc                                                                                                                                                        <  1:80881/250‑1 (MQ=255)                                                                                           tCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGAtt                                                                                                                                                           <  1:393776/247‑1 (MQ=255)                                                                                                ccAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGc                                                                                                                                                  <  2:3278/251‑1 (MQ=255)                                                                                                         gCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGc                                                                                                                                         <  1:432822/251‑1 (MQ=255)                                                                                                                  ttCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGAtt                                                                                                                                 >  1:55963/1‑250 (MQ=255)                                                                                                                  ttCATCGGTATGGTGAGAGATGCAAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGAtt                                                                                                                                 <  1:284515/250‑1 (MQ=255)                                                                                                                       cGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCt                                                                                                                           <  1:114796/251‑1 (MQ=255)                                                                                                                                   agaTGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGc                                                                                                               >  2:166672/1‑251 (MQ=255)                                                                                                                                   agaTGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCACAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGc                                                                                                               >  2:48197/1‑251 (MQ=255)                                                                                                                                       gCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGa                                                                                                           <  2:275525/251‑1 (MQ=255)                                                                                                                                        ccAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGAt                                                                                                          <  1:256327/251‑1 (MQ=255)                                                                                                                                          aaaaGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCGGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGt                                                                                                        <  1:230693/251‑1 (MQ=255)                                                                                                                                              gTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTtgt                                                                                                    >  2:132510/1‑251 (MQ=255)                                                                                                                                                catcaCGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTgg                                                                                                  <  1:132510/251‑1 (MQ=255)                                                                                                                                                          cGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAgggg                                                                                        >  1:133633/1‑251 (MQ=255)                                                                                                                                                          cGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAgggg                                                                                        <  1:331884/251‑1 (MQ=255)                                                                                                                                                            aaTACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGc                                                                                       >  1:262352/1‑250 (MQ=255)                                                                                                                                                             aTACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCAt                                                                                     <  1:93136/251‑1 (MQ=255)                                                                                                                                                             aTACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCAt                                                                                     <  2:153082/251‑1 (MQ=255)                                                                                                                                                             aTACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCAt                                                                                     >  2:385827/1‑251 (MQ=255)                                                                                                                                                              tACCGGCGTTGTCGATGCGTTCAATTTCACGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑GTCAACAAAGATCGGCCACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACAAAGCGCGCAGCGCAGGAGATTTTCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGACTCGAGTGGcgt                                                                                     >  1:321574/1‑248 (MQ=255)                                                                                                                                                               aCCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCGGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATa                                                                                    <  2:43748/250‑1 (MQ=255)                                                                                                                                                                       tGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTg                                                                           >  1:271167/1‑251 (MQ=255)                                                                                                                                                                              gCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAg                                                                     <  1:242403/250‑1 (MQ=255)                                                                                                                                                                               cGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCCCGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAgg                                                                    <  1:180899/250‑1 (MQ=255)                                                                                                                                                                                    aagttCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAggg                                                              <  2:99880/248‑1 (MQ=255)                                                                                                                                                                                    aaTTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAggg                                                              >  1:259356/1‑251 (MQ=255)                                                                                                                                                                                          gcgTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGttt                                                        <  1:21810/251‑1 (MQ=255)                                                                                                                                                                                          gcgTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGttt                                                        <  1:218989/251‑1 (MQ=255)                                                                                                                                                                                                   aTGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCCGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGCAGCAGGGTAGGGATGTTTTACGAAAGt                                               >  2:123891/1‑251 (MQ=255)                                                                                                                                                                                                     ggttCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTg                                             <  2:406449/249‑1 (MQ=255)                                                                                                                                                                                                     gTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTg                                             >  1:40832/1‑251 (MQ=255)                                                                                                                                                                                                        tctcTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTc                                           <  2:133633/250‑1 (MQ=255)                                                                                                                                                                                                             ggAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTcgcg                                        <  1:212070/248‑1 (MQ=255)                                                                                                                                                                                                              gAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGt                                    <  1:148749/251‑1 (MQ=255)                                                                                                                                                                                                                  gcgcATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGAc                                <  1:403454/251‑1 (MQ=255)                                                                                                                                                                                                                   cgcATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGAc                                >  2:245723/1‑250 (MQ=255)                                                                                                                                                                                                                    gcATCACGCCTGGCATGGCTTCAACCGGTTTGGAGGTGGT‑AATTTTTTCCTCCAACAA‑GAACCGCAAACCCAGGGCGGTTAGGGC‑AAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACaa                              <  2:249866/251‑1 (MQ=255)                                                                                                                                                                                                                      aTCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACaaga                            <  2:426549/251‑1 (MQ=255)                                                                                                                                                                                                                        cACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACaagaa                           <  1:306582/250‑1 (MQ=255)                                                                                                                                                                                                                        cACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACaagaa                           <  1:58287/250‑1 (MQ=255)                                                                                                                                                                                                                                          aaCGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGcc         >  1:188695/1‑250 (MQ=255)                                                                                                                                                                                                                                           aCGGCTTTGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAg       <  1:422448/251‑1 (MQ=255)                                                                                                                                                                                                                                                tttGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACAt  >  1:292495/1‑251 (MQ=255)                                                                                                                                                                                                                                                 ttGTAG‑TTGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACa   <  1:106709/249‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     GTGTATTCACAGAAGCAGGTGTCTGACATAACGATCATTTCTGGCACGGTCTGCTTGCAGATGCGCGACATACGCGCCACCAGTCCATCTTCCCGCCAGGCATCGCTGCCGGTTTCATCGGTATGGTGAGAGATGCCAAAAGTCATCACGGAACGAATACCGGCGTTGGCGATGCGTTCAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAG‑TCGTCAATTTCTTC‑TTCAACAAAGATCGGCAACACCAGGTCG‑TTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACAT  >  NC_000913/389354‑389843

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