Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,619,197	G→A	P323P (CCG→CCA)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,619,197	0	G	A	98.2%	171.5 / NA	57	P323P (CCG→CCA)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (29/26);  minor base C (1/0);  total (30/27)
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.84e-01

ACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGA  >  NC_000913/4618950‑4619441                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     acaggAATCAGGTTCTGGCCCTCAGTGCGGGTAGCGTGGTGTAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGGA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGTGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTAGGGTAGCGGCACAAAAAGGCccaacc                                                                                                                                                                                                                                                   <  1:360786/247‑1 (MQ=255) aaaTGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGGCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCccaacc                                                                                                                                                                                                                                                   <  2:148549/249‑1 (MQ=255)  caggAATCAGGTACAGGCCTCCAGGGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCATTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCGTTGGTCGTATGGTAGCGGCACAAAAAGGCccaaccg                                                                                                                                                                                                                                                  <  2:367642/248‑1 (MQ=255)  cATGAATCGGGTACGGGCCCCCAGTGCAGGTACCGCGGTTGAAGTTCGTGAATCGGTGCGGTTCCTGAGGGGTGA‑TACTCTTAGCCCGCGGCTGTTTGATGTCACGATGGCGCTGTGCGTGGGGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGTAGCGCGCGCGAAATGGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTGTGGTCGTATGGTAGCGGCACAAAAAGGCccaaccg                                                                                                                                                                                                                                                  <  2:63378/251‑1 (MQ=255)   aTGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGACGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGGAGCGGCACACAAAGGCccaaccg                                                                                                                                                                                                                                                  >  2:483374/1‑250 (MQ=255)              aCTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGAGAGTA‑TCGTAACACGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGACGCGCGCGAGACATTGCCAGCAGTGCTGGCCAACGGTAAAGCGGAAGAAGTCTATGGTCGCATGGTAACGGCACAAAAAGGCCCAACCGACTTCGTTgcg                                                                                                                                                                                                                                       >  2:19236/1‑248 (MQ=255)                 ggCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTACCCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACt                                                                                                                                                                                                                                   <  1:483374/251‑1 (MQ=255)                  gCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTGGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTa                                                                                                                                                                                                                                  <  1:78930/251‑1 (MQ=255)                     tCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTAcgc                                                                                                                                                                                                                               >  1:81775/1‑251 (MQ=255)                            cAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTAt                                                                                                                                                                                                                        <  1:38844/251‑1 (MQ=255)                            cAGGAAACGCGGTTGAAGTTCGTGAAGCGGGGCAGTTCCTGAAGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCAGTGAGTGGAGATGCTGAACTCCGGCACACTGGCGAAAGAAGACGCCGAAGCGCGAGCGAAATTGCAAGAAGTGCTGGACAAAGGTAAAGCGGCAGAAGTCTGTGGTCGTATGGTAGCGGCAAAAAAAGGACCAACCGTATTCGATGAGGACTACGCGaggt                                                                                                                                                                                                                          >  2:353377/1‑246 (MQ=255)                             cgggaaCGCGGTTGAGGTGCGTGAAGCGGTGCAGTTCCTGAGGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTGCGTTGAGAACTACGCGAAGTATc                                                                                                                                                                                                                       <  2:275808/247‑1 (MQ=255)                                           aaGTTCGTGACGCGGTGCAGTTCCGGACGGGTGAGTC‑TCGAAACCCGCGTCTGTTTGATGTCTCGATGGCGCTGTGCGTGGAGATGCTTATCGCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACATCGATg                                                                                                                                                                                                         <  2:406497/251‑1 (MQ=255)                                             gTTCGTGAGGCGGTGCAGTTCCGGCCGGGGGAGTA‑TCGTAACCCGCGGCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGAGGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCt                                                                                                                                                                                                       <  2:134121/251‑1 (MQ=255)                                             gTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCCCAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGcc                                                                                                                                                                                                                   >  2:35095/1‑239 (MQ=255)                                             gTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGcc                                                                                                                                                                                                                   <  1:35095/239‑1 (MQ=255)                                                            cAGTTCCTGACGGGAGAGTA‑TCGTAACCCGCGTATGTTTGATGACACGATGGAGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCACAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCt                                                                                                                                                                                         >  1:167956/1‑250 (MQ=255)                                                                     aCGGGTGAGTA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGAGGCACAAAAAGGCCCAACCGACTTAGTTGGAAACTACGCAAAGTCTCTGCAGACAGAGATGCTGACGAAAGCAGGTTATGCTGATAc                                                                                                                                                                               >  2:274465/1‑251 (MQ=255)                                                                           gagta‑tCGTCACCCGCGTCTGTTTGATGTCAAGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGACCCAACCGACTTCGTTGAGCACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAgg                                                                                                                                                                         >  1:369989/4‑251 (MQ=255)                                                                                     aaCCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGACGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAGGCAGTCTATGCTGATACCGAAGGTTTTGTc                                                                                                                                                                  >  2:135193/1‑249 (MQ=255)                                                                                         cgcgTCTGTTTGATGTCACGCTGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAACGCGGCAGAAGTCTTTGGTCGATTGGCAGCGGACCACAAAGGCCCAACCGACTTCGTTGGGAACTACGCGAAGTATCTGCCGCCCGCGATGCTGACGAAAGCAGTCTATGCTGCTACCGTAGGGTTTTCCAGTGa                                                                                                                                                             >  2:43828/1‑250 (MQ=255)                                                                                             tCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGCACTACGCGAAGTATCTGCCGACAGCGATGCTCACGACAGCAGTCTATGCTGATACCGAGGGATTTGTCAGTGAAATgg                                                                                                                                                        >  2:181474/1‑251 (MQ=255)                                                                                             tCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATgg                                                                                                                                                        >  1:65141/1‑251 (MQ=255)                                                                                               tGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAt                                                                                                                                                      >  1:22137/1‑251 (MQ=255)                                                                                                gTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAt                                                                                                                                                      >  1:100305/1‑250 (MQ=255)                                                                                                          aCGATGGCGCTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCt                                                                                                                                           <  1:54370/251‑1 (MQ=255)                                                                                                              gggCGCTGGGGGTGGAGATGCTTATCTCCGGAAAACTGGCGAAAGATTACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGGATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGGCTATGCTGATACCTAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTtggg                                                                                                                                       <  2:234539/250‑1 (MQ=255)                                                                                                                 cgcTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGACGCGCGCGAGAAATTGCAAGCAGTGCTGGACCACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAACGGCCCAACCGACTTCGTTGAGCACTACGCGAAGTATCTGCCGACAACCATGCTGACGAAAGCAGTCTAGGCTGATAACGACGGATTTGAAAGTGAAATGGCTAcagacccgctggggt                                                                                                                                      >  1:44852/1‑235 (MQ=255)                                                                                                                  gcTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATgg                                                                                                                                   >  1:86294/1‑251 (MQ=255)                                                                                                                  gcTGTGCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATgg                                                                                                                                   <  2:214338/251‑1 (MQ=255)                                                                                                                    ggtgCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCa                                                                                                                                 <  1:319973/250‑1 (MQ=255)                                                                                                                     gtgtgtgtagaTGCTGATCTCCGGCAAACTGGCAGAAGATTACGCCGAAGCGCGCGCTAAATTGCAATCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGTTCGTATGGTAGCGGCCCAAAAAGGCCCAAACGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAACGCAGTCTTTGCTGAGACCGATGGTTTTGTCAGTGAAAGGGATACCCGCTCGCTGGGGATGGc                                                                                                                                  <  2:324667/241‑1 (MQ=255)                                                                                                                     gtgCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTACAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGc                                                                                                                                  <  1:449344/249‑1 (MQ=255)                                                                                                                     gggCGTGGAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCa                                                                                                                                 <  2:66729/248‑1 (MQ=255)                                                                                                                           ggAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTg                                                                                                                          >  1:356684/1‑251 (MQ=255)                                                                                                                           ggAGATGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGGCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTGCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGtt                                                                                                                           <  2:213759/250‑1 (MQ=255)                                                                                                                               aTGCTTATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAAt                                                                                                                      <  1:242599/251‑1 (MQ=255)                                                                                                                                    taTCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGggc                                                                                                                  <  1:420623/249‑1 (MQ=255)                                                                                                                                    taTCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAAGTGCAATCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCGTTGGTCGTATGGTATCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTGTTGTCAGTGAAAGGGATACCCTCGCGCTGGGGATGGCAGTGGTTGCAATGggc                                                                                                                  <  2:294957/249‑1 (MQ=255)                                                                                                                                      tctcCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTCAAGCGGCAGAAGTCTTTGGTCGTATGGAAGCGGCACAACAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGATTGGTCAGTGCAATGGATACCAGCGCGCTGGGGATGGCGGTGGTTGCAATGggcggc                                                                                                               >  1:259812/1‑251 (MQ=255)                                                                                                                                              aaaCGGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGcggcggtcgtcg                                                                                                       <  2:65141/251‑7 (MQ=255)                                                                                                                                                aCTGGCGAAAGATGACGCCGCGGCGCGCGAGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCCAAAGGTAGCGGCACAAAACGGCCCATCCGACTTCGTTGAGAAATACGCGAAGTATCTGCCGACAGCGATGATGACGAAAGCAGTCTATGCTGATACCGAAGGGTTTGTCCGTGAAATGGGTACCCGCGCGATGGGGACGGCAGTGGTTGCAAGGGGgggcggtcgtccgc                                                                                                     >  2:490118/1‑243 (MQ=255)                                                                                                                                                 cTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGgtcgtcgtc                                                                                                     >  1:205469/1‑250 (MQ=255)                                                                                                                                                      gAAAGATGACGCCGAAGCGCGCGCGAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGGCTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCAt                                                                                               >  2:379803/1‑251 (MQ=255)                                                                                                                                                              aCGCCGAAGCGCGCGCGAACTTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGCCTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGCCGCAAGCCGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGCTGGCGGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGGATCTGACAcc                                                                                       >  2:496009/1‑251 (MQ=255)                                                                                                                                                                             cgAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTc                                                                        >  2:479376/1‑251 (MQ=255)                                                                                                                                                                              gAAATTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCg                                                                       <  1:228064/251‑1 (MQ=255)                                                                                                                                                                                 aTTGCAAGCAGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAAACGACTGCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCt                                                                    <  2:372800/251‑1 (MQ=255)                                                                                                                                                                                      gcgcagTGCTGGACAACGGTCAAGCGGCGGAAGTCTTTGGTCGTATGGTAGCGGAACAACAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCAGACAGCGCTGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGTTCCCCGCGCGCTGGGGGTGGCAGTGGTTGCAATGGGCGGCGGTCGTCGGCAGGCCTCTGACCCCCTTGATTACaacggagcctttctt                                                               >  2:103888/6‑237 (MQ=255)                                                                                                                                                                                      aagcagTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAgg                                                                                                                                                                         <  2:73679/140‑1 (MQ=255)                                                                                                                                                                                      aagcagTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAgg                                                                                                                                                                         >  1:73679/6‑145 (MQ=255)                                                                                                                                                                                                     aCGGAAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAACAGGCACAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGATGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGTAATGGATACCCGCCCGCTGGGGATGGGAGTGTTTGCAATGGGGGGCGGTCGGCGTCGAGCATCTGACACCATTGATTACAGCGACGGCTTTACTGATATGGCGCGtttg                                                >  2:355067/1‑248 (MQ=255)                                                                                                                                                                                                         tAAAGCGGCCGAAGGCGTTGGTCGAATGGGAACGGCCCAAAAAGGCCCCACCGACTTCGTTGAGAACGCAGCGAAGTATCAGCAGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACACGCGCGATGGGGATGGCAGTGGTTGCCATGGGCGGCGGTCTACTTCCGGCATCTGCCACCATTGTTTACAGCGTCGGCGTTACTGCTATGGCGCGTTTggggg                                            >  2:492891/1‑249 (MQ=255)                                                                                                                                                                                                                          ttGGCCGTCTGGTACCGGTACAAAAGTGCCCGACCGACTTCATGGAGACCTCCGCGAATCACCTGACGTCCGCCATGCTGACGAACGCAGTCTATGCTCATACCGAAGGTTTTGTCAGGGAATCGGTTACCCGCGCGCTGGGGATGGCAGTGGTTGCAAGGGGCGGCGGTCTTCGTCAGGCATCTGACGCCATTGATTACAGCGTCGGATTTACTGATATGGCGCGTTTGGGCGagcattgagccggtca                            <  2:62730/250‑16 (MQ=255)                                                                                                                                                                                                                           tGGTCGTATGGTAGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAg                           <  1:274465/250‑1 (MQ=255)                                                                                                                                                                                                                                       aGCGGCACAAAAAGGCCCAACCGGCTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCg               <  1:379803/250‑1 (MQ=255)                                                                                                                                                                                                                                       aGCGGCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGTCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCg               <  1:77781/250‑1 (MQ=255)                                                                                                                                                                                                                                          ggCACAAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTa           <  1:19236/251‑1 (MQ=255)                                                                                                                                                                                                                                              caAAAAGGCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACACGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGGCGCCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCATTACTGATATGGCGCGTTTGGGCGGCGAGGTAGACGGGCAGCGTCCGCTGGCGGTTaaccc       >  2:193677/1‑247 (MQ=255)                                                                                                                                                                                                                                                   aGGCCCAACCGACTTCGTTGCGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGAAACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGGCAGGCATCTGACACCATTGATTACAGCGGCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGGCAGCGTGCGCTGGGGGGTATctcaccgg  >  2:47189/1‑244 (MQ=255)                                                                                                                                                                                                                                                    ggCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCAcgcgc  >  1:476055/1‑249 (MQ=255)                                                                                                                                                                                                                                                    ggCCCAACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGTCGTCGTCAGGCATCTGACACCATTGATTACAGCGTCGGCTTTACTGATATGGCGCGTTTGGGCGACCAGGGAGACGGTCAGCGTCCGCTGGCGGGTATCCAcgcgg  >  2:437241/1‑249 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     ACATGAATCAGGTACTGGCCTCCAGTGCAGGTAACGCGGTTGAAGTTCGTGAAGCGGTGCAGTTCCTGACGGGTGAATA‑TCGTAACCCGCGTCTGTTTGATGTCACGATGGCGCTGTGCGTGGAGATGCTGATCTCCGGCAAACTGGCGAAAGATGACGCCGAAGCGCGCGCGAAATTGCAGGCGGTGCTGGACAACGGTAAAGCGGCAGAAGTCTTTGGTCGTATGGTAGCGGCACAAAAAGGCCCGACCGACTTCGTTGAGAACTACGCGAAGTATCTGCCGACAGCGATGCTGACGAAAGCAGTCTATGCTGATACCGAAGGTTTTGTCAGTGAAATGGATACCCGCGCGCTGGGGATGGCAGTGGTTGCAATGGGCGGCGGACGCCGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGA  >  NC_000913/4618950‑4619441

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