Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,621,321	A→G	I147V (ATT→GTT)	deoD →	purine nucleoside phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,621,321	0	A	G	95.5%	126.0 / ‑2.0	45	I147V (ATT→GTT)	deoD	purine nucleoside phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (19/23);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-01

GGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACG  >  NC_000913/4621077‑4621552                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        gggTATCCCGTCCTGCTCCATCTACACCAAAGAACAGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCTCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCTCCGATTCCAAAGTTAACCGTATCCGTTTTAAATACCATGACGTTGCCGCTATCGCTGATTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATg                                                                                                                                                                                                                                   <  1:230529/251‑1 (MQ=255)                cTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGGGTTGATGCTCGCGTGGGTAAcc                                                                                                                                                                                                                    >  1:147850/1‑251 (MQ=255)                cTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGGGTTGATGCTCGCGTGGGTAAcc                                                                                                                                                                                                                    >  1:304555/1‑251 (MQ=255)                   cATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTTTTGATGCTCGCGTGGGTAACCTGt                                                                                                                                                                                                                 <  1:347029/251‑1 (MQ=255)                                        gacgagatcgaccatcaaatcggcaaTCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTAct                                                                                                                                                                                            <  1:355799/226‑1 (MQ=255)                                        cACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTAct                                                                                                                                                                                            >  2:420575/1‑251 (MQ=255)                                             aTTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCg                                                                                                                                                                                       >  1:381747/1‑251 (MQ=255)                                                  ggCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGGTCTCCGCTGGCCTGTTCTACTCTCCGGACgg                                                                                                                                                                                  >  2:143762/1‑251 (MQ=255)                                                             aaTTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCg                                                                                                                                                                       >  1:380352/1‑251 (MQ=255)                                                             aaTTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCg                                                                                                                                                                       <  1:420575/251‑1 (MQ=255)                                                                tATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGAc                                                                                                                                                                     <  2:381747/250‑1 (MQ=255)                                                                    cgcgTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGAt                                                                                                                                                                <  1:143762/251‑1 (MQ=255)                                                                           gTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGaaaaa                                                                                                                                                         <  1:161775/251‑1 (MQ=255)                                                                              ccTGGGGCGCAGTTCTGCCGCACGAAAAACGGCGCGACGTCGTTATCGGTATGGGGGCATGCACCGATTCCAAAGTTAACCGCAGCCGTGTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACAGGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGGGGGTAACCTGTTCTCCGCTGAACTGTTCTAATCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATa                                                                                                                                                       <  2:155509/250‑1 (MQ=255)                                                                                 gtgGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGc                                                                                                                                                   >  1:23434/1‑251 (MQ=255)                                                                                               ccGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTgg                                                                                                                                      >  1:164849/1‑250 (MQ=255)                                                                                                    cGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATgg                                                                                                                                <  1:201000/251‑1 (MQ=255)                                                                                                       aaaaCTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGCTTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGGAATGTTCGCCGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGagg                                                                                                                             >  2:470667/1‑249 (MQ=255)                                                                                                        aaaCTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCCAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGCCGGCGAAATGTTCGACGTGATGGAAAAATACGGCATACTCGGCGTGGAAATTGaa                                                                                                                              >  2:95056/1‑249 (MQ=255)                                                                                                        aaaCTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGc                                                                                                                            >  1:228044/1‑251 (MQ=255)                                                                                                        aaaCGGCGCGACGTCGTTAGCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGc                                                                                                                            <  2:228044/251‑1 (MQ=255)                                                                                                           cggcgcGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTAGCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGc                                                                                                                         <  2:36583/249‑1 (MQ=255)                                                                                                             gcgcGACGTCTTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGTGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGc                                                                                                                         <  2:304555/249‑1 (MQ=255)                                                                                                              cgctaCGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTgg                                                                                                                      <  2:23434/247‑1 (MQ=255)                                                                                                                   cgtcgtTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCt                                                                                                                 >  1:186352/1‑251 (MQ=255)                                                                                                                         tATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGc                                                                                                            <  2:186352/250‑1 (MQ=255)                                                                                                                          aTCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCg                                                                                                           <  1:405699/250‑1 (MQ=255)                                                                                                                              gTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCCTCCTCGGCGGGGAAATGGAAGCGGCTGGGATATACGGCGTCGCt                                                                                                      >  2:307685/1‑251 (MQ=255)                                                                                                                                  gggTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCa                                                                                                   <  1:76856/250‑1 (MQ=255)                                                                                                                                          gCACCGATTCCAAAGTTAACCTTATCCGTTTTAAAGTCCATGACTTTGTCGCTATCGCTGACTTTGTCATGGTGCGTAACGCGGGAGTCGCGGCTAAAGCACTGGGTGTTTATGCTCGCGTGGGTAACCTGGTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTTATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGgc                                                                                          <  1:229477/251‑1 (MQ=255)                                                                                                                                            aCCGATTCCAAAGTTAACCGTATCCGTTTTAACGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGCAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGGATCTACGGCGGCGCGGCCGAGTTTGgcgc                                                                                        >  2:270502/1‑251 (MQ=255)                                                                                                                                            aCCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGgcgc                                                                                        >  1:429131/1‑251 (MQ=255)                                                                                                                                                     aaaGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCt                                                                               >  2:201175/1‑251 (MQ=255)                                                                                                                                                                        ttAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCt                                                            >  1:128516/1‑251 (MQ=255)                                                                                                                                                                          aaaGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGa                                                          <  1:201175/251‑1 (MQ=255)                                                                                                                                                                             gACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACca                                                       <  1:307685/251‑1 (MQ=255)                                                                                                                                                                                cATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACAt                                                    <  1:423575/251‑1 (MQ=255)                                                                                                                                                                                  tGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGGAAGCCCTGACCATCTGCACCGTTTCTGACCACATcc                                                  >  1:40757/1‑251 (MQ=255)                                                                                                                                                                                                         ttgACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTTGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTg                             <  1:311499/247‑1 (MQ=255)                                                                                                                                                                                                          ggACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGcc                           <  2:128516/249‑1 (MQ=255)                                                                                                                                                                                                                    cGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCa                <  1:227909/251‑1 (MQ=255)                                                                                                                                                                                                                       aacgcgggagacgctgCTAAAGTACTGGGTGTTGATGCTCGCGTGGGTAACCTTTTCTCCGCTGACCTGTTCTACTCGCCGGACGGCAAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGTAGACCACTGCCGCTGAGCGTCAGa              <  2:347927/235‑1 (MQ=255)                                                                                                                                                                                                                       aaCGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGAc             <  1:3327/251‑1 (MQ=255)                                                                                                                                                                                                                          gcggTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGAccac          <  2:164849/248‑4 (MQ=255)                                                                                                                                                                                                                                  cgcggCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACCACCTTCAACg  >  1:254209/5‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        GGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACTACCTTCAACG  >  NC_000913/4621077‑4621552

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