Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,298,401	G→A	L322L (CTC→CTT)	fdhF ←	formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,298,401	0	G	A	94.9%	224.4 / 2.1	79	L322L (CTC→CTT)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (34/41);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

GCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGCGACGAACGCTTT  >  NC_000913/4298161‑4298643                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   gCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTAACCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTGTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCgg                                                                                                                                                                                                                                          >  1:191741/1‑251 (MQ=255)   ccATAATGTACGCGGCACGCATTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTGTCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTc                                                                                                                                                                                                                                        <  1:120961/251‑1 (MQ=255)   ccATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTc                                                                                                                                                                                                                                        >  1:424484/1‑251 (MQ=255)        aTGTACGCGGCACGCACTGCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCAt                                                                                                                                                                                                                                   <  2:141137/251‑1 (MQ=255)        aTGTACGCAGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCACATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGATTACCAAGGTTACCGGTCAGCAt                                                                                                                                                                                                                                   >  1:333981/1‑251 (MQ=255)          gTACGCTGCACGCACGTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGTATCTTTCACGTACTGGTAGCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGGCAGCATc                                                                                                                                                                                                                                  <  2:485833/250‑1 (MQ=255)               cgGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAgg                                                                                                                                                                                                                            <  1:118040/251‑1 (MQ=255)                        cTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCCCAGGCCCCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGCGGCTGGGCagg                                                                                                                                                                                                                   >  2:163119/1‑250 (MQ=255)                               atgcgtggctcgGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCTGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGGTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCa                                                                                                                                                                                                            <  2:173680/241‑1 (MQ=255)                                    cggcgcGGTGCGGCAGCTCGCTGAGGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGAGACGTATCCGGCAGCGCGCCCATAGCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGtt                                                                                                                                                                                                       <  2:424484/249‑1 (MQ=255)                                      gcgcGGTGCGGCAGCGCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCGTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGtt                                                                                                                                                                                                       <  1:466437/249‑1 (MQ=255)                                       cgcgGTGCTGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTcc                                                                                                                                                                                                    <  1:425507/251‑1 (MQ=255)                                       cgcgGGGCGGCTGCTCGCTGATGCGAGAGCCTGTATGCTCTGGGAGGCTGTCCACGCCCCAGGCTTTGGCGAATGACTCGCGGTTAGCCGGATCTTACACGTACGGGTATCCCGGATACGTATCCGGCAGCGCGCCCATAGCGCAGGCACCCTGAACGTGGTGCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGGTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTcc                                                                                                                                                                                                    <  2:152639/251‑1 (MQ=255)                                          ggTGCTGCAGCTCGCTGATGCGATGGCCGGTATGCGCTGGCAGGCTTACCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCGCCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGTCTTACCAAGGTTACCGGTCAGCATCGCTAGGCTGGTCAGAGAACGCACGGTTTCCACg                                                                                                                                                                                                 <  1:357243/251‑1 (MQ=255)                                           gTGCAGCAGCTCGATGATGCGATAGCCGGCATGCGCTAGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGCTCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGATCGTTGTTCTGACCCCGACCCGGGTTAACACCCGCAAGCGGCTTACCACGGTTACCGGTCAGCATCGCGAGCATCGTCAGAGAACGCCCGCTTTCCCCGc                                                                                                                                                                                                >  2:242265/1‑251 (MQ=255)                                              ggggcgcccgcggtttcGATGGCCGGGAGGCGCGGGCAGTCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCt                                                                                                                                                                                             <  2:313644/235‑1 (MQ=255)                                                gcaccccGCTGATGGGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCGCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTgg                                                                                                                                                                                           <  2:381388/245‑1 (MQ=255)                                                 ccgCTCGCTGATGCGATAGCCGGTATGCGATGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACACGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGCGACTGGTCAGAGAACGCCCGGTTCCCCCGCCCTGGt                                                                                                                                                                                          >  2:200306/3‑251 (MQ=255)                                                   gCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGt                                                                                                                                                                                          >  1:284466/1‑249 (MQ=255)                                                             gCGATAGCCGGTATGCGCTGGCAGTCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTa                                                                                                                                                                              <  1:135704/251‑1 (MQ=255)                                                              cGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGGTCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCCCGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCCCGGTTTCCACGCCCTGGTAGTACTGGGttcc                                                                                                                                                                             >  2:217486/1‑249 (MQ=255)                                                                   gccgacaTGCTCTGGGAGGCTTTCCACGCCCCATGCCTTGGCTAATTTCTCGCGGTGAGCCGGATCTTTCACGTACTGGTATCCCGGAGACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCGTACCAAGGTTACCGGGCAGCATCTCGAGGCTGGTCAGAGAGCGCACGGTTGCCACGCCATGGTAGAATTGGATTACTCCCa                                                                                                                                                                        <  2:464717/245‑1 (MQ=255)                                                                         aTGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGcccc                                                                                                                                                                  >  1:358110/1‑251 (MQ=255)                                                                                ggCAGGGTTTCCACGCACCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGAt                                                                                                                                                           <  1:71977/251‑1 (MQ=255)                                                                                ggCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAGCCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGAt                                                                                                                                                           <  1:474596/251‑1 (MQ=255)                                                                                 gCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATg                                                                                                                                                          >  1:198020/1‑251 (MQ=255)                                                                                     gCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATggcgg                                                                                                                                                      <  1:374351/251‑1 (MQ=255)                                                                                               ccccAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGCTCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGc                                                                                                                                            >  1:183653/1‑251 (MQ=255)                                                                                                 ccttgCTTTTGCGAATTTCGCGCGGTTAGGCGGAGCCTTCACGTCCCGGCATCCCGGAGACGTAACAGGCAGCGCGCCCAGAGCGCATGCACCCTGAACGGTGTTCTGACCACGAACCTGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGATAACGCACGGTTTCCTCGCCCTGGTAGAATTGTGTTACACCCATGCCCCACAGGATTGCGGCGCTTTTCGCCt                                                                                                                                          <  2:12394/247‑1 (MQ=255)                                                                                                 ccAGGCTTTGGCGAATTTCTCGCGGGTAGCCGGATCTCTCACGTACTGTGATACCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCCCCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGACCTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTGTCGCCt                                                                                                                                          <  2:129713/251‑1 (MQ=255)                                                                                                   aGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTgg                                                                                                                                        >  1:66614/1‑251 (MQ=255)                                                                                                   aGGCGTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCATCGCGCCCATATCGCAGGCACCCTGAACGTTGGTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTgg                                                                                                                                        <  1:436225/251‑1 (MQ=255)                                                                                                          ggCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATAc                                                                                                                                  >  1:296420/1‑250 (MQ=255)                                                                                                              aaTTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCg                                                                                                                             <  1:275439/251‑1 (MQ=255)                                                                                                              aaTTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCg                                                                                                                             >  1:435043/1‑251 (MQ=255)                                                                                                               aTTTCTCGCGGTTAGCCGGCTCTTTCACGGACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCCCGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGt                                                                                                                            <  2:254833/251‑1 (MQ=255)                                                                                                                   ctcGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACt                                                                                                                                                                                    >  1:352547/1‑191 (MQ=255)                                                                                                                   ctcGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACt                                                                                                                                                                                    <  2:352547/191‑1 (MQ=255)                                                                                                                   ctcGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTgccc                                                                                                                        >  2:164771/1‑250 (MQ=255)                                                                                                                       cgGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTg                                                                                                                    >  1:270051/1‑251 (MQ=255)                                                                                                                       cgGTTAGCCGGATATTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTg                                                                                                                    <  1:267920/251‑1 (MQ=255)                                                                                                                           tAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACg                                                                                                                 >  1:156397/1‑250 (MQ=255)                                                                                                                              ccGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAAtc                                                                                                             <  1:348476/251‑1 (MQ=255)                                                                                                                                gtaTCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAAtctc                                                                                                           >  1:258760/3‑251 (MQ=255)                                                                                                                                 gATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTAGACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAAtctc                                                                                                           >  1:59064/1‑250 (MQ=255)                                                                                                                                   tCTTTCACGTACTGGTATCC‑GGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTg                                                                                                       >  1:69887/1‑251 (MQ=255)                                                                                                                                      ttCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAgaga                                                                                                                                                                                                                 >  2:15038/1‑143 (MQ=255)                                                                                                                                      ttCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAgaga                                                                                                                                                                                                                 <  1:15038/143‑1 (MQ=255)                                                                                                                                              cTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGCCGACTCTCACTGGCGCTGACGc                                                                                             >  2:34809/1‑251 (MQ=255)                                                                                                                                                 gTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCgg                                                                                          <  1:152969/251‑1 (MQ=255)                                                                                                                                                    tcccgacccgtcgcccgtccgatctcccATAGCGCAGGCACCCCGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTGCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCTGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTg                                                                                        <  2:199837/225‑1 (MQ=255)                                                                                                                                                    tCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTAACCACGACCCGCGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACAATTTCCACGCCCTGGAAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCTTCCAACCGTGGCGCCTGACGAATCTCACTGGCGCTGCCGCCGGTgt                                                                                       >  2:21194/1‑250 (MQ=255)                                                                                                                                                              gTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAAc                                                                             >  1:343165/1‑251 (MQ=255)                                                                                                                                                                aTCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAAcc                                                                            <  2:435043/250‑1 (MQ=255)                                                                                                                                                                  ccGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGAc                                                                         >  1:236856/1‑251 (MQ=255)                                                                                                                                                                   cGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACt                                                                        <  1:8645/251‑1 (MQ=255)                                                                                                                                                                    ggCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGCCTc                                                                       >  1:142187/1‑251 (MQ=255)                                                                                                                                                                     gCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTcc                                                                      <  1:335519/251‑1 (MQ=255)                                                                                                                                                                      cAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCg                                                                     <  2:156397/251‑1 (MQ=255)                                                                                                                                                                       aGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAGCGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCgg                                                                    >  1:406730/1‑251 (MQ=255)                                                                                                                                                                             ccttatcgctggcccccttatagtttgttC‑GACCACGAACCGGGTTAACCCCGGCATGCGGCTTCCCAAGGTGACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTTCCCGCCCGGGTAGAACTGGGTTACACCCATGCCCCATAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCCCCTGACGAATCTCACTGGCTCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAg                                                             <  2:191741/227‑1 (MQ=255)                                                                                                                                                                             cccATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCCGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTg                                                                                        >  1:199837/1‑225 (MQ=255)                                                                                                                                                                                          cACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGAtt                                                 <  2:343165/251‑1 (MQ=255)                                                                                                                                                                                          cACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGAtt                                                 <  1:378755/251‑1 (MQ=255)                                                                                                                                                                                                  aCGTTGGTCTGCCCACGAACCGGGTTAGCTCCCGCAGGCGGCTTACCAAGGTTACCGGTCAGCATCGCGATGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGTCGCTGACGCCGGTGAGATCTTCAACCGACGCCGGCGTGTAGCCTTCAACGAGTTTACGATa                                         <  2:321098/251‑1 (MQ=255)                                                                                                                                                                                                   cGTTGTTCTGACCACGAACCGGGTTAACAACCGCATGCGCCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGATTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTTCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGCTTTTACGATAc                                        >  1:132119/1‑251 (MQ=255)                                                                                                                                                                                                       gttCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCGATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATActct                                     <  1:32429/250‑1 (MQ=255)                                                                                                                                                                                                            gACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAg                               <  1:164771/251‑1 (MQ=255)                                                                                                                                                                                                            gACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTAGCCTGGGCATACATCAGTGCCGCCTGACGTATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAAAGATTTTACGATTCTCTTCaa                                 >  2:473758/1‑249 (MQ=255)                                                                                                                                                                                                                      cGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTCTGTAc                     >  1:466002/1‑251 (MQ=255)                                                                                                                                                                                                                                ccctcATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTTACGAATCTCACTGGCTCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTCTGTACTTGAAGcgac           <  1:443072/247‑1 (MQ=255)                                                                                                                                                                                                                                cccGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGcgac           <  1:424309/251‑1 (MQ=255)                                                                                                                                                                                                                                cccGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTGCCACGCCCTGGTAGATCTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTTACGCCGGTGATAGCTTCAACCGACTCCGGCGTGTAGCCTTCAACGAATTTACGATTCTCTTCAAAGCTTTCTGTACGTGAAGcgac           <  2:69638/251‑1 (MQ=255)                                                                                                                                                                                                                                   gCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGcgacga         >  2:77502/1‑250 (MQ=255)                                                                                                                                                                                                                                    caggCGGATTACCAAGGTTCCCGGTCAGCATCGCGAGGGTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTGCTGTACGTGAAGCGACGAAc       <  2:94611/248‑1 (MQ=255)                                                                                                                                                                                                                                    cATGCGGCTTCCCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGTCGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATAATATTCAAAGCCTTCTGTACGTGAAGcgccaaa        >  2:353538/1‑245 (MQ=255)                                                                                                                                                                                                                                       gCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGCGACGAACGCt    <  1:118344/251‑1 (MQ=255)                                                                                                                                                                                                                                         ggCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACCCCCATGCCCCACAGGATGGCGGCGCTGTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGCGACGAACGCttt  <  2:358400/251‑1 (MQ=255)                                                                                                                                                                                                                                         ggCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGACGTGTAGCCTTCAACGCTGTTACGATTCTCTTCAAAGCCTTCTGTACGTGAAGCGACGAACGCttt  >  2:301730/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   GCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAACGATTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGCGACGAACGCTTT  >  NC_000913/4298161‑4298643

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