Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,559,230	T→G	S99S (TCA→TCC)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,559,230	0	T	G	100.0%	252.6 / NA	74	S99S (TCA→TCC)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (40/34);  total (40/34)

AAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCA  >  NC_000913/4559001‑4559451                                                                                                                                                                                                                                      |                                                                                                                                                                                                                              aaaCGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACggacgt                                                                                                                                                                                                                          >  1:37300/1‑235 (MQ=255)      tCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACggacgtgacg                                                                                                                                                                                                                      >  1:397591/1‑234 (MQ=255)       ccGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGcc                                                                                                                                                                                                                    >  2:353166/1‑235 (MQ=255)               gCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGGGCCCAGCAGACCAACAACGGCGGTGACGCCCGCTTCCGt                                                                                                                                                                                                           >  2:58246/1‑236 (MQ=255)                       gTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGc                                                                                                                                                                                                      <  2:456806/233‑1 (MQ=255)                        tGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGGGCCCAGCAGACCAACCACGGACGGGACGCCCGCTTCCGTACGGCGAc                                                                                                                                                                                                   >  1:232696/1‑235 (MQ=255)                             agagATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGGGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGCCTTAga                                                                                                                                                                                              >  2:444791/1‑234 (MQ=255)                               agaTGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGAATTAgcga                                                                                                                                                                                            >  2:621181/1‑234 (MQ=255)                                  tGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACt                                                                                                                                                                                        >  1:80668/1‑236 (MQ=255)                                  tGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACt                                                                                                                                                                                        >  1:381969/1‑236 (MQ=255)                                       cATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCgg                                                                                                                                                                                   >  1:457677/1‑236 (MQ=255)                                           tCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGt                                                                                                                                                                                <  2:381969/235‑1 (MQ=255)                                           tCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTg                                                                                                                                                                               <  1:539650/236‑1 (MQ=255)                                           tCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTg                                                                                                                                                                               >  2:680337/1‑236 (MQ=255)                                            cACGCCAAGCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTg                                                                                                                                                                               <  2:80668/235‑1 (MQ=255)                                              cggcAATCACACGATCGATAATCGCCCAGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTgcgc                                                                                                                                                                            <  1:354306/233‑1 (MQ=255)                                              cgcgagtCACACGAGTGATAAGCGCCAAGTCTTTTTCCACGGAACCCGTAATGGTGCGGCAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTgcgc                                                                                                                                                                            <  2:491451/230‑1 (MQ=255)                                                             cgatAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGACTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCTGCGTTCGCGTCGCGGGACCTGc                                                                                                                                                              >  1:537209/1‑235 (MQ=255)                                                               atAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGGGCGCGTCGTGGGACCTGCTTc                                                                                                                                                           >  1:293592/1‑236 (MQ=255)                                                                      ccACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCgccgcc                                                                                                                                                     >  1:26237/1‑235 (MQ=255)                                                                               tctCCGCGGAACCCTGAATGGGGCGGGAAGGGACATGATAAGCGCCGGTCAGCAGCCAGGCGCTGATGCCTTCTTCATTGAGCGCAGGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCCCGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATc                                                                                                                                           <  2:457677/234‑1 (MQ=255)                                                                                 tCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCaa                                                                                                                                         <  1:680337/236‑1 (MQ=255)                                                                                      ggAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCCAATGGa                                                                                                                                    >  1:194779/1‑236 (MQ=255)                                                                                        acccggTAATGGTGGGGGACGGGACATGATAGGCCCTGGTCCGCAGCCGGGCCTTAATGCCTTCTTGATGGAGCGCCCGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTCCCCAGCAGCCCACCCTCGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCCCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACg                                                                                                                                  <  2:232696/231‑1 (MQ=255)                                                                                           ccGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGCCCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACgtg                                                                                                                                <  1:183284/235‑1 (MQ=255)                                                                                             gtaagtgTGCGGGAAGGGACCTGAAAAGCCGCGGTCGACATCCAGGCGGTGATGCCTTCTTCATTGACCCCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGCCCACCCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTg                                                                                                                             <  2:53391/230‑1 (MQ=255)                                                                                              tAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGCCGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGa                                                                                                                            <  1:353166/236‑1 (MQ=255)                                                                                              tAATGGTGCGGGAAGGGACATGATAAGCCCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGa                                                                                                                            <  1:111080/236‑1 (MQ=255)                                                                                                    tGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGGGTCGTGTGACCTGCTTCGCCGCCACCGCCAATCAACTGGACGTGGTGATCAAt                                                                                                                       >  2:426677/1‑235 (MQ=255)                                                                                                    tGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTggg                                                                                                                                                                    >  2:291375/1‑190 (MQ=255)                                                                                                    tGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTggg                                                                                                                                                                    <  1:291375/190‑1 (MQ=255)                                                                                                    tGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACCGGTTGTTCAAt                                                                                                                       >  1:367425/1‑235 (MQ=255)                                                                                                      cGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGCGTTGATCAATa                                                                                                                      >  1:711172/1‑234 (MQ=255)                                                                                                        ggAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCACTa                                                                                                                      >  1:281347/1‑232 (MQ=255)                                                                                                          aaGGGACATGATAAGCCCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAAcc                                                                                                                 <  1:289691/235‑1 (MQ=255)                                                                                                              gACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGc                                                                                                            >  2:490796/1‑236 (MQ=255)                                                                                                              gACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGGGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGGGCCGCTCGTGGGACCTGCTTTGCCGCCACCGCCAATCAAATGGACTGGTTGATCAATAAAACCtctgc                                                                                                            >  2:536305/1‑232 (MQ=255)                                                                                                               aCATGATAAGCGCCGGTCAGCATCCAGGCGGTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGc                                                                                                            <  1:423327/235‑1 (MQ=255)                                                                                                               aCATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCa                                                                                                           >  2:695080/1‑236 (MQ=255)                                                                                                                   gATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCACCCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAgg                                                                                                       <  1:279998/236‑1 (MQ=255)                                                                                                                       ggcgcCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTggg                                                                                                             <  1:58246/225‑1 (MQ=255)                                                                                                                         cgcCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGa                                                                                                                            >  2:169772/1‑209 (MQ=255)                                                                                                                         cgcCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGa                                                                                                                            <  1:169772/209‑1 (MQ=255)                                                                                                                             ggTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCGCGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCa                                                                                              <  1:546575/235‑1 (MQ=255)                                                                                                                              gTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAgc                                                                                                                                                                                              <  1:708037/138‑1 (MQ=255)                                                                                                                              gTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAgc                                                                                                                                                                                              >  2:708037/1‑138 (MQ=255)                                                                                                                               tccacATCCAGGCCCTTGTGCCTTCTTCATGGAGCACACGGTTCTTGGGCAGCCGGGATTCCGGGTGGCGAGGGTTAGCGTCGGTCCCCAGGAGACCAACCAGGGACGTGACGCCCGTTCCCGTCAGGCGACTTAGCGCCACTTCCGGGTGGCGGGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTg                                                                                           <  2:748899/232‑1 (MQ=255)                                                                                                                                             cTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGGGTTGATCAATAAAACCTGGGCAGAGGATATGCCCACTGAGATCGACAAACgg                                                                             >  2:54696/1‑235 (MQ=255)                                                                                                                                              tGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAAc                                                                                >  1:124375/1‑232 (MQ=255)                                                                                                                                                aTGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCCAATGGCCGTGTTGATCAATAAAACCTGGGCAGCGCACCTGCCCACTGAGATCTACAACCTTGc                                                                           >  1:722862/1‑235 (MQ=255)                                                                                                                                                aTGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGCGTTGACCAATAAAACCTGGGCCTAGGATCTGCCCACTGAGACCGATCACCGTGCa                                                                          >  1:173561/1‑236 (MQ=255)                                                                                                                                                    tttctgcATTGAGCGCACGGGTCTGGGCGAGCAGGGATCCCGGGTGGCGAGAGATAGAGTCGGTCACCAGCAGCCCAACCCCGGACGTGACGCCCGTTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGtt                                                                       <  2:367425/234‑1 (MQ=255)                                                                                                                                                     ttcttcATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTc                                                                      >  1:128779/1‑235 (MQ=255)                                                                                                                                                             tGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATg                                                             >  2:550586/1‑236 (MQ=255)                                                                                                                                                                 cgcACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAg                                                         >  2:447017/1‑236 (MQ=255)                                                                                                                                                                  caccggtttggggCGAGCAGGGATCCCGGGTGGCGAGAGAGAGGGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGGGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCAGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAg                                                         <  2:704061/224‑1 (MQ=255)                                                                                                                                                                        tgttgggCGAGCAGGTATTCCGGGTGGCGAGAGATAGAGGCGGTGCCCAGCAGACCAACCACGGACGGGACGCCCGCTCCCGTCAGGCGACTTAGCCCCACTTCCGGCGTGCGGGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGat                                                  <  2:194779/231‑1 (MQ=255)                                                                                                                                                                            tGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAgggct                                                  >  1:726518/1‑230 (MQ=255)                                                                                                                                                                                  gCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCa                                        >  2:22714/1‑236 (MQ=255)                                                                                                                                                                                  gCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCa                                        >  2:386566/1‑236 (MQ=255)                                                                                                                                                                                  gCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGATCGGTACAATGTCCGAAGGGAAATTGCTGGCa                                        >  2:411749/1‑236 (MQ=255)                                                                                                                                                                                   cggggATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGCCCAACCCGGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCa                                        <  1:570246/233‑1 (MQ=255)                                                                                                                                                                                   cAGGAATTCGGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGCCCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGGGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCa                                        <  1:695264/235‑1 (MQ=255)                                                                                                                                                                                          tCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGa                                 <  1:386566/235‑1 (MQ=255)                                                                                                                                                                                              gttGGCAAGAGATAGCGTAGGTGCCTAGCAGCCCAACCACGGACGTGAAGCCGGCTTCCGTCAGGCGTCTTAGCTCCACTTCGGACGTGCGCGTGGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAAtt                            <  2:78518/234‑1 (MQ=255)                                                                                                                                                                                                  gCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCaaa                                                                                                                                        <  2:409915/124‑1 (MQ=255)                                                                                                                                                                                                  gCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCaaa                                                                                                                                        >  1:409915/1‑124 (MQ=255)                                                                                                                                                                                                  gCGAGAGAGAGAGGCGGTGCCAAGCAGACCAGCCAGGGAGGTGACGCCCGGTTCCGTAAGGCGACTTGGCGCCCCTTCCGGCGTGCGCGGCGTGGGACCTGCTGCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTg                         <  2:722862/235‑1 (MQ=255)                                                                                                                                                                                                         aTAGAGTCGGTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGc                 >  2:98713/1‑236 (MQ=255)                                                                                                                                                                                                                 ggTGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAgg          >  1:134881/1‑235 (MQ=255)                                                                                                                                                                                                                  gTGCCCACCGCACCACCCACGGACGTGACGCCCGTTTCGTTCAGGCGACTTAGCGCCCCTTCCGGCGTGCGGGTCGGGGGACCTGCTTGGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACGATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGa         <  2:711172/235‑1 (MQ=255)                                                                                                                                                                                                                   tGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTa                                                                  <  1:262049/177‑1 (MQ=255)                                                                                                                                                                                                                   tGCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTa                                                                  >  2:262049/1‑177 (MQ=255)                                                                                                                                                                                                                    gCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGa                                                        >  2:125414/1‑186 (MQ=255)                                                                                                                                                                                                                    gCCCAGCAGACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGa                                                        <  1:125414/186‑1 (MQ=255)                                                                                                                                                                                                                        agcagACCAACCACGGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCCTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGGGCAGTTCGGTACAATGTCCGAAGGGCTATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGAGGCGGACATCGCa  >  2:62489/1‑236 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                              AAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCTTCATTGAGCGCACGGGTCTTGGCGAGCAGGGATTCCGGGTGGCGAGAGATAGAGTCGGTGCCCAGCAGACCAACCACTGACGTGACGCCCGCTTCCGTCAGGCGACTTAGCGCCACTTCCGGCGTGCGCGTCGTGGGACCTGCTTCGCCGCCACCGCCAATCAAATGGACGTGTTGATCAATAAAACCTGGGCAGAGGATCTGCCCACTGAGATCGACAACCGTGCAGTTCGGTACAATGTCAGAAGGGATATTGCTGGCAACGGCGATAATTTTGCCGTTAGCGACGAGGACATCGCA  >  NC_000913/4559001‑4559451

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