Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	281,778	G→A	S69F (TCC→TTC)	yagA ←	integrase core domain‑containing protein YagA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	281,778	0	G	A	100.0%	120.7 / NA	36	S69F (TCC→TTC)	yagA	integrase core domain‑containing protein YagA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/20);  total (16/20)

GTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACC  >  NC_000913/281493‑282062                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                             gTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATg                                                                                                                                                                                                                                                                                >  2:96703/1‑300 (MQ=255)                    cACCGCCAAAATGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGAc                                                                                                                                                                                                                                                             <  1:217163/300‑1 (MQ=255)                                       gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCTGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGc                                                                                                                                                                                                                                          <  1:582938/300‑1 (MQ=255)                                       gCCCTGAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGTTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGGTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGc                                                                                                                                                                                                                                          <  2:166428/300‑1 (MQ=255)                                              aaaTCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAgcgc                                                                                                                                                                                                                                   <  1:403376/300‑1 (MQ=255)                                                                               tgtTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATAAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCAGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCAAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATg                                                                                                                                                                                                  >  2:85889/1‑300 (MQ=255)                                                                                                        tGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGa                                                                                                                                                                         <  1:421050/300‑1 (MQ=255)                                                                                                                 aaGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGccc                                                                                                                                                                <  1:503063/300‑1 (MQ=255)                                                                                                                           cagcagGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTNATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGANTGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGa                                                                                                                                                      <  2:524706/300‑1 (MQ=255)                                                                                                                                                   gTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCCGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAAAAAACTCAGTACg                                                                                                                              >  1:143944/1‑300 (MQ=255)                                                                                                                                                       tgtaCGGTGCTTAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCTTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAAt                                                                                                                          <  2:134308/297‑1 (MQ=255)                                                                                                                                                       tGGATGGTGATGAAGGCGGTCATGGTGTTCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCGTGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCCGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATATTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACTTAAt                                                                                                                          <  2:134756/300‑1 (MQ=255)                                                                                                                                                         gACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGa                                                                                                                        <  1:504794/300‑1 (MQ=255)                                                                                                                                                             gTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATg                                                                                                                    >  1:29101/1‑300 (MQ=255)                                                                                                                                                                 tGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCNGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTAt                                                                                                                <  1:483606/300‑1 (MQ=255)                                                                                                                                                                              ggTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCa                                                                                                   >  2:254419/1‑300 (MQ=255)                                                                                                                                                                              ggTGTGCACATGGTCCTCGAGCCAGCGCTTAATCTTGCAGGCTCCCCAGCGTTCATGAAGGTCGTGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAGTGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTGCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGAAGGCAGAGGGAACGGATGTTCGACCAGTCCTGTGAGGCGAAAAAAACAAACTCAGTACGTACTGACCTGGTactttccgcatcaca                                                                                                   >  2:255282/1‑286 (MQ=255)                                                                                                                                                                                    ccccTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTATCGCATCCCAGGGCAt                                                                                             >  1:337263/1‑300 (MQ=255)                                                                                                                                                                                       cTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAg                                                                                          <  1:304214/300‑1 (MQ=255)                                                                                                                                                                                                     aGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACg                                                                            <  1:667447/300‑1 (MQ=255)                                                                                                                                                                                                      gcgcTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCCCGCATCACAGGGCATAAGCGACTACATAAACgg                                                                           >  2:205366/1‑300 (MQ=255)                                                                                                                                                                                                                         ccccAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGt                                                        >  2:202452/1‑300 (MQ=255)                                                                                                                                                                                                                            cAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAg                                                     <  2:408680/300‑1 (MQ=255)                                                                                                                                                                                                                              gCGTTCATGACGGTCATGGGCCATACGCAGCNGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCNGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTACATAAACGGGTTCTTATGCATTAGTTGTACgtg                                                   >  2:476405/1‑300 (MQ=255)                                                                                                                                                                                                                                   cATGACGGTCATGGGCCATACGCAGCAGAGCCATGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAATACCGGCTGCACCTTCCTGAGCCCAGCGCTGGAGCCACGTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAGCTCAGTACGTAATGACATGTTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGCGTCTAc                                              <  2:403964/300‑1 (MQ=255)                                                                                                                                                                                                                                          gTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGNCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTcc                                       >  1:482603/1‑300 (MQ=255)                                                                                                                                                                                                                                          gTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTACTGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTcc                                       >  2:407699/1‑300 (MQ=255)                                                                                                                                                                                                                                                ggCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAAc                                 <  1:342986/300‑1 (MQ=255)                                                                                                                                                                                                                                                    aTACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAgt                             >  1:649325/1‑300 (MQ=255)                                                                                                                                                                                                                                                     tACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAgtg                            <  1:627873/300‑1 (MQ=255)                                                                                                                                                                                                                                                           gcagGGCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACt                      <  1:594944/300‑1 (MQ=255)                                                                                                                                                                                                                                                                gCCGTGATGTCGTCAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTc                 >  1:67032/1‑300 (MQ=255)                                                                                                                                                                                                                                                                       tGTCGTCAGATGAGCGGTTCTGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTACCCGGAc          <  1:203443/300‑1 (MQ=255)                                                                                                                                                                                                                                                                          cgtcAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGt       >  1:390970/1‑300 (MQ=255)                                                                                                                                                                                                                                                                          cgtcAGATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCAGTACTGAAGACCAGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGACAGGTGAAATGCCGAAGCGACGACAGGGGGAACGGGTGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACCTGGTATCCCTCGAATCCCAGGGCATAAGCGACTCCATAAACGGGGTCTTATGTATTAGTTGTAAGTGTCTACCATGTAACCGAAAAAGTGTTCACTATGtacaagggccgt       >  2:99251/1‑289 (MQ=255)                                                                                                                                                                                                                                                                               gATGAGCGGTTCGGGAAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACAcc  <  1:34215/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                             GTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCATGTCCCCGAACAAGTGTTCACTATGTCCCCGGACCGTACACC  >  NC_000913/281493‑282062

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