Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,411	A→T	G135G (GGA→GGT)	creB →	DNA‑binding transcriptional regulator CreB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,636,411	0	A	T	100.0%	198.6 / NA	66	G135G (GGA→GGT)	creB	DNA‑binding transcriptional regulator CreB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (40/26);  total (40/26)

GGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTA  >  NC_000913/4636172‑4636655                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      ggtttGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTTGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGGATCTTACTGCGTCGGGTGAAGAAGTGCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTttgaat                                                                                                                                                                                                                                             <  1:254314/246‑1 (MQ=255) ggtttGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTttgaatt                                                                                                                                                                                                                                            <  1:135494/247‑1 (MQ=255) ggtttGCCTGATATCAGCGGCTTTGAATTATGACGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGACGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTttgaatt                                                                                                                                                                                                                                            >  1:529579/5‑251 (MQ=255)     tGCCTGATATCAGCGGCTTTGAGTTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATAGCCTGCTTGGGCTGGTAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCACGGTGCGCACATTACTGCGGCGGGAGACGAAGGTCTCGTCGCCGTCCGCCGTCATCCGTATTGGTCATTttgaattgaa                                                                                                                                                                                                                                         >  2:380348/1‑250 (MQ=255)     tGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTttgaattgaa                                                                                                                                                                                                                                         <  1:7405/250‑1 (MQ=255)     tGCCTGATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTAACGGCCCGAAGTGCAGAGGTCGATCGCCTGCTTGGGCTGGCAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCCCCTTACTACGTCGGGTGAAGAAGTTCTAGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTtgcattgaat                                                                                                                                                                                                                                        >  1:610973/1‑251 (MQ=255)      gcctgATATCAGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATtgaatg                                                                                                                                                                                                                                       >  2:478890/5‑251 (MQ=255)      gcctgATATCAGCGGCTGTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAAGTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAATAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTtgaattgaat                                                                                                                                                                                                                                        <  2:373643/246‑1 (MQ=255)       cctgATATCAGCGGCTTTGAATTATGCCGCCAGTTCCTGGCGCTCCATCCGGCGTTCCCTGTACTGTTCCTGACGGCCAGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATtgaatga                                                                                                                                                                                                                                      <  2:331832/248‑1 (MQ=255)           atatCAGCGGTTTTGAATTATGCCGCCAGTTACTGGCGCTCCAGCCGGCGTTACCTGTACTGTGCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAATCCGGTTTCACCCCTCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCGCCCGTCAGCCGTATGGGTCATTTTCAATTGAATGAAccc                                                                                                                                                                                                                                  <  2:536295/251‑1 (MQ=255)               caGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCgcgg                                                                                                                                                                                                                              >  1:220967/2‑251 (MQ=255)                aGCGGCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGCGGTCGATCGCCTGCTTGGGCTGGAACTTGGTGCTGACGCCTATGTGGCTAAACCGTTTTCACCCCGCGCAGTGTGCGCCAGGGTGAGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGGTGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATAAAcccacgac                                                                                                                                                                                                                             >  2:337207/1‑246 (MQ=255)                   ggCTTTGAATTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGTCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGgcgc                                                                                                                                                                                                                           >  1:578310/1‑250 (MQ=255)                          aaTTATGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGc                                                                                                                                                                                                                   <  2:376213/251‑1 (MQ=255)                                 ccgccAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTg                                                                                                                                                                                                           >  1:218253/1‑251 (MQ=255)                                      aGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGACGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACg                                                                                                                                                                                                      >  1:290044/1‑251 (MQ=255)                                            tGGCGCTCCATCCGGCGTTACCAGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGCTCAGCTGG‑TTTGACACGCCAtt                                                                                                                                                                                                 >  2:536207/1‑250 (MQ=255)                                             ggCGCTCCACCCGGCGTTCCCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGAGCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCAtt                                                                                                                                                                                                 <  1:589089/249‑1 (MQ=255)                                              gcgcTCCATCCGGCGTTACCTGTCCTGTTCCTGCCGGCCCGAAGTGAAGAGGTCGATCGCCTGTTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTTTGCGCCAGGGTGCGCACCTTACTGCGTCGGTTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATGGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTaac                                                                                                                                                                                              <  1:363218/251‑3 (MQ=255)                                                gcTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGCATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAacgc                                                                                                                                                                                            >  1:58268/1‑251 (MQ=255)                                                gcTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCCTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAacgc                                                                                                                                                                                            >  2:119306/1‑251 (MQ=255)                                                     aTCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGGGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGAcg                                                                                                                                                                                       >  2:131973/1‑250 (MQ=255)                                                       ccGGCGTTACCTGTACTGTTACTGACGGCCCGGAGTGAAGAGGTCGCTCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGCGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAGCCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAAAGCTGACTCg                                                                                                                                                                                     >  2:56928/1‑251 (MQ=255)                                                         ggCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCggg                                                                                                                                                                                   >  2:350513/1‑250 (MQ=255)                                                          gCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTa                                                                                                                                                                                  <  1:512392/251‑1 (MQ=255)                                                                 cggTACTGTTCTTGACGGCCCGAATTGAAGAGGTCGATCTCTTGCTTGGGCTGGAAATTGGTGCTGACGACTAGGTGGCTAAACCGTTTTCACCCCGCGAAGGGTGCTCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCTTCTCCCTTCATCCTTATTGGTCATTTTGAATTGAATGAACCCGTGGCGCAGATCAACGGG‑TGTGACACGCCATTAACGCTGACGCGGTATGAGttt                                                                                                                                                                           <  2:394870/249‑1 (MQ=255)                                                                     aCTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTAt                                                                                                                                                                       <  1:323248/251‑1 (MQ=255)                                                                          tgcTGACGGCCCGAAGTGAATAGTTCGATCGCCTGCTTGGGCTGTAGATTGGTGCTGACGACGATGTGGCTAAACCGTTTTTCCCCCGCGAGGTGTGCGCCTGGGTGCGCCCTTTATTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCGGTCATCCCTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTAttgct                                                                                                                                                                   <  2:439896/248‑3 (MQ=255)                                                                                    ccGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTACTGACGACTATGGGGCTAAACCGTTTTCACCCCGCGAAGTGAGCGCCAGGGTGCGCACCTTACTGCGTCGGGCGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGCATGGGTTTATATTGCTGAGGACGTTAc                                                                                                                                                        >  2:492588/1‑251 (MQ=255)                                                                                           gAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGt                                                                                                                                                  <  1:577571/250‑1 (MQ=255)                                                                                              gagGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTcgcc                                                                                                                                              >  1:12238/1‑248 (MQ=255)                                                                                                ggTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTcgccag                                                                                                                                            <  2:7311/251‑6 (MQ=255)                                                                                                 gTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTCAACCGTTTTCACACCGCGAAGTGTGCGCCAGGGGGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCCTCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCCTTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTcgcccgg                                                                                                                                           >  2:627592/1‑245 (MQ=255)                                                                                                   cgagcGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGcc                                                                                                                                         <  1:350398/247‑1 (MQ=255)                                                                                                    gATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCcg                                                                                                                                        >  2:421100/1‑251 (MQ=255)                                                                                                        cccTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTc                                                                                                                                    <  1:511383/250‑1 (MQ=255)                                                                                                         ccTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTc                                                                                                                                    >  2:511383/1‑250 (MQ=255)                                                                                                          cTGCTTGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTg                                                                                                                                  >  1:204090/1‑251 (MQ=255)                                                                                                              ttGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCCCCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTc                                                                                                                               <  2:117364/250‑1 (MQ=255)                                                                                                               tGGGCTGGAAATTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCCGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTccc                                                                                                                             >  2:303407/1‑251 (MQ=255)                                                                                                                gggCTGGAAATTGGTGCTTACGACTATGTGGCTAACCCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATGGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTccc                                                                                                                             <  1:615555/250‑1 (MQ=255)                                                                                                                       aaaTTGGTGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCa                                                                                                                      >  1:312264/1‑250 (MQ=255)                                                                                                                                       ctatgTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGtgt                                                                                                     >  2:149188/5‑249 (MQ=255)                                                                                                                                           gTGGCTAAACCGTTTTCCCCCCGCTAAGTGGGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCGCGTCGCCGTCTCCCGACTTCCGTAGTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGAcgcg                                                                                                                                                                                     <  1:589709/167‑4 (MQ=255)                                                                                                                                           gTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTCCTGCGTAGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCCTCCGTATTGGTCATTTTGACTTGGATGAACCCGCGGCGCAAATCAGCTGG‑TTTGACACGCCCTTAACGCTGAcgcg                                                                                                                                                                                     >  2:589709/1‑164 (MQ=255)                                                                                                                                                  aaCCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGCCGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATgcg                                                                                          >  1:329279/1‑251 (MQ=255)                                                                                                                                                         ttCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTCATGGATAGCGTGTGGGAAGATGCGCAGGac                                                                                    >  1:328963/1‑250 (MQ=255)                                                                                                                                                              cccGCGAAGTGTGCGCCAGTGTGCCCCCTTTCCTGCGTCGGGTAAAGAAGTTCTCGTCGCCGTCTCCCTTCATCCGTATTGGTCATTTTGAATGGAATGACCCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTGCTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTAc                                                                              <  1:270189/251‑1 (MQ=255)                                                                                                                                                                    aagtggcccccaggtggcgcCCCTTACTGGGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGCCACCCGATTTGGTATTTTTAAATTGAATGACCCCGGGGCCCAGATCAGTTGG‑TTTGACACGCCATTAACGCTGTCGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGc                                                                        <  1:460151/235‑1 (MQ=255)                                                                                                                                                                     aGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCCTTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACg                                                                             >  2:90094/1‑245 (MQ=255)                                                                                                                                                                          gcgcCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTCATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGtt                                                                  >  2:283526/1‑250 (MQ=255)                                                                                                                                                                                    gcgcACCTTACTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCcac                                                         >  2:24130/1‑250 (MQ=255)                                                                                                                                                                                     cgcACCTTACTGCGTCGGGTGAAGACGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGACTGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGcc                                                                                                                          >  2:604425/1‑184 (MQ=255)                                                                                                                                                                                              cTGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACgtt                                              >  2:166902/1‑249 (MQ=255)                                                                                                                                                                                               tGCGTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACgttg                                             >  1:95485/1‑248 (MQ=255)                                                                                                                                                                                                 cgtgggttgaagaatgtcgcgtcGCCGTCCCCCGTCATCCGTATTGGTTATTTTGATTTAATTGAACCCGCGGCGCAAATCAGTTGT‑TTTGCCACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGTTTCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGc                                            <  1:417476/228‑1 (MQ=255)                                                                                                                                                                                                  gTCGGGTGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATAAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGt                                          >  2:391575/1‑251 (MQ=255)                                                                                                                                                                                                        tGAAGAAGTTCTCGTCGCCGTCTCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATgg                                                                                                              >  1:587863/1‑177 (MQ=255)                                                                                                                                                                                                            gaatttCTCGTCGCCGTCTCCCTTCATCCTTATTGGTCATTTTAATTTGAATGAACCCGCGGCGCAGATCAGCGGG‑TTTAACACGCCATTACCGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTgc                                <  1:397570/247‑1 (MQ=255)                                                                                                                                                                                                            gaagtgccctcccccgcctcCCTTCATCCTTATCGGTTATTTTGAATGGAATGACCCCGCTGCGCAGAACA‑CTGGGTTTCACCCGCCATTAACGCTGACGCGGTATGAGTTCTTATTGCTGAAGACGTTACCCAAGCCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATAACCACATTAAAACGTTGCGTGCCAAGCTgc                                <  1:338264/234‑1 (MQ=255)                                                                                                                                                                                                                        gCCGTCGCCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCCCAGATCAGCTGG‑TTTGACACGCCATTACCGCTGCCGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCAACATTAAAACGTTGCGTGCCAAGCTGCGCGCAATCAAcc                    >  1:257142/1‑251 (MQ=255)                                                                                                                                                                                                                            tctcCCGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCg                 >  1:89430/1‑250 (MQ=255)                                                                                                                                                                                                                               cccGTCATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGctct             >  1:249341/1‑251 (MQ=255)                                                                                                                                                                                                                               cccGTCATCCGTATGGGTCATTTTGAATTGAATGACCCCGCGGCGAAGATCAGTTGG‑TTTGACACGCCATTAACGCTGACGCGGTAGGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGGTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCt             <  1:480694/251‑1 (MQ=255)                                                                                                                                                                                                                                    cATCCGTATTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCAc        >  1:92914/1‑251 (MQ=255)                                                                                                                                                                                                                                           aTTGGTCATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGTCACGCCATTAACGCTGACGCGGTATGAGTTTTTATTGCTGAAGACGTTACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGCTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTa  >  2:143138/1‑250 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      GGTCTGCCGGATATTAGCGGCTTTGAATTGTGCCGCCAGTTACTGGCGCTCCATCCGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTACTGCGTCGGGTGAAGAAGTTCTCGACGCCGTCTCCCGTCATCCGTATTGGACATTTTGAATTGAATGAACCCGCGGCGCAGATCAGCTGG‑TTTGACACGCCATTAGCGCTGACTCGGTATGAGTTTTTATTGTTGAAGACGTTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTA  >  NC_000913/4636172‑4636655

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