Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,244,908	C→T	L61L (CTG→TTG)	ymgE →	PF04226 family protein YmgE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,244,908	0	C	T	100.0%	198.2 / NA	61	L61L (CTG→TTG)	ymgE	PF04226 family protein YmgE
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/33);  total (28/33)

GTTCTTAAATCAGCCGACATACGCCCAGCATCGATCCCGCCCTTGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATCTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCG  >  NC_000913/1244619‑1245202                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                       gttCTTAAATCAGCCGACATACGCCCAGCATCGATCCCGCCCTTGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCtt                                                                                                                                                                                                                                                                                              >  1:161572/1‑300 (MQ=255)       aaaTCAGCCGACATACGCCCAGCATCGATCCCGCCCTTGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCtggt                                                                                                                                                                                                                                                                                        <  1:217501/300‑1 (MQ=255)        aaTCAGCCGACATACGCCCAGCATCAATCCCGCCCTTGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGTCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTCCCtggtg                                                                                                                                                                                                                                                                                       <  1:196560/300‑1 (MQ=255)         aTCAGCCGACATACGCCCAGCATCGATCCCGCCCTTGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGANTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCtggtgg                                                                                                                                                                                                                                                                                      <  1:492919/300‑1 (MQ=255)          tCAGCCGACATACGCCCAGCATCGATCCCGCCCTTGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGc                                                                                                                                                                                                                                                                                     <  1:115075/300‑1 (MQ=255)             gCCGACATACGCCCAGCATCGATCCCGCCCTTGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCggt                                                                                                                                                                                                                                                                                  >  1:203443/1‑300 (MQ=255)             gCCGACATACACCCAGCATCGATCCCGCCCTTGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGCCCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTACTGGTGGCggt                                                                                                                                                                                                                                                                                  >  1:534271/1‑300 (MQ=255)                  cATACGCCCAGCATCGATCCCGCCCTTGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCggtggt                                                                                                                                                                                                                                                                               >  1:362684/1‑298 (MQ=255)                  cATACGCCCAGCATCGATCCCGCCCTTGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCATGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCggtggt                                                                                                                                                                                                                                                                               >  1:103450/1‑298 (MQ=255)                       gCCCAGCATCGATCCCGCCCTTGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCt                                                                                                                                                                                                                                                                        >  2:483226/1‑300 (MQ=255)                            gCATCGATCCCGCCCTTGCCTATACTTAGAGCGTTGCTGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATtct                                                                                                                                                                                                                                                                   >  1:692807/1‑300 (MQ=255)                            gCATCGATCCCGCCCTTGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATtct                                                                                                                                                                                                                                                                   >  1:345542/1‑300 (MQ=255)                                 gATCCCGCCCTTGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGtt                                                                                                                                                                                                                                                               <  1:561326/299‑1 (MQ=255)                                 gATCCCGCCCTGGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCTGCATTATCGCCAAGCTNATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTc                                                                                                                                                                                                                                                              <  2:216680/300‑1 (MQ=255)                                  aTCCCGCCCTTGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCt                                                                                                                                                                                                                                                             >  2:56438/1‑300 (MQ=255)                                            tGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGTTTGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTc                                                                                                                                                                                                                                                   <  2:631140/300‑1 (MQ=255)                                            tGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTc                                                                                                                                                                                                                                                   <  2:431462/300‑1 (MQ=255)                                            tGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTc                                                                                                                                                                                                                                                   >  1:51197/1‑300 (MQ=255)                                              ccTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTcc                                                                                                                                                                                                                                                  >  1:378439/1‑299 (MQ=255)                                                                         tttCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCa                                                                                                                                                                                                                      >  1:466926/1‑300 (MQ=255)                                                                                 gTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCTGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGc                                                                                                                                                                                                                <  2:670323/298‑1 (MQ=255)                                                                                        tCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGa                                                                                                                                                                                                       <  1:615319/300‑1 (MQ=255)                                                                                           tCGAAACGAGTGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGTCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATAc                                                                                                                                                                                                    <  2:163703/300‑1 (MQ=255)                                                                                             gAAACGAGGGCTTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGTTGTTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAg                                                                                                                                                                                                  <  1:478298/300‑1 (MQ=255)                                                                                             gAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAg                                                                                                                                                                                                  <  1:480270/300‑1 (MQ=255)                                                                                                aCGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATg                                                                                                                                                                                               <  1:5385/300‑1 (MQ=255)                                                                                                    gggCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCGGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCg                                                                                                                                                                                            <  2:363823/299‑1 (MQ=255)                                                                                                     ggCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCt                                                                                                                                                                                          >  1:384277/1‑300 (MQ=255)                                                                                                                  gAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTccca                                                                                                                                                                             >  2:587112/1‑297 (MQ=255)                                                                                                                               gATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATTCCTGTGCTTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCTGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGTGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATAc                                                                                                                                                                <  1:685956/300‑1 (MQ=255)                                                                                                                               gATTATTTTCGACCTGATAGCCGGCATTATCGCCAAGCGAATCATGCCTGTGCGGGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCTAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATAc                                                                                                                                                                <  2:205366/300‑1 (MQ=255)                                                                                                                                attattTTCGACCTGACAGCCGGCATTATCGCCAAGCTAACCATGCCTGGGCGTGATGGTTGTGGATTTTTCCTGACCTGTATTCTCTGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCTACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACt                                                                                                                                                               <  2:362684/300‑1 (MQ=255)                                                                                                                                 ttattTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTc                                                                                                                                                              <  1:5848/300‑1 (MQ=255)                                                                                                                                    tttttGACCTGATAGCCTGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGac                                                                                                                                                           <  1:337198/300‑1 (MQ=255)                                                                                                                                    tttttGACCTGATAGCCNGCATTATCGCCAAGCTAANCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATNAGAACGTATACTCGac                                                                                                                                                           <  1:577684/300‑1 (MQ=255)                                                                                                                                                 aGCCGGCATTATCGCCAATCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACgg                                                                                                                                              <  2:399662/300‑1 (MQ=255)                                                                                                                                                  gCCGACATTATCGCCAAGCTAAGCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACggg                                                                                                                                             <  2:307442/300‑1 (MQ=255)                                                                                                                                                     tgCATTATCTCCAAACTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCGGTATTCTCGGGATAGTCGGTGCGTTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTAc                                                                                                                                                <  2:179641/293‑1 (MQ=255)                                                                                                                                                                 aaGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGAACTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTACTTTGCCCAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCAGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGAGTGAAAAAAATCaa                                                                                                                                >  1:218226/1‑298 (MQ=255)                                                                                                                                                                      aaTCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGNTACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCa                                                                                                                         >  1:505319/1‑300 (MQ=255)                                                                                                                                                                      aaTCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCa                                                                                                                         >  2:423486/1‑300 (MQ=255)                                                                                                                                                                      aaTCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTACTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCa                                                                                                                         >  2:331755/1‑300 (MQ=255)                                                                                                                                                                       accATGCCGGGGCGTGATGGTGGTGGATTTTTCCGGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCGGGCTTGCGACCAGGTTTGGCATTGGCGGCTCCATCAGGGGTTTTAATTTGCACAGCTTCCTGGTGGCGGGGGTGGGAGCTATTCTCGTTCTGGGCATATTCCTCCTCCTGCGAAGAGAATAAGATTTTCAGAGGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTAcctgtctcttatacacatctcgcag                                                                                                                        <  1:488354/298‑25 (MQ=255)                                                                                                                                                                                  gCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATcc                                                                                                              >  1:227273/1‑299 (MQ=255)                                                                                                                                                                                               ggATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTACTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTc                                                                                                >  1:377276/1‑300 (MQ=255)                                                                                                                                                                                               ggATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGGGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGTTGGCGGGGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAATGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCGTTACGGGTGAAAAAAATCAATATAGCACTTCGGTATTCGATAAAGTCGGTTc                                                                                                <  2:568169/300‑1 (MQ=255)                                                                                                                                                                                                gATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGAGGTGGTGGGAACTATTCTCGTTCTGGGCATATTCCGCCTACTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATCTAGCACTTCGCTATCCGATAAAGTCGGTTc                                                                                                >  2:265587/1‑299 (MQ=255)                                                                                                                                                                                                          gACCTTTACTCTCGGGATATTCAGCGCGGTGGTCGGCAGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGGGGGAGCGATTCTCGTTTTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATAc                                                                                                                                                                <  2:203302/225‑1 (MQ=255)                                                                                                                                                                                                               gTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAacaacaaca                                                                                 >  2:686953/1‑299 (MQ=255)                                                                                                                                                                                                                 aTTCTCGGGATAGTCGGTGCGGTGGTCGGCGGTTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAg                                                                                <  2:345542/298‑1 (MQ=255)                                                                                                                                                                                                                 aTTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATANGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGCCGGTTCCCGGTAACAACAACAGaa                                                                              >  2:222572/1‑300 (MQ=255)                                                                                                                                                                                                                                                   ggctggcGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGTTGGCGGTGGTGGGAGCTATTCTCGTTCTGTGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAATGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAAGAGAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCt                                                                          <  2:161572/270‑1 (MQ=255)                                                                                                                                                                                                                                                                    ggCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATANAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAaaataaat                           >  1:512241/1‑300 (MQ=255)                                                                                                                                                                                                                                                                     gCATTGGCGGCTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAaataaata                          >  2:653256/1‑300 (MQ=255)                                                                                                                                                                                                                                                                               cTCCATCAGTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACa                 >  2:374750/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                  cATCAGTGGTTTTAATNTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTNTGGTGTGATTCAGAAAATNAATATTTTAGACAGTa              >  1:589822/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                      aGTGGTTTTAATTTGCACATCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGTGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGTTGTGATTCAGAAAATAAATATTTTAGACAGTAGcc           <  2:149212/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                       gttgTTTTAATTTGCACAGCTTCCTGTTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCTAAGAGAATAAGATTTTCATAAGGCGGATAGCGACACAGATGCCGCTATCCGCTTTCATATCAGAACGTATACTCGACACCTGCTTGACNGGTGAAAAAAATCAATATAGCACTGCGATATCNGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTGCTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGt        <  2:534271/297‑1 (MQ=255)                                                                                                                                                                                                                                                                                       gTGGTTTTAATTTGCACAGCTTCCTTGTGTCGGTGGTGGGAACGATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGt        <  2:203443/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                       gTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTCGACAGTAGCCAGt        >  1:273290/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                       gTGGTTTTAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGt        >  1:272978/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                             ttAATTTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCg  <  1:56438/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                       GTTCTTAAATCAGCCGACATACGCCCAGCATCGATCCCGCCCTTGCCTATACTTAGAGCGTTGATGTAAGCATTTCTTGCGTCCGATTCATCGAAACGAGGGCCTGAACATGGGAATTATTGCCTGGATTATTTTTGACCTGATAGCCGGCATTATCGCCAAGCTAATCATGCCGGGGCGTGATGGTGGTGGATTTTTCCTGACCTGTATTCTCGGGATAGTCGGTGCGGTGGTCGGCGGCTGGCTGGCGACCATGTTTGGCATTGGCGGCTCCATCAGTGGTTTTAATCTGCACAGCTTCCTGGTGGCGGTGGTGGGAGCTATTCTCGTTCTGGGCATATTCCGCCTCCTGCGAAGAGAATAAGATTTTCATAAGGCGGATAGCGATACAGATGCCGCTATCCGCTTTCACATCAGAACGTATACTCGACACCTGCTTTACGGGTGAAAAAAATCAATATAGCACTTCGATATCCGATAAAGTCGGTTCCCGGTAACAACAACAGAATCCTTTTTTAATTATTGTTTCGTTGTTTTTGGTGTGATTCAGAAAATAAATATTTTAGACAGTAGCCAGTCAGGCG  >  NC_000913/1244619‑1245202

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