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.9 / NA	45	S139S (AGC→AGT)	mrr	Type IV methyl‑directed restriction enzyme EcoKMrr
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/22);  total (23/22)

AGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGC  >  NC_000913/4587125‑4587607                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   catccagaacaaTGACCGTTTAAACCGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGAGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGACGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATc                                                                                                                                                                                                                                          >  2:52061/12‑251 (MQ=255)                        aaCGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGCAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAag                                                                                                                                                                                                                     >  2:237003/1‑249 (MQ=255)                         aCGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGc                                                                                                                                                                                                                  >  1:250120/1‑251 (MQ=255)                                                ggTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTGTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAg                                                                                                                                                                                           <  2:199347/250‑1 (MQ=255)                                                      gTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTg                                                                                                                                                                                     >  1:111207/1‑251 (MQ=255)                                                      gTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTg                                                                                                                                                                                     <  2:137082/251‑1 (MQ=255)                                                      gTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTg                                                                                                                                                                                     <  1:65390/251‑1 (MQ=255)                                                            aaTGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAAc                                                                                                                                                                               <  2:81381/251‑1 (MQ=255)                                                                   ccTGACTCCTGCGGGTTTTGACTGGGTTGCGTCGCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCGTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGc                                                                                                                                                                        <  1:70863/251‑1 (MQ=255)                                                                             gcgggtttttaccggactgcacctcATCCCCAGCCAATGACGGCGCAGGAGACGAACCATATGACCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATACGAAAGCCGACTCTCACGATCATGAAGAACTTGCACAGAGCAGTCCGGCAGCTTGGTTCGACCAAGAGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAGAACTTATTGAAGGTTtctac                                                                                                                                                              >  2:386328/23‑249 (MQ=255)                                                                                           gtgtgcgtttcgcccccAGCCAATGACGGAGCAGGAGACGAACCATCTGGCTGTCGCTTTTGTGAATGTCAAACTTAAGGCACGGCCGGATGCCGTCGATTTAGAACCGAAAGCCGACTCTCCCGATCATGAAGAACGTGCAAAGAGCTGTCCGGACGATCGGTTAGAGCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAg                                                                                                                                                <  1:308127/238‑1 (MQ=255)                                                                                            gTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGCTCAGGCGCTAAAAGCGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCtttaacttctcgcc                                                                                                                                                >  2:336022/1‑238 (MQ=255)                                                                                               gcggctcATCCCCATCCAATGACGGAGCAGGAGACGAACCATCTGGCCTGCGCTTTTGTGAATGTCAAACTTAAGTCCCGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCTGTGGCTGATGAGGTTCTGGAAATCTTAGTGCAGGTTTCTCCTTCGCGCTTTGAAGTCa                                                                                                                                             <  2:44368/246‑1 (MQ=255)                                                                                                 gTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGGCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTg                                                                                                                                          <  2:39681/251‑1 (MQ=255)                                                                                                  tctcctgctccgtcAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACGTATTGCAGGTTTCTCCTTCGCGCt                                                                                                                                                      <  2:224829/225‑1 (MQ=255)                                                                                                  tctcctgctccgtcAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGATAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTATGGAAAACTTATTGCAGGTTTCTCCTTCGCGCt                                                                                                                                                      >  1:224829/14‑238 (MQ=255)                                                                                                     cATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGt                                                                                                                                         <  2:433733/248‑1 (MQ=255)                                                                                                          ccAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATg                                                                                                                                 <  1:52061/251‑1 (MQ=255)                                                                                                          ccAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCCGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCACGATCTTGATGCACTTGCATAGAGCAGTCCGGACGCTCGGTTAGATCAGGCGCTACAAGCGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGCATg                                                                                                                                 >  1:228230/1‑251 (MQ=255)                                                                                                                        aGCAGGAGACGAACCATCTGGCCTTCGCTGTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGACGTTTTGCTTCGCCTg                                                                                                                   <  2:449159/251‑1 (MQ=255)                                                                                                                              agaCGAACCCTCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTAt                                                                                                             >  1:83575/1‑251 (MQ=255)                                                                                                                                   aaCCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATggcgg                                                                                                        <  1:398876/251‑1 (MQ=255)                                                                                                                                             ccTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGCAAGCCGACTCTACCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCt                                                                                                                    >  2:72751/1‑229 (MQ=255)                                                                                                                                             ccTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCt                                                                                                                    <  1:72751/229‑1 (MQ=255)                                                                                                                                                tCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTgatgat                                                                                           >  1:348140/1‑251 (MQ=255)                                                                                                                                                 cGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCCCCGTGATGAtt                                                                                          >  1:317385/1‑251 (MQ=255)                                                                                                                                                     tttGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCCGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGCTGATTTGc                                                                                       >  1:450792/1‑250 (MQ=255)                                                                                                                                                            aTGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGtt                                                                               >  1:306549/1‑251 (MQ=255)                                                                                                                                                                cAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTGGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCAGCGCCTGGGGTATGGCGGCCTCCGTGATGATTGGCAGCGTGTTGGc                                                                            <  2:80773/250‑1 (MQ=255)                                                                                                                                                                 aaaCGTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCgg                                                                          <  2:30014/251‑1 (MQ=255)                                                                                                                                                                     ttAAGTCACGGCCGGCTGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGCCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGCTGATTTGCAGAGTGTTGGCGGtc                                                                        >  2:386885/1‑248 (MQ=255)                                                                                                                                                                            aCGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGAGGATTTTCAGCGTGTTGGCGGGACTGgaggtg                                                               >  2:371261/1‑248 (MQ=255)                                                                                                                                                                               gCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGAtggtg                                                            <  1:260420/251‑1 (MQ=255)                                                                                                                                                                                ccGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGAtggtgg                                                           <  1:264016/251‑1 (MQ=255)                                                                                                                                                                                  ggATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGTCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCa                                                         <  1:386885/251‑1 (MQ=255)                                                                                                                                                                                      gCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGAGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGGGGTGGTGGCATCgt                                                     >  1:427860/1‑250 (MQ=255)                                                                                                                                                                                         gTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATgg                                                  >  1:236622/1‑251 (MQ=255)                                                                                                                                                                                                 aGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGatat                                          <  2:376696/251‑1 (MQ=255)                                                                                                                                                                                                         aaGCCGACTCTACCGATCATGAAGAACATGCAAAGAGCAGTCCGGACGAGCGGTTAGATCAGGCGCTAAAAGAGCTGCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTAGTGCAGGTTTCTCCTGCGCGCGTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCGGGGGTGGGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCTAGGGTGTGATATCGCTTGAt                                  <  2:318412/251‑1 (MQ=255)                                                                                                                                                                                                             cGACTCTCCCGCACATGAAGACCTTGCAAAGAGCAGTCAGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGCTGGGGTTCTGGAACACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGCTGATTTGCAGCGTGTTTGCGGTAGTGGGGTTGGGGGAATTGATGGGGTGATATCGGTTggtcaa                               >  2:243539/1‑245 (MQ=255)                                                                                                                                                                                                                    cccGATCATGAAGAACTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTgg                          <  1:402370/248‑1 (MQ=255)                                                                                                                                                                                                                          cATGAAGAACTTGCAAAGAGCAGTCCGGACGATGGGTTAGATAAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGGTGGTGGCATCGATGGTGTGATTTCGCTTGGTAAACTTGGGCTGgaga                  >  2:434467/1‑250 (MQ=255)                                                                                                                                                                                                                            tGAAGAACTTGCAAAGAGCCGTCCGGACGATCGGTTCGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCATTCGCGCTTTGAAGTCATTGTTATGGATGTTTTGCATCGCCTGGGGTATGGCGGCCAACGTGATGATTTGCAGCGTGTTGGCGGTCCTGGAGATGGAGGCATCGTTGGGGTGATATCGCTTGATAAAATTGTCCTGGAGaaat               >  2:185123/1‑250 (MQ=255)                                                                                                                                                                                                                                   cTTGCAAAGAGCAGTCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGGACTGGAGATGGTGGCATCGGTGGTGTGGTATCGCTTGATAAACTTGGCCTGGGGAAAGTTTATGt        >  2:320251/1‑251 (MQ=255)                                                                                                                                                                                                                                         aaGAGCAGTCCGGCCGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGCTGATTTGCAGCGTGTTGGCGGGACTGGTGATGGGGGCATCGAGGGTGTGTTATCGCTTGGTAGACTTGGCCTGGGGAAAGTTTAttgtccggc  >  2:161612/1‑243 (MQ=255)                                                                                                                                                                                                                                            agcagtCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGAt                                                    <  2:123105/198‑1 (MQ=255)                                                                                                                                                                                                                                            agcagtCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGAt                                                    >  1:123105/1‑198 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   AGGCTGGGCGCATGACCGTTTAAAACGTGCCGGGTTGTCGCAAAGTTTGTCGCGTGGCAAATGGTGCCTGACTCCTGCGGGTTTTGACTGGGTTGCGTCTCATCCCCAGCCAATGACGGAGCAGGAGACGAACCATCTGGCCTTCGCTTTTGTGAATGTCAAACTTAAGTCACGGCCGGATGCCGTCGATTTAGATCCGAAAGCCGACTCTCCCGATCATGAAGAACTTGCAAAGAGCAGCCCGGACGATCGGTTAGATCAGGCGCTAAAAGAGCTTCGTGATGCGGTGGCTGATGAGGTTCTGGAAAACTTATTGCAGGTTTCTCCTTCGCGCTTTGAAGTCATTGTTCTGGATGTTTTGCATCGCCTGGGGTATGGCGGCCACCGTGATGATTTGCAGCGTGTTGGCGGTACTGGAGATGGTGGCATCGATGGTGTGATATCGCTTGATAAACTTGGCCTGGAGAAAGTTTATGTTCAGGC  >  NC_000913/4587125‑4587607

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