Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,587,365	C→T	S139S (AGC→AGT)	mrr →	Type IV methyl‑directed restriction enzyme EcoKMrr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,587,365	0	C	T	100.0%	145.3 / NA	48	S139S (AGC→AGT)	mrr	Type IV methyl‑directed restriction enzyme EcoKMrr
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/29);  total (19/29)

TCGTGCAGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGC  >  NC_000913/4587119‑4587607                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   tCGTGCAGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCgg                                                                                                                                                                                                                                                >  1:201968/1‑251 (MQ=255)       aGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATc                                                                                                                                                                                                                                          <  1:407132/251‑1 (MQ=255)       aGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATc                                                                                                                                                                                                                                          >  1:260704/1‑251 (MQ=255)            gtgcgcATGACCCTTTAAAGCGGGTCGGGTTTTCGCAAAGTGTGTCGCGCGGCAAATGGTTCCTTATTCCTGCGGGTTTTGACTGGTTTTCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCCCGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAATAGCAGTCCGGACGATCGGt                                                                                                                                                                                                                                       <  2:524632/247‑1 (MQ=255)            ggtcgcttGACGGTTTAAAACGGGCCGGGTGTTCGCAAAGTTTGTCGCTGGGCAAATGGTGCCTGACTCCTGCGGGGTTTTTCCGGGTTGCGT‑CTCACCTCCAGCCAAGGACTGCGCTGGCGACGAACCATCTTGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGAGGCCGTCGATTTAGATCCGAAAGCAGACCCTTCCGATCATGAAGAACTTGCAAAGAGTAGTCCGGACGATCGGtg                                                                                                                                                                                                                                      <  2:112967/243‑2 (MQ=255)                       cGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGTTCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGTACGATCGGTTAGATCAGgcggt                                                                                                                                                                                                                          >  2:426480/1‑249 (MQ=255)                                  ggCCGGGTTGGCGCAAAGTTTGTCTCGTGGCAAATGGTGCCTGACACCTGCGGGTTTTGACTGGGTTGCGT‑CTCACCCCCAGCCAATGACGGAGCAGGAGACGAACCATCGGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATGCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTc                                                                                                                                                                                                               <  1:531344/250‑1 (MQ=255)                                        tttGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTGTTGTGCATGGCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCTTGATg                                                                                                                                                                                                         >  2:483678/2‑251 (MQ=255)                                              gCAAAGTTTGTCGCGTGGCAAATTGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTgg                                                                                                                                                                                                   <  1:55825/251‑1 (MQ=255)                                                  aGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCtga                                                                                                                                                                                               <  1:483678/251‑1 (MQ=255)                                                              ggcaaggggtgcctatctctctcGGGTTTTGACTGGGTTGCGT‑CTCATCCCCAGCCTATGACGGAGCAGGAGACGGACCATTTTGCCTTCGGTTTTGTGAATGTCAATCTTAAGTCACGGCCGGATGCCGTTGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCTGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCTTGTTGCGGTGGCCGATGAGGTTCTGGa                                                                                                                                                                                   <  2:249284/229‑1 (MQ=255)                                                               gCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGaa                                                                                                                                                                                  <  1:386950/251‑1 (MQ=255)                                                                     ggTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTa                                                                                                                                                                            <  1:395696/251‑1 (MQ=255)                                                                      gggCCTGACTCCTGCGGGTTTTGACTGGGTTGCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTAt                                                                                                                                                                           <  2:292379/249‑1 (MQ=255)                                                                      gTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGCTGAGGTTCTGGAAAACTTAt                                                                                                                                                                           >  2:116326/1‑251 (MQ=255)                                                                            gACTCCTGCGGGTTTTGACTGGGTTGCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAgg                                                                                                                                                                     <  1:347898/251‑1 (MQ=255)                                                                                ccggCGGTTTTTGCCTGGGTTG‑GTGCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCTGCTGGCCTTTGCGTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGAGGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTtc                                                                                                                                                                 <  2:340587/248‑1 (MQ=255)                                                                                                gggTTGCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGaa                                                                                                                                                 <  1:95159/251‑1 (MQ=255)                                                                                                            cATCCCCAGCCAATGACTGAGCAGTAGACGAACCATATGGCCTTCTCTTTTGTGAATGTCCAACTTAAGTCACGGCCGGATGCCGTCGATTTAGAGCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTATAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGTGCTTTGAAGTCATTGt                                                                                                                                         <  2:378543/248‑1 (MQ=255)                                                                                                                 ccAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATg                                                                                                                                 >  1:302813/1‑251 (MQ=255)                                                                                                                          gACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGTTGTTTTGCATc                                                                                                                        <  2:137262/251‑1 (MQ=255)                                                                                                                          gACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCAt                                                                                                                         <  2:302813/250‑1 (MQ=255)                                                                                                                              gAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCt                                                                                                                    >  1:256232/1‑251 (MQ=255)                                                                                                                                 cAGGAGACGAACCATGTGGCCTTCGCTGTTGTGAATGGCAATCTTAGGTCACGGCAGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCGCGATCAGGAAGAACTGGAAAAGAGCAGTCCGGACTATCGGTTAGAACAGGCGCTAAAGGAGCTTCGGGATGCGGTTGCTGATGAGGGTCTTGAAAACTTATTGCAGGTTTGTCCGTCGCGCTTTGAAGGCATTGTTCTGTATGTTTTGCATCGCCTggg                                                                                                                 <  1:39772/251‑1 (MQ=255)                                                                                                                                        cGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATggc                                                                                                          <  1:116326/251‑1 (MQ=255)                                                                                                                                         gAACCATCTGGCCTTCGCTTTTGTGACTGTCAAACTTAAGTCAAGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATggcg                                                                                                         >  1:223989/1‑251 (MQ=255)                                                                                                                                                                                 tCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGCAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGGAATGgagg                                                                 >  1:550287/1‑250 (MQ=255)                                                                                                                                                                                  cACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTgg                                                                    >  1:134576/1‑247 (MQ=255)                                                                                                                                                                                      gCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGAtggtg                                                            >  1:272193/1‑251 (MQ=255)                                                                                                                                                                                           aagCCGTCGATTTAGAGCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATc                                                       <  2:201968/249‑1 (MQ=255)                                                                                                                                                                                                  cGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGgtg                                                <  1:388313/251‑1 (MQ=255)                                                                                                                                                                                                   gATTTAGATACGACAGCCGACTCTCCCGATCCTGAAGAACTTGCAAAGAGCAGTCCGGACGCTCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGAATGGCGGCCACCGAGCTTATTTGCAGAGTTGTGGCGGTTCAGGGGGTGGTGGCacggtaggtgg                                               >  2:445291/1‑241 (MQ=255)                                                                                                                                                                                                      ttAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCg                                                                                   <  2:324384/212‑1 (MQ=255)                                                                                                                                                                                                      ttAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCg                                                                                   >  1:324384/1‑212 (MQ=255)                                                                                                                                                                                                             cGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGGGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCtt                                     >  1:359247/1‑251 (MQ=255)                                                                                                                                                                                                                          tcCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTa                                                                                                              <  2:193116/165‑1 (MQ=255)                                                                                                                                                                                                                          tcCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTa                                                                                                              >  1:193116/1‑165 (MQ=255)                                                                                                                                                                                                                                 cATGAAGAACTTGCAAAGAGCAGTTCGGACGAGCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGgaga                  <  2:30867/250‑1 (MQ=255)                                                                                                                                                                                                                                 cATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGga                                                                   <  1:145928/201‑1 (MQ=255)                                                                                                                                                                                                                                 cATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGAGGCATCGGTGGGGTGATATCGCTTGATAAAATTGGCCTGgaga                  >  1:271764/1‑250 (MQ=255)                                                                                                                                                                                                                                    gaagaaCTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGt              <  1:426480/251‑1 (MQ=255)                                                                                                                                                                                                                                            tGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGt        <  1:554980/249‑1 (MQ=255)                                                                                                                                                                                                                                            tGCAAAGAGCAGTCAGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACGTATTGCAGGTTTCTCCTTCGCGCTTGTAAGCCATTGTGCTGGATGATTTGCAGCGACTGGGGTAGGGCGGCGACCTTGACGATGTGAAGCTTGTTGGCGGTACTGGAGATGCTGGCATCGATGCTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGtt       <  1:59125/250‑1 (MQ=255)                                                                                                                                                                                                                                              cAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCt                                      <  2:411273/217‑1 (MQ=255)                                                                                                                                                                                                                                              cAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGGGTGATATCGCt                                      >  1:411273/1‑217 (MQ=255)                                                                                                                                                                                                                                                aaGAGGAGTCCGGACGAGTGGTTAGATCAGGCACTAAAAGAGCGTCGTGAATCGGTGGATGATGAGGTTCTTGAAAACTTATTGCAGGTTTCTCCTTCGCGATGTGAAGTCATTGTTCTGGATTTTTTGCATCGCCTTGGGTATGGCGGCCACCGTGATGATTTGCATCGTTTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCTCTTGATAAACTTGGCCTGGACAAAGTTTATGTTCAGGc  <  2:244790/251‑1 (MQ=255)                                                                                                                                                                                                                                                aaGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGGGTGATATCGCTTGATAAACTTGGCCTGGAGGAAGTTTATGTTCAggg  >  1:190311/1‑250 (MQ=255)                                                                                                                                                                                                                                                 agagCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGGGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGGTAAACTTGGCCTGGGGAAAGTTTATGTTCAGGc  >  1:427095/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   TCGTGCAGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGT‑CTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGC  >  NC_000913/4587119‑4587607

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