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%	165.8 / NA	53	G161D (GGC→GAC)	recA	DNA recombination/repair protein RecA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/24);  total (29/24)

CGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGT‑CGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATC  >  NC_000913/2823046‑2823530                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   cGCTCCCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGT‑CGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCg                                                                                                                                                                                                                                            <  1:227306/251‑1 (MQ=255)  gCTACCCACCACGTTTTCGCCCTCTTTCCCCGCGCCGATACGACGGAGGT‑CGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGc                                                                                                                                                                                                                                           <  1:385997/251‑1 (MQ=255)  gCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGT‑CGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGACGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGCTTTCGc                                                                                                                                                                                                                                           >  1:204681/1‑251 (MQ=255)              gTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGT‑CGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGAttt                                                                                                                                                                                                                                              <  1:441337/236‑1 (MQ=255)              gTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGT‑CGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGAttt                                                                                                                                                                                                                                              >  2:441337/1‑236 (MQ=255)                       ctctTTCACAGCGCCGATACGACGGATGT‑CGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCACTTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCACATGTGAGCGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCg                                                                                                                                                                                                                      >  2:216964/1‑251 (MQ=255)                       cccttgcACCGCGCCGATACGACG‑TTGTTCGAGACGAACAGGGGCGTAGAAGTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACAAATTTTCATACGGATTGGGGTGATGAAGATCAGCAGCGTGTGGGACTGCTTCAGGTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCg                                                                                                                                                                                                                      <  2:201778/245‑1 (MQ=255)                           ttCACCGCGCCGATACGACGGATGT‑CGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTc                                                                                                                                                                                                                    >  1:6366/1‑249 (MQ=255)                            tCACCGCGCCGATACGACGGATGT‑CGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGt                                                                                                                                                                                                                 >  1:291975/1‑251 (MQ=255)                                    ccGATACGACGGATGT‑CGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCAc                                                                                                                                                                                                         <  1:237132/251‑1 (MQ=255)                                         ccgacgGATGT‑CGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAg                                                                                                                                                                                                     <  1:410083/249‑1 (MQ=255)                                            acgGATGT‑CGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCa                                                                                                                                                                                                  >  1:486985/1‑250 (MQ=255)                                                            aaCAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCACTTTTCATACGGATCTGGTTGCTGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGCCGATAACGTCt                                                                                                                                                                                   >  2:307716/1‑250 (MQ=255)                                                             aCAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGCGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTAc                                                                                                                                                                                 >  2:58707/1‑251 (MQ=255)                                                                  ggCGTAGCATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTCCTgcgc                                                                                                                                                                            >  2:53159/1‑251 (MQ=255)                                                                            ttCAGCGCGTTACCACCTGTAGTGGTTTCCGGGTTACCGAACATCACTCCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGTATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGGCGTCGATTTCGCTTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAAcgcg                                                                                                                                                                   <  1:51340/250‑1 (MQ=255)                                                                                     ttACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTc                                                                                                                                                         >  1:500079/1‑251 (MQ=255)                                                                                      ttccaccGGTATTGGTTCCCGGTTAACCGAAAATCCCCCCATTTTCCATACGACTCTGTTTGATGAAGATCAGCAGCGTGTTGGTCTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTc                                                                                                                                                         <  2:284623/248‑1 (MQ=255)                                                                                        ccaccGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTcac                                                                                                                                                       >  1:240637/1‑250 (MQ=255)                                                                                         caccGGTAGTGGTTTCCGGGTTACCGACCATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGCCGATAACGTCTACTGCGCCCGAACGCGCCAGGGCGTcac                                                                                                                                                       >  2:70705/1‑249 (MQ=255)                                                                                                                ccGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCggg                                                                                                                              >  1:492384/1‑250 (MQ=255)                                                                                                                           cccAATTTTCATACGGATCTGGTTTATGAAGATCAGCAGCGTGTTGGCCTGCTTCAGGTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTGTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTggg                                                                                                                   <  2:498079/250‑1 (MQ=255)                                                                                                                                ttttCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGCGcac                                                                                                              >  1:433582/1‑251 (MQ=255)                                                                                                                                    cagaCGGATCTGGTTGAGGAAGATCAGCAGCGTGTGGGACTGCTTCAGTTTACCCGCCAGCTTACGCAGCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGTGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCAcagca                                                                                                          <  1:219547/248‑1 (MQ=255)                                                                                                                                        cGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGtt                                                                                                      >  1:445151/1‑251 (MQ=255)                                                                                                                                         ggATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTg                                                                                                     >  1:49013/1‑251 (MQ=255)                                                                                                                                               gttgtataaaGACCAGCAGCGGGTGGGTCGGCTGCAGGTTACCCTCCAGCTGACGCATCGCCAGGCTCATTATACGTGCCGCAAGGCCCATGTTATAGTCGTCGATTTCGCCTTCGACTTCCGCTTTCGGCGTCAGTGCCGCCACGTAGTCAACGACGATAACGTCTACAGCGCTAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTACGGCTGGGAGCACAGCAGGGTGTCGata                                                                                               <  2:498424/243‑1 (MQ=255)                                                                                                                                               ggTTGAGGAAGATACGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCAGCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCg                                                                                                  <  2:354423/248‑1 (MQ=255)                                                                                                                                                          aTCAGCAGCGTGTGGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAg                                                                                    <  1:218081/251‑1 (MQ=255)                                                                                                                                                             cgcagcGTGTTGCACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGTCTCATCATACCTGCCGCAATGCCCATGTGAGAGTCGTCGCTTTCGCCTTCGCTTTCCGCTTTCGTCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGACCGCTCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGTAGTCCGGCTGGGAGCCCAGCCGGTTGTCGTTATCGCCGCCCCGttt                                                                                 >  2:191122/2‑251 (MQ=255)                                                                                                                                                               cagcGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTAAGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATCACGTCTACTGCGCCAGCACGCGCCAGGGCGTCACAGATTTCCAGTGACTGCTCGCCAGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGCCTCCCAGTTTAc                                                                               >  1:355358/1‑251 (MQ=255)                                                                                                                                                                 gcGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGt                                                                             <  1:442690/251‑1 (MQ=255)                                                                                                                                                                 gcGTGTTGGACTGCTGCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGt                                                                             <  1:58707/251‑1 (MQ=255)                                                                                                                                                                    tgtTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGGGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGGCGATATCGGCGCCCAGTTTACGTGCg                                                                          >  2:211006/1‑251 (MQ=255)                                                                                                                                                                        ggactctgtCAGGTTCCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAg                                                                       <  2:6366/242‑1 (MQ=255)                                                                                                                                                                            tGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCCGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGTTTgg                                                                  >  2:307922/1‑251 (MQ=255)                                                                                                                                                                            tGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCt                                                                                                                      >  1:371/1‑199 (MQ=255)                                                                                                                                                                              cTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGt                                                                >  1:297726/1‑251 (MQ=255)                                                                                                                                                                                   ggTTCCCCGCCAGCTTACGCATCGCCTGGCTCAGCATACGGGCCGCAAGGCCCATGTGAGAGTCTTCGATTTCGCCTTCGATGTCCGCTGTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGGCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAgc                                                           <  2:229641/251‑1 (MQ=255)                                                                                                                                                                                   gctcaCCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATTTGAGAGTCGTCGATTTCGCCGTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAgc                                                           <  1:111976/247‑1 (MQ=255)                                                                                                                                                                                                                 tcatcaTACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGtt                             >  1:336948/1‑251 (MQ=255)                                                                                                                                                                                                                       tACGTGCCGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCt                       >  1:235741/1‑251 (MQ=255)                                                                                                                                                                                                                              cGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGTGTAGATTGGGTCCAgcg                                                          <  1:331155/209‑1 (MQ=255)                                                                                                                                                                                                                              cGCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcgct                <  1:163923/251‑1 (MQ=255)                                                                                                                                                                                                                               gCAAGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcgctg               <  1:276162/251‑1 (MQ=255)                                                                                                                                                                                                                                 aaGGCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcgctgcg             >  2:31212/1‑251 (MQ=255)                                                                                                                                                                                                                                    gCCCATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCAcgctgcgctg          <  2:299596/251‑1 (MQ=255)                                                                                                                                                                                                                                     cccATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACTCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCTTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACgctgcgctgc         <  1:174622/251‑1 (MQ=255)                                                                                                                                                                                                                                     cccATGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCCCCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCAGGCTGGGAGCACAGCAGGTTGTAGATATCGACGCCCAGTTTACGTGCGTCGATTGGGTACAGCGCGTGTTCAGCATCGATACACGAACAGGTTTTACCTTCAAgctgcgctgc         >  2:286335/1‑251 (MQ=255)                                                                                                                                                                                                                                        aTGTGAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATCACGTCTACTGCGCCAGAACGAGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTgcggc      >  1:288415/1‑251 (MQ=255)                                                                                                                                                                                                                                          gtgAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCgtg                                                      >  1:200477/1‑201 (MQ=255)                                                                                                                                                                                                                                          gtgAGAGTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGGCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCgtg                                                      <  2:200477/201‑1 (MQ=255)                                                                                                                                                                                                                                             agagTCGTCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATc  >  2:71823/1‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   CGCTACCCACCACGTTTTCGCCCTCTTTCACCGCGCCGATACGACGGATGT‑CGAGACGAACAGAGGCGTAGAATTTCAGCGCGTTACCACCGGTAGTGGTTTCCGGGTTACCGAACATCACACCAATTTTCATACGGATCTGGTTGATGAAGATCAGCAGCGTGTTGGACTGCTTCAGGTTACCCGCCAGCTTACGCATCGCCTGGCTCATCATACGTGCCGCAAGGCCCATGTGAGAGTCGCCGATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATC  >  NC_000913/2823046‑2823530

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