Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	414,163	G→A	R531* (CGA→TGA)	sbcC ←	ATP‑dependent structure‑specific DNA nuclease ‑ SbcC subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	414,163	0	G	A	100.0%	174.4 / NA	51	R531* (CGA→TGA)	sbcC	ATP‑dependent structure‑specific DNA nuclease ‑ SbcC subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/18);  total (33/18)

AGCTGGCGTTCGTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCGAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGA  >  NC_000913/413919‑414399                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             aGCTGGCGTTCGTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTg                                                                                                                                                                                                                                        >  1:324787/1‑251 (MQ=255)     ggCGTTCGTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTgctt                                                                                                                                                                                                                                     >  1:392419/1‑247 (MQ=255)          tCGTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGATGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGcc                                                                                                                                                                                                                               >  1:91622/1‑251 (MQ=255)            gTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAg                                                                                                                                                                                                                             <  2:413491/251‑1 (MQ=255)             tGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAgg                                                                                                                                                                                                                            >  1:449161/1‑251 (MQ=255)                      tgtgCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGTCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCGTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTATCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAg                                                                                                                                                                                                                      <  2:506844/248‑1 (MQ=255)                                  cACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTc                                                                                                                                                                                                       <  1:415506/251‑1 (MQ=255)                                            tatCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCggg                                                                                                                                                                                             >  1:24311/1‑251 (MQ=255)                                             atCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGGCGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGt                                                                                                                                                                                            <  2:351565/251‑1 (MQ=255)                                                     gTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGATGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCtggtg                                                                                                                                                                                    <  2:91622/251‑1 (MQ=255)                                                           gCAAGGTGATATTGAGGCTGGACGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGAGGCTGGTGGAAcc                                                                                                                                                                               >  1:111050/1‑250 (MQ=255)                                                            cAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACcac                                                                                                                                                                             <  1:14112/251‑1 (MQ=255)                                                                      ttGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTCAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACCAAGTGGGc                                                                                                                                                                    >  2:129387/1‑250 (MQ=255)                                                                        gAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGACCCACAAAGTGGGCa                                                                                                                                                                   >  1:372711/1‑249 (MQ=255)                                                                              ggCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGAc                                                                                                                                                           <  2:444495/251‑1 (MQ=255)                                                                                    gACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATCCGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGcc                                                                                                                                                      <  1:146383/250‑1 (MQ=255)                                                                                    gACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCt                                                                                                                                                     >  1:42505/1‑251 (MQ=255)                                                                                         cTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAAc                                                                                                                                                >  2:386142/1‑251 (MQ=255)                                                                                           tGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCACCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTGGCCCCCGTAGCG‑GCGCACCTTCTCCACCGAGCTTTTTAACT‑TCGTTTCCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCCCCAGCGCCTGATACGCCTCGCCCGCCGGGTGGCTGGTGAAACCAAAAAGTGGGCAAGGCTGCCCCGCCTGTAAc                                                                                                                                                <  1:386142/249‑1 (MQ=255)                                                                                               aTTGTTGTGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGCTTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGcc                                                                                                                                          >  1:365251/1‑250 (MQ=255)                                                                                                     gattaggggctggctCTTCTTGTGGGAGGCTTGGCGCTTCGTTTTCATCCCGCTGAAGCTGCTTTGTTATGGCGTCCAGTGGCCCACGTAGCG‑CCGCCCTTTCTCCCCCGATCTTTTTAACT‑TCGTTTCCCAGCGCCAGTAATCAAGCCTGATAACCGCCAGGCTCCAGCGCCTGATGCGCCTCGCCCGCCGGGTGGCGGGGGAAACCACAAAGGGGGAAGGGCTGCCCCGCCTGTAACTGTGCACGTTGa                                                                                                                                    <  1:212988/239‑1 (MQ=255)                                                                                                           gTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTAGGCTTcc                                                                                                                              >  2:405540/1‑251 (MQ=255)                                                                                                             gcttgctCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTACCT‑TCGTTTTCCCGCGCCAGTAATCAAGACTGATTAACGCCCGGCCCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTcc                                                                                                                              >  1:208782/1‑249 (MQ=255)                                                                                                                ggctcgtcgtGTCGGCGGCTTTCCTCTTCGTTTTAGTCACGCTGAAACTGCTTTGTTATGGCGTACAGTTGCACACGTAGCG‑TCGCACCTTCTTCACCCAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGACTCGACCGCCGGGTGGCTGGCGGACCCACAAAGTGGGCAAGGCTGACCAGCCTGTACCTGTGCAATTTGAGCTTCCAGCGt                                                                                                                         >  2:49559/10‑251 (MQ=255)                                                                                                                  ctCATCTTGTCGGCGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGttt                                                                                                                       >  2:254254/1‑251 (MQ=255)                                                                                                                   cctttttGTCGGAGGCTTTGCGTTTCGTTTTCAACCCGCTGAAGTTGCTTGTTTATGGCTTCCATTTTCCCCCGTAGCG‑TCCCCCCTCCTCCACCGCGTTTTTTAACT‑TCGTTTCCCAGCGCCAGTAATCAAGATTGATAACCGCCGGGCTCCAGCGCCTGATACGCCTCGCCCGCCGGGTGGCGGGGGGACCCACAAGGTGGGCAAGGCGTACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGtt                                                                                                                        <  2:382879/244‑1 (MQ=255)                                                                                                                      tCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGCAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTACCTGTGCACGTTGAGCTTCCAGCGTTTTGAt                                                                                                                   >  2:352664/1‑251 (MQ=255)                                                                                                                            cGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATgcgcgc                                                                                                             >  2:141447/1‑251 (MQ=255)                                                                                                                                        gcTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGc                                                                                                  >  1:126157/1‑250 (MQ=255)                                                                                                                                                   atcaCGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCAAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAACGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCCGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTc                                                                                       >  2:156606/1‑250 (MQ=255)                                                                                                                                                   atcaCGCTGAAGCTGCTTTATTATGGCGTCCAGTTGCCCACGTAGCG‑TAGCACCTTCTTCACCGAGCTTTTTAACT‑TCGGTTTCAAGCGCCAGTAATCAAGACTGATTAACGCCAGGGTCCAGCGCCTGATACGCCTCGACCGCCGGGGGGCTGGTGGAACCACAAGGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGGTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTc                                                                                       >  1:468753/1‑250 (MQ=255)                                                                                                                                                     ccccggtaagggtctTTTTTATGGCCTCCATTCGCCCGCGTAGCGCTC‑CCCCTTCTTCCCCGAG‑TTTTTTACTCTCGTTTTCCAGCGCCATTAAAAAAGACTGTTTACCGCCAGTCTCCGGCGCCTGATACGCCTCGCCCGCGGGGTGGTTGGTGGACCAACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTcaca                                                                                    <  2:18669/238‑1 (MQ=255)                                                                                                                                                        gCTGAAGATGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGCTTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGGGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGCGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTAACATCg                                                                                 >  2:32520/1‑251 (MQ=255)                                                                                                                                                          tGAAGGTGCTTTGATGTGGCGTCAAGTTGGCCACGTAGCG‑TCGCACCTTATTCACCGAGCTTTTTAACT‑TCGTTTTCAAGCGACAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATAAGCATCGAACGCCGGATGGCTGGTGGAACCACAATGTGGGCAAGGCCGACCCGCCTGTACCTGTGCCGGTTGCGCTTCCAGCGTTATGCTGCGAGCTTCCTTCTCGCAACAGGTTTTAACAtagac                                                                               >  2:311514/1‑247 (MQ=255)                                                                                                                                                          agAAGCTGTTTTGTTAGGGCTTCCCTTTGCCCCCGTAGCG‑TCGCCCCTTCTTCCTCGATCTTTTTAACT‑TCGTTTTCCCGCGCCAGTAACCAAGCCTGATTACCGCCAGGCCCCAGCTCCTGATACGCCTCGACCGCGGGGTGGCGGGTGGACCCACAAAGGGGGCAAGGCTGCCCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGACGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGc                                                                               <  2:57247/250‑1 (MQ=255)                                                                                                                                                                        tATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGtct                                                                 >  2:339290/1‑251 (MQ=255)                                                                                                                                                                         aTGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGtctt                                                                >  2:26837/1‑251 (MQ=255)                                                                                                                                                                          gggCGCCCAGGTGCCCACGTAGGG‑TCCCCCCTTCTTCCCCGAGTTTTTTATCT‑TCGTTTTCCAGCGCCATTAATCAAGACGGATAACCGCCAGGCTCCGGCCCCTGATACGCCTCGACCGCCGGGTGGTTGGTGTACCCACAAAGGGGGCAAGGCGGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGtcttt                                                               <  1:402521/250‑1 (MQ=255)                                                                                                                                                                           ggCGTCCCGTTGCCAACGTAGAG‑TCGCACATCCTTCACCGAGATTTTTAACT‑TCGTTTTCCAGTGCCAGTACTCACGACTGATTAACGCCAGGCTCCAGCGCCTGATCCGCCTCGCCCGCCGGGTGGCTAGTGGAACCACAAGATGGGCAAGGCTGACCCGCCTGTACCTGTGCACGTTGAGCTTCCCGCGCTTTGCTGCGCGCTTCCTGCTCGCAAATGGTTTTCACATAGGACAGTTGTTGAGtctttt                                                              >  2:77948/1‑251 (MQ=255)                                                                                                                                                                              gTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGtcttttctt                                                           >  1:138115/1‑251 (MQ=255)                                                                                                                                                                                         cACGTAGCG‑TCCCCCCTTCTCCACCGAGTTTTTTACCT‑TCGTTTCCCAGCGCCAGAAATCAAGACGGATAAACCCCAGGCTCCAGCGCCGGATACGCCCCGCCCGCCGGGGGGCTGGTGGACCCAAAAAGTGGGAAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTgg                                                <  1:149150/251‑1 (MQ=255)                                                                                                                                                                                           cGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGgcg                                              >  1:418459/1‑251 (MQ=255)                                                                                                                                                                                                 g‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTc                                        >  1:281209/1‑251 (MQ=255)                                                                                                                                                                                                             tcccccaatcTTTTTAACT‑TCGTTTCCCAGCTCCAGAAATAAAGACTTATTACCGCCAGGCTCCAGCGCCTGATGCGCCTCGACCGCGGGGTGGCTGGGGGACCCACAAAGTGGGCATGGCTGCCCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGTCGCATTTCGTTAAGTgcg                              <  1:439680/241‑1 (MQ=255)                                                                                                                                                                                                              accccGAGTTTTTTAACT‑TCGTTTTCCAGCGCCATCAATCAAGCCTGATTAACGCCAGGCTCCGGCCCCTGAAGCCCCTCGCCCGCCGGGTGGCGGGTGGAACCACAAGGTGGCCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTgcgg                             <  1:339290/247‑1 (MQ=255)                                                                                                                                                                                                                   gAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTAc                       >  2:185960/1‑251 (MQ=255)                                                                                                                                                                                                                        ttttAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAACGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGc                  >  1:225034/1‑251 (MQ=255)                                                                                                                                                                                                                          ttAACT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGCGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGt                >  2:57056/1‑251 (MQ=255)                                                                                                                                                                                                                              cT‑TCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAACTGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGt            >  2:420102/1‑251 (MQ=255)                                                                                                                                                                                                                                    ttttCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTg       >  1:177951/1‑251 (MQ=255)                                                                                                                                                                                                                                         cggcccaaGTAATCAAGACTGATTAACGCCAGGCGCCAGCCCCTGATACCCCGCGACCGCCGGGGGGCGGGTGGACCCACAAGGTGGGAAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGa  <  1:141447/244‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             AGCTGGCGTTCGTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCG‑TCGCACCTTCTTCACCGAGCTTTTTAACT‑TCGTTTTCCAGCGCCAGTAATCGAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCGTCTGTTCTTGCGTGA  >  NC_000913/413919‑414399

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