Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,235,651	T→C	intergenic (+244/‑255)	pgi → / → yjbE	glucose‑6‑phosphate isomerase/uncharacterized protein YjbE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,235,651	0	T	C	100.0%	262.6 / NA	69	intergenic (+244/‑255)	pgi/yjbE	glucose‑6‑phosphate isomerase/uncharacterized protein YjbE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/46);  total (23/46)

AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCA  >  NC_000913/4235406‑4235889                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               aaTCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTaa                                                                                                                                                                                                                                           <  1:534047/251‑1 (MQ=255) aaTCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTaa                                                                                                                                                                                                                                           <  1:137959/251‑1 (MQ=255)     atcGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAAt                                                                                                                                                                                                                                        >  1:64640/1‑250 (MQ=255)        gttgATATGTAGGCCGGATAAGGCCGTCACGCCGCATCCGGCCACATATGCCTGATGCGACTCGGTCGCGTCTTATCAGGCTTACAGGTCGATGCCTATATGTACATCGTATTCGGCAATTAATACATAGCTCGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGAt                                                                                                                                                                                                                                    <  1:440767/248‑1 (MQ=255)          cgctatttatgCCGGATAAGGCTTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTTGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATATCGTGAGGTTAACCCTAAAATAGAtt                                                                                                                                                                                                                                   <  2:835093/240‑1 (MQ=255)          cgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGAttt                                                                                                                                                                                                                                  <  2:414244/251‑1 (MQ=255)          cgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGAttt                                                                                                                                                                                                                                  <  1:97362/251‑1 (MQ=255)          cgATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGAtt                                                                                                                                                                                                                                   <  1:420957/250‑1 (MQ=255)            atatGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGAttttt                                                                                                                                                                                                                                <  2:378254/251‑1 (MQ=255)                    gCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCgttg                                                                                                                                                                                                                        <  2:466290/251‑1 (MQ=255)                     ccGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCgttgt                                                                                                                                                                                                                       >  1:507575/1‑251 (MQ=255)                     ccGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCgttg                                                                                                                                                                                                                        <  2:467736/250‑1 (MQ=255)                     ccGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCgttg                                                                                                                                                                                                                        >  1:467736/1‑250 (MQ=255)                     ccGGATAAGGCGTTCACGCCGCATACGGCAACCGATGCCTGATGCGAAGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCAGATATGTACATCGTATTCGGCAATTAATACATAGACCGATTGATTAAATAACATTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAACTCATACCGTGAGGTTAACCCTAAAATAGCTTTTTAATCgttgt                                                                                                                                                                                                                       >  2:761406/1‑251 (MQ=255)                                   cACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAt                                                                                                                                                                                                         <  1:479153/251‑1 (MQ=255)                                     cgccgcATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATAc                                                                                                                                                                                                       <  1:750632/251‑1 (MQ=255)                                      gaccattCCGGCACCCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAAGAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACg                                                                                                                                                                                                      <  2:567962/245‑1 (MQ=255)                                       ccgcATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGc                                                                                                                                                                                                     >  1:760142/1‑251 (MQ=255)                                        cgcATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGCTATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATaagc                                                                                                                                                                                                     >  2:795803/1‑247 (MQ=255)                                           aTCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCGACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGAGTTGTAATCGTTGTTTATGTCGGAAAATACGCAGAt                                                                                                                                                                                                 <  1:258033/251‑1 (MQ=255)                                               ggCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTGTATTTCGGAAAATACGCAGATTAAt                                                                                                                                                                                             <  1:625631/251‑1 (MQ=255)                                                 cAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTg                                                                                                                                                                                           <  2:64640/251‑1 (MQ=255)                                                 cAACCGATGCCTGATGCGAAGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCCTTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGATAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTg                                                                                                                                                                                           >  2:643839/1‑251 (MQ=255)                                                  aaCCGATGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGc                                                                                                                                                                                          >  2:55621/1‑251 (MQ=255)                                                         gCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGt                                                                                                                                                                                    <  1:568591/250‑1 (MQ=255)                                                          ccTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGttt                                                                                                                                                                                  <  1:660270/251‑1 (MQ=255)                                                                gCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTttttattt                                                                                                                                                                            >  1:6757/1‑251 (MQ=255)                                                                  gACGCGGTCCGGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACTTCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTttttatttta                                                                                                                                                                          <  2:540567/251‑1 (MQ=255)                                                                    cgcgGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTttttatttta                                                                                                                                                                          >  1:492156/1‑249 (MQ=255)                                                                     gcgGTCGCGGCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGAAGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGt                                                                                                                                                                       <  2:481228/251‑1 (MQ=255)                                                                         gtacctCTTATCATGACTACAGTGCGACGCCGATATGTACATCGTATTCGGCAATTAATTCATAGCACGATTTATTAAATAACCTTAAAATCAATGCCGACGTTATGTCGGCAGTTTGTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCATGAGGTTAACCCTAAAATAGATGTTTAATCGTTTTTTATTTCGTAAATTGCGCAGATTAATTGCTGTTGTTTTTTTTTTAAGTTTAt                                                                                                                                                                   <  2:447917/246‑1 (MQ=255)                                                                            cgTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGAGTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGAt                                                                                                                                                                <  2:563165/251‑1 (MQ=255)                                                                                           aCAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTa                                                                                                                                                 <  2:507575/251‑1 (MQ=255)                                                                                           aCAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTa                                                                                                                                                 <  1:24157/251‑1 (MQ=255)                                                                                            cAGGTCGATGCCGATATGTACATCGTATGCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTaa                                                                                                                                                <  2:873947/251‑1 (MQ=255)                                                                                              ggTCGATGCCGATATTTACATCGTATTCGGCAATTAATACATATCACGATTGATTAAATACCCGTAATAACAATGCCGACGTTATTTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATGCCGTGAGGTTAACCCTAAAATAGTTTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAat                                                                                                                                              <  2:825188/251‑1 (MQ=255)                                                                                                     gCCGATATGTACATCGTAATCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGtt                                                                                                                                       <  2:649521/251‑1 (MQ=255)                                                                                                       cGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGAt                                                                                                                                                                >  1:440108/1‑224 (MQ=255)                                                                                                       cGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGAt                                                                                                                                                                <  2:440108/224‑1 (MQ=255)                                                                                                             gTACATCGTATTCGGCAATTAATACATAACACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACtg                                                                                                                                <  2:59614/250‑2 (MQ=255)                                                                                                                 aTCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATAtt                                                                                                                           <  1:159898/251‑1 (MQ=255)                                                                                                                     tATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGat                                                                                                                       >  1:845987/1‑251 (MQ=255)                                                                                                                         cGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTCTTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATtttc                                                                                                                   >  2:389922/1‑251 (MQ=255)                                                                                                                                          cACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTaaa                                                                                                  <  2:882262/251‑1 (MQ=255)                                                                                                                                                       aTAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGt                                                                                                                                                                                    >  2:810712/1‑156 (MQ=255)                                                                                                                                                       aTAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGt                                                                                                                                                                                    <  1:810712/156‑1 (MQ=255)                                                                                                                                                               aataaCAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCcgc                                                                             <  2:225198/251‑1 (MQ=255)                                                                                                                                                                  aaCAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGcc                                                                          <  1:466617/251‑1 (MQ=255)                                                                                                                                                                   aaaatgtcgaggttatacgggCATTTTTTTATCAGTTCAATCCCTTTGTCGGTATTTTAACGGTGATCATACCGTGAGGATAACCCTAATATAGTTTTTTAATCGTTATTTATTTCGGAAAATACACAGCTTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATATGTTGATACTTATATTTGTTATTCATTCCACTTATCCTAAAACGTCATCGCTAATTCCCCGCGCCg                                                                         >  2:513922/20‑251 (MQ=255)                                                                                                                                                                   aCAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCg                                                                         <  2:404070/251‑1 (MQ=255)                                                                                                                                                                      aTGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTACCGCAAATCATACCGTGAGGTTACCCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCCGATTCATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTCTTGTTATTTAAATATAAGTTGACACTTATATTTGCTATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGtc                                                                       >  1:295117/1‑249 (MQ=255)                                                                                                                                                                             cGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACAGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCAt                                                               <  2:492156/251‑1 (MQ=255)                                                                                                                                                                                 aTGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCt                                                            <  2:10111/250‑1 (MQ=255)                                                                                                                                                                                         aTTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGa                                                                                                                                     >  2:646153/1‑169 (MQ=255)                                                                                                                                                                                         aTTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGa                                                                                                                                     <  1:646153/169‑1 (MQ=255)                                                                                                                                                                                                   cAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAACATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGAt                                         >  2:279419/1‑251 (MQ=255)                                                                                                                                                                                                    aGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATc                                        <  1:795803/251‑1 (MQ=255)                                                                                                                                                                                                      aTAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGATTTTGTTTTTATTTTCAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATCTTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCCTGCTTTAGTTGTGTATACTCGATccc                                      >  2:322309/1‑251 (MQ=255)                                                                                                                                                                                                        aaaTCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTtgt                                                    >  2:458934/1‑235 (MQ=255)                                                                                                                                                                                                        aaaTCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTtgt                                                    <  1:458934/235‑1 (MQ=255)                                                                                                                                                                                                              ccTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGaa                               >  2:37839/1‑250 (MQ=255)                                                                                                                                                                                                                 tGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAAt                             <  1:82311/249‑1 (MQ=255)                                                                                                                                                                                                                      gTAATTTAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGttt                         <  1:339715/248‑1 (MQ=255)                                                                                                                                                                                                                        aaTTTAACGGGAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCCGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAGGTTTAAGCTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCACGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAATTGATTTTggg                    >  2:6831/1‑251 (MQ=255)                                                                                                                                                                                                                         aTTTAACGGAAATCATAACGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTGTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATTTTCATTCCAATTATCCTAAATCGTCATCGCTAATCCCCCGCGCCGTAATTCTAATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGt                   <  2:875576/251‑1 (MQ=255)                                                                                                                                                                                                                           ttAACGGAAATCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTGTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTa                  <  2:173588/250‑1 (MQ=255)                                                                                                                                                                                                                                    aTCATACCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCcatt        >  2:256590/1‑251 (MQ=255)                                                                                                                                                                                                                                         aCCGTGAGGTTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCcattcatt    <  1:55621/250‑1 (MQ=255)                                                                                                                                                                                                                                          ccGTGAGGGTAACCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTATAAGTTTATGATTTTTATTGTTATGTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCAGCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATTCTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCAttcattca  <  2:558610/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATTCCAATTATCCTAAAACGCCATCGCTAATTCCCCGCGCCGTAATTCGCATGCTTTAGTTGTGTATACTCGATCCCGCCCGAAATGTTTTTGGGTAAATCTCCATTCATTCA  >  NC_000913/4235406‑4235889

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