Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,132,292	T→C	intergenic (+25/‑324)	metL → / → metF	fused aspartate kinase/homoserine dehydrogenase 2/5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,132,292	0	T	C	100.0%	104.4 / NA	33	intergenic (+25/‑324)	metL/metF	fused aspartate kinase/homoserine dehydrogenase 2/5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/21);  total (12/21)

GCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTC  >  NC_000913/4132054‑4132534                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   gcgTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTc                                                                                                                                                                                                                                        <  2:265316/251‑1 (MQ=255)   gTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCgg                                                                                                                                                                                                                                      >  1:25726/1‑251 (MQ=255)         aTGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTAt                                                                                                                                                                                                                                <  1:86065/251‑1 (MQ=255)           gCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTTCCGTGCGATAACGTCTTGGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATgg                                                                                                                                                                                                                              <  1:241377/251‑1 (MQ=255)           gCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATTTTTTCGGGTTTATg                                                                                                                                                                                                                               <  1:113280/250‑1 (MQ=255)            ccAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGt                                                                                                                                                                                                                             <  1:137642/251‑1 (MQ=255)              aaCGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGttt                                                                                                                                                                                                                           >  1:101070/1‑251 (MQ=255)                    aaaGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCAGCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGTGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCATAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATg                                                                                                                                                                                                                     <  2:44857/251‑1 (MQ=255)                          ggtgtAGGCGTGGAAGGGGGGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAatat                                                                                                                                                                                                               <  2:330691/250‑1 (MQ=255)                                           ggTGTGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGTTCGCGACGTCACCGCCTGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGa                                                                                                                                                                                              <  2:329078/251‑1 (MQ=255)                                               cGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGTTAACGTCATGGACATAGAAAGCCGCTGGTATCGCGATCAACCTCTGGCGATACGCGGACCTGGCGCTGGGCGTGACGTCACAGTCGGGGCGATTCCGTCGGATATCAACCGGATGGCGCCGTTGTTGTAATTTCAGAAATTTAATAATGCCAGGCAATCATGTTTTCGGGTTTATGGTTTCTAGTGAACTATATTGCATTATCATAGGAtt                                                                                                                                                                                            >  1:49542/1‑249 (MQ=255)                                                   aaGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGGCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCGGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTTAATTATCATAGGATTAGGCCg                                                                                                                                                                                      <  2:217364/251‑1 (MQ=255)                                                    aGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCgg                                                                                                                                                                                     >  1:319506/1‑251 (MQ=255)                                                                            gcCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCg                                                                                                                                                              <  1:317725/250‑1 (MQ=255)                                                                                      aaCGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATACCCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACt                                                                                                                                                                                                                                                 <  2:251678/157‑1 (MQ=255)                                                                                      aaCGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACt                                                                                                                                                                                                                                                 >  1:251678/1‑157 (MQ=255)                                                                                             gtGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCAGGATTAAGCGTTTCCGACGAATCCGGCAAGAAGCAATAAGt                                                                                                                                              <  2:279921/248‑1 (MQ=255)                                                                                                                                            ggACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttt                                                                                             >  1:360274/1‑251 (MQ=255)                                                                                                                                                  ggCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttcttt                                                                                        <  1:3828/250‑1 (MQ=255)                                                                                                                                                        gggCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAc                                                                                 <  1:152232/251‑1 (MQ=255)                                                                                                                                                          gcgcGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCa                                                                               <  1:235824/251‑1 (MQ=255)                                                                                                                                                              gACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCt                                                                           <  1:164892/251‑1 (MQ=255)                                                                                                                                                                 gTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGaa                                                                        <  1:108921/251‑1 (MQ=255)                                                                                                                                                                              cGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCtt                                                           <  2:328742/251‑1 (MQ=255)                                                                                                                                                                                               aaCCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCAc                                          >  1:306003/1‑251 (MQ=255)                                                                                                                                                                                                aCCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCAcc                                         <  1:293112/251‑1 (MQ=255)                                                                                                                                                                                                   ggCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTGATGGTTTCTAATGAAATATATTGAATGATCATAGGATTAGGCCTGATTAAGCGTTTACGACGAATCCTGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGTTTTGCATACCTGATTTTCTTTTTCTTACCAGCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACcatc                                      <  1:222590/251‑1 (MQ=255)                                                                                                                                                                                                       ggCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTg                                  >  1:12345/1‑251 (MQ=255)                                                                                                                                                                                                       ggCACAGTTGTTCTAATTTCAGACATTTAATAATGCGCGGCACTCCTGTTTACGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATATTCCCTGATTTTTCCTCACCATCATTg                                  >  2:363391/1‑251 (MQ=255)                                                                                                                                                                                                                        ttCAGAAATTTAATAATGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGAc                 >  1:363088/1‑251 (MQ=255)                                                                                                                                                                                                                                       atGCCCGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTc  >  1:85704/1‑251 (MQ=255)                                                                                                                                                                                                                                           ccGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCgg                                                                                                                                                                                     >  2:161165/1‑68 (MQ=255)                                                                                                                                                                                                                                           ccGGCACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCgg                                                                                                                                                                                     <  1:161165/68‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   GCGTTTCGATGCCAACGGTAAAGCGCGTGTAGGCGTGGAAGCGGTGCGTGAAGATCATCCGTTGGCATCACTGCTGCCGTGCGATAACGTCTTTGCCATCGAAAGCCGCTGGTATCGCGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTTTTCATCTTCCCTGATTTTTCCTCACCATCATTGGTCATTTTTCGGTTGACGCCCTTCGGCTTTTC  >  NC_000913/4132054‑4132534

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