Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,325,837	C→T	intergenic (‑96/+562)	yjdN ← / ← yjdM	PF06983 family protein YjdN/zinc ribbon domain‑containing protein YjdM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,325,837	0	C	T	100.0%	216.8 / NA	73	intergenic (‑96/+562)	yjdN/yjdM	PF06983 family protein YjdN/zinc ribbon domain‑containing protein YjdM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/36);  total (37/36)

GGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCT  >  NC_000913/4325592‑4326082                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      ggAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCAGGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGTGAGCGTCGCATCAGGTAGTCg                                                                                                                                                                                                                                                  <  1:90528/251‑1 (MQ=255)       gCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTg                                                                                                                                                                                                                                            >  2:10269/1‑251 (MQ=255)         aGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGtt                                                                                                                                                                                                                                          >  1:390195/1‑251 (MQ=255)          gTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTg                                                                                                                                                                                                                                         <  2:384819/251‑1 (MQ=255)               cGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCgg                                                                                                                                                                                                                                     <  1:137528/250‑1 (MQ=255)                  cgcTGTCCTGCGCTGAGTTTGGCATGTCGCCGAAGGTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTGGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGAGAATACGCGCGAGCGTAGCATCAGGTAGTCGGCACTGTTGCCGGATgc                                                                                                                                                                                                                                 <  2:280882/251‑1 (MQ=255)                             gcTGATTTTGGCATTTCGCCGAAGGTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTAGACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGc                                                                                                                                                                                                                      <  2:419749/251‑1 (MQ=255)                                 aTTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGAAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCTCGAGCGTCGCATCAGGTAGGCGGCACTGTTTCCAGATGCGGCGTAAACTCCtt                                                                                                                                                                                                                   >  2:16725/1‑250 (MQ=255)                                    ttGGCATTTCGCCGAAGATGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGGGAGGCGGCCTGATAAGACGCGCGAGCGTAGCATCAGGTAGTCGGCACTGTAGCCGGAGGCGGCGTAAACGCCTTAtac                                                                                                                                                                                                               >  2:151743/1‑249 (MQ=255)                                        cATTTCGCCGAAGTTGAGTTTATAGAGCAGCTCCGCGCCCAACGGACGTTGATAATAGGCAATCGCGGCGTAACAGTTACCGGCAAAAGAGAGGTAGTGACGTAACGGCAGGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTCCGAAATTTATACATATGAGAGGCGGTCTGATAAGTCGCTCGAACGGCGCATCAGGTAGTCGGCACTGTTGCCGTATGCGTCTTAAACGCCTTATCCggac                                                                                                                                                                                                           <  2:437641/251‑3 (MQ=255)                                           ttCGCCGAAGTTGATTTTATAGAGCAGTTCCTCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                                                                                         <  2:425497/250‑1 (MQ=255)                                            tCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                                                                                                        <  1:122533/250‑1 (MQ=255)                                             cGCCGAAGTTGATTTTATAGAGCAGTTTCGCGCCCAACGTACGTTGAGAATAGGCCATCGCGTCGGAAATGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGATAAGCCAGTTAAGTGTAGTGCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGTGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGGTTTAACGCCTTATCCGGCCTACg                                                                                                                                                                                                       <  2:410514/250‑1 (MQ=255)                                             cGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACgg                                                                                                                                                                                                      >  2:350923/1‑251 (MQ=255)                                                       gATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCgg                                                                                                                                                                                             <  1:357739/250‑1 (MQ=255)                                                             ataGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAttt                                                                                                                                                                                       <  1:350923/250‑1 (MQ=255)                                                                    aGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGCCCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTCCGGCGGGTGCGGGAATTTGTAGGCCt                                                                                                                                                                               >  1:150345/1‑251 (MQ=255)                                                                      ttCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTg                                                                                                                                                                              >  1:438251/1‑250 (MQ=255)                                                                       tCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCGGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGa                                                                                                                                                                             <  2:183117/250‑1 (MQ=255)                                                                       tCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGATACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTAGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAt                                                                                                                                                                            >  1:191383/1‑251 (MQ=255)                                                                            gcCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAc                                                                                                                                                                       >  1:430116/1‑251 (MQ=255)                                                                               cAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcgcg                                                                                                                                                                    <  1:43201/251‑1 (MQ=255)                                                                               cAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTGAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGGAGTTCGAAATTTATACAGAGGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTAGCCGGTCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcgcg                                                                                                                                                                    <  2:150345/251‑1 (MQ=255)                                                                                   gtacgtTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTGTAGGCCTGATAAGACACGCa                                                                                                                                                                  >  1:287371/1‑249 (MQ=255)                                                                                    tacgtTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGc                                                                                                                                                               >  1:227891/1‑251 (MQ=255)                                                                                     acgtTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCg                                                                                                                                                              >  1:230253/1‑251 (MQ=255)                                                                                     acgtTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCg                                                                                                                                                              <  2:390195/251‑1 (MQ=255)                                                                                          gATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCCGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGaa                                                                                                                                                                                          >  1:383979/1‑218 (MQ=255)                                                                                          gATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGaa                                                                                                                                                                                          <  2:383979/218‑1 (MQ=255)                                                                                                 ggCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCa                                                                                                                                                  <  1:7973/251‑1 (MQ=255)                                                                                                  gCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGGCGCATCAGGc                                                                                                                                                   >  2:323641/1‑249 (MQ=255)                                                                                                   cAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGAGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGAATCAGGTCGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTATAGGCCTGATAAGACGAGCACGCGTCGCATGAGGCAgg                                                                                                                                                >  2:439677/1‑250 (MQ=37)                                                                                                    aaTCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATACGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTc                                                                                                                                               >  2:437706/1‑251 (MQ=255)                                                                                                     aTCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCg                                                                                                                                              <  1:457333/251‑1 (MQ=255)                                                                                                          gTCGAAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCCGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTATTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGGAGCCCTGATAAGCCGCGCAAGCGTCGCATCAGTCAGTCGGCAcc                                                                                                                                         >  2:254004/1‑251 (MQ=25)                                                                                                           tCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                                                                                                                                                                                           >  1:444468/1‑200 (MQ=255)                                                                                                           tCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGAGGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTGCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTATGGAGGGTGCGGGa                                                                                                                                                                                           <  2:444468/200‑1 (MQ=255)                                                                                                            cGGAACAGTTACCGGCAAAAGAGAGGGAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTAGACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCa                                                                                                                                        <  2:387751/250‑1 (MQ=255)                                                                                                              gAACAGTTACCGGCAAAAGCGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGAGTCGCATAAGGTAGTCGGCACTGTTGCCGGACGCGGCGTAAACGCCTTATCCGGCCTAAGGAGGGTGCGGGAATTTGTAGGCCTGATAAGCCGCGACAGCGTCGCATAAGGCAGCCGTCACCATTg                                                                                                                                     >  1:37279/1‑251 (MQ=25)                                                                                                              gAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTg                                                                                                                                     <  2:356500/251‑1 (MQ=255)                                                                                                                      aCCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATg                                                                                                                              <  1:10269/250‑1 (MQ=255)                                                                                                                            aaaaGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGGTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCATGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAAAGTCGCATCAGGCCGGCGGCACCATAGCCgggtgaggagtc                                                                                                                       >  2:246192/1‑241 (MQ=25)                                                                                                                                      taggtacgtAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGGAGTTCGAAATTTATAGAGATGAGAGGCAGCCTGATAAGAAGCGGGAGCGTCGCATAAGGTAGTCGGCACTGTTGACGGATGCGGCGTAAAAGCCTTATCCGGCCTCCGGAGGGTGCGGGAATTGATAGGCATGATAAGACGCGGAAACGGCGCAACAGGATGAAGGCACCATTGCCGGAAGAGGCGGAACCGCCtttt                                                                                                             >  2:261859/9‑249 (MQ=255)                                                                                                                                       aGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                            >  1:110777/1‑251 (MQ=255)                                                                                                                                          gACTTAACGGCATGATGGTGACCGCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTAGCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                         <  2:430116/251‑1 (MQ=255)                                                                                                                                               aaCGGCATGATGGTTACTTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAAGTTATACAGATGAGAGGCGGCCTTATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGGAGGGTGCGGGAATTTGTAGGCCGGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                    <  1:466516/251‑1 (MQ=18)                                                                                                                                               aaCGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                     <  1:323641/250‑1 (MQ=255)                                                                                                                                                   gCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGAGTAGTTCGAAATTTATCCAGATGAGAGGCGGCCTCCGAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTACACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTGTAGGACTGCTAAGACGCGCAAGCGTCGCATCCGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCATATCCGGCCTACgg                                                                                                  >  2:351125/1‑249 (MQ=255)                                                                                                                                                       gatgGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCCGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAggggg                                                                                            >  1:449282/1‑249 (MQ=255)                                                                                                                                                       gatgGTGACCTAAGGTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTAAACAGTTGAGAGGCGGCCTGATACGACGCGGGAGCGGCGCATCAGGTAGTCGGCAGTGTTGCCGGATGAGGCGTAAACGCCTTATCCGGACTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTg                                                                                            <  2:250593/251‑1 (MQ=18)                                                                                                                                                           gTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCggg                                                                                        >  1:235548/1‑251 (MQ=255)                                                                                                                                                               ccTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAtt                                                                                    <  1:14321/251‑1 (MQ=255)                                                                                                                                                                  cAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTg                                                                                  <  1:403089/250‑1 (MQ=255)                                                                                                                                                                  cAGTTAAGAGAAGCCAGTTAAGTGGAGTTCGAAATGTATAAAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATGTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAttt                                                                                   <  2:110777/249‑1 (MQ=255)                                                                                                                                                                   aGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTc                                                                                                                                               >  1:421394/1‑188 (MQ=255)                                                                                                                                                                   aGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTc                                                                                                                                               <  2:421394/188‑1 (MQ=255)                                                                                                                                                                      tAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGc                                                                             >  1:65700/1‑251 (MQ=255)                                                                                                                                                                        agagAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCt                                                                           <  1:320569/251‑1 (MQ=255)                                                                                                                                                                                   ttAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGAATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGTATTTGTAGACCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTCAACGCCTTATCCGGCCTACGGAGGGTGCGGGACTTTGTAGGCCTGATAAAAcgcg                                                                >  2:86698/1‑251 (MQ=17)                                                                                                                                                                                     aaGTGGAGTTAGAAAATTATACGGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCAGTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCa                                                              <  2:227891/251‑1 (MQ=255)                                                                                                                                                                                            gTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGTGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCCGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCg                                                       <  1:151743/251‑1 (MQ=255)                                                                                                                                                                                                gAAATTGATACAGATGAGAGGCGGCCGGATAAGACGCGCGAGCGGCGCATCAGGTAGTCGGCGCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCCCCATTGCCGGATGCGTCGTAAACGCCTTATCCGGCCTACGTAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                                                   <  1:407318/251‑1 (MQ=14)                                                                                                                                                                                                      tataCAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGAAGGCCTGATAAGACGCGCAAGCGTCGCTTCAGGc                                               >  2:245768/1‑249 (MQ=255)                                                                                                                                                                                                          cAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCgg                                         >  1:117545/1‑251 (MQ=255)                                                                                                                                                                                                          cAGATGAAAGGCGGCCTGATAAGACGCGCGAGCGTAGCATCCGGTAGTCGGCACTGTTGACGGATGCGGCGTCAACGCCTTATCCGGCCTACGGAGGGGGCGGGCATTTGTAGGCATGATAAGAAGCGCAAGAGTCGCATCAGGCAGTCGGCACCAGTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGAGGGCATTCGTTGGCCTGATAAGACGAGCCAGTGTCGACCCAGGCAGTCgg                                         >  2:235775/1‑251 (MQ=14)                                                                                                                                                                                                                agagGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGTAGTCGGAACTGTTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGAGGGTGCGGGGATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAATTGTAGGCCCGATAAGACGCGCAAGCGTCGCCTCAGGGAGTCGGCACCGt                                   >  2:85178/1‑251 (MQ=14)                                                                                                                                                                                                                      ggcCTGATAAGACGCGAGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCgg                             >  1:143050/1‑251 (MQ=255)                                                                                                                                                                                                                            aTAAGACGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCCCCGTTGCCGGATgcgg                       >  1:431737/1‑251 (MQ=255)                                                                                                                                                                                                                                 aCGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTaa                  <  1:180859/251‑1 (MQ=255)                                                                                                                                                                                                                                 aCGCGCGAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTa                   >  1:236102/1‑250 (MQ=255)                                                                                                                                                                                                                                  cgcgcgAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACgg                                                                                                  <  2:306150/170‑1 (MQ=38)                                                                                                                                                                                                                                  cgcgcgAGCGTCGCATCAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACgg                                                                                                  >  1:306150/1‑170 (MQ=38)                                                                                                                                                                                                                                                  cAGGTAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt  >  1:284388/1‑250 (MQ=39)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      GGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCT  >  NC_000913/4325592‑4326082

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