Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,937	A→G	P68P (CCA→CCG)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,937	0	A	G	98.8%	255.8 / ‑5.3	83	P68P (CCA→CCG)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (40/41);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

GGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTA  >  NC_000913/3696‑4177                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  ggCTGGATACGGCGGGCGCACGATTACTGGAAAACTAAAGGAAACTCTACAATCTGAAAGAGCACAATGAGCAGGGCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCGGTTTTTTCCGCACGACCTGCCGGAATGCAGCCTGACTGAAAGTGATTAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGAGCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATc                                                                                                                                                                                                                                          <  2:568463/251‑1 (MQ=255)           ggcgggcgCACGAGTACTGGAAAACTAAATGATACTCTACAATCTGAAAGATCACAAAGAGCCGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGAACGACCTGCCGGAATTCAGCCTGAATGAAATTGATGAGATGCTGAAGATGGATTTTGTCACCCGCAGGGCGAAGATCCTCTCGGCGTTTATGGGTGATCAAATCCCGCAGGCAATCATGGacaggg                                                                                                                                                                                                                                >  2:452767/1‑246 (MQ=255)           ggcgggcgCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAgc                                                                                                                                                                                                                                >  1:36548/1‑251 (MQ=255)              gggcgCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGcgcg                                                                                                                                                                                                                             <  1:134337/251‑1 (MQ=255)                        gTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGACGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCgcggc                                                                                                                                                                                                                   >  1:469803/1‑251 (MQ=255)                             ggAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGtt                                                                                                                                                                                                                <  1:442421/249‑1 (MQ=255)                                          aaaCTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTcc                                                                                                                                                                                                 >  1:120943/1‑251 (MQ=255)                                          aaaCTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTCTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTcc                                                                                                                                                                                                 <  1:582599/251‑1 (MQ=255)                                          aaaCTATACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGCTCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTcc                                                                                                                                                                                                 >  2:149845/1‑251 (MQ=255)                                                        aaaaGATCACAATGAGCAGGTCAGCGTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAAGCAGGGGCTGTTTTTGCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCGGGAGTTTGTCACCCGCAGTGCGAAGATCCGCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTg                                                                                                                                                                                   <  1:289172/250‑1 (MQ=255)                                                         aaaGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCAGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAGGATCCTCTCGGCGTTTATTGGCGATGAAATCCCCCAGGAAATCCTGGAAGCGCGCGTGCGCGCGGCGTTTGCCTTCCAGGCTCCGGTCGCCAATGTTg                                                                                                                                                                                   >  2:238623/1‑250 (MQ=255)                                                               cacaATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTg                                                                                                                                                                                   <  2:591959/244‑1 (MQ=255)                                                               cacaATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTg                                                                                                                                                                                   >  1:591959/1‑244 (MQ=255)                                                               cacaATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCg                                                                                                                                                                             >  2:286167/1‑250 (MQ=255)                                                                acaATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGAt                                                                                                                                                                           <  1:173170/251‑1 (MQ=255)                                                                         aGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGt                                                                                                                                                                                                                            >  1:380983/1‑193 (MQ=255)                                                                         aGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAAGCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGt                                                                                                                                                                                                                            <  2:380983/193‑1 (MQ=255)                                                                             ccgtgtgctgcAAGCGGTACCCCAGTGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTTACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCATCCGCAGTGCGAAGATCCTTGCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCAAATGTTGAAAGCGATGTCGGTTGTCTgg                                                                                                                                                              <  2:582811/242‑1 (MQ=255)                                                                                cTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAAt                                                                                                                                                           >  1:350878/1‑251 (MQ=255)                                                                                   tgggcAAGCAGTATCCCAGGGGTTGGGCAAAAATCAGGGTCTGTTTTTTCCGCACGCCCTGACGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTGTTTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGGGCGCGCGGCGTTTTCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTg                                                                                                                                                         <  2:365009/247‑1 (MQ=255)                                                                                    gcgcAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGtt                                                                                                                                                       <  1:232839/251‑1 (MQ=255)                                                                                     cgcAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTc                                                                                                                                                      >  1:229175/1‑251 (MQ=255)                                                                                       cAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGaaa                                                                                                                                                                                <  1:312132/223‑1 (MQ=255)                                                                                       cAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGaaa                                                                                                                                                                                >  2:312132/1‑223 (MQ=255)                                                                                               aaCCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGCATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCa                                                                                                                                             >  2:194429/1‑250 (MQ=255)                                                                                                 cccAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAAAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACg                                                                                                                                          <  1:130919/251‑1 (MQ=255)                                                                                                    aGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGGATTGTTCCCCGGGCCAACGCTg                                                                                                                                       >  2:190235/1‑251 (MQ=255)                                                                                                    aGGGGTTGGGCAAAAATCAGGGGCTGGTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGCTGCTGAAGCTGGATTTTGTCACCAGCAGTGCGAAGATCCACTCGGCGTTTATTGGTGATGAAATCCCGCAGGCAATCCTGGAAACACGCGTGCGCGCGGAGTTTGCCTTCCCCGCACCGGTCGCCAATGTTGGAAGCGTTGTCGGTTTTCTGGTATTGTTTCACGGGCCACCGCtt                                                                                                                                       >  2:114863/1‑250 (MQ=255)                                                                                                     ggggTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAAt                                                                                                                                                                                       >  2:327737/1‑202 (MQ=255)                                                                                                     ggggTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAAt                                                                                                                                                                                       <  1:327737/202‑1 (MQ=255)                                                                                                       ggTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCa                                                                                                                                    <  1:35974/251‑1 (MQ=255)                                                                                                                agaagttggggcggtgttgtCCGCACGACCTGCCGGAATTCAGCCTGACTGAAAGGGATGAGACGCTGAGGCGGGATTTTGTCACCCGCAGTGCGAAGAACCGCTCGGCGTCTATTGGTGATGAAATGCCGCTGGCAAGCCTGGAAGAGCGCGCGCGCGCCGCGTTTGCCTTCCCGGCTCGGGTCGCCAATGTTGAAAGCGATGGCTGTTGTCTGGACCTGGTCCACGGGCGACTGCTGCCTTTTAAAGAt                                                                                                                           <  2:133575/232‑1 (MQ=255)                                                                                                                aaaaTAAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGTGATGCGGAAGCTGGATTATGTCACCCGCAGTGCGACGATCCTCTCGGCGTTTCTTGGTGATGAAATCCCGCAGTAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGCAGTAGGTCGTATAGAATTGTTCCCAGGGACCACGCTGGCAtttta                                                                                                                               >  2:360297/1‑245 (MQ=255)                                                                                                                     cAGGGGCTGTTTTTTCCGCCCGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGACAGCGATTTCGGTTGTCTGGAATTGTTCCACGTGCCAAAGCTGGCATTTCAAGAttacgt                                                                                                                      >  2:467670/1‑247 (MQ=255)                                                                                                                          gCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTc                                                                                                                 <  1:185754/251‑1 (MQ=255)                                                                                                                          gCTGTTTGTTCCGCACGTCCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATGTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAATTCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCTTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTc                                                                                                                 <  2:350878/251‑1 (MQ=255)                                                                                                                            tGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGc                                                                                                                                                                                                    >  1:120700/1‑166 (MQ=255)                                                                                                                            tGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATGGATGAGATGCTGAAGCTGGATGTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGc                                                                                                                                                                                                    <  2:120700/166‑1 (MQ=255)                                                                                                                              ttttttCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCtt                                                                                                             >  1:602442/1‑251 (MQ=255)                                                                                                                                  ttCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATg                                                                                                         >  1:29177/1‑251 (MQ=255)                                                                                                                                   tCCGCACGACCTGCCGGAATTCAGTCTGACTGAAATTGATGAGATGCTGAAGCTGGAGTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTTTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATg                                                                                                         <  1:672127/250‑1 (MQ=255)                                                                                                                                     cGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCCATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCCCGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGGTTTATGGc                                                                                                       >  2:259999/1‑250 (MQ=255)                                                                                                                                       cACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCgg                                                                                                                                                                      <  1:300673/185‑1 (MQ=255)                                                                                                                                       aaCGACATGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCATCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCATGGACGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCgg                                                                                                                                                                      >  2:300673/2‑185 (MQ=255)                                                                                                                                          gACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGATTTGTTCCACGGGCCAACGCTGGCATTTAACGATTTCGGAGGTCGGTTTTTGGccaaaa                                                                                                  >  2:647825/1‑245 (MQ=255)                                                                                                                                              tGCCGGAATTCTGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTGTGGTACCCGAAGTGCGAATGTCCTGCCGGCGCTTATTGGGGACGAATTCCCGTAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGGTTGCCTTGCCGTCTCTGGTCTCCAATGTTGAAAGCGATGTCGGGTGTCTGGACTTGTTCCACGGGCCCACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATTTCACAAATGCt                                                                                              <  2:170117/250‑1 (MQ=255)                                                                                                                                                ccGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGCAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGTGCCAACGCTGGCCTTTAAAGATTTCGGCGGTCGCtttttgtgccaaatgatgac                                                                                           >  2:245566/1‑234 (MQ=255)                                                                                                                                                     aTTCAGCCTGACTGAAATTGAGTAGATGCTGAAGCTGGATTTTGTCTCCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCCGGAAATCCTGGAAGAGCGCGTGCGCTCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATTTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCata                                                                                     <  1:467866/251‑1 (MQ=255)                                                                                                                                                         aGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGc                                                                                 <  2:333529/251‑1 (MQ=255)                                                                                                                                                              gactgaAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCgggcg                                                                            >  1:508450/1‑249 (MQ=255)                                                                                                                                                                 tgaAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGAt                                                                          >  1:122206/1‑250 (MQ=255)                                                                                                                                                                 tgaAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGAAGTGCGACGATCCTCTCGGCGTTTATTGGTGCTGAAATCCCGCAGGAAATCCTGGAAGCGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATCGTTCCACGGGCCAACGCTGGCATTTACAGATTTCGGCGGTCGCTTTTTTGCACAAATGCT‑GACACATCATGCGGGCGAt                                                                          >  1:252535/1‑250 (MQ=255)                                                                                                                                                                  gaAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATa                                                                         <  1:149845/250‑1 (MQ=255)                                                                                                                                                                  gaAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGCCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACTGGCCAACGCTGGCATTTAAAGATTTCGGCTGCGGCTTTATGTCCCAAAAGCT‑CACCCATATTGccggacata                                                                         >  2:33812/1‑242 (MQ=255)                                                                                                                                                                   aaaTTGATGAGATGCGTAAGCTGGATGTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTGATTGGTGATGAAACCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAg                                                                       <  2:300655/251‑1 (MQ=255)                                                                                                                                                                    aaTTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGACGCTTTATGGCACAAATGTT‑GTCCCATATTTCGGGCGAtt                                                                         >  2:561051/1‑247 (MQ=255)                                                                                                                                                                            agaTGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGGTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATGTAAAGATTTCGGCGGTCGGTTTTTGGGCCAACAGCT‑GACACAAATTGTGGGCGATAAGCCAGGGAcc                                                              >  2:313348/1‑251 (MQ=255)                                                                                                                                                                              aTGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCAt                                                            <  1:531250/251‑1 (MQ=255)                                                                                                                                                                                gCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCAtt                                                           <  1:431656/250‑1 (MQ=255)                                                                                                                                                                                 cTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGCATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTTTGGCCAAAATACT‑GACCCCTATTGCGGGGGATAAACCAGTGaacctatt                                                         >  2:344184/1‑244 (MQ=255)                                                                                                                                                                                     aGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACCAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGAcc                                                     >  1:456319/1‑251 (MQ=255)                                                                                                                                                                                      gCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACcg                                                    <  1:232735/251‑1 (MQ=255)                                                                                                                                                                                       cTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGACTTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCCTATTGCGGGCGATAAGCCCGTGACCATTATGACcgc                                                   >  1:78186/1‑251 (MQ=255)                                                                                                                                                                                           aTTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGACGCTTTATGGCACAAAAGCTTGACC‑ATATTGCGGGCGGTAAGcccgagcacattctgcacccgccc                                               >  2:98832/1‑229 (MQ=255)                                                                                                                                                                                                  cACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTTACCGCGCCCTCCGGTg                                        >  1:588811/1‑251 (MQ=255)                                                                                                                                                                                                         aGTGCGAAGATCCTCTCGGCGTTTATTGGTGAAGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCgg                                 <  2:343832/251‑1 (MQ=255)                                                                                                                                                                                                          gTGCGAAGATCCTTTCGGAGTTTATGGGTGATGAAACCCCGCAGGAAATCCTGGAAGGGCGCGTGCGCGCGGCGGTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGAGGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCAGATTGCGGGCGATAAGCCAGTGTCCATTCTGACCGCGACCTCCGGTGATACCGGa                                <  2:602442/251‑1 (MQ=255)                                                                                                                                                                                                               aaGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGCAAGCGATGTCGGATGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGGGACCATTCTGACCGCCACCTCCGGCGATAACGGAGCGGc                           >  1:342248/1‑251 (MQ=255)                                                                                                                                                                                                                     ctctCGGCGTTTATTGGTGATGAAATCCCGCAAGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTAAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTAGCTTTATGGCGCAAATCCT‑GACCCATATTGCGGGCGATAAACCAGTGGCCATTATGACCGCGAGCTCCGGTGaatatcgggagagcacccg                     >  1:350777/1‑233 (MQ=255)                                                                                                                                                                                                                         cGGCGTTTATTTGTGATGAAATCCCGCAGGAAATCCTGGAAGGGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATTGCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCAt                 <  1:339337/251‑1 (MQ=255)                                                                                                                                                                                                                          ggCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGACCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGTGATAAGCCAGTAACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATg                <  1:144349/251‑1 (MQ=255)                                                                                                                                                                                                                           gCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCCATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGa                                >  1:437188/1‑234 (MQ=255)                                                                                                                                                                                                                           gCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGGGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTGCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCTACCTCCGGTGATACCGGa                                <  2:437188/234‑1 (MQ=255)                                                                                                                                                                                                                                tATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCCCAAATGCT‑GCACCATCTTGCgggcga                                                                           >  1:274255/1‑183 (MQ=255)                                                                                                                                                                                                                                tATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATGTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGTgggcga                                                                           <  2:274255/186‑4 (MQ=255)                                                                                                                                                                                                                                 aTTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGGGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTa         <  2:508450/251‑1 (MQ=255)                                                                                                                                                                                                                                  ttGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGTCGTTTGCCGTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACg                                                                                                                                          <  2:428558/122‑1 (MQ=255)                                                                                                                                                                                                                                  ttGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACg                                                                                                                                          >  1:428558/1‑122 (MQ=255)                                                                                                                                                                                                                                   tGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTAc        <  1:116731/250‑1 (MQ=255)                                                                                                                                                                                                                                     gTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGt     <  1:647825/251‑1 (MQ=255)                                                                                                                                                                                                                                        atgaAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTa  <  1:561051/251‑1 (MQ=255)                                                                                                                                                                                                                                         ggaAATCCCGAAGGAAACCCTGGAAGAGCGCGGGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATGGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTa  <  2:177311/249‑1 (MQ=255)                                                                                                                                                                                                                                          gaAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGt                                                                  <  1:155235/185‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  GGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTA  >  NC_000913/3696‑4177

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