Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,586,096	A→G	R222R (CGT→CGC)	hsdR ←	type I restriction enzyme EcoKI endonuclease subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,586,096	0	A	G	100.0%	178.3 / NA	53	R222R (CGT→CGC)	hsdR	type I restriction enzyme EcoKI endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/28);  total (25/28)

CCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGG  >  NC_000913/4585851‑4586339                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    ccGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTc                                                                                                                                                                                                                                                >  2:136767/1‑251 (MQ=255)  cGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTcc                                                                                                                                                                                                                                               <  1:150681/251‑1 (MQ=255)          aGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGt                                                                                                                                                                                                                                       >  2:270678/1‑251 (MQ=255)           gCACATAATCCGCAAGGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTc                                                                                                                                                                                                                                      <  1:306890/251‑1 (MQ=255)                 aaTCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGtttt                                                                                                                                                                                                                                <  1:33200/251‑1 (MQ=255)                 aaTCCGCAAAGCCCTGATTACCCGTTTCAGCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGGTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGtttt                                                                                                                                                                                                                                <  2:251313/251‑1 (MQ=255)                    ccGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGc                                                                                                                                                                                                                              <  2:67935/250‑1 (MQ=255)                    ccGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTGGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGGCAGCCTGCTTTACGCAGTTGCGCAGCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGcc                                                                                                                                                                                                                             <  1:221495/251‑1 (MQ=255)                          aGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTCCGGTCTGTTTTGCCAGTTcc                                                                                                                                                                                                                       >  1:334230/1‑251 (MQ=255)                           gCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCcg                                                                                                                                                                                                                      >  1:29348/1‑251 (MQ=255)                                 aTTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGttct                                                                                                                                                                                                                >  2:100782/1‑251 (MQ=255)                                              ttgtttCCGGTCGGCCATTCGGTAATGGCTTGCGTGACGCCGGGTTACGGACGTGCGTCTTTGGATAAGCGCAGGGTTTTGCTGTCTGCCTGCCAGCCTGCGTTACGCAGGTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGATTGTGCGCTTGATGGCCTTATCGGTAATTTCTTTGTTGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCGTCTCGGCGAGCt                                                                                                                                                                                                    <  2:316681/247‑1 (MQ=255)                                                ttttCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTACTTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCtgtg                                                                                                                                                                                                 >  2:107016/1‑251 (MQ=255)                                                      ggTCGGCCAGTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCa                                                                                                                                                                                           <  2:332181/251‑1 (MQ=255)                                                         cGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTg                                                                                                                                                                                        >  1:251800/1‑251 (MQ=255)                                                                     aaTGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCg                                                                                                                                                                              >  2:348617/1‑249 (MQ=255)                                                                      aTGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGc                                                                                                                                                                           >  1:120808/1‑251 (MQ=255)                                                                        ggCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCt                                                                                                                                                                         <  1:103380/251‑1 (MQ=255)                                                                                      cGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGGCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCg                                                                                                                                                           <  1:118362/251‑1 (MQ=255)                                                                                         gTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGtt                                                                                                                                                        <  2:254441/251‑1 (MQ=255)                                                                                                 cGTGCGTCTTTGGCGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTATGGTGGAGTGTGCGCTGGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCATTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGTTTGCctgctg                                                                                                                                                <  1:321015/251‑1 (MQ=255)                                                                                                  gTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGt                                                                                                                                               <  1:225040/251‑1 (MQ=255)                                                                                                     cgccttgtgcgaAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAAATTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTg                                                                                                                                            <  2:261293/241‑1 (MQ=255)                                                                                                     cgcCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAAACAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAAATTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTTTGCTTCCAGTTCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGttcg                                                                                                                                            <  2:261316/251‑3 (MQ=255)                                                                                                     cgcCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGGTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTg                                                                                                                                            >  2:369617/1‑251 (MQ=255)                                                                                                       ccTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGcc                                                                                                                                          <  2:318207/251‑1 (MQ=255)                                                                                                             gagaAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTcc                                                                                                                                      >  2:346039/1‑249 (MQ=255)                                                                                                                aaGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCGTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAAt                                                                                                                                 <  2:32719/251‑1 (MQ=255)                                                                                                                 aGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCa                                                                                                                                                                                   <  2:324085/200‑1 (MQ=255)                                                                                                                 aGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCa                                                                                                                                                                                   >  1:324085/1‑200 (MQ=255)                                                                                                                           tttGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAg                                                                                                                       <  1:348617/250‑1 (MQ=255)                                                                                                                            ttGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGcc                                                                                                                     >  2:337028/1‑251 (MQ=255)                                                                                                                            gggCTGTCGGCCTGCCAGCCTGCTTTCCGCAGTTGCGCATCATTCAGGAAGCGTCTCTCTTCTTCGCTACGGTTGAGTGTGCGCTTGATGGCCTGATTGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTGCTTCTCGGCGTTCTTTGCGTCCAGTTCCGCTGGGGTAGTTTGGGTTTGCGCTTCGGTTTCCTGCTTTTTTCCTTCCAGAATTGCGATATAGc                                                                                                                      <  2:163239/248‑1 (MQ=255)                                                                                                                                                tGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCTCTAATGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGATTGGCGATATAGCCGTTCAGGGCAACAAGCTTct                                                                                                 >  2:250799/1‑251 (MQ=255)                                                                                                                                                       cGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTg                                                                                          <  2:302226/251‑1 (MQ=255)                                                                                                                                                          aGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTgcgc                                                                                       >  2:180850/1‑251 (MQ=255)                                                                                                                                                            ttGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCtt                                                                                     >  2:35766/1‑251 (MQ=255)                                                                                                                                                                    caatcaGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTc                                                                              <  1:65263/250‑1 (MQ=255)                                                                                                                                                                      atcaGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGc                                                                           >  2:10558/1‑251 (MQ=255)                                                                                                                                                                                gACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAAACAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTgc                                                                 >  2:221289/1‑251 (MQ=255)                                                                                                                                                                                  ctctctTATTCACTACGGTTGAGTGCGCGCTTGATGGCCTGATCGGTAATTTACTTGTGGTAAGCCTTGCGTTCCTATTCGGTCTGTTTTGCCCGTTCCGCGTTCTTCTCGGCGAGCTGCGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTgcgc                                                               >  2:186322/1‑251 (MQ=255)                                                                                                                                                                                            cGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGc                                                     <  2:125981/251‑1 (MQ=255)                                                                                                                                                                                               tAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGCTCTTTCTTTGTGGGACGACTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTATTCTCGGCGTGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAACGGCGATATAGCCGTTCAGTGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCGGCGTGAGTCTGCGCTTTTTCTCGCAcc                                                  >  2:51149/1‑251 (MQ=255)                                                                                                                                                                                                     tGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGt                                            >  2:73969/1‑251 (MQ=255)                                                                                                                                                                                                      gAGTGTGCTCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCTCGTTCTTCTCGGCGAGCTGTGCTTCCATTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCTTTCAGGGCAACCAGCTTCTGCTGTTGCTCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGtt                                           <  1:185076/251‑1 (MQ=255)                                                                                                                                                                                                       aGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTc                                          >  2:149992/1‑251 (MQ=255)                                                                                                                                                                                                       aGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTc                                          <  1:10558/251‑1 (MQ=255)                                                                                                                                                                                                                   aTGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGttt                              <  1:270678/251‑1 (MQ=255)                                                                                                                                                                                                                            tCGGTAATTTCTTTGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCACGCTGTTGTTTTAGCGTc                       >  2:138147/1‑249 (MQ=255)                                                                                                                                                                                                                                        ttGTGGTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTg         >  1:47964/1‑251 (MQ=255)                                                                                                                                                                                                                                        ttGTGGGAAGCCTTGCGTTCCTGTCCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGt          <  2:334230/250‑1 (MQ=255)                                                                                                                                                                                                                                             gTAAGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGATTGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGAt       >  2:251515/1‑248 (MQ=255)                                                                                                                                                                                                                                               aaGCCTTGCGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAgg  <  1:35766/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    CCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTGTGCGCTTGATGGCCTGATCGGTAATTTCTTTGTGGTAAGCCTTACGTTCCTGTTCGGTCTGTTTTGCCAGTTCCGCGTTCTTCTCGGCGAGCTGTGCTTCCAGTGCCGCAAGGCGAGCCTGGGTTTGCGCTTCGGTTTCCTGCTGTTTGCCTTCCAGAATGGCGATATAGCCGTTCAGGGCAACCAGCTTCTGCTGTTGCGCTTCGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGG  >  NC_000913/4585851‑4586339

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