Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,823,288	C→T	G161D (GGC→GAC)	recA ←	DNA recombination/repair protein RecA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,823,288	0	C	T	100.0%	183.2 / NA	58	G161D (GGC→GAC)	recA	DNA recombination/repair protein RecA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/26);  total (32/26)

GTTTCGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCG  >  NC_000913/2823042‑2823527                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                gTTTCGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGAt                                                                                                                                                                                                                                             >  1:245302/1‑251 (MQ=255)    tCGCTACCCACCACGTTTTCGCCCTCTTTCCCCGCGCCGATCCGACGGATGTCTAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTTAGAGTCGTCGAttt                                                                                                                                                                                                                                           <  1:260874/250‑1 (MQ=255)     cGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTCGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCg                                                                                                                                                                                                                                         >  2:138036/1‑251 (MQ=255)                  gTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGCTTTCCg                                                                                                                                                                                                                            >  2:215713/1‑251 (MQ=255)                   ttttCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGTACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGGGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGc                                                                                                                                                                                                                           >  2:326845/1‑251 (MQ=255)                     ttCGCCCTCTTTCCCCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTCCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCt                                                                                                                                                                                                                          <  2:28533/250‑1 (MQ=255)                             ctTTCCCCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCg                                                                                                                                                                                                                   <  1:214495/249‑1 (MQ=255)                                tCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGt                                                                                                                                                                                                              >  1:197198/1‑251 (MQ=255)                                 cACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTg                                                                                                                                                                                                             >  1:215467/1‑251 (MQ=255)                                  aCCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTgc                                                                                                                                                                                                            >  2:219006/1‑251 (MQ=255)                                   ccGCGCCGATACGACGGATGTCGTGACGAACAGAGGCGTAGAATTTCAGCTCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCTTACGTATCTGGTTGATGAAGTTCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCTCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTgcc                                                                                                                                                                                                           <  1:205376/251‑1 (MQ=255)                                       gcCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTAAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCCTGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTgccgcc                                                                                                                                                                                                        >  2:218559/1‑250 (MQ=255)                                             acgacgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGt                                                                                                                                                                                                 >  1:46896/1‑251 (MQ=255)                                             acgacgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGt                                                                                                                                                                                                 >  2:73999/1‑251 (MQ=255)                                                 cgGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTTCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAac                                                                                                                                                                                             <  1:55198/251‑1 (MQ=255)                                                              gAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCt                                                                                                                                                                                >  1:72715/1‑251 (MQ=255)                                                               aaCAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCt                                                                                                                                                                                >  1:125057/1‑250 (MQ=255)                                                                   gtggCGTATAATTTCAGCTCGTGACCCCCGGTAGTGGTTTCCGGGTTACTGAACATCACACCAATTTTCATACGGATCTGGGTGATGAAGATCAGCTGGGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTTGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTgc                                                                                                                                                                           <  2:264591/249‑1 (MQ=255)                                                                               ttCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCGCACCAATTTGCATACGGATCGGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAAcgcg                                                                                                                                                                <  1:138036/250‑1 (MQ=255)                                                                               ttCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGACCgcgc                                                                                                                                                               >  2:110727/1‑251 (MQ=255)                                                                                      cgTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCg                                                                                                                                                        >  1:81050/1‑251 (MQ=255)                                                                                       gTTACCACCGGTAGTGGTTTCCGGGTTCCCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGATTACGCATCGCCTGGCTCATCATACGTGCCGACAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCAAGGAGTCAACGACGAAAACGTCTACTGCGCCAGACCGCGCCAGGGCGt                                                                                                                                                       >  1:241832/1‑251 (MQ=255)                                                                                               cGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGAtt                                                                                                                                               <  1:25962/251‑1 (MQ=255)                                                                                                 gTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTc                                                                                                                                             >  2:298199/1‑251 (MQ=255)                                                                                                            cGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCt                                                                                                                                  >  1:86078/1‑251 (MQ=255)                                                                                                                   ccGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGgt                                                                                                                           >  1:199150/1‑251 (MQ=255)                                                                                                                           cacaCCAATTTTCATCCGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCt                                                                                                                   <  1:215713/251‑1 (MQ=255)                                                                                                                                   ttttCATACGGATCTGGTTGATGAAGATCAGCAGCGGGTTGGACTGCTTCAG‑GTTACCCGCTAGCTTACGCATCGCCTGGCTCATCATACGAGCCGCAAGGCCCATGTGAGTGTCGTCGATTTTGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCTACGGAGACAAAGACGATAACATCTACTGCGCGAGAAAGAGACAAGGAGTCACAGACATCCAGAGACTGATCGCCggagttctactgggaaca                                                                                                            >  1:322795/1‑234 (MQ=255)                                                                                                                                      tCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGATGGGAGCAcagc                                                                                                        >  2:313885/1‑251 (MQ=255)                                                                                                                                        aTACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCAcagcag                                                                                                      <  1:298199/251‑1 (MQ=255)                                                                                                                                          aCGGATCTGGTTGATGACGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGACGCAAGGCCCATGTGCGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGGGCACAGCAGGt                                                                                                    >  1:325199/1‑251 (MQ=255)                                                                                                                                           cggttCTGGTTGATGAAGATTCTCAGCGTGTTGGACTGCTTCAG‑GTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCTTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTTCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGTCGGGGCCCGGCTGGGAGCACAGCATGtt                                                                                                   <  2:343098/247‑1 (MQ=255)                                                                                                                                                  ggTTGAAGACGATCAGCAGCGTGTTGGACTGCTTAAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGCGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCCGTGCCGCCACGGCGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGat                                                                                             >  2:312189/1‑250 (MQ=255)                                                                                                                                                   gTTGATGAAGAGCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGatat                                                                                           <  2:6865/251‑1 (MQ=255)                                                                                                                                                       atgaAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTAAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGAc                                                                                       >  1:346797/1‑251 (MQ=255)                                                                                                                                                         gaAGATCAGCAGCGTGTTGGCCTGCTTCAG‑GTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACg                                                                                      <  2:283077/250‑1 (MQ=255)                                                                                                                                                           aGATCAGCAGCGTGTTGGCCTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGccc                                                                                   <  2:197198/251‑1 (MQ=255)                                                                                                                                                             aTCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGTCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAAcg                                                                                                                                                                  >  1:172434/1‑170 (MQ=255)                                                                                                                                                             aTCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGTCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAAcg                                                                                                                                                                  <  2:172434/170‑1 (MQ=255)                                                                                                                                                                     cGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTg                                                                         >  1:109200/1‑251 (MQ=255)                                                                                                                                                                      gtgtTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGc                                                                        <  1:219006/251‑1 (MQ=255)                                                                                                                                                                        gtTGGACTGCTCCAGTGTT‑CCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGAGTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGt                                                                      <  2:131673/251‑1 (MQ=255)                                                                                                                                                                          tGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCGCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGAGTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTa                                                                     <  1:50000/250‑1 (MQ=255)                                                                                                                                                                                gCTTCAG‑GTTACCCGCCAGCTTACTCATCGCCTGGCTCAGCATACGTGCCGCGAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTTCGTAGATTgg                                                               <  2:51052/250‑1 (MQ=255)                                                                                                                                                                                 gctccagtttACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTg                                                                <  2:72715/241‑1 (MQ=255)                                                                                                                                                                                          tACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGcgcg                                                     >  1:337926/1‑251 (MQ=255)                                                                                                                                                                                                       aCGCATCGCCTGGCTCAGCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCg                                         <  2:125057/250‑1 (MQ=255)                                                                                                                                                                                                               ccTGGCTCATCATACGTGCCGCAAGGCCCATATGAGAGTCGTCGCTTTCGCCTTCGATTTCCGCTATGGGCGTCAGTGCCGCCACGGAGTCAAAGACGATCACGTCTACTGCGCCAGAACGCGCCCGGGAGTCACAGATTTCCAGTGCCTGCTGGCCGGTGACTGGCAGGGAGCACAGCAGGGTGTGGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTACAGCGCGTGTTCGGCATCGATAAACGc                                 >  1:211274/1‑250 (MQ=255)                                                                                                                                                                                                                               ggCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAc                <  2:314382/250‑1 (MQ=255)                                                                                                                                                                                                                                 ccGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcgc              >  1:97587/1‑251 (MQ=255)                                                                                                                                                                                                                                  cGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTAAGTGCCGCCACGGAGTCAACGACGCTAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAagct             >  2:314181/1‑251 (MQ=255)                                                                                                                                                                                                                                    cAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcgctgc           >  1:224720/1‑251 (MQ=255)                                                                                                                                                                                                                                       ggCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTc                  >  2:58551/1‑241 (MQ=255)                                                                                                                                                                                                                                       ggCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTc                  <  1:58551/241‑1 (MQ=255)                                                                                                                                                                                                                                          ccATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTGCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCTCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGctgcgctgcg     <  1:291423/251‑1 (MQ=255)                                                                                                                                                                                                                                            aTGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTgcggc   <  2:97587/251‑1 (MQ=255)                                                                                                                                                                                                                                             tgtgAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTgcggcg  >  2:73127/1‑251 (MQ=255)                                                                                                                                                                                                                                              gtgAGAGTCGTCGATTTCGCCTTCGATTGCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGGTTTACCTTCACGCTGCGCTgcggcg  <  2:283009/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                GTTTCGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGTCGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAG‑GTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCG  >  NC_000913/2823042‑2823527

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