Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,571	A→G	E323E (GAA→GAG)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,571	0	A	G	96.3%	254.9 / ‑3.0	83	E323E (GAA→GAG)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (35/44);  total (36/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAA  >  NC_000913/4620331‑4620813                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   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                                                                                                                                                                     >  2:564375/1‑134 (MQ=255)                                                                                                                                                                                            cGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCTGCGGGTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGtt                                                                                                                                                                     <  1:564375/134‑1 (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)                                                                                                                                                                                                                                       gcgGGTCTGGAGCTTTTCGACCGCCGTCTGCCGTATCTGATGTCTCTGCTGCGCGATGACGATGTCCTGATCCTCACTGCTTACCACGGTTGTGATCCGTCCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGGGAAACCTTCGCGGATATTGGTCAGACGCGGGCAAAATATTGTGGTACCGCTGATATGGAATATgg     <  2:329055/251‑1 (MQ=255)                                                                                                                                                                                                                                           gTCTGGAGCTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGACAAGGAATATGGCaa  >  2:170604/1‑250 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   GCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAAGCGACTGGCCTGGACGCGCTGTTTGACGCCACCATCAAAGAGATGAAAGAAGCG‑GGTGATAACACCATCGTCTTCACCAACTTCGTTGACTTCGACTCTTCCTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAA  >  NC_000913/4620331‑4620813

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: ‑