Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,694,217	C→T	A443T (GCA→ACA)	purL ←	phosphoribosylformylglycinamide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,694,217	0	C	T	100.0%	222.4 / NA	62	A443T (GCA→ACA)	purL	phosphoribosylformylglycinamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/28);  total (34/28)

GCCCGGTTCGTCGCTTAGAATCTCGCGCAGTTCAAATTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGT  >  NC_000913/2693936‑2694506                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                  gCCCGGTTCGTCGCTTAGAATCTCGCGCAGTTCAAATTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGaa                                                                                                                                                                                                                                                                                 >  1:184305/1‑299 (MQ=255)            cgCTTAGAATCTCGCGCAGTTCAAATTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTTGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCAc                                                                                                                                                                                                                                                                      <  2:49385/300‑1 (MQ=255)            cgCTTAGAATCTCGCGCAGTTCAAATTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCAc                                                                                                                                                                                                                                                                      <  1:229791/300‑1 (MQ=255)            cgCTTAGAATCTCGCGCAGTTCAAATTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCAc                                                                                                                                                                                                                                                                      <  2:352195/300‑1 (MQ=255)              cTTAGAATCTCGCGCAGTTCAAATTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCg                                                                                                                                                                                                                                                                    >  1:439812/1‑300 (MQ=255)                                                     ccgTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATg                                                                                                                                                                                                                             <  1:586563/300‑1 (MQ=255)                                                     ccgTCGCTCACCAGGTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACATGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGGACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCAATGTTCATTTTCGGGCCACCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATg                                                                                                                                                                                                                             <  2:351411/300‑1 (MQ=255)                                                              aCCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCNACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATc                                                                                                                                                                                                                    <  1:523226/300‑1 (MQ=255)                                                                             gCGTTAGGAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTACGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCACGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATg                                                                                                                                                                                                     >  1:549051/1‑300 (MQ=255)                                                                             gCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAATAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCAAAGAACGACCAGCTTCGCACCGACGTTGATCTCACCTTTTTGTACGTGATCGGCGCGAACGTTGCAGATCCAGCCCGACAGCATg                                                                                                                                                                                                     >  2:278090/1‑300 (MQ=255)                                                                             gCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCNCCAAGCTACCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATNTCGCATTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATg                                                                                                                                                                                                     >  2:695038/1‑300 (MQ=255)                                                                             gCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATg                                                                                                                                                                                                     >  1:668029/1‑300 (MQ=255)                                                                                          cGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGAt                                                                                                                                                                                        <  1:616254/300‑1 (MQ=255)                                                                                                cgcCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGNCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGAACCGACGTTGATCTCGCCTTTTTGTACGTGGTCGGCGCGGATGTTGCCGATCCCGCCCGCCAGCATGATTGGTTTGTGATGACcgc                                                                                                                                                                                  >  2:165260/1‑300 (MQ=255)                                                                                                cgcCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACcgc                                                                                                                                                                                  >  1:242290/1‑300 (MQ=255)                                                                                                cgcCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTACTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATAAGACTGACCAGACGCCATAGAAGACGCTGCACCAACGNCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCACCTTTTTGTACGTGATCGGTGCGAATGTTGCCGATTCCGCCCGCCAGCATGATCGGTTTATGCTAACcgc                                                                                                                                                                                  >  2:162297/1‑300 (MQ=255)                                                                                                       gtcgtGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGctc                                                                                                                                                                           <  2:675548/300‑1 (MQ=255)                                                                                                       gtagtGGATAAACAGGCTTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTAAGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCACGACCGATGTTCATTGTCGGGCCGCCGAGAACGACTAGCTTCGCACCGACGTTGATCTAGCCTTTTTGTACGTGATCGTCGAGAACGTCGCCGATCCTGCCCGCCCGCATGATCGGTTTGCGGCAACCgcgcgactc                                                                                                                                                                           >  2:41203/4‑295 (MQ=255)                                                                                                                      gATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCAGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCt                                                                                                                                                            <  1:659808/300‑1 (MQ=255)                                                                                                                          gggTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTACGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACNGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTc                                                                                                                                                        >  1:522816/1‑300 (MQ=255)                                                                                                                            gTTGGCATCACCAAGCTGCCAGCAACGTTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCAc                                                                                                                                                      <  1:304019/300‑1 (MQ=255)                                                                                                                            gTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCAc                                                                                                                                                      <  1:305863/300‑1 (MQ=255)                                                                                                                               ggCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTACGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCAGTTGTGGCTGTTCACttt                                                                                                                                                   >  2:559916/1‑300 (MQ=255)                                                                                                                                               cAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTACGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGAC‑AGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGa                                                                                                                                  >  2:3049/1‑300 (MQ=255)                                                                                                                                                 gCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAGCGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTNATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGaa                                                                                                                                 >  2:645704/1‑300 (MQ=255)                                                                                                                                                               cTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGNGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTNCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTg                                                                                                                   >  2:695517/1‑300 (MQ=255)                                                                                                                                                               cTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTg                                                                                                                   >  1:244831/1‑300 (MQ=255)                                                                                                                                                                ttCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCNAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGc                                                                                                                  <  2:652742/300‑1 (MQ=255)                                                                                                                                                                     ggCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACNGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGgac                                                                                                             <  2:637585/300‑1 (MQ=255)                                                                                                                                                                     ggCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTTTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGgac                                                                                                             <  1:202993/300‑1 (MQ=255)                                                                                                                                                                     ggCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGgac                                                                                                             <  1:697127/300‑1 (MQ=255)                                                                                                                                                                     ggCAGCGACGCTCCATCTCCGGGTTGTCACGCGGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCATACGCCGTAGGAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCACAGCTCTTCGCCGTTGCGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCATTCAGTGCCGgac                                                                                                             <  2:127572/300‑1 (MQ=255)                                                                                                                                                                                    tctcCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCgttgt                                                                                              >  1:376950/1‑300 (MQ=255)                                                                                                                                                                                         gggTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAAc                                                                                         <  1:314313/300‑1 (MQ=255)                                                                                                                                                                                                 gcgcTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAACCGTTCAGTGCCGGACGACCAAATTCGTTGTTACgcgccgcgcc                                                                                 >  2:563004/1‑300 (MQ=255)                                                                                                                                                                                                     tGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGccc                                                                             <  1:67291/300‑1 (MQ=255)                                                                                                                                                                                                      gTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGTCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCa                                                                            >  1:752883/1‑300 (MQ=255)                                                                                                                                                                                                      gTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATNTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTNGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTAGTTGTTACACGCCGCGCCGCCCa                                                                            >  2:590656/1‑300 (MQ=255)                                                                                                                                                                                                      gTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCa                                                                            >  1:574503/1‑300 (MQ=255)                                                                                                                                                                                                      gTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCAGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGAACACATTCGTTGTTGAACgcagcgccgcaca                                                                            >  2:328781/1‑297 (MQ=255)                                                                                                                                                                                                      gTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCAATGTTCATTGTCGGGCCGCCGAGAACAACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCa                                                                            >  1:682580/1‑300 (MQ=255)                                                                                                                                                                                                          ggAAGCAAAGTCGAGTTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCgg                                                                        <  2:402723/300‑1 (MQ=255)                                                                                                                                                                                                                   gTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGt                                                               >  1:721577/1‑300 (MQ=255)                                                                                                                                                                                                                         gTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGAACAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCatgat                                                         >  1:326704/1‑300 (MQ=255)                                                                                                                                                                                                                          tCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCatgatg                                                        <  1:313045/300‑1 (MQ=255)                                                                                                                                                                                                                             gCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTcc                                                     >  1:76657/1‑300 (MQ=255)                                                                                                                                                                                                                             gCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCACCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTcc                                                     <  2:435664/300‑1 (MQ=255)                                                                                                                                                                                                                                  aGACTGACCAGACGCCATAGAAGACGCTGCCCCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGAACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGANTGTTGCCGATCCCGCCCGCCAGCATGNTTGGTTTGTGGTAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGACGCCCAGCGGGCCTTCGGTCATGATGTCCAgcgc                                                >  2:496273/1‑300 (MQ=255)                                                                                                                                                                                                                                  aGACTAACCATACGCTATAGAGGACGCTGCACCACCGCCAATGCCGATGTTCATTGTCGGGCCACCGAGAACGACCAGCTTCACACCGACGTTGACCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCNATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAgcgc                                                <  2:164341/300‑1 (MQ=255)                                                                                                                                                                                                                                       gACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTNACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGa                                           >  1:484705/1‑300 (MQ=255)                                                                                                                                                                                                                                       gACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGa                                           >  1:465107/1‑300 (MQ=255)                                                                                                                                                                                                                                            gACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATg                                      <  1:560193/300‑1 (MQ=255)                                                                                                                                                                                                                                                ccATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTAATAAGTACGGAAGTAGCCGTTCAGTGCCGGATGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATgcact                                  >  2:153730/1‑297 (MQ=255)                                                                                                                                                                                                                                                     gaagaCGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGNCCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCNCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGc                             <  2:678595/300‑1 (MQ=255)                                                                                                                                                                                                                                                     gaagaCGCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGc                             >  1:447125/1‑300 (MQ=255)                                                                                                                                                                                                                                                           gCTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGGTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCg                       <  1:385732/300‑1 (MQ=255)                                                                                                                                                                                                                                                            cTGCACCACCGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGa                      >  1:608024/1‑300 (MQ=255)                                                                                                                                                                                                                                                                     cGCCAAGACCGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATcttctt             >  1:274042/1‑300 (MQ=255)                                                                                                                                                                                                                                                                             ccGATTTTCATTGTCTGACCGCCGAAAACGGCCCGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCGGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCGAATTCGTCGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCCTGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGc     <  2:456667/300‑1 (MQ=255)                                                                                                                                                                                                                                                                             ccGATGTTCATTGTCGTGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGc     <  2:621454/300‑1 (MQ=255)                                                                                                                                                                                                                                                                             ccGATGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTCACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGc     <  2:417717/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                aTGTTCATTGTCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGt  >  1:82027/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                  GCCCGGTTCGTCGCTTAGAATCTCGCGCAGTTCAAATTTACCGCCGCGCCCGCCGTCGCTCACCAGTTCCGGCATGGCGTTAGAAAGACCGCCAGCGCCAACGTCGTGGATAAACAGGATTGGGTTGGCATCACCAAGCTGCCAGCAACGGTCGATCACTTCCTGGCAGCGACGCTCCATCTCCGGGTTGTCGCGCTGTACGGAAGCAAAGTCGAGGTCGGCATCAGACTGACCAGACGCCATAGAAGACGCTGCACCACCGCCAAGACCGATGTTCATTGCCGGGCCGCCGAGAACGACCAGCTTCGCACCGACGTTGATCTCGCCTTTTTGTACGTGATCGGCGCGAATGTTGCCGATCCCGCCCGCCAGCATGATCGGTTTGTGATAACCGCGCAGCTCTTCGCCGTTGTGGCTGTTCACTTTTTCTTCATAAGTACGGAAGTAGCCGTTCAGTGCCGGACGACCAAATTCGTTGTTAAACGCCGCGCCGCCCAGCGGGCCTTCGGTCATGATGTCCAGCGCGGTGACAATGCGCTCAGGCTTACCGAAATCTTCTTCCCACGGCTGT  >  NC_000913/2693936‑2694506

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