Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	5,642	A→C	intergenic (+112/+41)	yaaX → / ← yaaA	DUF2502 domain‑containing protein YaaX/DNA binding and peroxide stress response protein YaaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	5,642	0	A	C	100.0%	234.6 / NA	72	intergenic (+112/+41)	yaaX/yaaA	DUF2502 domain‑containing protein YaaX/DNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (35/37);  total (35/37)

TGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATGATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACC  >  NC_000913/5402‑5887                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      tggAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑CTTCAACGCCCGATGCGACGCTGGGGCGTCTTATC‑AGGCCTACGTGAATTCTGCCATATATTTAATCTGCATGCTTTTGTAGGCCGGATAAGGCg                                                                                                                                                                                                                                             >  1:138324/1‑251 (MQ=255) gggAAACAACATTATGAAGGGCGAGGCAACCGCTGCCACCCACCCGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTGATC‑AGGCTTACGGGAAGTCGGCAATATATTGAATCTGCATGCTTGTGTAGGCCGGAGAAGGCg                                                                                                                                                                                                                                             <  2:93906/250‑1 (MQ=255)    aaaCACCATTATGAAGGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                                                                                                          <  1:179172/251‑1 (MQ=255)               tGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGGCATGGGCCAAGCAAAAAACACAGCTAAATGACAAAAGTCGGGCAACAATCCGGCCATT‑AGCCCCTGATGCGACGCTGGCGCGGCTTATCCAGGCTAAGGGCAATTCTG‑AATATATTGAATCTGCATGGTTTTGTAGGACGGGTAAGGAGTTCACGCCGCacac                                                                                                                                                                                                                               >  2:176286/1‑248 (MQ=255)                gAATGGCCAGGCAACCGCTGCCACCCACCCGGCCCGCCGCCCCCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                                                                                                              <  2:550629/251‑1 (MQ=39)                   tGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTTAATCTGCATGCTTTTTTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGc                                                                                                                                                                                                                            >  1:149492/1‑250 (MQ=255)                       gAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCCGGCAAACATCACCGCTAAATGACAAAAGCCGGGGAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGGATTTGCATGCTTTTTTAGGCCGGATAAGGCGGTCCCGCCGCCTCCGGCAttt                                                                                                                                                                                                                        >  1:182571/1‑249 (MQ=255)                        aGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGAc                                                                                                                                                                                                                      <  1:286175/251‑1 (MQ=255)                           cAATGGCTGGCCCCCACCCGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCCCGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGc                                                                                                                                                                                                                   <  1:478938/251‑1 (MQ=38)                                gCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGCTCATCACGGCGGACCTGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGGAACAATCCGGC‑ATTAAGCGCCTGCTGCCACGCTGGCGCGTCTTATC‑AAGCCTCCGTGAATTCTTCCAAATATTGAATCTGCATGCTTTTGTTGGCCGGATAAGGCGTTTACGCCCCCTCCGGCCTTGACTGCCAACt                                                                                                                                                                                                              >  2:283693/1‑251 (MQ=255)                                   ggCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTa                                                                                                                                                                                                            >  1:43710/1‑250 (MQ=255)                                   gcgtccccaccggcccgccgCCCCCGCCCCGCCACCATAAGAAAGCTCCTCATGAGCAACACGTCGGTCATGGTCCAGGCAAACATCACCGATAAACGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGTCTGAGGCTACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCGTCAATATATGGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCATTGACTGCAAACTTaa                                                                                                                                                                                                           <  2:219553/237‑1 (MQ=25)                                        ccacacGGACCGCCGCCACCCCCTCGCCACCATAAGAACGCTCCTCATGATCATCACGGCGGCCATGGTCCAGTCAAACATCGCCGATAAATGCCAAATGCCGGGTAAAAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTAATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACgctg                                                                                                                                                                                                      <  2:346485/249‑1 (MQ=38)                                               gacgcccgccACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGTTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAACGCCTGATTCGACGCTGGCGCGTCTGATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGTTTTTGTAGGCCGTATAAGGCGTTCACGCCTCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAg                                                                                                                                                                                               <  2:89584/246‑1 (MQ=255)                                               gACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAAAGCCGGGGAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGGGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGGTCACGACGCATCCGGCATTTACTGCAAAATTTAAGCTGCTCgg                                                                                                                                                                                                 >  1:604917/1‑248 (MQ=255)                                                aCCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATTCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCACTATATTGAATCTGCATGCCTTTTTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACACTGCTCGTAGc                                                                                                                                                                                              >  1:219502/1‑251 (MQ=255)                                                    ccgccACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGttt                                                                                                                                                                                          <  1:411276/251‑1 (MQ=255)                                                      gccACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCCAACATCAACGCTAAATGAAAAAAGACGGGAAACAATCCGGC‑AATCAGCGCCACATGCGACGCTGTCGCCGCTTTTA‑AGGCCTACCTGAATTCTGCAATATTTTGAATCTGCCTGCTTTTTTAGGCCTGATAAGGCGTTCCCGCCACATCCGGCATTTACTACaaaactcacgcttctcctaacgtttcc                                                                                                                                                                                        >  1:356292/1‑227 (MQ=255)                                                         aCCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAaca                                                                                                                                                                                     <  1:283693/251‑1 (MQ=255)                                                             ccgcGCAACCATAAGAAAGCTCCACATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAg                                                                                                                                                                                 <  2:411235/251‑1 (MQ=255)                                                                     ccaTAAAAAAGCCCCCCATGATCATCACGGCGGTCATGGCCCAGGCAAACATCACCGCTAAATGACAAATGCCGGTTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCAt                                                                                                                                                                         <  2:130876/251‑1 (MQ=39)                                                                       aTAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTg                                                                                                                                                                       <  2:67416/251‑1 (MQ=255)                                                                        tAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCATGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAACGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTTTAGGCCGGGTAAGGCGGTCACGCCGCATACGGCCTTGACTGCAAAATTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGc                                                                                                                                                                      >  2:546843/1‑251 (MQ=255)                                                                          aGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTg                                                                                                                                                                    <  1:529777/251‑1 (MQ=255)                                                                             aaGCTCCTCATGATCATCACGGCGGTCAGGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAAGCCGGC‑ATCCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTggag                                                                                                                                                                 <  1:162828/251‑1 (MQ=255)                                                                                  cctcATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATGCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGaag                                                                                                                                                             <  1:147303/250‑2 (MQ=255)                                                                                        gagcatcaCGGCGGTCATGGTCCAGGCAAACATCCCCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCTACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTtcatc                                                                                                                                                      <  2:227347/248‑1 (MQ=255)                                                                                                    ggTCATGGTCCAGGCAAACATCACCGCTAAATTACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATGGCTGGAGTAAGCTTCATCAAAGAAGTAAcc                                                                                                                                          <  1:109328/251‑1 (MQ=255)                                                                                                     ggcATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCt                                                                                                                                         <  1:525338/249‑1 (MQ=255)                                                                                                        aTGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCg                                                                                                                                      <  2:138324/251‑1 (MQ=255)                                                                                                             ccAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTAtt                                                                                                                                 >  2:343330/1‑251 (MQ=255)                                                                                                             ccAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTAtt                                                                                                                                 <  1:366763/251‑1 (MQ=255)                                                                                                               aGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTaa                                                                                                                               <  2:181074/251‑1 (MQ=255)                                                                                                                   aaaCATCACCGCTAAATAACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAAcc                                                                                                                           >  1:443451/1‑251 (MQ=255)                                                                                                                    aaCATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAAcc                                                                                                                           >  1:513330/1‑250 (MQ=255)                                                                                                                       aTCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAACGAAGTAACCTTCGCTATTAAAACcagtc                                                                                                                       >  2:128946/1‑251 (MQ=255)                                                                                                                          cccGCTAAAGGACAAATGCCGGGTAACATTCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCTTACGTGAATTCTGCAATATATTGAATCTGCATGCTGTTGTAGGCCGGATAAGGCTTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCTTTTAAACTCCAGTTCGCCATTGCTGGATGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACctgtcagt                                                                                                                    <  2:601037/250‑1 (MQ=255)                                                                                                                                aaaTGACAAATGCCGGGTAACAATCCGAC‑ATGCGGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATGGAAGCTGCATGCTTTTGTAGGCCGGATAAGGCTGTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATGGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGctc                                                                                                               <  2:559210/250‑1 (MQ=255)                                                                                                                                       aaaTGCCGGGTAACAATCCGGC‑ATTCAGCGCATGATGCGAAGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAATCCAGTCAGTTGCTCTGGTTTgg                                                                                                       >  2:77020/1‑251 (MQ=39)                                                                                                                                        aaTGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGCATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAACCCAGTCAGTTGCTCTGGTTTGGt                                                                                                      >  2:101231/1‑251 (MQ=255)                                                                                                                                         aTGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGTCGCTGGCGCGTCTTATC‑AGGCCTACGTTAATTCTGCAATATATGGAAGCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTGCGCTATTAAAACCAGTCAGTTGCGCTGTTTTGGTc                                                                                                     <  2:568926/251‑1 (MQ=255)                                                                                                                                          ggCCGGTTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATTCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCAGCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTTTTAAAACCAGTCAGTTGCTCTGGGGTGGTCa                                                                                                    <  2:386087/250‑1 (MQ=255)                                                                                                                                           gCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCa                                                                                                    <  2:182571/250‑1 (MQ=255)                                                                                                                                                  aaCAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGAttt                                                                                            >  2:439937/1‑251 (MQ=255)                                                                                                                                                         cGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCaata                                                                                      >  1:588982/1‑250 (MQ=255)                                                                                                                                                          ggC‑ATTCAGAGCCTGATGCGTCGCTGGGGGGTCTTATC‑AGGCCTCGGTGCATTCTGCAATATATTGAATCTGCATGCTTTTGTAGCCCGGATAAGGCGTTAACGCAGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCCGTTCGCCATTGCTGGAGGAAGCTTAAACAAAGTAGTAACCTTAGCTATTAACACGAGTCTGTTGCTCTGGTTTGGTCATCCGATCTTCtataa                                                                                     >  2:30322/1‑250 (MQ=255)                                                                                                                                                                 cAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACg                                                                              >  1:444429/1‑251 (MQ=255)                                                                                                                                                                      ccTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCATTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACtca                                                                         <  1:499542/251‑1 (MQ=255)                                                                                                                                                                         gaggCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCTTTAACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGGTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACtcatca                                                                      <  2:581382/248‑1 (MQ=255)                                                                                                                                                                             cGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGAcc                                                                  >  2:563882/1‑251 (MQ=255)                                                                                                                                                                              gACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACcg                                                                 >  1:121295/1‑251 (MQ=255)                                                                                                                                                                                cGCTGGCGCGTCTTATC‑CGGCCTACGTGAATTACGCCATATATTGCATCTGCATGCTTTTGTAGGCCGGATACGGCGTTCACGCCGCATCCGGTATTGACTGCAAACTTAACGCTGCTCGTCGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGACGTCACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGCCTCATCAGACagcg                                                               >  2:30890/1‑251 (MQ=255)                                                                                                                                                                                cGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACcgcg                                                               <  1:439937/251‑1 (MQ=255)                                                                                                                                                                                 gCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAACACGACTCATCAGACcgcat                                                              >  1:371235/1‑249 (MQ=255)                                                                                                                                                                                            ttATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGc                                                   >  1:23864/1‑251 (MQ=255)                                                                                                                                                                                              aTC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATACTAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGt                                                 >  2:400483/1‑251 (MQ=255)                                                                                                                                                                                                  aGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCCGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGt                                                 >  1:343904/1‑248 (MQ=255)                                                                                                                                                                                                     ccTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGc                                           >  2:325508/1‑251 (MQ=255)                                                                                                                                                                                                           tgaaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCtgatt                                      >  2:488606/3‑249 (MQ=255)                                                                                                                                                                                                            gaaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTgg                                                                                                       >  2:533205/2‑184 (MQ=255)                                                                                                                                                                                                            gaaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTGCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTgg                                                                                                       <  1:533205/183‑1 (MQ=255)                                                                                                                                                                                                                      atatatTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTg                          <  2:170605/251‑1 (MQ=255)                                                                                                                                                                                                                          atTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCg                       <  1:344895/250‑1 (MQ=255)                                                                                                                                                                                                                                 cTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTtctc               >  1:545457/1‑251 (MQ=255)                                                                                                                                                                                                                                   gCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATACAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCAt             >  2:533227/1‑251 (MQ=255)                                                                                                                                                                                                                                   gCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCAt             >  2:502037/1‑251 (MQ=255)                                                                                                                                                                                                                                   gCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCAt             <  1:343330/251‑1 (MQ=255)                                                                                                                                                                                                                                   gCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTAGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTtctc               <  1:458069/249‑1 (MQ=255)                                                                                                                                                                                                                                    cATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATc            <  1:344232/251‑1 (MQ=255)                                                                                                                                                                                                                                    cATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGAAAACTTAACGCTGATCGTAGCGTTTAAACACAAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGCTTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATc            >  1:577858/1‑251 (MQ=255)                                                                                                                                                                                                                                            tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAGGAac     <  2:543858/250‑1 (MQ=255)                                                                                                                                                                                                                                              tAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAGGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCGGCCGTTTTTCAATAATAAAACGCCTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTATTCGAGGaccacc  >  2:136384/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      TGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATGATCTTAAATTTGCCGTTCTTCTCATCGAGGAACACC  >  NC_000913/5402‑5887

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