Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,956,003	G→A	D225N (GAT→AAT)	ilvA →	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,956,003	0	G	A	100.0%	155.0 / NA	48	D225N (GAT→AAT)	ilvA	threonine deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/32);  total (16/32)

ACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAA  >  NC_000913/3955762‑3956241                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               aCTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTc                                                                                                                                                                                                                                       <  2:268147/251‑1 (MQ=255)        cagcagcagGGGTTCACCGGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCATGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCgtt                                                                                                                                                                                                                                <  2:99352/251‑3 (MQ=255)         agcagcagGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGtt                                                                                                                                                                                                                               <  1:94467/251‑1 (MQ=255)         agcagcagGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGGCATCCGGtt                                                                                                                                                                                                                               >  1:158568/1‑251 (MQ=255)                ggggTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGc                                                                                                                                                                                                                        <  1:85055/251‑1 (MQ=255)                 gggTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAACGGAAGACTCCCCCCGCCTTAAAGCAGCGCTGAAAGCGGGGCATCCGGGTGGTCTGcc                                                                                                                                                                                                                       >  2:274959/1‑251 (MQ=255)                      cACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCcgcgcg                                                                                                                                                                                                                  <  2:262746/251‑1 (MQ=255)                         cTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGGTCAAACAACTGATGCCGCCAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGGCATCCGGTTGATCTGCCGCGCGtcg                                                                                                                                                                                                               >  2:73373/1‑249 (MQ=255)                                        cGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGaa                                                                                                                                                                                                 <  1:274959/250‑1 (MQ=255)                                        cGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGGGATCGCCGGAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGGTGGGCTCTTTGCTGa                                                                                                                                                                                                  >  1:212719/1‑249 (MQ=255)                                        cGACCAGCCGATGGGGATTGCCGGGCTAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATGTGTGCCAGACGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGCCTCCGCGTGCCTGAAAGCAGCGCGGAATGCGGGTCATTCGGTTGATCTGCCGCGGGTATGGCTATTTGCTGAAg                                                                                                                                                                                                <  2:33287/251‑1 (MQ=255)                                         gACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCCAAGGGGTCGCCGGAGAAGGGGAAAACTCCGCCCCCCTGAAAGCCGCGCCGAATGCGGGTCCTCCGGTTGATCTGCCGCGCGGAGGGCGATTTGCTGGAgg                                                                                                                                                                                               >  2:36178/1‑251 (MQ=255)                                          aCCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGGGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAggg                                                                                                                                                                                              >  1:159078/1‑250 (MQ=255)                                               ccGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCGGAAGGCGTAg                                                                                                                                                                                          <  2:279330/250‑1 (MQ=255)                                                    ggTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTaa                                                                                                                                                                                    <  1:196479/251‑1 (MQ=255)                                                                            ggcgctggAACTGCTCCAGCAGGACGCCCATCGCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAGCTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGTGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTTAAGGCGTAGCGGTAAAACGCATCCGTGACGAAACCTTc                                                                                                                                                             <  1:139996/250‑1 (MQ=255)                                                                            ggcgcgggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTcc                                                                                                                                                            <  2:27345/251‑1 (MQ=255)                                                                                  ggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTa                                                                                                                                                       <  2:78912/250‑1 (MQ=255)                                                                                  ggAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGGGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTAt                                                                                                                                                      <  1:196309/251‑1 (MQ=255)                                                                                    aaCTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAAATGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCCGAAAACAGCCCTGAATTCGGGGCATCCGGTTGATCTGCCCCCCGTAGGGCTCTTTTCTGAAGGCGTAGCGGTAAAACGCATCGGTGGCCGAACCTTCCGTTTAttc                                                                                                                                                    >  2:217981/1‑249 (MQ=255)                                                                                        gCTCCAGCAGGACGCGCATCTCGACCGCGTAGTTGTGCCAGGCGGCGGCGGCGTTCGGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCAGCCGGTTGATCTGCCGCGCGTAGGTCTATTTGCTTAAGGCGTAGCTGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCa                                                                                                                                                  <  1:200582/249‑1 (MQ=255)                                                                                           ccAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGt                                                                                                                                             <  2:197575/251‑1 (MQ=255)                                                                                             agcagGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAACCGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTAt                                                                                                                                           <  1:25045/251‑1 (MQ=255)                                                                                               cagGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTAtct                                                                                                                                         >  1:248467/1‑251 (MQ=255)                                                                                                 ggACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTcg                                                                                                                                       <  2:6536/251‑1 (MQ=255)                                                                                                                         tgtgCCGGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgc                                                                                                                <  2:7671/250‑1 (MQ=255)                                                                                                                               agtgggcggcggcggTCCGGCTGGTGGCGTGGCGGTGCTGATCAAACAACTGAGGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCGGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTGGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTGATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAgcgatgcga                                                                                                         <  2:272784/247‑1 (MQ=255)                                                                                                                                         ggcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTgcggc                                                                                               <  1:289312/251‑1 (MQ=255)                                                                                                                                          gcggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTgcggcg                                                                                              <  2:242592/251‑1 (MQ=255)                                                                                                                                           cggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACCACTGGTGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCCGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGaaaa                                                                                                                                                                 >  2:300007/1‑182 (MQ=255)                                                                                                                                           cggTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAAc                                                                                                                                                                 <  1:300007/183‑1 (MQ=255)                                                                                                                                            ggTCTGGATGCTGGCGTGGCGGTGATGAGCAAACAACTGATGCCGCAAAACAAAGTGATCGCCGTAGAAGAGGGAGACTCCGCCGGCCTGAAAGCAGCGCTGAATGCGGGTCAGCCGGTTGAGCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGGAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGATATCGGCACCGTTTATATCGATACGATCTGTGCGGCGAt                                                                                            <  2:285769/251‑1 (MQ=255)                                                                                                                                                 ggCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAgg                                                                                       <  2:212424/251‑1 (MQ=255)                                                                                                                                                    tgctgGCGTGGCGGTGCTGATCAAACAACTGATGCCGCCAATCCAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCCTCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGAtt                                                                                    >  1:30918/1‑251 (MQ=255)                                                                                                                                                           gtggcgGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTAATTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGa                                                                             <  2:212719/251‑1 (MQ=255)                                                                                                                                                            tggcgGTGCTGATCAAACAACTGATGCCGCAACTCAAAGTGATCGACGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGAAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGa                                                                             >  1:174254/1‑250 (MQ=255)                                                                                                                                                                       aTCAAACAAATGATGCCGCAACTCAAAGGGATCGCCGTAGAAGCGGAAGACTACGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATAAGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGCAGGATTTATTCGAAGATGTgcgcgc                                                                 >  2:21443/1‑251 (MQ=255)                                                                                                                                                                           aacaacTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTg                                                             <  1:170949/251‑1 (MQ=255)                                                                                                                                                                                                aaGTGATCGCCGTAGACGCGGAAGACTCCGCCTGCCTGAAAGCAGAGCTGAATGCGGGTCATCCGGTTGATCTGCCGGGCGTAGGGCTATTTGCTGAATGCGTAGCGGTAAAACGCATCGGTGATGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGACCCCTCTGGCgcgctg                                        >  2:212516/1‑251 (MQ=255)                                                                                                                                                                                                   tGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTg                                                >  1:10251/1‑240 (MQ=255)                                                                                                                                                                                                   tGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTg                                                <  2:10251/240‑1 (MQ=255)                                                                                                                                                                                                   tGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggcg                                     >  1:184083/1‑251 (MQ=255)                                                                                                                                                                                                     aTCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggcgct                                   <  1:127684/251‑1 (MQ=255)                                                                                                                                                                                                       cGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCgcgctggcgctg                                  <  2:288813/250‑1 (MQ=255)                                                                                                                                                                                                          cGTAGAAGCGGCAGACTCCGCCTGCCTGAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGCTGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACACTCTGGCGCGCTGGCGCTGGCgg                              >  1:13161/1‑251 (MQ=255)                                                                                                                                                                                                                                 gcctgAAAGCAGCGCTGAATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGccc         <  2:248467/249‑1 (MQ=255)                                                                                                                                                                                                                                 gcctgAAAGCAGCGCTGAAGGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGTAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGccc         <  2:106682/249‑1 (MQ=255)                                                                                                                                                                                                                                      aaaGCAGCGGTGAATGCGGGGCATAAGGTTGAGATGACGCGCGTAGGGCTAGTTGCTGAAGGCGTAGCGGTAAAACGCAACGGTGACGAAACCTTACGTGTATGCCAGGAGTAGCTCGACGACAGCATCACCGTCGATAGCGATGCGATCTGTGCGGCGAGGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCGGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCacaa  <  2:91149/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               ACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGACCGCGTATTTGTGCCAGTCGGCGGCGGCGGTCTGGCTGCTGGCGTGGCGGTGCTGATCAAACAACTGATGCCGCAAATCAAAGTGATCGCCGTAGAAGCGGAAGACTCCGCCTGCCTGAAAGCAGCGCTGGATGCGGGTCATCCGGTTGATCTGCCGCGCGTAGGGCTATTTGCTGAAGGCGTAGCGGTAAAACGCATCGGTGACGAAACCTTCCGTTTATGCCAGGAGTATCTCGACGACATCATCACCGTCGATAGCGATGCGATCTGTGCGGCGATGAAGGATTTATTCGAAGATGTGCGCGCGGTGGCGGAACCCTCTGGCGCGCTGGCGCTGGCGGGAATGAAAAAATATATCGCCCTGCACAA  >  NC_000913/3955762‑3956241

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