Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	410,666	C→T	P175S (CCT→TCT)	mak →	fructokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	410,666	0	C	T	95.7%	124.2 / ‑3.4	47	P175S (CCT→TCT)	mak	fructokinase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (18/26);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

CGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCCCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGA  >  NC_000913/410422‑410903                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              cgcgAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATTGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGGTGGACGAGGACGGACTGCGTTATCGCGAGGAAGTCTCtttgt                                                                                                                                                                                                                                         >  1:145278/1‑248 (MQ=255)            aGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGAGGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTaaa                                                                                                                                                                                                                              <  2:250842/251‑1 (MQ=255)                   aGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGtt                                                                                                                                                                                                                       <  1:294035/251‑1 (MQ=255)                          gtGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGa                                                                                                                                                                                                                <  2:105821/250‑1 (MQ=255)                             gCAAATGACGGTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAAc                                                                                                                                                                                                             <  1:265432/251‑1 (MQ=255)                                                      tttCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGGGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGGGGGGGAACAATCCGCTGCCGTGGGTGGAAGAAGACGGACGTCGTGATCTCGCGGGAGGGTCTTGGTATTTGGGTAAGCAAGTGTGTAGTGAAACCTTTATTTCGGGTACGGTATTcc                                                                                                                                                                                      >  2:130523/1‑248 (MQ=255)                                                          aGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTACCGGCACGGGATGCGGCGCGGGCGTGGCATTCACTGGGCGGGCGCATATCGGCGGCAATGGCCCGGCAGGAGGGTGGGGAACCCATCCGCTACCGTGGATGGGCGCAGACGCACTGCGTTATCGCGAGGAAGGCTCTTGTTATTGCGGTAAAAAAGGCTGTATTGAAACCTTTATTTCAGGCACGGGATTCGCGGTgg                                                                                                                                                                                >  1:201412/1‑251 (MQ=255)                                                                aGTAGATGGCGCGGCAGCGGGTGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTAtc                                                                                                                                                                          <  1:141475/251‑1 (MQ=255)                                                                    gATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGAGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCCGGGGGGGGGGGGCACAATCCGACACCGTGGATGGCCGGAGAAGAACAGCGGTTTCGCGGAGAGGTCTCTTGTTAATTCGGGAAAAAAGGGTGTATTTAAACCTTTTTTTCGGGCCGGGGGTTTGCGATGGAtttctgtcg                                                                                                                                                                      >  2:264267/1‑244 (MQ=255)                                                                       ggCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGttt                                                                                                                                                                   >  1:435539/1‑251 (MQ=255)                                                                            ggCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCg                                                                                                                                                              >  1:382642/1‑251 (MQ=255)                                                                             gCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACAAAATCCGCTACCGAGGATGGAAGAAGACGAACTGCGTTATCGCGGGGAAGGCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGCTTCGCGATGGATTATCTTCGTTTGGGCgg                                                                                                                                                             >  1:293263/1‑251 (MQ=255)                                                                                 ggggAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACAt                                                                                                                                                         <  1:251649/250‑1 (MQ=255)                                                                                      gcgcAGACGGGATTTGCCGTGATTTTCGGAACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCt                                                                                                                                                    <  2:376741/251‑1 (MQ=255)                                                                                              ggTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACGATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCCCTTGTTATTGCGGTAAACAAGGtcg                                                                                                                                                                                                                      <  1:184243/177‑3 (MQ=255)                                                                                              ggTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGGGGGTGGGGACACGATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCCCTTGTTATTTCGGTAAACAAGGtcg                                                                                                                                                                                                                      >  2:184243/1‑175 (MQ=255)                                                                                                 gtttGCCGTGGTTATCGGCACGGGCTGCGGCGCGGGCGTGGAACTCAAGGGGCGGGCGCGTATCGGTGGGAATGGCACGGCAGGTGAGTGGGGACACAAGCCGCTGCCGTAGATGGACGAAGACGGACGGCGTTACCGCGATGAAGTCGTGTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTGTATTTCGGGCACGGGATTCGCGATGGATTATCGGCGTTTGAGCGGACATGCACTGAAAGGCAGt                                                                                                                                          <  2:193583/249‑1 (MQ=255)                                                                                                   ttGCCGTGATTATCGGCACGGGAGGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGa                                                                                                                                        <  2:439925/250‑1 (MQ=255)                                                                                                                   cACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTg                                                                                                                       >  1:109899/1‑251 (MQ=255)                                                                                                                      gggATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGAAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGGCGTTTGAGCGGACATGAGCTGAAAGGCAGGGAAATTATCCGCCTGGTTga                                                                                                                      >  1:502403/1‑249 (MQ=255)                                                                                                                           gcggcgCGGGCGTGGCATTTAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTTGGGACGCAATCCGGTATCGTGGATTGACGAAGTCGAACTGCGTTATGGCGAGGAAGTCTCTTGTTATTGCGGTAAACATGGTTGTAGTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTAGCGTCGTTTGAGCGGACATGCGCGGAAAGGCAGTGAAATTATCCGCATGGTTgaagaacg                                                                                                                <  2:337355/250‑3 (MQ=255)                                                                                                                            gggcgCGGGCGTGGCATTCAATGGGCGGGCGCATAGCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTTCCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGAGTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGACATTATCCGCCTGGTTTAAGAAAGCg                                                                                                              <  2:365401/250‑1 (MQ=255)                                                                                                                               ggggggggtggactgcAATGGGCGGGGGCATACGGGCGGCATGGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGATCTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATc                                                                                                           <  2:195981/236‑1 (MQ=255)                                                                                                                                   ggCGTGGCATGCAAGGGGCGGGCGCATATCGGGGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGAGCCGGt                                                                                                       <  2:332054/251‑1 (MQ=255)                                                                                                                                                            atatCGGCGGCAATGGCACGGCAGGTGAGTGGGGTCACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTGCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTAc                                                                              <  2:301117/251‑1 (MQ=255)                                                                                                                                                             tatCGGCAGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCtccg                                                                             >  1:379202/1‑248 (MQ=255)                                                                                                                                                                    ggaaagggcccggtaGGTGAGTGGGGACACAATCCGCTTCCGTGGAGGGACGAAGACGAACTGCGTTATCGCGAGGAAATCTCTTGTTATTGCGGTAAACAAGGTTGGATTGAAACCGTTATTGCGGGCACGGGATTCGCGATGGATTAGCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCg                                                                      <  2:293263/237‑1 (MQ=255)                                                                                                                                                                               ggCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAAt                                                           <  1:45574/251‑1 (MQ=255)                                                                                                                                                                                gCAGGTGAGTGGGGACACAATCCGCGACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTACGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAATTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATc                                                          <  1:364025/251‑1 (MQ=255)                                                                                                                                                                                        gTGGGGACACAATCCGCTACAGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGAGATGGATTATCGTCGTTTGCGCGGACATGCGATGAACGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGAATTGCGTCGCTACGCGCTGCGGCAGGCAAAATCGCTGGcac                                                  >  2:499873/1‑251 (MQ=255)                                                                                                                                                                                        gTGGGGACACAATCCGCCACCGTGGATGGACGAAGACGAACTGCGGTAGCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGTGACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACAAGCGCTGAAAGGCAGAGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGcac                                                  <  2:288082/251‑1 (MQ=255)                                                                                                                                                                                               cacaAGCCGCTACCGGGGATGGACGAAGACGAAGGGCGTTATCTGGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTGTATTTCGGGCACGGGAGACGCGATGGATGATCGTCGTTTGAGCGGACAGGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCAGTGCGGCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGt                                           <  2:201412/251‑1 (MQ=255)                                                                                                                                                                                                   aTCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAat                                       <  1:267103/251‑1 (MQ=255)                                                                                                                                                                                                    cccGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAata                                      <  1:133254/250‑1 (MQ=255)                                                                                                                                                                                                     ccGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGt                                              >  2:232601/1‑242 (MQ=255)                                                                                                                                                                                                     ccGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGt                                              <  1:232601/242‑1 (MQ=255)                                                                                                                                                                                                       gCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAACGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCACACTCGCTGGCCCATTTCGTGAATATtc                                   >  2:29448/1‑251 (MQ=255)                                                                                                                                                                                                       gCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATtc                                   <  1:369285/251‑1 (MQ=255)                                                                                                                                                                                                          agcGTGTATGGACGAGGACGAACTGCGTTTTCGCGAGTAAGTCTCTTGTTATTGCGGTAAACAAGTTTGTATTCATACCTTTATTTCGGGCACGGGATTCGCGATTGATTATCGTCGTTTGAGCGGACATGCTCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCTCTCCGAGCTGCGGCTGGCAAAATCGCTGGCACCTGTCGTGAATATtctc                                 <  2:452899/248‑1 (MQ=255)                                                                                                                                                                                                            cGTGGATGGACGAAGACGAACTGCGTTATCGCTAGGAAGGCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGATCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGAt                              <  2:145278/251‑1 (MQ=255)                                                                                                                                                                                                                     acgaagacgaaCTGCGTTATCGCGAGGAAGTCTCTTGTTATTTCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGTCATGCGCTGAAAGGCAGTGGAATTATCCGCCTGGTTGAAGAAAGCGATCCGTTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGATGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGTTCCggttttg                     >  2:49090/1‑246 (MQ=255)                                                                                                                                                                                                                     acgaagacgaaCTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGAtgtg                     >  1:96192/1‑251 (MQ=255)                                                                                                                                                                                                                     acgaagacgaaCTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGAtgtg                     >  1:423824/1‑251 (MQ=255)                                                                                                                                                                                                                             gaaCTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCt             >  1:289728/1‑251 (MQ=255)                                                                                                                                                                                                                                cTGCGTTATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCgggg          >  2:459427/1‑251 (MQ=255)                                                                                                                                                                                                                                      tATCGCGAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCg     <  1:264267/250‑1 (MQ=255)                                                                                                                                                                                                                                         cgcgAGGAAGTCTCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGGTTGTCCTGGGGGGCGGGa  <  1:73786/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              CGCGAGGTTGCAGCGGGAAGTGCGGCTGGCAAATGACGCTAACTGTCTGGCGGTTTCAGAAGCAGTAGATGGCGCGGCAGCGGGAGCGCAGACGGTATTTGCCGTGATTATCGGCACGGGATGCGGCGCGGGCGTGGCATTCAATGGGCGGGCGCATATCGGCGGCAATGGCACGGCAGGTGAGTGGGGACACAATCCGCTACCGTGGATGGACGAAGACGAACTGCGTTATCGCGAGGAAGTCCCTTGTTATTGCGGTAAACAAGGTTGTATTGAAACCTTTATTTCGGGCACGGGATTCGCGATGGATTATCGTCGTTTGAGCGGACATGCGCTGAAAGGCAGTGAAATTATCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGA  >  NC_000913/410422‑410903

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