Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,193,820	C→T	G129R (GGA→AGA)	thiE ←	thiamine phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,193,820	0	C	T	100.0%	222.5 / NA	61	G129R (GGA→AGA)	thiE	thiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/29);  total (32/29)

GGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCC  >  NC_000913/4193537‑4194111                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                    ggATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAagag                                                                                                                                                                                                                                                                                     <  1:575526/300‑1 (MQ=255) ggATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCAGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAaga                                                                                                                                                                                                                                                                                      >  1:5835/1‑299 (MQ=255) ggATTTGGCGGCTATAACGCATAAAGGCACGGTCATTCATCGCCAACTCCTGCAGTTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCTGCTTGAGTAGTGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCACCAATCGCCACAGTTGGATAATCCACCAGTCACTCAACATGCCGTGCCAGCCGTTCCAGTCCCTGTGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAagag                                                                                                                                                                                                                                                                                     <  2:60441/300‑1 (MQ=255) ggATTTGGAGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAagag                                                                                                                                                                                                                                                                                     <  1:196914/300‑1 (MQ=255)     ttGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGc                                                                                                                                                                                                                                                                                 <  1:352464/300‑1 (MQ=255)           gcTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGtgct                                                                                                                                                                                                                                                                            >  2:690260/1‑299 (MQ=255)           gcTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTAGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAACCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGtgct                                                                                                                                                                                                                                                                            >  2:439290/1‑299 (MQ=255)             tataACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGcc                                                                                                                                                                                                                                                                         <  1:469028/300‑1 (MQ=255)               taACGCATAAATTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTTGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTATTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGtgcttcc                                                                                                                                                                                                                                                                         <  2:587745/298‑4 (MQ=255)                        aaGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTc                                                                                                                                                                                                                                                              >  1:280502/1‑300 (MQ=255)                              cGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCAAGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTAAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTc                                                                                                                                                                                                                                                        >  2:553465/1‑300 (MQ=255)                                gTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCa                                                                                                                                                                                                                                                      >  1:208300/1‑300 (MQ=255)                                        tCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGGTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGTCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCa                                                                                                                                                                                                                                              <  2:170407/300‑1 (MQ=255)                                          gCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCAtg                                                                                                                                                                                                                                            <  1:627554/300‑1 (MQ=255)                                                   tGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGaa                                                                                                                                                                                                                                     <  1:833556/298‑1 (MQ=255)                                                             aGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGcc                                                                                                                                                                                                                           <  1:550476/298‑1 (MQ=255)                                                                 aCTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTAAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGc                                                                                                                                                                                                                     >  1:345938/1‑300 (MQ=255)                                                                  cTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGcc                                                                                                                                                                                                                    >  1:401201/1‑300 (MQ=255)                                                                    gtgCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGANACCGCCAATCGCCACGGTGGGATAATCCGCCANTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCt                                                                                                                                                                                                                   >  2:637196/1‑299 (MQ=255)                                                                         gTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCcg                                                                                                                                                                                                               <  1:439378/298‑1 (MQ=255)                                                                         gTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTAAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCAGCAGGCCTGCcgc                                                                                                                                                                                                              >  2:382253/1‑299 (MQ=255)                                                                                        tCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTgaga                                                                                                                                                                                              <  1:216154/300‑1 (MQ=255)                                                                                             ggCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTNTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGt                                                                                                                                                                                         <  2:615012/300‑1 (MQ=255)                                                                                              gCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGt                                                                                                                                                                                         >  2:825580/1‑299 (MQ=255)                                                                                                      aaTGGCGCTGACGGCGGCGATACTGCCGACACCCGTTGCTATCACCGCATGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTTGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATTGGCGGCTTGCaa                                                                                                                                                                                <  1:394506/300‑1 (MQ=255)                                                                                                                  gacgGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAATGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGccc                                                                                                                                                                    <  1:96055/300‑1 (MQ=255)                                                                                                                                     cccGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGcc                                                                                                                                                 >  1:662440/1‑300 (MQ=255)                                                                                                                                      ccGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCt                                                                                                                                                >  2:676774/1‑300 (MQ=255)                                                                                                                                             tattatcGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCt                                                                                                                                          <  2:345938/293‑1 (MQ=255)                                                                                                                                                              gTGCCAGACTGATACCGCCAATCGCCAAGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGt                                                                                                                        >  2:150239/1‑300 (MQ=255)                                                                                                                                                                    gACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGt                                                                                                                  <  2:805484/300‑1 (MQ=255)                                                                                                                                                                    gACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGt                                                                                                                  >  1:765057/1‑300 (MQ=255)                                                                                                                                                                                      ccACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGt                                                                                                >  1:416648/1‑300 (MQ=255)                                                                                                                                                                                         cGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTc                                                                                                                                                                                                                                                        >  2:828999/1‑145 (MQ=255)                                                                                                                                                                                         cGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTc                                                                                                                                                                                                                                                        <  1:828999/145‑1 (MQ=255)                                                                                                                                                                                           gTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACNCATGTCTCAGCGCGATATAAGATGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAg                                                                                           <  2:662440/300‑1 (MQ=255)                                                                                                                                                                                            tGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTNCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAgc                                                                                          <  1:161462/300‑1 (MQ=255)                                                                                                                                                                                             gggATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCACCAGTAATCGCTGATAAACAATCGCGCGTtata                                                                                            >  2:65134/1‑297 (MQ=255)                                                                                                                                                                                                    tCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGccc                                                                                  <  1:694760/300‑1 (MQ=255)                                                                                                                                                                                                        ccAGTCGCTCAACATGCCGTGCAAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCAGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTACTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCATCCGCCAGTACTCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCa                                                                                 >  2:681226/1‑297 (MQ=255)                                                                                                                                                                                                                tCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTg                                                                        <  2:335378/298‑1 (MQ=255)                                                                                                                                                                                                                 cAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATGGCGGCGGCCCAGCGCAATTgcc                                                                      >  2:562814/1‑299 (MQ=255)                                                                                                                                                                                                                 cAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTgcc                                                                      >  1:136431/1‑299 (MQ=255)                                                                                                                                                                                                                               aGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACATATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGc                                                        <  2:759679/299‑1 (MQ=255)                                                                                                                                                                                                                                   gTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGTGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTACTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGCTAAACAATCGCGCGTTATAGTGGCGGCCCAGCGCAATTGCCGACACGACCTCGACttaca                                                   >  2:162557/1‑297 (MQ=255)                                                                                                                                                                                                                                   gTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCa                                                   >  1:281909/1‑300 (MQ=255)                                                                                                                                                                                                                                   gTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCATGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCAAGACATCGGCTTCCa                                                   >  1:282230/1‑300 (MQ=255)                                                                                                                                                                                                                                   gTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTAAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTACCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCACCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCACAGCGCAATTGCCGCCAAGACATCGGCTTCCa                                                   >  2:162607/1‑300 (MQ=255)                                                                                                                                                                                                                                    ttCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCAc                                                  <  1:825580/300‑1 (MQ=255)                                                                                                                                                                                                                                     tCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCAcc                                                 <  1:741201/300‑1 (MQ=255)                                                                                                                                                                                                                                        aGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACctct                                              >  1:440825/1‑300 (MQ=255)                                                                                                                                                                                                                                                  gTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGcc                                    >  1:3501/1‑300 (MQ=255)                                                                                                                                                                                                                                                  gTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGACAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCgcgc                                     >  2:225340/1‑299 (MQ=255)                                                                                                                                                                                                                                                         aGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCATGATGCTTGATCGCCAGCCGACAGTAGTCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAGTTGCAGCCACGACATCGGCTTACACCTCTTCATCGCGCCGGTCtt                              >  2:196672/1‑299 (MQ=255)                                                                                                                                                                                                                                                             aTCTGTTTGGTTTGCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACTTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATg                         <  1:676774/300‑1 (MQ=255)                                                                                                                                                                                                                                                                tGTTTGGTTTGCGTCGGGAACACATGTCTAAGCGCGATATAAGAGGGGCGTCCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGt                      >  1:403828/1‑300 (MQ=255)                                                                                                                                                                                                                                                                 gtttggtttgCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGNTAAACAATCGCGCGTTATAGCGGCGGCACAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTa                     >  2:664171/1‑300 (MQ=255)                                                                                                                                                                                                                                                                          gCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGgaga             >  2:555256/1‑299 (MQ=255)                                                                                                                                                                                                                                                                          gCGTCGGGAACACATGTCTCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGagag            <  1:67985/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                  aaCACATGTCTCAGCGCGATATAATAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCATGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGtacgtacg    <  2:208300/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                     acaTGTCTCAGCGCGATATAAGATGGGCGTGCTGCCAGCGCGGCGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACTTACGcc  <  2:401201/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                    GGATTTGGCGGCTATAACGCATAAAGTCACGGTCATTCATCGCCAACTCCTGCAATTTCCAGCAACTGTGCCGTTGCCAAACGCCAGTCTGCGGCTTGAGTAATGGCGCTGACGACGGCGATACTGCCGACACCCGTTGCTATCACCGCAGGCGCGCGTGCCAGACTGATACCGCCAATCGCCACGGTGGGATAATCCGCCAGTCGCTCAACATGCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATATCGTCATGTGTCGAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACGCCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCC  >  NC_000913/4193537‑4194111

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