Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,809,487	T→C	Y316Y (TAT→TAC)	waaA →	KDO transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,809,487	0	T	C	100.0%	217.4 / NA	62	Y316Y (TAT→TAC)	waaA	KDO transferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/29);  total (33/29)

TATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGC  >  NC_000913/3809240‑3809732                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      tATGGATTGCCACCAGCACTCACGAAGGCGCAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCGCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGCTACGATGGGCGGGTTGATGTTACTGTACGGc                                                                                                                                                                                                                                                    >  2:80493/1‑251 (MQ=255)    ggATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCAtt                                                                                                                                                                                                                                                 <  1:424813/251‑1 (MQ=255)        tGCCACCAGCAATCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGGGGGCGATACGATGGGCGTGTTGATGTTACTGTACGGCCTTGCCg                                                                                                                                                                                                                                             >  2:393255/1‑251 (MQ=255)          ccaccaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGAt                                                                                                                                                                                                                                           >  1:159359/1‑251 (MQ=255)                gCACTCACGAAGGCGAAGATAGTGTGGTGATCGCCGCACATCAGGCATTGCTACAGCAACTCCCGAATTTATTGCTCATCCTGGTACCCCGGCAACCAGAACGCTGCCTGGATGCGATTAAACGTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCTCCAGCACGCAGGTGGTGGTGGGCGATACGATGGGCGAGTTGATATTACTGTACGTCATTGCCGATCTCGcc                                                                                                                                                                                                                                     <  2:152495/251‑1 (MQ=255)                 cACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCt                                                                                                                                                                                                                                    >  1:314245/1‑251 (MQ=255)                      aCGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTTGTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGCACGGCATTGCCGATCTCGCCTttgtt                                                                                                                                                                                                                               <  1:242195/251‑1 (MQ=255)                                  agagTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTAGTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTGGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCGTTGTTGGCGGTTCACTg                                                                                                                                                                                                                   <  2:179341/251‑1 (MQ=255)                                      tgtgGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTg                                                                                                                                                                                                               >  2:125539/1‑251 (MQ=255)                                          gtgATCGCCGCACATCAGGCATTGTGACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCTCCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACg                                                                                                                                                                                                           <  2:256107/251‑1 (MQ=255)                                            gATCGCCGCACATCAGGCATTGTTCCAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACgtg                                                                                                                                                                                                         >  2:296621/1‑251 (MQ=255)                                                 ccgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTggg                                                                                                                                                                                                    >  1:258740/1‑251 (MQ=255)                                                  cgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGc                                                                                                                                                                                                   >  2:221749/1‑251 (MQ=255)                                                  cgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTCGTAACCCGTCATCCAGAACGCTTCCCGGCTGCGATTCATCTTGTCAGCCAGGCTGGACTAAGCTATATCACCAGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGGTGAACGTGGTGGGc                                                                                                                                                                                                   >  1:415969/1‑251 (MQ=255)                                                             cATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGc                                                                                                                                                                                                   <  1:131547/240‑1 (MQ=255)                                                             cATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGc                                                                                                                                                                                                   >  2:131547/1‑240 (MQ=255)                                                                  ttACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGc                                                                                                                                                                                   >  1:267692/1‑251 (MQ=255)                                                                    aCAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTg                                                                                                                                                                                 >  1:259727/1‑251 (MQ=255)                                                                        cAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGc                                                                                                                                                                             >  2:252935/1‑251 (MQ=255)                                                                            tCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGcaca                                                                                                                                                                          <  2:159359/250‑1 (MQ=255)                                                                             cccGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCATCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCTCTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCacac                                                                                                                                                                         <  2:144006/250‑1 (MQ=255)                                                                                    ttATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTcc                                                                                                                                                                 >  1:68380/1‑251 (MQ=255)                                                                                     tATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCg                                                                                                                                                                >  1:219828/1‑251 (MQ=255)                                                                                     tATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAAGCTTGTCCGCCAGGCTGGACTAAGCTATATCTCACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCTCACGCTATTcc                                                                                                                                                                 <  2:190323/250‑1 (MQ=255)                                                                                      aTTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCgg                                                                                                                                                               >  1:375838/1‑251 (MQ=255)                                                                                           tCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCt                                                                                                                                                                                  >  2:373584/1‑227 (MQ=255)                                                                                           tCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCt                                                                                                                                                                                  <  1:373584/227‑1 (MQ=255)                                                                                                 gggTACCCCGTCATCCAGAAAGCTTCCCGGATTCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCCCTCCAGGGGAAGGCCCCTCGACCAGCACGCAGGTTGTGGTGGGCGATTCGATGGGCGAGTTGATGTTACTGGTCGGCATGTCCGATGCCGTCTTTGTTGGCGATTCTCTGGTTGAGCGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACATGCTATTCCGGGATTGATgggg                                                                                                                                                    <  2:415969/250‑1 (MQ=255)                                                                                                     aCCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGc                                                                                                                                                >  1:277970/1‑251 (MQ=255)                                                                                                       cccGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCAt                                                                                                                                              <  2:348174/251‑1 (MQ=255)                                                                                                        ccGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATa                                                                                                                                             <  2:430504/251‑1 (MQ=255)                                                                                                                cagAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAAc                                                                                                                                     <  2:258740/249‑1 (MQ=255)                                                                                                                  gAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACtt                                                                                                                                   >  1:131849/1‑251 (MQ=255)                                                                                                                         tCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGCCTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGGGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGATGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATCCTTTTAACTTTAAAGa                                                                                                                             <  2:40555/250‑1 (MQ=255)                                                                                                                                    aTTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGcgcgcg                                                                                                                 <  2:122361/251‑1 (MQ=255)                                                                                                                                     ttAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCgg                                                                                                                >  2:8376/1‑251 (MQ=255)                                                                                                                                          cTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCt                                                                                                              >  2:241334/1‑248 (MQ=255)                                                                                                                                          cTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCt                                                                                                              <  1:241334/248‑1 (MQ=255)                                                                                                                                                    aGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGCAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCAAACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGc                                                                                                 >  1:55225/1‑251 (MQ=255)                                                                                                                                                          gACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTg                                                                                           >  1:437273/1‑251 (MQ=255)                                                                                                                                                               aGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCgg                                                                                                                >  2:426253/1‑225 (MQ=255)                                                                                                                                                               aGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCgg                                                                                                                <  1:426253/225‑1 (MQ=255)                                                                                                                                                                 cTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGAttaccg                                                                                    >  1:67943/1‑251 (MQ=255)                                                                                                                                                                     atCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGAttaccgttac                                                                                >  1:235467/1‑251 (MQ=255)                                                                                                                                                                       cacaCGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGGCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGAttaccgttacc                                                                               <  1:169485/250‑1 (MQ=255)                                                                                                                                                                          acGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATg                                                                           <  1:167191/251‑1 (MQ=255)                                                                                                                                                                           cgcgctTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGCTGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGACCGTGGTGGGCATAATCCGCTGGAAGCTGCCGAACACGCTATTCCGGTATTGGTGGGGCCGCATACTTTTACCTTTAAAGCAATTTGCGCGCGGCTGGCGCAGGAAAGCGGGCTGGTTACCGTTACCGaggc                                                                          >  2:52437/5‑248 (MQ=255)                                                                                                                                                                                ttCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCAcgt                                                                     >  2:414311/1‑249 (MQ=255)                                                                                                                                                                                 tCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCAcg                                                                      <  1:177895/249‑2 (MQ=255)                                                                                                                                                                                 tCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGACGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCAc                                                                       <  1:331572/248‑1 (MQ=255)                                                                                                                                                                                    ggggAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCAcgtcgct                                                                 <  1:289548/251‑7 (MQ=255)                                                                                                                                                                                    ggggAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCAcgtcgct                                                                 >  1:113130/1‑245 (MQ=255)                                                                                                                                                                                                   cccaGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACGTCGCTTGCAAAAGAGGTTTc                                                  <  1:181337/250‑1 (MQ=255)                                                                                                                                                                                                           gCAGGTTGTGGTGGGAGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGGACGTGGTGGGGATAATCCGCTGGACGCTGCCGCACACGCTATTCCGGCATTGAGGGGGCCGCAGACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGGGCAGGCAAGCGGGCTGATTAACGTTACCGATGCAACGTCGCTTGACAAAGAGGTTTCCTCTTTa                                           >  2:358669/1‑250 (MQ=255)                                                                                                                                                                                                             aGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACGTCGCTTGCAAAAGAGGTTTCCTCTTTACTc                                        >  1:17515/1‑251 (MQ=255)                                                                                                                                                                                                               gTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACGTCGCTTGCAAAAGAGGTTTCCTCTTTACTCAc                                      <  2:131849/251‑1 (MQ=255)                                                                                                                                                                                                                    ggtgggCGATACGATGGGCGAGTTGATGTTACTGTACGGTATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACGTCGCTTGCAAAAGAGGTTTCCTCTTTACTCTCCGACg                                 <  2:163872/247‑1 (MQ=255)                                                                                                                                                                                                                                gATGTGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGTCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTACTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACGTCGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTa                     <  1:71778/251‑1 (MQ=255)                                                                                                                                                                                                                                       gAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACGTCGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTa              <  1:414311/251‑1 (MQ=255)                                                                                                                                                                                                                                               gTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACGTCGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGCTTACCGTAGTTTCTATGGCCGTc      >  2:159404/1‑251 (MQ=255)                                                                                                                                                                                                                                               gTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACGTCGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTc      <  1:8376/251‑1 (MQ=255)                                                                                                                                                                                                                                                   cTGTACGGCATTGCCGATCGCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACGTCGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGc  >  2:157622/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      TATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCACCGACGCCGATTACCGTAGTTTCTATGGCCGTCATGC  >  NC_000913/3809240‑3809732

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