Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,862,227	C→T	G82G (GGG→GGA)	glvB ←	putative PTS enzyme II component GlvB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,227	0	C	T	100.0%	243.9 / NA	72	G82G (GGG→GGA)	glvB	putative PTS enzyme II component GlvB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/38);  total (34/38)

TGGTCATTATACTGCCTCCGTAATGGCAACATTTTCTGCTGATTGATGAGAATTAATTAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACCCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAAC  >  NC_000913/3861981‑3862470                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    tGGTCATTATACTGCCTCCGTAATGGCAACATTTTCTGCTGATTGATGAGAATTAATTAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCg                                                                                                                                                                                                                                                  <  2:566868/250‑1 (MQ=255) tGGTCATTATACTGCCTCCGTAATGGCAACATTTTCTGCTGATTGATGAGAATTAATTAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCg                                                                                                                                                                                                                                                  <  1:499262/250‑1 (MQ=255)     cATTATACTGCCTCCGTAATGGCAACATTTTCTGCTGATTGATGAGAATTAATTAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGccc                                                                                                                                                                                                                                               >  1:589848/1‑249 (MQ=255)       ttATACTGCCTCCGTAATGGCAACATTTTCTGCTAATTGATGAGAATTAATTAAGCTATCGAGCTGTTCAAGCAGCTGAGATACATGCAGACCAATTATCAGCTGAATGGCATCGCCACTACGGCAGACGCCGTGCGATACCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCCGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGACGACTCCGCCCAggg                                                                                                                                                                                                                                           >  1:253870/1‑251 (MQ=255)         ataCTGCCTCCGTAATGGCAACATTTTCTGCTGATTGATGAGAATTAATTAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCt                                                                                                                                                                                                                                         <  2:354159/251‑1 (MQ=255)               ccTCCGTAATGGCAACATTTTCTGCTGATTGATGAGAATTAATTAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCTAGCTTTGTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGTATCCGCCCAGGGCTTGCAgg                                                                                                                                                                                                                                   <  1:462784/251‑1 (MQ=255)                     tAATGCCACCATTTTCTACTGATTGCTGAGAATTAATTAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCACCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACCTGTCATGAAGTGCAATACGTAAACGCGTCACGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGAtcccc                                                                                                                                                                                                                             >  2:624385/1‑247 (MQ=255)                      aaTGGCAACATTTTCTGCTGATTGATGAGAATTAATTAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCgg                                                                                                                                                                                                                            >  2:513190/1‑251 (MQ=255)                       aTGGCAACATTTTCTGCCGGTTGGTGACAATTAATTAAGCTATCGAGCTGTTCACACAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCCTCGCCACTACCGAAGACGCCGTGCGCTCCCACCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACCTGTCATGCAGTGCAATACGTAAACGCGTCGCCCACTTGTTAATGCTGGAGATATTGCCGACTCCGCACAGTGCTTGCAGGATACCGGc                                                                                                                                                                                                                           >  2:258054/1‑251 (MQ=255)                         ggCAACATTTTCTGCTGATTGATGAGAATTAATTAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAg                                                                                                                                                                                                                         <  2:375514/251‑1 (MQ=255)                           cAACATTTTCTGCTGATTGATGAGAATTAATTAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCCTCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGTCTTGCAGGATACCGGCAGcc                                                                                                                                                                                                                       <  2:451057/251‑1 (MQ=255)                              cATNTTCTGCTGATTGATGAGAATTAATTAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGTCTTGCAGGATACCGGCAGCCTGa                                                                                                                                                                                                                    <  2:589848/251‑1 (MQ=255)                                  ttCTGCTGATTGATGAGAATTAATTAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATcc                                                                                                                                                                                                                 >  2:187258/1‑250 (MQ=255)                                              atgaGAATTAATTAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTg                                                                                                                                                                                                    >  2:92141/1‑251 (MQ=255)                                                  gaATTAATTAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGctc                                                                                                                                                                                                >  2:466670/1‑251 (MQ=255)                                                   aattaattAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGctct                                                                                                                                                                                               <  2:102060/251‑1 (MQ=255)                                                       aattaaGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCcg                                                                                                                                                                                           >  2:40923/1‑251 (MQ=255)                                                       aagtaaGCTATCTAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAAGACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGTCTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCATCTCTTTTTTTGGCTCTGCcg                                                                                                                                                                                           <  2:398436/248‑1 (MQ=255)                                                                  cGAGCTGTTCACGCAGCTGAGATACATGCATACCAATGATCACCTGAATGGCCTCGCCACTTCGTAAGACGCCGTGCGCTCCCAGCTTTTTAAATACTTCGTCATCCAGCGTTTGTGACATGTCTTGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGTATACCGGCAGCCTGATCCAGCTCTTTTTTTGGATCTGCCGCTGTGGTTTg                                                                                                                                                                                 <  2:378416/250‑1 (MQ=255)                                                                      cTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTc                                                                                                                                                                            >  2:610740/1‑251 (MQ=255)                                                                          tCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTcgcg                                                                                                                                                                         >  2:211255/1‑250 (MQ=255)                                                                              gCAGCTGAGATACATGCAGACCAATGATCACCTGAATGTCATCGCCACTACGGAAGACGCCGTGCGCGCCCAGCTTGTGAAAGACTTCATCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAAGGCTGGTGATATTGCCGACTCCGccc                                                                                                                                                                                                                                               <  2:378461/176‑1 (MQ=255)                                                                              gCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCt                                                                                                                                                                    >  2:557275/1‑251 (MQ=255)                                                                              gCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCt                                                                                                                                                                    <  1:211255/251‑1 (MQ=255)                                                                              gCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGAGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGccc                                                                                                                                                                                                                                               >  1:378461/1‑176 (MQ=255)                                                                                agttGAGATACATGCAGACCAATCATCTCCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATCCCGGCAGCCTGATCCATCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCttt                                                                                                                                                                  <  1:584615/248‑1 (MQ=255)                                                                                    gagaTACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGctc                                                                                                                                                                                                            >  1:353265/1‑205 (MQ=255)                                                                                    gagaTACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGctc                                                                                                                                                                                                            <  2:353265/205‑1 (MQ=255)                                                                                     agaTACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTAtt                                                                                                                                                             <  1:268527/251‑1 (MQ=255)                                                                                      gaTACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTc                                                                                                                                                            >  2:476569/1‑251 (MQ=255)                                                                                        tACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCgg                                                                                                                                                          >  1:304184/1‑251 (MQ=255)                                                                                           aTGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCtt                                                                                                                                                       <  2:578443/251‑1 (MQ=255)                                                                                                          cACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCa                                                                                                                                        <  1:432955/251‑1 (MQ=255)                                                                                                          cACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCGGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTGGTATTCGGCTTTTGAGTAGAGTTTc                                                                                                                                         <  2:100501/250‑1 (MQ=255)                                                                                                                            aCTACGGCAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCTTCCAGCGTTTGAGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCGCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTc                                                                                                                      <  2:660245/251‑1 (MQ=255)                                                                                                                              tACGGAAGAGGCCGTGAGCTCACCGGTTTTTAAAGACTCCGTCCTCAAGCGTTTGTGACATGTCATGCAGTGCACTACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTCTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTGTCACTTCAGCATCTTCACGTccc                                                                                                                    >  2:79423/1‑251 (MQ=255)                                                                                                                               aCGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACGTCGTCATCCAGCTTTTGTGACATGTCATTCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGTCCTCGCGAGTCTTGGTATTCGTCTTTTGAGTAGAGTTTCTCTTCCGCATCTTCACGTccc                                                                                                                    <  1:120296/250‑1 (MQ=255)                                                                                                                                cGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCgg                                                                                                                  >  1:432182/1‑251 (MQ=255)                                                                                                                                cGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGACATCCAGCGTTTATGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGACTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCgg                                                                                                                  >  2:660876/1‑251 (MQ=255)                                                                                                                                    aGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTg                                                                                                              <  1:573333/251‑1 (MQ=255)                                                                                                                                    aGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTGGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACGTCCGCATCTTCACGTCCCGGCGTg                                                                                                              <  2:364042/251‑1 (MQ=255)                                                                                                                                     gACGCTGTGCCCTCCCAGCTTTTTAATGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTAATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCTGCGTGc                                                                                                             >  2:90657/1‑251 (MQ=255)                                                                                                                                       cGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGcac                                                                                                           >  1:508004/1‑251 (MQ=255)                                                                                                                                         ccGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACAt                                                                                                         <  1:258773/251‑1 (MQ=255)                                                                                                                                           gTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATGCGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGt                                                                                                       <  1:602316/251‑1 (MQ=255)                                                                                                                                                 tCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGcc                                                                                                                                                                                                                       >  2:167228/1‑133 (MQ=255)                                                                                                                                                 tCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGcc                                                                                                                                                                                                                       <  1:167228/133‑1 (MQ=255)                                                                                                                                                  cccAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTGATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTg                                                                                                <  2:618221/251‑1 (MQ=255)                                                                                                                                                  cccAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAAc                                                                                                  <  2:55431/249‑1 (MQ=255)                                                                                                                                                  cccAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTg                                                                                                <  1:203474/251‑1 (MQ=255)                                                                                                                                                    cAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGAAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGc                                                                                               >  1:503247/1‑250 (MQ=255)                                                                                                                                                                     cGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCATGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGaaa                                                                             <  2:405070/251‑1 (MQ=255)                                                                                                                                                                     cGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGaa                                                                              <  1:92141/250‑1 (MQ=255)                                                                                                                                                                           ccAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACg                                                                       <  2:107933/251‑1 (MQ=255)                                                                                                                                                                             aGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGa                                                                      >  1:272449/1‑250 (MQ=255)                                                                                                                                                                               cGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCAGGCGTGCACATGTTGAACTGCAGAATCAGCGGGCGGAAAACCACGAAGt                                                                   >  1:383223/1‑251 (MQ=255)                                                                                                                                                                                gTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGt                                                                   >  1:5399/1‑250 (MQ=255)                                                                                                                                                                                    gtgACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTAcagc                                                              <  1:476569/251‑1 (MQ=255)                                                                                                                                                                                      gACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTAcagc                                                              >  1:37102/1‑249 (MQ=255)                                                                                                                                                                                           gTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCGTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTaa                                                       <  2:213259/251‑1 (MQ=255)                                                                                                                                                                                              aTGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGAGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAACACCACGAAGTACAGCAGGGTAAAGc                                                    >  2:326009/1‑251 (MQ=255)                                                                                                                                                                                                  aGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCCGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGcaca                                                 >  2:534000/1‑250 (MQ=255)                                                                                                                                                                                                          aCGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATgg                                        <  1:466670/251‑1 (MQ=255)                                                                                                                                                                                                           cGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGc                                       <  1:187258/251‑1 (MQ=255)                                                                                                                                                                                                                         gccaTTGTTAATGCTGGATATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCGCATGTTGAACTGCAGAATCAGTGTTCGGAAAACCACGGAGTACAGCAGGGTAAAGCACAGCCCAATGGCGATTTGGGTCTGca                         <  2:490610/248‑1 (MQ=255)                                                                                                                                                                                                                                     gCTGGAGATATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCcgcgg             >  1:547234/1‑250 (MQ=255)                                                                                                                                                                                                                                          agaTATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTg        >  1:364516/1‑251 (MQ=255)                                                                                                                                                                                                                                           gaTATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGc       <  1:534000/251‑1 (MQ=255)                                                                                                                                                                                                                                           gaTATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGc       >  2:229971/1‑251 (MQ=255)                                                                                                                                                                                                                                           gaTATTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACACGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGCACAGCAGGGTAAAGCACAACCCAATGGCCATTTGGGTCAGCATCCTAGCCGCGAGGttcc       >  2:337918/1‑249 (MQ=255)                                                                                                                                                                                                                                              atTGCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTg     <  2:531022/250‑1 (MQ=255)                                                                                                                                                                                                                                                 gCCGACTCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAAc  >  1:270506/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    TGGTCATTATACTGCCTCCGTAATGGCAACATTTTCTGCTGATTGATGAGAATTAATTAAGCTATCGAGCTGTTCACGCAGCTGAGATACATGCAGACCAATGATCACCTGAATGGCATCGCCACTACGGAAGACGCCGTGCGCTCCCAGCTTTTTAAAGACTTCGTCATCCAGCGTTTGTGACATGTCATGCAGTGCAATACGTAAACGCGTCGCGCAATTGTTAATGCTGGAGATATTGCCGACCCCGCCCAGGGCTTGCAGGATACCGGCAGCCTGATCCAGCTCTTTTTTTGGCTCTGCCGCGGTGGTTTGGCCTCGCGAGGCTTTGTATTCGGCTTTTGAGTAGAGTTTCACTTCCGCATCTTCACGTCCCGGCGTGCACATGTTGAACTGCAGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAAC  >  NC_000913/3861981‑3862470

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