Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,292	G→C	R159R (CGC→CGG)	pyrE ←	orotate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,292	0	G	C	100.0%	211.9 / NA	62	R159R (CGC→CGG)	pyrE	orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (31/31);  total (31/31)

TGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTG  >  NC_000913/3815045‑3815535                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    tgatTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCCGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGt                                                                                                                                                                                                                                                  >  2:94954/1‑251 (MQ=255)  gatTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGt                                                                                                                                                                                                                                                  <  1:329448/250‑1 (MQ=255)   atTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTc                                                                                                                                                                                                                                                >  2:416608/1‑251 (MQ=255)        gACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATAGGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCGCGCCGCGCCCCCGTTCCTGAc                                                                                                                                                                                                                                           <  1:89646/251‑1 (MQ=255)           ccccccgggCCGCTAATGGCGGCACAGTTGCATTAATATGACGCCGGATGACTTTTCATCCGGCTAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAAGCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAAGCACGCTCAACTTCCTGAATCGCCGAAAGCTCGCCGCGCCCCCGTTCCTGACGAt                                                                                                                                                                                                                                        <  2:408637/244‑1 (MQ=255)               gCTGGCCGTTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACCCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATcaag                                                                                                                                                                                                                                    <  1:434363/251‑4 (MQ=255)                  ggCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGc                                                                                                                                                                                                                                   <  2:223301/249‑1 (MQ=255)                               ggCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCa                                                                                                                                                                                                                    <  2:324606/251‑1 (MQ=255)                                 cacaGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGc                                                                                                                                                                                                                  >  2:95575/1‑251 (MQ=255)                                  acaGTTGCAGTAATATGACGCCGGATGCCTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCa                                                                                                                                                                                                                 <  1:342714/251‑1 (MQ=255)                                           gTAATATGGCGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAATCGCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTcgcg                                                                                                                                                                                                        <  1:174601/251‑1 (MQ=255)                                                       cGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTg                                                                                                                                                                                             >  1:429248/1‑250 (MQ=255)                                                        ggATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATca                                                                                                                                                                                                                                      <  2:338462/208‑2 (MQ=255)                                                        ggATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATca                                                                                                                                                                                                                                      >  1:338462/1‑207 (MQ=255)                                                                  tCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCATTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCAt                                                                                                                                                                                 >  2:377257/1‑251 (MQ=255)                                                                      ccggctatttgcTTTAAACGCCAACCTCGTCGCGATAGGCTTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCTTCGCGCCATTGGCCTGAAAAATCTCCATCGAc                                                                                                                                                                             <  2:28090/240‑1 (MQ=255)                                                                       cGGCGAGTTTCTTTAAACGCAAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCGTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACt                                                                                                                                                                            <  1:70307/251‑1 (MQ=255)                                                                        ggCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTc                                                                                                                                                                           >  1:181226/1‑251 (MQ=255)                                                                         gCGAGTTTCTTGAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTcg                                                                                                                                                                          <  1:18341/251‑1 (MQ=255)                                                                            aGTTTCTTTACACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTCATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTcgcg                                                                                                                                                                        >  2:419750/1‑250 (MQ=255)                                                                              tttCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGCCTCGCGAAt                                                                                                                                                                     >  2:375656/1‑251 (MQ=255)                                                                                tCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGACTAATCTCCATCGACTCGCGAATCg                                                                                                                                                                   >  2:332958/1‑251 (MQ=255)                                                                                       aCGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTcc                                                                                                                                                            <  2:149585/251‑1 (MQ=255)                                                                                               ctctTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCcggcggggg                                                                                                                                                    >  2:47352/1‑248 (MQ=255)                                                                                                      cgATAGACCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGTTAAcatc                                                                                                                                             >  2:291167/1‑251 (MQ=255)                                                                                                       gATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCCGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGCTCAAGCGAAATCAATACGCCAGCAAGCTACGCGCCATTGGCCTGAATAATCTCCATCCACTCGCGCATCGCCGTTCCGGCGGCGGTCAcatcc                                                                                                                                            >  2:113307/1‑250 (MQ=255)                                                                                                          aGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGTTCACATCATCt                                                                                                                                         >  1:95839/1‑251 (MQ=255)                                                                                                           ggCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGCTCACATCctct                                                                                                                                         >  2:363140/1‑250 (MQ=255)                                                                                                                tAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAg                                                                                                                                    <  2:341093/250‑1 (MQ=255)                                                                                                                  aCCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCAt                                                                                                                                 <  1:37471/251‑1 (MQ=255)                                                                                                                                           ccGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTCCCAGCATTACGCGTCCTTGTAACGCGCTaccc                                                                                                        >  1:256790/1‑250 (MQ=255)                                                                                                                                           ccGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAAGCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTacc                                                                                                         <  1:285501/250‑1 (MQ=255)                                                                                                                                             ggCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTaccaac                                                                                                      <  1:416608/251‑1 (MQ=255)                                                                                                                                               cTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATc                                                                                                                                                                                                                                       <  1:167222/120‑1 (MQ=255)                                                                                                                                               cTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATc                                                                                                                                                                                                                                       >  2:167222/1‑120 (MQ=255)                                                                                                                                                            tAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGaata                                                                                                                                                                                         >  2:152189/1‑153 (MQ=255)                                                                                                                                                            tAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGaata                                                                                                                                                                                         <  1:152189/153‑1 (MQ=255)                                                                                                                                                             aaGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTcc                                                                                      >  2:118146/1‑250 (MQ=255)                                                                                                                                                                    cAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGACATTGACCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCACCCAGATTGACGCCTTCACCGTGGt                                                                               >  1:165499/1‑251 (MQ=255)                                                                                                                                                                    cAGGTCTTTCAGGGTGATGATAGAGATCACGTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGt                                                                               <  1:375656/251‑1 (MQ=255)                                                                                                                                                                      ggTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGCGCTCACATCATCTACCAGCCTTACGCGTCCTTGTAACGCGCCTCCACCCAGATTGCCGCCTTCACCGTGGTCt                                                                             >  2:220175/1‑251 (MQ=255)                                                                                                                                                                             cAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGctt                                                                      >  2:159503/1‑251 (MQ=255)                                                                                                                                                                              aGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAAGCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCAGCTACCAGCATTACGCGTCCTTGTAGCGCGCTACCAACCAGATTTCCGCCTTCTCCCTGTTCTTTTGcttc                                                                     <  2:60905/251‑1 (MQ=255)                                                                                                                                                                              aGGGTGATGATAGAGATCACTTTGCAGGTGTAATCACGCGCAACTTCCTGAATCGCCGAAATCTCGCCGCGTCCCCGTGCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAAGCGCCGTTCCGGCGGTGATCACATCTTCTACCAGCATTACGCGTCCTTGTAACGCGCTGCCAACCAGATTGCCGCCTTCACCGTGGTCTTCTGcttc                                                                     <  2:167443/251‑1 (MQ=255)                                                                                                                                                                                           agaTCAATTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGCGGTCTTTTGCTTCTTTGCGGTTaa                                                          >  2:280198/1‑249 (MQ=255)                                                                                                                                                                                                  tttGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACg                                                 >  1:136267/1‑251 (MQ=255)                                                                                                                                                                                                   ttGCATTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTTACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTGGCTTCTTTGCGTTTAAAGCAGTACg                                                 <  1:173224/250‑1 (MQ=255)                                                                                                                                                                                                   ttGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACg                                                 <  2:95839/250‑1 (MQ=255)                                                                                                                                                                                                   ttGCAGATGTAATCAAGCTCAACTTCCTGAATCGCCGAAATCTCGCAGCGCCCCCGTTCATGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCATTGTTACGCGCTACCAACAAAATTGCCGCCTTCACAGTGGGCTTTTGTTTTTTTGCGGGTAAAGAAGTaagg                                                >  2:175760/1‑248 (MQ=255)                                                                                                                                                                                                       aGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGcagg                                            <  1:44011/251‑1 (MQ=255)                                                                                                                                                                                                           gTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGcaggtcca                                        <  2:161649/251‑1 (MQ=255)                                                                                                                                                                                                           gTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGcaggtcca                                        <  1:278269/251‑1 (MQ=255)                                                                                                                                                                                                            tAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCCTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGcaggttctg                                       >  1:239176/1‑247 (MQ=255)                                                                                                                                                                                                            tAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACAACATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGGCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGccgggtccg                                       >  2:129500/1‑243 (MQ=255)                                                                                                                                                                                                               tCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGccggtccagggc                                    >  1:271684/1‑249 (MQ=255)                                                                                                                                                                                                                  cGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTg                                 <  1:332958/251‑1 (MQ=255)                                                                                                                                                                                                                       aaCTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCCTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTAAGGcagggccatgatcgggttgc                            >  2:366916/1‑239 (MQ=255)                                                                                                                                                                                                                          ttCCTGACTCGCCGCAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGAACGGcagggcccgggcggtgtgctccc                         >  2:61709/1‑232 (MQ=255)                                                                                                                                                                                                                           tCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTcc                          <  1:47352/249‑1 (MQ=255)                                                                                                                                                                                                                                      cGAAATCTCGCCGCGCCCCCGTTCCTGTCGATCAAGAGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCAAATCATCTACCAGGATTACGCGTCCTTGTAACGCACTACCAACCAGATTGCCGCCTTCACCGGGGGCTTTAGCTTCTTTTCGGTTAAAGAACTACGGCAGGCTCAGGGCGTTATACTCCGCCatgcaccccg             >  1:237039/1‑242 (MQ=255)                                                                                                                                                                                                                                              cgccgcGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTtgtg     <  1:361218/251‑1 (MQ=255)                                                                                                                                                                                                                                                 cgcgCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGGCGCGCCATTGGCCGGAATAATCTCCACCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGGCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCCCCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTgtggtg  <  1:113307/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    TGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTG  >  NC_000913/3815045‑3815535

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