Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,489,016	C→T	E350K (GAA→AAA)	yjgR ←	DUF853 domain‑containing protein YjgR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,489,016	0	C	T	100.0%	147.6 / NA	45	E350K (GAA→AAA)	yjgR	DUF853 domain‑containing protein YjgR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/22);  total (23/22)

GAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCTGCGGG  >  NC_000913/4488788‑4489239                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                gAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGcc                                                                                                                                                                                                           >  1:197225/1‑251 (MQ=255)      ccTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCAt                                                                                                                                                                                                      >  1:111256/1‑251 (MQ=255)            tGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTgg                                                                                                                                                                                                             <  1:371814/238‑1 (MQ=255)            tGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTgg                                                                                                                                                                                                             >  2:371814/1‑238 (MQ=255)            tGCAACATCTAATAGGCGGATTCCCGGTCCACCTCATCCCCATTTTTGCCATACACCGGAGAGTGATTAACCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGttt                                                                                                                                                                                                 <  2:151090/250‑1 (MQ=255)             gCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGCTCACCATCGCACGCTCCACCCCAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGACTCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTgc                                                                                                                                                                                               >  2:308085/1‑251 (MQ=255)             gCAACATCTCATAGGCGGAGTCCCGGTCCCCCTCATCCTCAGATTGGCCATACCCCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTgc                                                                                                                                                                                               <  2:111256/251‑1 (MQ=255)                     tcATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGc                                                                                                                                                                                       >  1:234816/1‑251 (MQ=255)                      cATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTg                                                                                                                                                                                        <  1:124950/249‑1 (MQ=255)                      cATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGcc                                                                                                                                                                                      <  1:174653/251‑1 (MQ=255)                       aTAGGCGGATTCCCGGTCCCCCTCACCCTCATATTTGCCATACCCCGGAGAGTGATTAATCAAGCCATTACTCTCATCTTCCGTCACCGGCCCCACCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCt                                                                                                                                                                                     <  2:234816/251‑1 (MQ=255)                                      gTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGAt                                                                                                                                                                      >  1:335591/1‑251 (MQ=255)                                       tCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATAGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATc                                                                                                                                                                     >  1:176183/1‑251 (MQ=255)                                              catcctcatATTTGCCATACACCGGAGAGTGATTAACCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTg                                                                                                                                                                <  2:4951/249‑1 (MQ=255)                                                             catccacCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCaaa                                                                                                                                               <  1:299133/247‑1 (MQ=255)                                                                      gagaGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGACCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGATTCCTTCACTGCTTTCTGATCTTTTGGCGTAAACGCCCGCAAAGCATGTTg                                                                                                                                       >  2:47945/1‑250 (MQ=255)                                                                         agTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGACGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGTCGTAAAAGCCCGCAAAGCATGTTGAAcg                                                                                                                                   >  1:75301/1‑251 (MQ=255)                                                                         agTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCAAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGACGGATTGGCCCGAATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAAcg                                                                                                                                   >  1:391496/1‑251 (MQ=255)                                                                          gTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAAcgc                                                                                                                                  >  1:197121/1‑251 (MQ=255)                                                                             atgaaTCAAGCCACTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGAt                                                                                                                               <  1:89079/248‑1 (MQ=255)                                                                             aTTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGCTTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATTTTGAACGCGAt                                                                                                                               >  2:89079/1‑251 (MQ=255)                                                                                 aTCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGCTCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGTTCTTTTGACGTAAAAGCCCGCAAAGCATGTGGAACGCGTTTAcc                                                                                                                           >  1:201914/1‑251 (MQ=255)                                                                                      gCCATGACGCTCATCTGCCGTCACCGGCCCCACCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCGTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCGTTCTGACCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCt                                                                                                                      <  2:71329/251‑1 (MQ=255)                                                                                       ccATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCt                                                                                                                      >  1:470947/1‑250 (MQ=255)                                                                                       ccATTACGCTCATCTTCCGTCACCGGCCCCAGCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTg                                                                                                                     <  1:326498/251‑1 (MQ=255)                                                                                       ccATTACGCTCACCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAATCAGGTTGAACGCGATTACCGAGCt                                                                                                                      <  2:470947/250‑1 (MQ=255)                                                                                               cTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGc                                                                                                              >  1:144844/1‑250 (MQ=255)                                                                                                     tccagtccccgtaccctcccGCGAACAAGGCGCGATCACAATCGCCCGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCTTGAATCGCCGTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCTCGGCTTTCTGATCTTGTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCACGAGCACAtt                                                                                                         <  2:437409/231‑1 (MQ=255)                                                                                                                   ccATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGAc                                                                                         >  1:191441/1‑251 (MQ=255)                                                                                                                         gcgAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACTCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATc                                                                                                                                                                                                                           <  1:416559/115‑1 (MQ=255)                                                                                                                         gcgAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATc                                                                                                                                                                                                                           >  2:416559/1‑115 (MQ=255)                                                                                                                               aaGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGa                                                                              >  1:358291/1‑250 (MQ=255)                                                                                                                                 ggCGCGCTCACCATCGCCCCCTCCCCCCCCGAAGGGCTTCCTTTCGCCTCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGCCGGGTTTTGCGaacc                                                                           <  2:7857/251‑3 (MQ=255)                                                                                                                                  gcgcGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACg                                                                          <  1:284864/251‑1 (MQ=255)                                                                                                                                      gATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAAcc                                                                      <  1:217531/251‑1 (MQ=255)                                                                                                                                         caccacCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAAGCGCCGTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGa                                                                   <  2:304889/251‑1 (MQ=255)                                                                                                                                                   cGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCCTGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCAACGcc                                                         >  2:242476/1‑251 (MQ=255)                                                                                                                                                    gCTCCACCACAGAAGGGCTTCCTTTCTCATCCAGAAAAGAAACCAACTCCTCGCCGTTGCCCAGTTCCTTAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCt                                                        <  2:287783/251‑1 (MQ=255)                                                                                                                                                                           ttCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATc                                   <  1:478034/249‑1 (MQ=255)                                                                                                                                                                           tgcGCATTCAGAAAAGAGATCAACGCCTCGCCGGTGCCAAGTTCCTGAATCGCCTTTTTGGTATCAAGCGCTGGATTGGCCCGCATGGTGTGCGCCGCTGACTTCACTGATTTCTGATCGTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTGCCGAGCTGCCCGGGCACATTATCCGGAATATCAGACGGGGCTTGCGAAACGAACCAGACGCCTACGCCTTTTGAGCGAATAAGCCTTAtta                                  <  1:398456/248‑3 (MQ=255)                                                                                                                                                                                  ccAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTc                           >  1:143161/1‑250 (MQ=255)                                                                                                                                                                                  ccAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTc                           <  1:308085/250‑1 (MQ=255)                                                                                                                                                                                              tCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCcagc              >  1:424998/1‑251 (MQ=255)                                                                                                                                                                                                          cGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCTGCggg  >  2:368103/1‑251 (MQ=255)                                                                                                                                                                                                          cGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACATGCggg  >  2:393931/1‑251 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                GAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCTGCGGG  >  NC_000913/4488788‑4489239

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: ‑