Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,598,984	G→A	P153L (CCG→CTG)	opgB ←	phosphoglycerol transferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,598,984	0	G	A	100.0%	182.3 / NA	49	P153L (CCG→CTG)	opgB	phosphoglycerol transferase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/21);  total (28/21)

AGGCGGAGGCGTTGCCTTCAAAGGGGGCAAACAGCGGTATGCCGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACGGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGAC  >  NC_000913/4598739‑4599225                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  aggcggaggcgTTGCCTTCAAAGGGGGCAAACAGCGGTATGCCGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGctct                                                                                                                                                                                                                                              <  2:481696/251‑1 (MQ=255)  ggcggaggcgTTGCCTTCAAAGGGGGCAAACAGCGGTATGCCGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCtt                                                                                                                                                                                                                                             <  2:406958/251‑1 (MQ=255)     ggaggcgTTGCCTTCAAAGGGGGCAAACAGCGGTATGCCGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTtat                                                                                                                                                                                                                                          >  1:232860/1‑251 (MQ=255)      gaggcgTTGCCTTCACAGGGGGCAAACAGCGGTATGCCGCACTGAGACGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTtat                                                                                                                                                                                                                                          >  2:283457/1‑250 (MQ=255)        ggcgTTGCCTTCAAAGGGGGCAAACAGCGGTATGCCGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTataat                                                                                                                                                                                                                                       <  2:254752/251‑1 (MQ=255)           gTTGCCTTCAAAGGGGGCAAACAGCGGTATGCCGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTataata                                                                                                                                                                                                                                      >  1:504913/1‑249 (MQ=255)              gCCTTCAAAGGGGGCAAACAGCGGTATGCCGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCg                                                                                                                                                                                                                                 >  1:116366/1‑251 (MQ=255)                     aaGGGGGCAAACAGCGGTATGCCGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTc                                                                                                                                                                                                                          >  2:49420/1‑251 (MQ=255)                      aGGGGGCAAACAGCGGTATGCCGCACTGAGCAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTcc                                                                                                                                                                                                                         >  1:529137/1‑251 (MQ=255)                        ggggCAAACAGCGGTATGCCGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCgg                                                                                                                                                                                                                       >  1:101758/1‑251 (MQ=255)                         gggCAAACAGCGGTATGCCGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCggg                                                                                                                                                                                                                      >  1:136640/1‑251 (MQ=255)                                cAGCGGTATGCCGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGt                                                                                                                                                                                                               <  2:116702/251‑1 (MQ=255)                                    ggTATGCCGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTcgcg                                                                                                                                                                                                           >  1:243935/1‑251 (MQ=255)                                    ggTATGCCGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTcgcg                                                                                                                                                                                                           >  1:306193/1‑251 (MQ=255)                                       aTGCCGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTg                                                                                                                                                                                                         <  1:20308/250‑1 (MQ=255)                                           cGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTgg                                                                                                                                                                                                    >  1:465036/1‑251 (MQ=255)                                                 gagaAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTc                                                                                                                                                                                              <  2:275664/251‑1 (MQ=255)                                                     aGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCAcc                                                                                                                                                                                           >  2:105350/1‑250 (MQ=255)                                                        caccaTGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGtt                                                                                                                                                                                       <  2:571497/251‑1 (MQ=255)                                                                       cGTGTAATCTGTTCCCGGCAGCTGCTGCGTGGGGCTGAAATCCAGGCCTTCCTTTTTCAACGCGCCCAGTTCAGGCGCGAGATACGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCCCAAAGTCCGGGTCGGCGCGGCGCGACTGGGATTCCACCAGTTCCGTTACCTGACGa                                                                                                                                                                         >  2:451782/1‑250 (MQ=255)                                                                          gtAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCAGTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAAt                                                                                                                                                                      <  2:29333/250‑1 (MQ=255)                                                                                 ggtCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCggg                                                                                                                                                                <  2:540382/247‑1 (MQ=255)                                                                                 gTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTGATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCggg                                                                                                                                                                <  2:276872/249‑1 (MQ=255)                                                                                 gTTCCCGGCAGATGCTGAGAGTGACTGAAATCCAGGCCTTCATTTTTCCACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAACTAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGCCTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGcc                                                                                                                                                              >  2:579547/1‑250 (MQ=255)                                                                                   tCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTg                                                                                                                                                             <  2:386402/250‑1 (MQ=255)                                                                                       ggCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTc                                                                                                                                                        >  2:41720/1‑251 (MQ=255)                                                                                        gCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTcc                                                                                                                                                       <  1:306488/251‑1 (MQ=255)                                                                                           gctgctGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCAcc                                                                                                                                                    >  1:123326/1‑251 (MQ=255)                                                                                                         cTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGcc                                                                                                                                        >  2:44742/1‑249 (MQ=255)                                                                                                           gAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGt                                                                                                                                     >  1:260052/1‑250 (MQ=255)                                                                                                                ccAGGCCTTCATTTTGCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAgc                                                                                                                               <  2:494688/251‑1 (MQ=255)                                                                                                                        tCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCcagc                                                                                                                         >  1:127635/1‑249 (MQ=255)                                                                                                                              ttCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTCTCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGt                                                                                                                  >  2:138602/1‑250 (MQ=255)                                                                                                                               tCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTa                                                                                                                 <  1:283457/250‑1 (MQ=255)                                                                                                                                      gcCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAAt                                                                                                         >  2:159639/1‑251 (MQ=255)                                                                                                                                       cccAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGCAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATg                                                                                                        >  2:564918/1‑251 (MQ=255)                                                                                                                                         cAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATGAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGt                                                                                                      >  1:454535/1‑251 (MQ=255)                                                                                                                                               aGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGgat                                                                                                >  2:444154/1‑251 (MQ=255)                                                                                                                                                          ccGGGAAAGCCTTGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTGCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGcgacg                                                                                      <  2:529137/250‑1 (MQ=255)                                                                                                                                                                             aaaaTAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATccagccc                                                                   >  2:412404/1‑250 (MQ=255)                                                                                                                                                                                aTAGTTCCGCTCGAGATTTTCGCCGTAGATATAAACCAGGTTGAGTTGCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATGTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAgcg                                                               <  2:490991/251‑1 (MQ=255)                                                                                                                                                                                              gaCTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCt                                                 >  2:331340/1‑251 (MQ=255)                                                                                                                                                                                                              ataAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACaa                                 <  1:49420/251‑1 (MQ=255)                                                                                                                                                                                                              ataAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACaa                                 <  1:207410/251‑1 (MQ=255)                                                                                                                                                                                                                aaaCCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAAt                                <  2:84733/250‑1 (MQ=255)                                                                                                                                                                                                                           aGTTTCGGGTCAGGGATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCgg                    >  1:308539/1‑251 (MQ=255)                                                                                                                                                                                                                                       gggATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTAttt        <  2:116366/251‑1 (MQ=255)                                                                                                                                                                                                                                       gggATAGTTTTCGACAGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTAttt        >  1:566778/1‑251 (MQ=255)                                                                                                                                                                                                                                             gTTTTCGACAGCGCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGAc  <  1:321687/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  AGGCGGAGGCGTTGCCTTCAAAGGGGGCAAACAGCGGTATGCCGCACTGAGAAGCCACCATGCCCGCAATCGTGTAATCTGTTCCCGGCAGCTGCTGCGTGTGGCTGAAATCCAGGCCTTCATTTTTCAACGCGCCCAGTTCAGGCGTGAGATCCGGGAAAGCCTCGTTATCAAAATAGGTCCGCTCGAGACTTTCGCCGTAGATATAAACCAGGTTGAGTTTCGGGTCAGGGATAGTTTTCGACGGCTCTTTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGAC  >  NC_000913/4598739‑4599225

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