Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,457,392	C→T	E158E (GAG→GAA)	yjgA ←	putative ribosome biogenesis factor YjgA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,457,392	0	C	T	100.0%	244.4 / NA	72	E158E (GAG→GAA)	yjgA	putative ribosome biogenesis factor YjgA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/37);  total (35/37)

AAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTCTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCT  >  NC_000913/4457145‑4457638                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       aaGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTtctttttct                                                                                                                                                                                                                                                     <  1:515688/251‑1 (MQ=255)    ccGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGAGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGAGCGGATTTAGGCGGCTTATTCCCttctttttctttt                                                                                                                                                                                                                                                  >  2:280027/1‑250 (MQ=255)    ccGCGCAGCAATTCGGCCAGAGCCAGATATAGCATGCGCAGGCCGGTTAAGTCGTTCACGCCGAATCCGGCATTTGCGCCTAGTGCCTGATGCGACGCTTGATTCGTCTTATCAGGCCTACAGCTGCAACGTATACGGAAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGATAACACGCGCAGATCCTGGCCAATCTGTCGGGCTGAATTCGGCTGATTATTCCCTtatttttctttt                                                                                                                                                                                                                                                  >  2:64825/1‑250 (MQ=255)      gcgcAGTAATTCGGCCTTAGCTAGAGATAGTATACGTATGCCGGATAAGGCTTTCACGCCGCATCCGTCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTTCATGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCTGCTTATTCCCTTCTTTTtctttctt                                                                                                                                                                                                                                                <  1:153315/251‑1 (MQ=255)       cgcAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTttcttc                                                                                                                                                                                                                                               <  1:426802/251‑1 (MQ=255)        gcAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTcg                                                                                                                                                                                                                                              <  1:326414/251‑1 (MQ=255)                tCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCgg                                                                                                                                                                                                                                      >  2:249759/1‑251 (MQ=255)                 cGGCCTTAGCCAGATATAGTAGGGGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAAGTACCCGTCTTTCTCCGGTAACTCGCTCAGATACTGGAAAAGCTGGCGTGCGGATTGCGGCGGCTTATTCCCTTCTTTTTCTGTCTTCGCGTTGCGGa                                                                                                                                                                                                                                     <  2:303154/251‑1 (MQ=38)                  ggCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGAt                                                                                                                                                                                                                                    >  1:413696/1‑251 (MQ=255)                    ccTTGGCCTGATATTGTATGCGGAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCGCATCAACAACGAAAATTTCCCTTCGTGCTCCGCTAACGCGCGCAGATGCTGGAAAATCTGGCGTGCGGATTGCGGCGGCTTATTCGCTTGTGTTTCTTTCTTCTCGTTGCGGATCa                                                                                                                                                                                                                                  <  2:262956/251‑1 (MQ=37)                    ccTTAGCCAGATAGAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCa                                                                                                                                                                                                                                  <  1:206409/251‑1 (MQ=255)                         gCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGt                                                                                                                                                                                                                              >  1:210714/1‑250 (MQ=255)                           cAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGt                                                                                                                                                                                                                              >  1:239708/1‑248 (MQ=255)                           cAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGt                                                                                                                                                                                                                              <  2:239708/248‑1 (MQ=255)                              atatAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGc                                                                                                                                                                                                                          >  1:65780/1‑249 (MQ=255)                                ataGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCAGCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGt                                                                                                                                                                                                                       <  1:40440/250‑1 (MQ=255)                                  aGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGc                                                                                                                                                                                                                    >  1:347021/1‑251 (MQ=255)                                    tATGCGTAGGCCGGAGAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTTATGCGTCTTATCAGGCTTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTGCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGtggtt                                                                                                                                                                                                                   <  2:434320/250‑5 (MQ=255)                                     aTGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTg                                                                                                                                                                                                                  <  1:568761/250‑1 (MQ=255)                                        cGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATACCGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACg                                                                                                                                                                                                               >  1:238200/1‑250 (MQ=255)                                           aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAACTTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGATGATCg                                                                                                                                                                                                           >  2:438038/1‑251 (MQ=255)                                           aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTGGCTGACGATc                                                                                                                                                                                                            <  2:562314/250‑1 (MQ=255)                                              ccGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGTTGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGTCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGGGCTGATTTCGGCGGCTTATTCCCTTCTTTTTCGTTCTTCGCGTGGCGGATCAGAGTACTCAGTTGCTGTCTATCGGCa                                                                                                                                                                                                        <  2:79790/251‑1 (MQ=255)                                               cGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGCTGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCAt                                                                                                                                                                                                       >  1:81793/1‑251 (MQ=255)                                                ggATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATc                                                                                                                                                                                                      >  1:249451/1‑251 (MQ=255)                                                 gATAAGGCGTTCACGCCGCATCCGGAATTTGCGCATAGTGCCTGATGCGACGATTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCg                                                                                                                                                                                                     >  2:333870/1‑251 (MQ=255)                                                   tAAGGCGTACACGCCGCATCCGGCATTTCAGCATCGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTCCAGGGGCCCCGTATCCGGCAAACCATACTCACATCCACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGAGGCTTATTCCCTTATTTTTCTTTCTTCGCGTTGCGGCTCAGAGTACGCAGTTGCTGACGTTCTGCATCggg                                                                                                                                                                                                   >  2:375246/1‑251 (MQ=255)                                                    aaGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGGCGCTTTATGCGTATTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACACTTCGTTCTCCGCTAACTCGCGCAGCTACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCTAGTTGCGGATCAGAGTACGCAGTTGCTGCAGATAGGCCTCGGGc                                                                                                                                                                                                  >  2:330325/1‑251 (MQ=255)                                                        cGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCcaca                                                                                                                                                                                              <  1:397546/251‑1 (MQ=255)                                                         cttcacgagccctCCGGTATTTGAGCATAGCACCTGATGCGCCGCTTGATGCGTCTTATCAGGCCTACAGGTGCGCCGACCCCCGCAAACCATACTCCCATCACCAACGCCCATTCCCCTTCGTTCTCCGCTAACTCGCGCAGATCCTGGCCAATCTGGCGTGCGGATTTCACCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGCGTCAGAGTGCGCCGTTGCTGACGAGCGGCATCTGGCcacac                                                                                                                                                                                             >  2:244391/13‑250 (MQ=255)                                                           tCACGCCGCATCCGTCATTTGAGCATAGTGCCTGATTCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCGCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGa                                                                                                                                                                                            <  2:318209/250‑1 (MQ=255)                                                                 cgcATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCCTCGGGCCACAGATTTAAt                                                                                                                                                                                      >  2:478578/1‑250 (MQ=255)                                                                  gcATCGGGTATTTGAGCATAGTGCCTGATGCGGCGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCATTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGAGCGGCATCGGGCCACAGATTTAATAc                                                                                                                                                                                    <  2:446689/251‑1 (MQ=255)                                                                  gcATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATAc                                                                                                                                                                                    <  1:201874/251‑1 (MQ=255)                                                                   cATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGCTCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATAcc                                                                                                                                                                                   >  1:117056/1‑251 (MQ=255)                                                                     tCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTc                                                                                                                                                                                 >  1:117086/1‑251 (MQ=255)                                                                                aGCATAGTGCCTGCTGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGTCGAt                                                                                                                                                                                                             >  2:429394/1‑212 (MQ=255)                                                                                aGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGAt                                                                                                                                                                                                             <  1:429394/212‑1 (MQ=255)                                                                                     aGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGtcatca                                                                                                                                                                 <  1:203070/251‑1 (MQ=255)                                                                                     aGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACAACAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGtcatc                                                                                                                                                                  <  2:104658/250‑1 (MQ=255)                                                                                      gTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGAGGTCGGCGGCTTAGTCCCTTCTTTTTCTTTCTTCGCGTTGCTGATCAGAGTACGCAGTTGCGGACGATCGGCATCGGGCCACAGATTTAAGACCTCAGCGAAGGCGCCATCAc                                                                                                                                                                <  2:63910/251‑1 (MQ=255)                                                                                         ccTGATGCGACGATTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCCAACCATACTCACATCAACAACGCAAATTACCCTTCGTTCTCAGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCTCGTTGCGGATAAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCCCAGCGATGGCGCCATCACCtt                                                                                                                                                             >  2:106061/1‑251 (MQ=255)                                                                                           tGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTCTTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGCGTACGCAGTTGCTGACGATCGGCATCGGGCAACAGATTTAATACCTCAGCGATGGCGCCATCACCTTg                                                                                                                                                            >  2:461385/1‑250 (MQ=255)                                                                                              tGCGACGCTTGATGCGTCTTATCAGGCCTACAGGCGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTgatcg                                                                                                                                                        >  2:220033/1‑251 (MQ=255)                                                                                                       tGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGAc                                                                                                                                                 <  1:348362/249‑1 (MQ=255)                                                                                                           gCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATcgc                                                                                                                                           >  1:251617/1‑251 (MQ=255)                                                                                                                    cAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGtttt                                                                                                                                  <  1:420332/251‑1 (MQ=255)                                                                                                                        ccTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAg                                                                                                                              <  2:347021/251‑1 (MQ=255)                                                                                                                                 gCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGaa                                                                                                                     <  1:371447/251‑1 (MQ=255)                                                                                                                                   cccGTATCCGGCAAACCATTCTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAAGCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGAGTTAATACCTCAGCGATGGCGTCATCTCCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACa                                                                                                                   <  1:366174/250‑1 (MQ=255)                                                                                                                                   aCCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACa                                                                                                                   <  1:434297/251‑1 (MQ=255)                                                                                                                                          ccGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGcaccac                                                                                                            <  1:165938/251‑1 (MQ=255)                                                                                                                                            ggCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACct                                                                                                          <  1:500237/251‑1 (MQ=255)                                                                                                                                                      aCTCACATCAACAACGAAACTTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTtgtg                                                                                                >  1:8813/1‑251 (MQ=255)                                                                                                                                                                 caacGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGCTACTGGCAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTcg                                                                                     >  2:53461/1‑250 (MQ=255)                                                                                                                                                                  aacGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATCCTGGAAAATCTGGCGTGCTGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACCGATTTAATACCTCAGCGGCGGCGTCATCACCTTGATCGCTCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTc                                                                                      >  2:42909/1‑249 (MQ=255)                                                                                                                                                                   acGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGa                                                                                                                      >  1:339806/1‑216 (MQ=255)                                                                                                                                                                   acGAAAATTACCCTTCATTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCTCCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGa                                                                                                                      <  2:339806/216‑1 (MQ=255)                                                                                                                                                                             cccTTCGGTCTCCGCTAACTCGCTCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTTATCGATCAGACGATCGCGCAGGTGTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGGGACGGTTTTGCAGCTTATCCAgcg                                                                         <  2:222571/251‑1 (MQ=255)                                                                                                                                                                              ccTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGAGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGGTGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAgcgc                                                                        >  1:251937/1‑251 (MQ=255)                                                                                                                                                                                   gTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGACCAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGAc                                                                   >  1:111500/1‑251 (MQ=255)                                                                                                                                                                                     tctcCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACg                                                                  >  2:203849/1‑250 (MQ=255)                                                                                                                                                                                                    gcAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCATTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATAGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTcc                                                  >  1:205500/1‑250 (MQ=255)                                                                                                                                                                                                      aGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCGTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACTAATAGGCTCGTCGTcgcg                                                <  1:229505/251‑1 (MQ=255)                                                                                                                                                                                                        aTACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCt                                              <  1:206691/251‑1 (MQ=255)                                                                                                                                                                                                               aaaaTCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCgcag                                       <  1:478578/251‑1 (MQ=255)                                                                                                                                                                                                               aaaaGCTGGCGTGCGGATTTCGGCGGCTTATGCCCATCTTTTTCTGTCTACGCGTTGCGGATCAGAGTACGCAGTTGTTGACGATCGGCGTCGGTCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTGCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGAAGGTGTTTCAGCTTATCCAGCGGCTGACGAATAGGCTCTACGTCGCGCTGGCgcag                                       <  2:339707/251‑1 (MQ=255)                                                                                                                                                                                                                        cGTGCGGCTTTCGGCGGCTTATTCCCTTATTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCa                                                                                                                                        >  2:583010/1‑145 (MQ=255)                                                                                                                                                                                                                        cGTGCGGATTTCGGCGGCTTATTCCCTTCTTGTTCTTTCTGCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCTATCAGACGATCGCGCa                                                                                                                                        <  1:583010/145‑1 (MQ=255)                                                                                                                                                                                                                                 ttCGGCGGCTTATTCCCTTCTTTTGCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGAGGGCGGCATCACCTTGAGCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCGGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCGGACGAATAGGCTCTACGTCGCGCTGGCGCATCATCTTACCGATGAGCt                      <  2:549957/250‑1 (MQ=255)                                                                                                                                                                                                                                      cggcTTTTTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTTCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTCCCGATGAGCTGCAACt                >  2:382199/1‑251 (MQ=255)                                                                                                                                                                                                                                      cggcTTATTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACt                >  1:364322/1‑251 (MQ=255)                                                                                                                                                                                                                                            aTTCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTggcggc          >  1:518774/1‑251 (MQ=255)                                                                                                                                                                                                                                            aTTCCCTTCGTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTggcgg           <  2:413696/250‑1 (MQ=255)                                                                                                                                                                                                                                             ttCCCTTCTTTTTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCCGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTggcggcg         >  1:85161/1‑251 (MQ=255)                                                                                                                                                                                                                                                    ctttttctttctTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCt  >  2:183081/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       AAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTCTCTTTCTTCGCGTTGCGGATCAGAGTACGCAGTTGCTGACGATCGGCATCGGGCCACAGATTTAATACCTCAGCGATGGCGTCATCACCTTGATCGATCAGACGATCGCGCAGGTTTTCCAGTTTATGGAACAGCACCACCTGCTGGTTGTGACGGTTTTTCAGCTTATCCAGCGCCTGACGAATAGGCTCTACGTCGCGCTGGCGCAGCATCTTACCGATGAGCTGCAACTGGCGGCGGCGACCT  >  NC_000913/4457145‑4457638

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