Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,559	C→T	A319A (GCC→GCT)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,559	0	C	T	95.2%	252.8 / ‑2.4	86	A319A (GCC→GCT)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (34/46);  total (37/49)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATAT  >  NC_000913/4620317‑4620784                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  gCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCt                                                                                                                                                                                                                           <  1:399864/251‑1 (MQ=255)   accgaccgTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCggg                                                                                                                                                                                                                              >  1:189881/1‑246 (MQ=255)   accgaccgTGCTGCAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTgg                                                                                                                                                                                                                         <  1:251500/251‑1 (MQ=255)   accgaccgTGCTGCAGAAACTGGGTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCggg                                                                                                                                                                                                                              <  2:189881/246‑1 (MQ=255)         cgTGCTGCAGCAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGAGATACCACCATCGTCTTCACCAACTTCGTTGACTACGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGagcttt                                                                                                                                                                                                                   >  2:505988/1‑246 (MQ=255)               gcAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAc                                                                                                                                                                                                              <  1:440735/250‑1 (MQ=255)               gcAGAAACTGGTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGAc                                                                                                                                                                                                              <  1:256456/250‑1 (MQ=255)                        ggTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTg                                                                                                                                                                                                     <  1:277375/250‑1 (MQ=255)                         gTTGATGAAAAACACGGTCAGGTGGTTTCTGTCGGTAAAAGTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGTCCGCCGTCTGcc                                                                                                                                                                                                   <  2:90776/251‑1 (MQ=255)                            gatgaAAAACACGGACAGGTGGATTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGTTAACACCATCGTCTTCACCAACTTCGTTGACTTCGACACTTCCTGGGGCCACCGTCGCTCCGTCGCCGGTTATGCTGCGGGCCTGGAGTTGTACCACCGCCGCCTGCCggg                                                                                                                                                                                                >  2:276264/1‑250 (MQ=255)                                aaaaaCACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTg                                                                                                                                                                                            >  1:408635/1‑251 (MQ=255)                                 aaaaCACGGTCAGGTGGTTTCTGTCGGTAAAAGTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCT‑GGTGATAACACCATCGTCTTCACCAACGTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGa                                                                                                                                                                                           <  1:98925/251‑1 (MQ=255)                                  aaaCACGGTCAGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGACAGTGAAAGCAACTGGCCTGGCCGAGCTGTTTGACGCCACCAACAAAGAGATGAAAGAAGCG‑GGTGTTAACACCATCGTCTTCCCCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCAGCCGCCGTCTGCCGGAGCTGAt                                                                                                                                                                                          >  2:601797/1‑251 (MQ=255)                                            aGGTGGTTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGGGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGCGCTGATGTCTctgctg                                                                                                                                                                                >  1:373386/1‑251 (MQ=255)                                             ggaggtTACTGTCGGTAAAATTGCAGACATCTACGACAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCGGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTTCGGGTCTGGAGCTGGTCGACCGAAGTCTGCCCGTGGTGATGtctttacttc                                                                                                                                                                               >  2:571722/4‑244 (MQ=255)                                              gtggtTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTgcg                                                                                                                                                                              <  1:675156/251‑1 (MQ=255)                                              gtggtTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCGGCCGGAGCTGATGTCTCTGCTgcg                                                                                                                                                                              >  1:329055/1‑251 (MQ=255)                                               tggtTTCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTgcgc                                                                                                                                                                             >  1:519052/1‑251 (MQ=255)                                                   ttCTGTAGGTAAAATTGCAGACATATACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGCCTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGCTCTGGAGCTGTTCGACCGCAGTCTGACGGAGCTGATGGCTCTGCTGCGCGATg                                                                                                                                                                         >  2:292990/1‑251 (MQ=255)                                                    tCTGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATg                                                                                                                                                                         <  1:311708/250‑1 (MQ=255)                                                      tGTCGGTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGCGCTGATGTCTCTGCTGCGCGATgacg                                                                                                                                                                      >  2:240837/1‑251 (MQ=255)                                                          ggTAAAATTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACAt                                                                                                                                                                  <  1:9619/251‑1 (MQ=255)                                                              aatttGCGGACATCTACGTCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCCAAGAGATGAAAGAAGCG‑GGTGATAACATCATCGTCTTCACCAACTCCGTTGACTTCGTCTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACGGTTCGACCGCCGTCTGCCGGAGCTGATGGCTCGGCCGCGCGATGACAGCATCCTg                                                                                                                                                              >  2:290173/4‑251 (MQ=255)                                                                aTTGCAGACATCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGAt                                                                                                                                                            >  1:342506/1‑251 (MQ=255)                                                                          tCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCt                                                                                                                                                  >  1:548329/1‑251 (MQ=255)                                                                          tCTACGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGGTTGACGCCACCATAAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTAGTTGACTTCGACTCTTCCTGGGGCCACTGTCGCGACGTCGCCGGTTATGCGGCGTGTCTGGAGCTGTTCGACCGCCTTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGTTCATCACTGCt                                                                                                                                                  >  2:535314/1‑251 (MQ=255)                                                                             aCGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCActgc                                                                                                                                                   <  1:678706/247‑4 (MQ=255)                                                                             aCGCCAACTGCGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCt                                                                                                                                                  <  1:398506/248‑1 (MQ=255)                                                                             aCGCCAACTGCGGTAGCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGAc                                                                                                                                               <  1:635946/251‑1 (MQ=255)                                                                             aCGCCAACGGCGGTATCGACAAGAAATAGAAAGAAACGGGCCTGGACGCGCTGTGTGACGCCACCAACAAAGAGATGAAAGAAGCG‑GGTGATAACCCCATCGTCTTCACCAACTTCGTTGACTTCGACGCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGTTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACAGCCTGATCCTctctgctg                                                                                                                                                 <  2:267888/249‑8 (MQ=255)                                                                                     cgCGGTATCACCAAGAAAGTGAAAGCCACTGGCCTGGACGCGCTGTTTGACTCCACCAGAAAAGAGATGAAAGAAGCG‑GTTGATACCACCATCGTCTTCACCATCGTCGTTGACTTCGACTCTTCCTGGGGCCCCCGTCGCGACGTCCCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTg                                                                                                                                       <  2:557871/250‑1 (MQ=255)                                                                                       cGGTATCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATAGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCAGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTgtg                                                                                                                                     >  1:402279/1‑249 (MQ=255)                                                                                            tCACCAAGAAAGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTc                                                                                                                                                                                                       >  1:482204/1‑180 (MQ=255)                                                                                              aCCAAGAAATTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGAc                                                                                                                              <  1:160310/251‑1 (MQ=255)                                                                                                    taagagAAAGCAACTGGCCTGGACGCGCTGTTTGACGCAACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACAGTCGCGACGTCGACGGGTTTGCTGCGGGTCCGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGATGCGCGATGACGACCTCCAGATCCACACTGCTAACCACGGTTgtgttcggccttgtgc                                                                                                                        >  2:157125/6‑236 (MQ=255)                                                                                                    aaaGTGAAAGCAACTGGCCTGGACGCGCTGTTTTACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGATTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATga                                                                                                                                                                        <  2:48107/203‑1 (MQ=255)                                                                                                    aaaGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATga                                                                                                                                                                        >  1:48107/1‑203 (MQ=255)                                                                                                      agggAAAGCAACTGGCCTGGAAGCGCTGGTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAAAACCAGCGTCTTCACCAACTTCGTTGATTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCGCTGCTGCGCGATGACGACATCATGAt                                                                                                                                                            <  2:52942/210‑1 (MQ=255)                                                                                                      aGTGAAAGCAACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGAt                                                                                                                                                            >  1:52942/1‑213 (MQ=255)                                                                                                              caaCTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGAc                                                                                                              <  2:345111/249‑1 (MQ=255)                                                                                                               aaCTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGTTACTGAcc                                                                                                             >  1:540559/2‑251 (MQ=255)                                                                                                                 cggtcctggaagttcTGTTCGACGCCACCATCAAAGAGATGAAAAAAGCG‑GGTGATAACATCAACGTCTTCACCAACTTCGTTGCCTTCGCCTCTTCCTGGGGCCACCTTCGCGACGTCGCCGTTTAGGTTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATTACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGAcctc                                                                                                           <  2:537280/237‑3 (MQ=255)                                                                                                                 cTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGCGATGAACGAAGCG‑GGTGATAACACCATCGGCTTCACCAACTTCGTTGACTTCGACTCTTCATGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGATCTGGAGCGGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGGGATcagaacagcacagggactgacccc                                                                                                           >  2:310731/1‑228 (MQ=255)                                                                                                                    gCCTGGACGCGCTGTTGGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACAcg                                                                                                        <  1:611779/251‑1 (MQ=255)                                                                                                                     ccTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGGCCtagatcggacgagcgtcgggtt                                                                                                       >  2:209451/1‑230 (MQ=255)                                                                                                                     ccTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCt                                                                                                                            <  1:209451/230‑1 (MQ=255)                                                                                                                           cgcgCTGTTTGACACCACAATCAAAGAGATGAAAGAAGAG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGCGCTGTTCGACCGCCGACTGCAGGAGCTGATGTCTCTGATGCGCGATGACGACATCATGCTCCTCAATGCAGACAACGGTTGTGATCCGACCTGGACCGGCAATGCCCACACGCGTAAac                                                                                                 >  2:619381/1‑251 (MQ=255)                                                                                                                              gcggattGACGCCCCCAGCAAAGAGATGAACGCAGCG‑GGTGATAACACCATCGTATTAACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCCACGTCCCCGGTTCGGCGGCAGGTCTGGAGCTGTTCGCCCGCCATCTACCGGAGCTGCTGTCTCTGCTGCCCGATTACGACATACTGATCCTCTCTGCTGCCCACGGTTGTGATCCGACCTCGACAGGTAATAACCACACACGttcacac                                                                                               >  2:315180/6‑244 (MQ=255)                                                                                                                                   ttGACGCCCCCATCAAAGAGATGAAGGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCTGCTCTTCCTGGGGCCACCGTCGCGACGGCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGGCCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTcc                                                                                          <  2:591204/250‑1 (MQ=255)                                                                                                                                   ttGACGCCACCAACAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACGAACTTCGTTGACTTCGACCCTGCCTGGGGCCACCGTCGCGATGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTGCTACCGCCGTCTGCCGGAGCTGATGTCTCTGCGGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATTCGACCTGGACCGGTACTGACCACACGCGTGAACTCATTcc                                                                                          <  2:465348/250‑1 (MQ=255)                                                                                                                                     gaagCCACCATCAAAGAGATGAAAGAAGTGTGGT‑ATAACCACATCGTCTTCACCAACTTCGTTGACTCCGCCGCTTCCTTGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGAGGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGt                                                                                       <  2:539849/248‑1 (MQ=255)                                                                                                                                       ccccaccaTCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCATCACCAACTTCGTTGACTTCGTCTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTTTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCGCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGGTCCGACCTGGCACGGTACTGACCACGCGCGTGAACACATTCCTGTAc                                                                                     <  2:478383/249‑1 (MQ=255)                                                                                                                                       cGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTAc                                                                                     >  1:159945/1‑251 (MQ=255)                                                                                                                                          caccaTCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTgg                                                                                  <  1:261468/251‑1 (MQ=255)                                                                                                                                           accaTCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGt                                                                                 <  2:519052/251‑1 (MQ=255)                                                                                                                                                cAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGtttacg                                                                            <  1:558897/251‑6 (MQ=255)                                                                                                                                                  aaGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGtttacggc                                                                          <  1:457009/251‑8 (MQ=255)                                                                                                                                                    gagaTGAAAGAAGCG‑GGAGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGCCGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGccc                                                                        >  2:180261/1‑251 (MQ=255)                                                                                                                                                         gAAAGAAGCG‑GGTGATAACACCAGCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACTGTCGCGACGCCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACAGACGTCTGCCGGAGCTGATGGCTCTGCTGCGCGATGACGACATCCTGATCCTAACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAg                                                                   <  1:657202/251‑1 (MQ=255)                                                                                                                                                          aaaGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGCAAGt                                                                  >  1:158440/1‑251 (MQ=255)                                                                                                                                                               agCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAAc                                                             <  1:293571/251‑1 (MQ=255)                                                                                                                                                                         tAACACCATCGTCTTCACCCACTTCGTTTACTTCGGCTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGCCCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCATGATCCTCACTGCTGACCACGGTTGTGATCTGACCTGGACCGGTACTGAACACACGCGTGAACACATTCAGGTACTGGTATACGGCCCGGAAGTAACCCCGGGCTCAc                                                    >  2:392731/1‑251 (MQ=255)                                                                                                                                                                          aaCACCATCGTCTTCACCAACTTCGTTGACTTCGACTCGTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACt                                                   <  2:315824/251‑1 (MQ=255)                                                                                                                                                                            cacCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGAt                                                                                                                                                            >  2:568826/1‑144 (MQ=255)                                                                                                                                                                            cacCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGAt                                                                                                                                                            <  1:568826/144‑1 (MQ=255)                                                                                                                                                                                   gTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCaccg                                          <  1:167342/251‑4 (MQ=255)                                                                                                                                                                                        cACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCCCGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGACCTGGTTTACGGCCCGAAAGTaaaccc                                                            >  2:146055/1‑225 (MQ=255)                                                                                                                                                                                        cACCAACTTCGGTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAAcc                                                            <  1:146055/228‑1 (MQ=255)                                                                                                                                                                                         aCCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAAc                                    >  1:511092/1‑251 (MQ=255)                                                                                                                                                                                             actgcGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTc                                <  1:292990/247‑1 (MQ=255)                                                                                                                                                                                                tatttGACTTGGACACGTCCTGGGTCCACCGTCGCGACGTCGCCGGTTAGGCCGCGGGTCTGGGGCTGGTCGTCCGCCGTCTGCCTGAGCTGATGGCTCTGCTGCGCGATGACGACATCCTGATCCTGACTGCTGACCACTGTTGTGATCCGACCGGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTcg                               <  2:648744/246‑1 (MQ=255)                                                                                                                                                                                                tCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTcgc                              <  1:371858/250‑1 (MQ=255)                                                                                                                                                                                                tCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTcg                               <  1:240837/249‑1 (MQ=255)                                                                                                                                                                                                 cGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTcgcg                             <  1:574922/250‑1 (MQ=255)                                                                                                                                                                                                  gTTGACTTCGCCTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGa                           <  1:245407/251‑1 (MQ=255)                                                                                                                                                                                                    tGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGat                          >  1:291151/1‑250 (MQ=255)                                                                                                                                                                                                       cTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGCTCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAATACCGGGCTCACTGGGTCCCCGGGAAACCTTCCCGGattttt                      >  2:575972/1‑247 (MQ=255)                                                                                                                                                                                                          cGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGtt                                                                                                                                        <  1:564375/134‑1 (MQ=255)                                                                                                                                                                                                          cGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGtt                                                                                                                                        >  2:564375/1‑134 (MQ=255)                                                                                                                                                                                                          cGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGatattggt                    <  1:507597/250‑5 (MQ=255)                                                                                                                                                                                                          cGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGcc                   <  1:290173/251‑1 (MQ=255)                                                                                                                                                                                                          cGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATACTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGcc                   <  1:332042/251‑1 (MQ=255)                                                                                                                                                                                                               ctTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGAct              >  1:652716/1‑251 (MQ=255)                                                                                                                                                                                                                  ccgggggCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACTCATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTg           <  2:342506/248‑1 (MQ=255)                                                                                                                                                                                                                          cACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGATCTGGTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAata   <  2:386476/251‑1 (MQ=255)                                                                                                                                                                                                                           aCCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGCCCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCCGCCTCTGGCAACatat  >  2:86529/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  GCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATAT  >  NC_000913/4620317‑4620784

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