Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	700,038	G→A	H97Y (CAT→TAT)	umpH ←	UMP phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	700,038	0	G	A	100.0%	239.0 / NA	66	H97Y (CAT→TAT)	umpH	UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/35);  total (31/35)

GGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTACGGCGACGTTATCGTG  >  NC_000913/699796‑700281                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    gggctTAACAACTTAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTTCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGt                                                                                                                                                                                                                                              <  1:131794/250‑1 (MQ=255)   gctTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGGAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTgcg                                                                                                                                                                                                                                           <  2:55463/251‑1 (MQ=255)       aCCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCtt                                                                                                                                                                                                                                       >  1:545713/1‑251 (MQ=255)       aCCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCATGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCCTTGTACAGTTCATAAATCAGTGCGCCtt                                                                                                                                                                                                                                       >  1:367237/1‑251 (MQ=255)                     gTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATAc                                                                                                                                                                                                                         >  1:42546/1‑251 (MQ=255)                             ccGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCtt                                                                                                                                                                                                                 >  1:380735/1‑251 (MQ=255)                             ccGGAGATTTTCTCAACCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCTGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGTAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCTCCACATACGCTTTCtt                                                                                                                                                                                                                 <  2:262555/251‑1 (MQ=255)                                          aaaccccTGCACACAACTCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATAATGTCCCAGTTGTAGGAACGCTTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGa                                                                                                                                                                                                    <  2:143657/247‑1 (MQ=255)                                                  gCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGa                                                                                                                                                                                             <  2:130860/250‑1 (MQ=255)                                                       caACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCgg                                                                                                                                                                                       <  1:140574/251‑1 (MQ=255)                                                            cgcCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTc                                                                                                                                                                                  <  1:399840/251‑1 (MQ=255)                                                                  aaGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATc                                                                                                                                                                            >  2:161721/1‑251 (MQ=255)                                                                  aaGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAAGCGGCAGTCGCCATc                                                                                                                                                                            <  1:71832/251‑1 (MQ=255)                                                                   aGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCg                                                                                                                                                                           >  2:432201/1‑251 (MQ=255)                                                                     cGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCa                                                                                                                                                                         <  1:25308/251‑1 (MQ=255)                                                                           aaaaaCCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGt                                                                                                                                                                    <  1:456251/250‑1 (MQ=255)                                                                                     cccggtggTGTCCGGATTGGTGTCGATAAAACGTTCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTTTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGc                                                                                                                                                         <  2:421238/245‑1 (MQ=255)                                                                                       cGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCAACGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTg                                                                                                                                                       >  1:479932/1‑251 (MQ=255)                                                                                         tGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTc                                                                                                                                                     >  1:432139/1‑251 (MQ=255)                                                                                                ccGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACa                                                                                                                                              <  2:458086/251‑1 (MQ=255)                                                                                                ccGGATGGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTAc                                                                                                                                               <  2:264966/250‑1 (MQ=255)                                                                                                      tggtggAGATAAAACGTGCACCGTTAGCGACGGAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGac                                                                                                                                         >  2:434895/1‑250 (MQ=255)                                                                                                                aaaaCGTGCCCCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGGCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTgg                                                                                                                              <  2:545713/251‑1 (MQ=255)                                                                                                                aaaaCGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTgg                                                                                                                              <  1:137456/251‑1 (MQ=255)                                                                                                                        cACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCgg                                                                                                                      >  1:467039/1‑251 (MQ=255)                                                                                                                                gCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCg                                                                                                                                                                                                     <  2:316781/164‑1 (MQ=255)                                                                                                                                gCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCg                                                                                                                                                                                                     >  1:316781/1‑164 (MQ=255)                                                                                                                                gCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACAACGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAg                                                                                                              <  2:326991/251‑1 (MQ=255)                                                                                                                                  gacgaAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGAt                                                                                                            <  2:232732/251‑1 (MQ=255)                                                                                                                                   acgaAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGcc                                                                                                                >  2:206456/1‑246 (MQ=255)                                                                                                                                     gaAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCtt                                                                                                         >  1:168661/1‑251 (MQ=255)                                                                                                                                           ggCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGccc                                                                                                     >  2:39043/1‑249 (MQ=255)                                                                                                                                            gctgctTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGt                                                                                                  <  1:432201/251‑1 (MQ=255)                                                                                                                                               gctTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCATTCTg                                                                                               <  1:141805/251‑1 (MQ=255)                                                                                                                                                  ttATGCATTATGTCCCCGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCCTGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTAGCCCAGTCTGc                                                                                              <  2:475052/249‑1 (MQ=255)                                                                                                                                                    aTGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCCGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCt                                                                                                >  2:70006/1‑245 (MQ=255)                                                                                                                                                    aTGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCt                                                                                                <  1:70006/245‑1 (MQ=255)                                                                                                                                                        atcatGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGAt                                                                                       >  1:346772/1‑250 (MQ=255)                                                                                                                                                           atGTCCCAGGTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGtt                                                                                   <  2:93789/251‑1 (MQ=255)                                                                                                                                                            tGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTg                                                                                  >  2:198092/1‑251 (MQ=255)                                                                                                                                                               cccAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGGCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTc                                                                               <  1:519716/251‑1 (MQ=255)                                                                                                                                                                 cAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTc                                                                               >  1:329163/1‑249 (MQ=255)                                                                                                                                                                             cgcgTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCgg                                                                 <  1:562828/251‑1 (MQ=255)                                                                                                                                                                               cgTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCCGTCTGCGAAGGATAGTTGGTCAGCACCACCAGCGGc                                                                >  2:182837/1‑250 (MQ=255)                                                                                                                                                                                    cGCCAACAATCACGAAATCAGGGTTCACAGCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCt                                                          <  2:464240/251‑1 (MQ=255)                                                                                                                                                                                       caacaaTCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCg                                                                                             >  2:99909/1‑213 (MQ=255)                                                                                                                                                                                       caacaaTCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCg                                                                                             <  1:99909/213‑1 (MQ=255)                                                                                                                                                                                              cACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATaa                                                <  1:28079/251‑1 (MQ=255)                                                                                                                                                                                                      cAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCCGCGGCAGGCCTTTATCCATAATCCCGTGc                                        >  2:377037/1‑251 (MQ=255)                                                                                                                                                                                                       aGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCa                                       <  1:198840/251‑1 (MQ=255)                                                                                                                                                                                                       aGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGATTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATATTGCCCAGTCTGCGAGGGGTAGTTTGTCAGCAACACCAGCGGCAGGCCTTTATCCTTAATccctggc                                        >  2:437876/1‑246 (MQ=255)                                                                                                                                                                                                        gggatCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGGTCTTGCCCAGTCTGCGAAGGATAGATGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCa                                       >  1:422598/5‑250 (MQ=255)                                                                                                                                                                                                             cacaTCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAAt                                  <  1:561430/250‑1 (MQ=255)                                                                                                                                                                                                             cacaTCGGTAATAGTGAAACCGGCTTTGGACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCGTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCGTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAAt                                  <  2:42546/250‑1 (MQ=255)                                                                                                                                                                                                                aTCGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCCGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAACCACGCTGTCAGGTACATCGACACCTGCGTGGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCACCACCAGCGGCCGGCCTTTCTCCAAAATCCCGTGCAAAAATTccc                              >  2:167987/1‑250 (MQ=255)                                                                                                                                                                                                                  cGGTAATAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGTAGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCAATAATCCCGTGCAAAAATTCCGCt                            >  1:344142/1‑251 (MQ=255)                                                                                                                                                                                                                        tAGTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAACTTCAGCTGCAccc                      >  2:153840/1‑251 (MQ=255)                                                                                                                                                                                                                          gTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCCGGCCTTTATCCATAATCCCGTGCAACAATTCCGCTGCACCCgg                    >  2:343768/1‑251 (MQ=255)                                                                                                                                                                                                                          gTGAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAA‑TTCCGCTGCACCCggg                   >  2:577795/1‑250 (MQ=255)                                                                                                                                                                                                                            gAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGt                   <  1:161721/250‑1 (MQ=255)                                                                                                                                                                                                                            gAAACCGGCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGt                   <  1:343768/250‑1 (MQ=255)                                                                                                                                                                                                                                  ggCTTTGTACAGTTCATAAATCAGAGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCACCACCAGCGGCAGGCCTTTATCAATAATCCCGTGCAAAAATTCCGCTGCCCCCGGAACGGctc            >  2:321774/1‑249 (MQ=255)                                                                                                                                                                                                                                   gCTTTGTACAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATTTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATTCATAATCCCGTGAAAAAATTCCGCTACACCCGGTACGGCGAc           >  2:314400/1‑251 (MQ=255)                                                                                                                                                                                                                                         caCAGTTCATAACTCAGTGCGCCTTCGTCCTCCACATACCCTTTCTTGCCCTCCTGTCGACGCAGGAAATCGGCAGTCGCAATCGCAGAGGTATAAAACACGCTGACAGGTACATCGACACCTGCGGTTGCAAAGCGGTTCGCCAGATCTTGCCTAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCTGTACGGCGTCGTTATc     <  2:126655/250‑1 (MQ=255)                                                                                                                                                                                                                                          aCAGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTACGGCGACGTTATCg    <  1:39043/251‑1 (MQ=255)                                                                                                                                                                                                                                            aGTTCATAAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAACGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAAGTCCGCAGACCCCGGTACGGCGGCGATATCgcg  >  2:354888/1‑249 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    GGGCTTACCAACATAGAACGGTTTGCGCCCGGAGATTTTCTCAATCCCTGCACACAACGCGCCACAAGCGGGATAAAAACCGCGCCCGTGGGTGTCCGGATTGGTGGCGATAAAACGTGCACCGTTAGCGACGAAATAGGCTGCTTTATGCATCATGTCCCAGTTGTAGGAACGCGTTTCGCCAACAATCACGAAATCAGGGTTCACATCGGTAATAGTGAAACCGGCTTTGTACAGTTCATGAATCAGTGCGCCTTCGCCCACCACATACGCTTTCTTGCCTTCCTGGCGACGCAGGAAATCGGCAGTCGCCATCGCAGAGGTATAAAACACGCTGTCAGGTACATCGACACCTGCGGTGGCAAAGCGGTTCGCCAGATCTTGCCCAGTCTGCGAAGGATAGTTGGTCAGCAACACCAGCGGCAGGCCTTTATCCATAATCCCGTGCAAAAATTCCGCTGCACCCGGTACGGCGACGTTATCGTG  >  NC_000913/699796‑700281

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