Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,306,035	C→T	G188D (GGC→GAC)	yjcS ←	linear primary‑alkylsulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,306,035	0	C	T	100.0%	153.9 / NA	50	G188D (GGC→GAC)	yjcS	linear primary‑alkylsulfatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/29);  total (21/29)

TAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAA‑CAATAATGCCTTTCTCGCCCTCAACG‑AACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGAC‑TTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACG  >  NC_000913/4305796‑4306276                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     tggttctcgtcgGGGAAATAATGCTCGTGGCGCCCGTTGCCAGCGTCCCGCCAAGGCCTTTACCCACATTGCCGTGCGCTTTGTGCGTGAGTAAAAGACCGTAAGAGTAGAGCGCAGGGCGGCTCATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCTCCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTg                                                                                                                                                                                                                                            <  1:248149/244‑1 (MQ=255)   gTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTggac                                                                                                                                                                                                                                         >  2:326977/1‑249 (MQ=255)     tttCGTCGGTGCAATAATGCTCGGGCCGCCCGTTGCCAGCGTCACGCCAAGGCCATTCCCCACATGGCCTTGCGCTTTGTGCGGCAGTAACAGCCCGTAAGAGTAGGGCGCAGGGGGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCg                                                                                                                                                                                                                                        <  1:359397/250‑1 (MQ=255)       tcgtcgGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGgtg                                                                                                                                                                                                                                     >  2:137926/1‑251 (MQ=255)          ccgGTGCAATAATGCTCGGGTCGCCCGTGGCCAGCGTCCCGCAAAGGCCATTCCCCACTTTGCCTTGCGGGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTgg                                                                                                                                                                                                                                  <  1:201456/250‑1 (MQ=255)             gTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCCCGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGCCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCCTAGTGGTCGGTGTGGCtg                                                                                                                                                                                                                               >  1:32637/1‑251 (MQ=255)                aaataatGCCCGGGTCGCCCGTTCCCGGCGCCCCGCAAGGGCCAT‑CCCCACATGGCTTTGCGGGTTGGGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAg                                                                                                                                                                                                                           <  1:388529/250‑1 (MQ=255)                                                  gCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCGAAAATCGGGTTTTGCGGAc                                                                                                                                                                                          >  2:153178/1‑251 (MQ=255)                                                          cATTCCCCACACTGCCGTGCGCGTGGGGCGGCAGCAACAGACCGGAAGAGTAGACCGCACGGCGGCTGATGATGTCACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCAGAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTgg                                                                                                                                                                                  <  2:160062/251‑1 (MQ=255)                                                                 cacaTTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAAAAATCGGTTTTTGCGGACGATGCTGGAAGTaa                                                                                                                                                                            >  1:239548/1‑250 (MQ=255)                                                                                  tgtgGGGCAGAACCAGACCGTAAGAGTTGAGCGCACGGCGGCTCATTATTTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCttt                                                                                                                                                           <  2:391532/250‑1 (MQ=255)                                                                                   gtgCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCCTGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGAGGTCGGTGTGGCTGTGAGTGTAGATAACGGCACCAATCGGTTTTTGCGGACGAAGCTGGAAGTAAAGGTCAAGGGCGGCttttt                                                                                                                                                         >  2:301998/1‑250 (MQ=255)                                                                                   gggcggtagaaaaaGACCGTAAGAGTAGAGCGCCCGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCtttt                                                                                                                                                          <  1:235276/237‑1 (MQ=255)                                                                                     gCGGCAGTAACAGACCGTAAGAGTAGAGCGCCGGGGGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCg                                                                                                                                                       <  1:275569/251‑1 (MQ=255)                                                                                      cGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACTATGCTGGAAGTAAAGGTCAAGGGCGGCCTTTGCgg                                                                                                                                                      >  2:11617/1‑251 (MQ=255)                                                                                              aCAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGCACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGGCGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTTCGGCAgggggg                                                                                                                                              >  2:461654/1‑247 (MQ=255)                                                                                              aCAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAggcggc                                                                                                                                              <  1:137926/251‑1 (MQ=255)                                                                                               aaGACCGTAAGAGTAGGGCGCACGGCGGCCTATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACCTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGgcggcg                                                                                                                                             <  1:326977/250‑1 (MQ=255)                                                                                                gcaCCGAAAGAGTAGAGCGCACGGCGGCCCATGATTTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGt                                                                                                                                            <  1:219767/249‑1 (MQ=255)                                                                                                   cgggaaGAGTAGATCGCGCGGGGTCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCGGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCTCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCAc                                                                                                                                         <  1:196118/247‑1 (MQ=255)                                                                                                        ggaggggggcgcccggtgctacATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCg                                                                                                                                    <  2:220977/230‑1 (MQ=255)                                                                                                          agTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCACTCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCCGCggg                                                                                                                                  >  1:289644/1‑251 (MQ=255)                                                                                                             agcgcgcACGGGTGCTCATGATGTCACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCg                                                                                                                               <  2:189291/248‑1 (MQ=255)                                                                                                             agagCGCACGGCGGCTCATGATGTTCCCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCg                                                                                                                               <  1:328968/251‑1 (MQ=255)                                                                                                             agagCCCAGGGCGGCTCATTATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTc                                                                                                                                <  1:279695/250‑1 (MQ=255)                                                                                                               agCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGAt                                                                                                                             <  1:345240/251‑1 (MQ=255)                                                                                                                gcgcACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGGTTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGGGAGGGTAGATATCGGAGGTCACCGGTTTTTTAGGACGATGCTGGAAGTAAAGGGCAGATGCGGCGTTTGCGGctgagtggtggcacaccgcgggggt                                                                                                                             >  2:99547/1‑226 (MQ=255)                                                                                                                         ggcTCCTGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCCGCGGGTGGAAAA‑Caaatatt                                                                                                                  >  1:384217/1‑246 (MQ=255)                                                                                                                             gagtattgccgccgccgGCACGTTTTCGCTGATGGCTTCGTCCAGAAAGCCTGCAGGGGCAATCCCCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCTGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCATGCGGCGTCACCAGCGGGTCGATAA‑CAATAATGcc                                                                                                                <  2:321782/238‑1 (MQ=255)                                                                                                                                gatgTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTAGGTAA‑CAATAAAGCCTTtct                                                                                                           >  1:399164/1‑251 (MQ=255)                                                                                                                                       gccGCCAGCACTTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTACCATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAAGCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAA‑CAATAATGCCTTTCTCGCCCt                                                                                                     <  2:83051/249‑1 (MQ=255)                                                                                                                                         cGCCAGCACGTTTTCGCTGATGGCTTCGTCCATACAGCCTGCAGGGGCAATAACCTGCACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATCGTGGTCGGTGTGGCTGTGAGAGTAGATAACGGCACAAATCGGATTTTGCgggcg                                                                                                                                                                                         >  2:295005/1‑162 (MQ=255)                                                                                                                                         cGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCACAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACg                                                                                                                                                                                         <  1:295005/165‑1 (MQ=255)                                                                                                                                                 cGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCCATAACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTGAAGGGCAAGTGCGGCTTTTGCGGCAGGCGGCGTCCCCAGCGGGTCGGTCA‑CAATAAGGCCCTTCTCGCCCTCACCG‑GACGtc                                                                                         >  1:189723/1‑250 (MQ=255)                                                                                                                                                  gTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGGCGATAA‑AAATAATGCCTTTCTCGccccccacggacgga                                                                                          >  1:571/1‑238 (MQ=255)                                                                                                                                                         cTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGCTGCTGGAAGTAAAGGTCAAGTGCg                                                                                                                                                                >  2:211135/1‑174 (MQ=255)                                                                                                                                                         cTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCg                                                                                                                                                                <  1:211135/174‑1 (MQ=255)                                                                                                                                                                      ccATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAA‑CAATAATGCCTTTCTCGCCCTCAACG‑GACGGAATGTTAGAGATATCCTGGCcg                                                                    >  1:92274/1‑251 (MQ=255)                                                                                                                                                                          aaaGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCCGCGGGTCGGTAA‑AAAAAATGCCTTTCTCGCCCTccac                                                                                                 >  1:57569/1‑216 (MQ=255)                                                                                                                                                                          aaaGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAA‑CAATAATGCCTTTCTCGCCCTCaac                                                                                                 <  2:57569/219‑1 (MQ=255)                                                                                                                                                                              ccTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAA‑CAATAATGCCTTTCTCGCCCTCAACG‑AACGTAATGTTAGAGATATCCTGGCCGCGCACCTg                                                            <  2:274379/251‑1 (MQ=255)                                                                                                                                                                              ccTGCAGGGGCAATCACCTGAACGTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAA‑CAATAATGCCTTTCTCGCCCTCAACG‑AACGTAATGTTAGAGATATCCTGGCCGCGCACCTg                                                            <  2:481047/251‑1 (MQ=255)                                                                                                                                                                                 gCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAA‑CAATAATGCCTTGCTCGCCCTCAACG‑AACGTAATGTTAGAGATATCCTGGCCGCGCACCTGAt                                                          <  1:153178/250‑1 (MQ=255)                                                                                                                                                                                  cAGGGGCAAGCACCGGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCTCGTCTCCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAA‑CAATAATGCCTTTCTCGCCCTCAACG‑AACGTAATGTTAGAGATATCCTGGCCGCGCTCCTGATAc                                                        <  2:290297/251‑1 (MQ=255)                                                                                                                                                                                     gggCAATCACCTGAACGTTGCCGGAGTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAA‑CAATAATGCCTTTCTCGCCCTCAACG‑AACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACAtt                                                     <  2:92274/251‑1 (MQ=255)                                                                                                                                                                                                   aCTTTGCCTGAGTTAACACCGGCTTCAGAGATAATGCCTTTCACGTCACAATAGTGGGCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGTTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAA‑CAATAATGCCTTTCTCGCCCTCAACT‑AACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGAC‑ttt                                      <  2:455893/251‑1 (MQ=255)                                                                                                                                                                                                    cTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGGCCCCAGCGGGGCGAACA‑CAATAATGCCTTTCTCGCCCTCCACT‑CACGGAATGTTTAAAAAATTCTGGCCGCGCACCCGATACAATTTATCggggccgttt                                      >  2:431765/1‑243 (MQ=255)                                                                                                                                                                                                                            cAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAA‑CAATAATGCCTTTCTCGCCCTCAACG‑AACGTAATGTTAGAGATAGCCTGGCCGCGCACCTGATACATTTTATCGGTGAC‑TTTGAACAGGCCAGAAATACCGTTGAt              <  1:308039/250‑1 (MQ=255)                                                                                                                                                                                                                             agagATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGGGGAAGGCGGCGTAACCAGCGGG‑CGGTAAACAATAATGCCTTCTTCGCCCT‑ACCGGAACGTAAGGTTGAAGGTATACTGGCCGGGCACCTGTTACATTTTTTCGGT‑ACATTTAACCAGGCCAGAAATAACGTGGAttt            >  1:337805/1‑251 (MQ=255)                                                                                                                                                                                                                                       ccTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAA‑CAATAATGCCTTTCTCGCCCTCAACG‑AACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGAC‑TTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACg  >  1:372413/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     TAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAA‑CAATAATGCCTTTCTCGCCCTCAACG‑AACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGAC‑TTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACG  >  NC_000913/4305796‑4306276

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