Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,160,157	G→A	L212L (CTG→TTG)	sthA ←	soluble pyridine nucleotide transhydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,160,157	0	G	A	100.0%	206.7 / NA	58	L212L (CTG→TTG)	sthA	soluble pyridine nucleotide transhydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/32);  total (26/32)

CAGCTGTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCG  >  NC_000913/4159912‑4160398                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  tagctgtaagctgCTGCCAGATTCAAGCCCAAAGATCTGTAACGCCTTCGAATCGGTATTACCGATTCGAAAGTTGGCATAGAGCAGGCAGTAAGCTTTCCGCTTTTTACCCGACTTCGGCTGCATGATAACACCCTCCACACAGCCTTCGATCTTCTCGTAATCCTCGTTGTGACGAATCACTACGCCACTGTTCCAGCAGTGATGGGAGAGAGAATCTGACATCTCTTGGTCGAGAAATGCCAACAGGc                                                                                                                                                                                                                                              >  2:287567/13‑251 (MQ=255)        ccGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGAgcgc                                                                                                                                                                                                                                        <  2:129631/250‑1 (MQ=255)        ccGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATcgc                                                                                                                                                                                                                                        <  2:382172/250‑1 (MQ=255)        ccGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGCATTACCGGTGCGACCGTTGGCATAGAGCAGGCACTCAGCTTTCAGTTTTCTACCCGACTTCAGCTGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGACTCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGCTcgcg                                                                                                                                                                                                                                       >  1:307176/1‑251 (MQ=255)           cgGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGgt                                                                                                                                                                                                                                    <  1:222593/251‑1 (MQ=255)                gTCAGTTTCTAGCCCAATTTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGAACTTCTCGTACTCATCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTTACATCTCTTGATCGAGAAATGCCAACATGCGATCGCGTGTGtt                                                                                                                                                                                                                                 <  2:301244/249‑1 (MQ=255)                 tCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGTTGTTGAt                                                                                                                                                                                                                              >  1:247853/1‑251 (MQ=255)                     tttCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGCCCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATc                                                                                                                                                                                                                             >  2:443333/1‑248 (MQ=255)                       tCTAGACCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGCGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATccaatc                                                                                                                                                                                                                        >  2:322600/1‑246 (MQ=255)                          aGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGTCCGCTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGTCTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCGTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATcc                                                                                                                                                                                                                       <  1:417772/249‑1 (MQ=255)                          aGCCCAATGTTCTGGAACGCCAGCGAATCGGTATTACTGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTCCCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGTCATCTCTTGATCGAGAAATGCAAACAGGCGATCGCGGGTGTTGATCAGATCCa                                                                                                                                                                                                                      <  2:275688/250‑1 (MQ=255)                                       gTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATcc                                                                                                                                                                                                           <  1:438962/248‑1 (MQ=255)                                       gTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCAt                                                                                                                                                                                                         <  1:368322/250‑1 (MQ=255)                                                    aTCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGa                                                                                                                                                                                                 >  2:100781/1‑245 (MQ=255)                                                              ccGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACGTTTACATCCATACCGCTGAAGATCGACGCATATTc                                                                                                                                                                                 <  1:177172/251‑1 (MQ=255)                                                                ggTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTcac                                                                                                                                                                               >  1:53806/1‑251 (MQ=255)                                                                   gCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGc                                                                                                                                                                            <  1:50115/251‑1 (MQ=255)                                                                   gCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCGTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAg                                                                                                                                                                             <  1:322600/250‑1 (MQ=255)                                                                         gTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATc                                                                                                                                                                       >  2:49755/1‑250 (MQ=255)                                                                            ggCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACt                                                                                                                                                                    <  1:326377/250‑1 (MQ=255)                                                                            ggCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTc                                                                                                                                                                   >  2:283032/1‑251 (MQ=255)                                                                             gCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTcc                                                                                                                                                                  >  1:31722/1‑251 (MQ=255)                                                                             gCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACGTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTcc                                                                                                                                                                  <  1:443333/251‑1 (MQ=255)                                                                                     caggcagTCAGATTTCAGTTTTTTACCCGACTTCAGATGCCTGATCACACCATCGTCCCAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGCGAAATGCCAACAGGCGATCGCGGGAGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCAt                                                                                                                                                          >  2:94412/1‑251 (MQ=255)                                                                                       ggcagTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTTACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAg                                                                                                                                                        <  2:274938/251‑1 (MQ=255)                                                                                        gcagTCAGCTTTCAGTGTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTTACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATTCCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGa                                                                                                                                                       <  1:414012/251‑1 (MQ=255)                                                                                             cAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGt                                                                                                                                                  <  2:31722/251‑1 (MQ=255)                                                                                               gCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTAc                                                                                                                                                <  1:69276/251‑1 (MQ=255)                                                                                                   tCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATg                                                                                                                                             <  1:100781/250‑1 (MQ=255)                                                                                                             ccccGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATa                                                                                                                                                     <  2:185500/231‑1 (MQ=255)                                                                                                             aCCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATa                                                                                                                                                     >  1:185500/1‑232 (MQ=255)                                                                                                             aCCCGACTTCAGATGCATGATCACACCATCGTCACAGCATTCGCTCTTCTCGTACTCTTCGTTGTGACGCATCACTACGCCACTGTTCCCGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCAACATAGCTAAGTACATGGCGCGGTTCgg                                                                                                                                  >  2:62488/1‑250 (MQ=255)                                                                                                                cGACTTCAGATGCATGATCACACCATCGTCACAGACTTCGATCTTCTCGTACTCTTCGTTGTGACGAATAACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATAGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGACGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCgtggg                                                                                                                               >  2:199085/1‑250 (MQ=255)                                                                                                                       aGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCt                                                                                                                        <  1:322369/251‑1 (MQ=255)                                                                                                                               gATCACACCAGCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAAttt                                                                                                                <  2:52952/251‑2 (MQ=255)                                                                                                                                   acacCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTc                                                                                                            <  1:444333/251‑1 (MQ=255)                                                                                                                                                aGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGCGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCCCTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCAACATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGGGAATTGAGACGCTGTCGTaattg                                                                                               >  1:405118/1‑248 (MQ=255)                                                                                                                                                 gCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGc                                                                                              <  1:74770/251‑1 (MQ=255)                                                                                                                                                   cTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCgg                                                                                            >  2:352865/1‑250 (MQ=255)                                                                                                                                                          cttctCGTACTCTTCGTTGTGACGAATCACTGCGCAACTGTTCCAGAAGGGATCGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGACAACAGGCGATCGCGGGTGTTGATCAGATACACTTTTACATCCATACCGCGGAAGATAGACGCATATTCACAGCCGATAACTCCAGCACCATCGATACGTACATGGCGCTGTTCGTGGTGCTTGCTGTGGCTTGAGTCACTGTCGTGAATGCGGGGCTggg                                                                                     >  1:240321/1‑251 (MQ=255)                                                                                                                                                             ctcGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTg                                                                                   >  1:353002/1‑250 (MQ=255)                                                                                                                                                               cGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTg                                                                                   <  1:283032/248‑1 (MQ=255)                                                                                                                                                                gTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAAt                                                                               >  1:307295/1‑251 (MQ=255)                                                                                                                                                                                  aaTCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAACGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGgatgat                                                             <  1:357792/251‑1 (MQ=255)                                                                                                                                                                                         aCGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGaa                                                       <  2:101917/250‑2 (MQ=255)                                                                                                                                                                                         aCGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGAc                                                       <  1:49755/250‑1 (MQ=255)                                                                                                                                                                                            ccACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAACTGCCAACAGGCGATCGAGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGCTCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACg                                                      >  1:233452/1‑248 (MQ=255)                                                                                                                                                                                                 gTTCCAGAAGTGATAGGAGAGAGAATCTGACAGCTCTTGAGCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCTCTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGcc                                                <  1:16658/249‑1 (MQ=255)                                                                                                                                                                                                             aTAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATTCCAACTGGCGGTCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGAGTGTCGCATATTCACAGCCGATCACTCCAGCACCATATATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAAGTGAGTGGCTGTCGTAAATGCGTTGATGGGTGAAATCAACATTTGTGGGATTATATGGACGATAGCCGcgcactataacaa                                  <  2:307176/251‑13 (MQ=255)                                                                                                                                                                                                                 gagagcgaATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAAttt                              >  1:139764/1‑251 (MQ=255)                                                                                                                                                                                                                      agaATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTc                           >  1:47371/1‑249 (MQ=255)                                                                                                                                                                                                                        aaTCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCACTAACAAATTTTTCAGCg                       >  1:127615/1‑251 (MQ=255)                                                                                                                                                                                                                        aaTCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCg                       <  1:438000/251‑1 (MQ=255)                                                                                                                                                                                                                           cTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGCTGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGtt                    >  1:20425/1‑251 (MQ=255)                                                                                                                                                                                                                                    ctTGATCGAGAAATGCCAACAGGCTATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGtt              <  1:18226/248‑1 (MQ=255)                                                                                                                                                                                                                                      tGATCGAGAAATGACAACAAGCGCTCGCAGGTGTTGATCAGATCCACTTTTACATCCATACAGCGGAAGATCGACGCATATTCACAGCGGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGAATGCGGCGAATTGAGTAGCTGCCGTAAATGCGTGGTTGGGTGAAATCAACAACTGTAGGATGATAGGGACAAGAGCCACAGGAACTAACAAAATTTTTAGAGATTAGCGTTTCAAc         >  2:250413/1‑251 (MQ=255)                                                                                                                                                                                                                                             gaAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCg  <  1:275024/251‑1 (MQ=255)                                                                                                                                                                                                                                             gaAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCg  >  1:47680/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  CAGCTGTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCG  >  NC_000913/4159912‑4160398

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