Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,276	G→A	A16A (GCG→GCA)	deoA →	thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,276	0	G	A	98.6%	203.4 / NA	72	A16A (GCG→GCA)	deoA	thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (35/34);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCG  >  NC_000913/4618029‑4618513                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              gtgcTTCCAGCCTGCTGGCAAGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTg                                                                                                                                                                                                                                            >  1:463604/3‑251 (MQ=255)              gctgGCAAGCCTGTTGAAAGCGCAGGGCCACGGTGACGGTACGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACACGAAATTATTAGTAAAAAACGTGATGGTAATGCACTGAGCGATgaagaa                                                                                                                                                                                                                                >  2:6262/1‑250 (MQ=255)                 ggCAAGCCTGTTGAAAGCGCTGGGCCTCGGTGACGTTAAGAGCGCCAGCAGCTACTAATTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAAtt                                                                                                                                                                                                                             <  1:624689/250‑1 (MQ=255)                    aaGCCTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGtt                                                                                                                                                                                                                         >  2:620202/1‑251 (MQ=255)                       ccTGTTGAAAGCGCTGGGCCACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCt                                                                                                                                                                                                                      >  1:166121/1‑251 (MQ=255)                                       ggccACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTcg                                                                                                                                                                                                      <  1:116612/248‑1 (MQ=255)                                          cACGGTGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTcgcg                                                                                                                                                                                                    >  1:309416/1‑250 (MQ=255)                                           aCGGTGACGATAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACACGAAGGCCGGCTAAGGCGTAGACGCATCCGACAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGATACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACATGATGGACATGCACTGAGAGATGAAGACATTAGATTCTTTATCCACGGCATAcgcgcc                                                                                                                                                                                                  >  2:681428/1‑249 (MQ=255)                                                     tAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTCGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGGGATGGTCATGCACTGAGCGCTGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTAtct                                                                                                                                                                                        >  2:546175/1‑251 (MQ=255)                                                      cagagCGCCAGCAGCTACTAAGAAAGATGCTTTACTCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCAGCCTCATACTCTTTTCCTCGGGCGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAAAGTGAGGGTCATGCACTGAGCGTTGAAGAAATTCGTTTCTTTATAAACGGCACTCGAGACAACACtctatc                                                                                                                                                                                       >  2:532936/2‑246 (MQ=255)                                                      aaGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTAtctc                                                                                                                                                                                       <  1:155270/251‑1 (MQ=255)                                                                 catcgtataagtaagATGCTTTACGCCTGATTCGCTGCGCTTATCAGG‑CCTACTGGACGTATTTACCCTTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGTCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGACGAAGAAATTCGTTTCTTTAGCAACGGCATTCGCGACAACACTATCGCCGAAGGGCAGa                                                                                                                                                                            <  2:643777/244‑1 (MQ=255)                                                                        taagtaagATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTgccgc                                                                                                                                                                     >  1:184547/1‑251 (MQ=255)                                                                         aagtaagATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTgccg                                                                                                                                                                      <  1:364265/249‑1 (MQ=255)                                                                            taagATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CATACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTc                                                                                                                                                                 >  1:178108/1‑251 (MQ=255)                                                                             aagATGCTTTACGCCTGATTCGTGGGGCTTATCAGG‑CCTTCGGGCCGTATCTTCCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGTCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTc                                                                                                                                                                 <  2:400001/250‑1 (MQ=255)                                                                                 tcctttaaggcTGATGCGCTGCGCTTAT‑AGGACCTAGGAAACGTATCTACCCGTCGGCCGGATAAGGCCTAGACGCATCAGGCAGAGGCCGCATCATACTCGTTTCATCGGTAGGGTCCCGTGGTTGTCGGACTAGAACTTAGTCGTAATAAACGTGATGGTCAAGCGCTGAGCGCGGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGGCAACAATGTCTCCGAAGGGCCGGTTGCAGCACTcgcaat                                                                                                                                                            >  2:713/11‑248 (MQ=255)                                                                                        cGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCAtt                                                                                                                                                     >  1:667498/1‑251 (MQ=255)                                                                                         gCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGCTGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGCCATTCGCGACAACACTATCTCCGCAGGGCAGATTGCCGCCCTCGCGATGACCAttt                                                                                                                                                    >  1:314432/1‑251 (MQ=255)                                                                                                         tatgctgg‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTATACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAAGTATTCGTAAAAAACGTGATGGTCATGCTCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACTATTTTCTTCCACGATATGa                                                                                                                                     <  2:416489/244‑1 (MQ=255)                                                                                                             cAGG‑CCTACGAGACGTATCTCCCAGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTCAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCACCGGAATTAGCGACAACACTATCTCCGAAGGGAAGATTGCCGCACTCGCGATGACCATTTTCTTCCCCGATATGACAATg                                                                                                                                >  2:680775/1‑251 (MQ=255)                                                                                                              aGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGc                                                                                                                               >  1:164723/1‑251 (MQ=255)                                                                                                                g‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCt                                                                                                                             <  1:6262/251‑1 (MQ=255)                                                                                                                      cGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACg                                                                                                                                                                                                             >  1:654172/1‑166 (MQ=255)                                                                                                                      cGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACg                                                                                                                                                                                                             <  2:654172/166‑1 (MQ=255)                                                                                                                      cGAGACGTATCTACCCGTAGGCCGGATAAGACGTCGTCCCATCCGGCCCACGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAACTTCGATTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGAGCAGATTGCCGCCCTCGCGTGGACCATTTTCTTCCACGATCTGACAATGCCTGGGCg                                                                                                                        >  1:602510/1‑251 (MQ=255)                                                                                                                       gagaCGTAACTACCCGGTGGCCGGATAAGGCGTAGACGCAGCCGGCAAAAGCCTCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCGCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCgt                                                                                                                       <  1:563166/251‑1 (MQ=255)                                                                                                                            gTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATCCTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCTCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGtctc                                                                                                                  <  1:365161/251‑1 (MQ=255)                                                                                                                                  aCCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGAc                                                                                                            <  1:24585/251‑1 (MQ=255)                                                                                                                                    ccGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGAcc                                                                                                           >  1:554177/1‑250 (MQ=255)                                                                                                                                      gTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATg                                                                                                        >  1:192080/1‑251 (MQ=255)                                                                                                                                             ggATTAGGCGTAGACTCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGc                                                                                                 <  2:597324/251‑1 (MQ=255)                                                                                                                                               aTAAGGCGTAGACGCATCCGGCAAAAGCAGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGAAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCg                                                                                                >  1:483595/1‑250 (MQ=255)                                                                                                                                                         gACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGa                                                                                      >  1:11353/1‑250 (MQ=255)                                                                                                                                                                    cAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTc                                                                                                                                                                 <  1:402146/164‑1 (MQ=255)                                                                                                                                                                    cAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGTTTCAGGACCCGTTCTCGAc                                                                          >  2:443023/1‑251 (MQ=255)                                                                                                                                                                         gccgccTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGaa                                                                     <  1:620202/251‑1 (MQ=255)                                                                                                                                                                           cgccTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCACCGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCTCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGaaa                                                                    >  2:547070/1‑250 (MQ=255)                                                                                                                                                                           cgccTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGaaa                                                                    >  1:600509/1‑250 (MQ=255)                                                                                                                                                                             ccTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGc                                                                                                      <  1:580515/214‑1 (MQ=255)                                                                                                                                                                             ccTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGc                                                                                                      >  2:580515/1‑214 (MQ=255)                                                                                                                                                                                   aCTCTTTTCATCGTGAGGTTTCCTTGTTTCTCGCACAAGAAATTAGTCGTAAAAAACGGGATTGTCAGGCACTGAGAGATGAAGAAATTCGTTTCTTTATCAACGGCAGTCGAGACAATACTATCTCCGAAGGGCAGGTTGCCGCCCTCGCGATGAACATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGATCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAGGCCTGCa                                                            <  2:597892/250‑1 (MQ=255)                                                                                                                                                                                           ccTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATg                                                    <  1:658722/250‑1 (MQ=255)                                                                                                                                                                                              cGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATTACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGccc                                                <  1:681428/251‑1 (MQ=255)                                                                                                                                                                                               gggAGGTTTCCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAAAGTGATGGTCATGCACTGAGCGATGAAGAAATTCTTTTCTTTATCAACGGCATGCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGGGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTGCTCGACTGGAAAAGCCTGCATCTGAATGGCCCg                                               <  2:45776/251‑1 (MQ=255)                                                                                                                                                                                                ggAGGTTACCTTGTTTCTCGCACAAGAAATTAGTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGTAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAATGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAATCTTTCTCGACTGGAAAAGCCTGCATCTGAATGGcc                                                 <  2:658054/248‑1 (MQ=255)                                                                                                                                                                                                     ttACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTcgcg                                                                                                                                                              >  2:375624/1‑134 (MQ=255)                                                                                                                                                                                                     ttACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTcgcg                                                                                                                                                              <  1:375624/134‑1 (MQ=255)                                                                                                                                                                                                       tccTTGTTGCCCGCACAAGAAAATATTCGTAAAAACCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCGCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGAttgttg                                        <  2:361568/249‑1 (MQ=255)                                                                                                                                                                                                        ccTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGACGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACAATATCTCCGAAGGGCAGATTGCCGACCTCGCGCTGACCATTTTCTTCCACGATCTGACAATGCCTGAGCGTGTCGCGATGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGa                                       >  2:673803/1‑250 (MQ=255)                                                                                                                                                                                                            gTTTCGCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTTATaaa                                   <  2:600509/250‑1 (MQ=255)                                                                                                                                                                                                              ttCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAaca                                 <  1:599107/250‑1 (MQ=255)                                                                                                                                                                                                                  cGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAacatt                               <  1:649906/248‑3 (MQ=255)                                                                                                                                                                                                                  cGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAacatt                               <  1:288046/248‑3 (MQ=255)                                                                                                                                                                                                                  cGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTcc                             <  2:309416/250‑1 (MQ=255)                                                                                                                                                                                                                    cacaAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTtctc                                                                             <  1:395056/200‑1 (MQ=255)                                                                                                                                                                                                                    cacaAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTtctc                                                                             >  2:395056/1‑200 (MQ=255)                                                                                                                                                                                                                     acaAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCActg                         >  1:346321/1‑249 (MQ=255)                                                                                                                                                                                                                     acaAGAAATTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATAGCCCGATTGGTGATAAAAATTCCActg                         >  2:465840/1‑249 (MQ=255)                                                                                                                                                                                                                         gAAATTATGCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTAGCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGCCCATGGCGATGCGAGATTCAGGAACCGTtc                                                                               <  1:70776/193‑1 (MQ=255)                                                                                                                                                                                                                           aaTTATTCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATGGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCg                   <  2:361341/251‑1 (MQ=255)                                                                                                                                                                                                                            attattCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGt                  <  1:581628/251‑1 (MQ=255)                                                                                                                                                                                                                               attCGTAACAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGTAGAGCAGATTGCCGCCCTCGCGATGCCCATTTTCTTCCCCGATATGACAATGCCTGAGCGTGTCTCGCTGACAATGGCGATGCGAGATTCAGGCACCGTTCTCGCCTGGCAAAGCCTGAATCTAGATGACCCGCTTGGTTATAAacattcacatggtggcgtagg               >  2:518031/1‑233 (MQ=255)                                                                                                                                                                                                                               attCGTAAAAAACGGGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCg                <  2:164723/250‑1 (MQ=255)                                                                                                                                                                                                                                ttCGTAACAAAAGTGAAGGTCATGCACTACGCGATGAAGAAATTCGTTTATTTATCAACGGCATTCGAGCCAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCATGTCGCGCTGACCATGGCGATACGAGATTCAGGAACGGTTCTAGACTTGAAAAGCATGCATGTGAATGGCACGCATGGTGATAAATATTCCACTGGTGGCGTCGGc              >  2:72791/1‑251 (MQ=255)                                                                                                                                                                                                                                ttCGTAAAAAACGTGATGGTCATGCCCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGAGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGGATGGCCCGCTTGTTGATAAACATTCCACTGGTGGCGTCGGc              >  1:492518/1‑251 (MQ=255)                                                                                                                                                                                                                                ttCGTAAAAAACGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGc              >  1:149719/1‑251 (MQ=255)                                                                                                                                                                                                                                 tCGTAAAAAACGGGATGGGCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGc              <  2:667498/250‑1 (MQ=255)                                                                                                                                                                                                                                      aaaaaCGTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGAtgtg        >  1:173055/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaaCGTTATGGTCATGCACTGAGCGATGAAGAAATTCGGTTCTTTATCAACGGCATTCGCGACAACTCTATCTCCGAAGGGCAGATTGCCGCCCTCGCGAGGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCACTGGTGGCGTCGGCGATGTGaa      <  2:466643/251‑2 (MQ=255)                                                                                                                                                                                                                                            gTGATGGTCATGCACTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGCATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACATTCCAATGGTGGCGTCGGCGATGTGACTTCg  >  2:567911/1‑251 (MQ=255)                                                                                                                                                                                                                                            gTGATGATCATGCACTGAGCGATGAAGCAATTCGTTTCTTTATCAACGGCATTCGAGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCCTTTTCTTCCACGATATGCCAATGCCTGAGCGTGTCTCGCTGCCCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGCTTGTTGATAAACATTCCACTGGGGGCGTCGTCGATGTGACTTCg  >  2:98009/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              GCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGG‑CCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTGAGCGTGTCTCGCTGACCATGGCGATGCGAGATTCAGGAACCGTTCTCGACTGGAAAAGCCTGCATCTGAATGGCCCGATTGTTGATAAACACTCCACCGGTGGCGTCGGCGATGTGACTTCG  >  NC_000913/4618029‑4618513

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