Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,813	Δ2 bp	coding (142‑143/1161 nt)	metB →	O‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,128,811	0	C	.	96.5%	124.3 / ‑3.3	29	coding (140/1161 nt)	metB	O‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (20/8);  total (20/9)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
*	NC_000913	4,128,812	0	G	.	96.5%	124.5 / ‑3.8	29	coding (141/1161 nt)	metB	O‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (20/8);  total (20/9)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAA  >  NC_000913/4128563‑4129043                                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                        cAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACtcgt                                                                                                                                                                                                                                        <  2:3889/251‑2 (MQ=255)          cGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAAc                                                                                                                                                                                                                             >  1:47785/1‑251 (MQ=255)                   gTTTAGATGTCCAGATGTATTGACGTCCATTAAAACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAAcgcgc                                                                                                                                                                                                                    >  1:107003/1‑251 (MQ=255)                      tAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGAt                                                                                                                                                                                                                 >  2:175947/1‑251 (MQ=255)                               ataTGTATTGACGTCCAGTAACACAATGTTTACTCTGGTGCCTGACAGTTCACCGACAAAGCCCAGTGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATTGTTGCGTTGTCCCACCGATCCATCTTGCCAGCACCTATAACTTTACCGGAGTTAATGAACCGCGCGCGCATGATTA‑‑CGCGTCGCGTCAACCCATCGCGTGAT‑TGGtttcag                                                                                                                                                                                                       <  2:146081/249‑5 (MQ=255)                                               aTTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAAc                                                                                                                                                                                        >  1:195760/1‑251 (MQ=255)                                                            ttACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGgtggtg                                                                                                                                                                            <  2:97093/250‑1 (MQ=255)                                                                tctGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGGATGGCTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTACTGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctggt                                                                                                                                                                       >  2:73763/1‑251 (MQ=255)                                                                 ctGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGtggtgctggtg                                                                                                                                                                      >  1:92862/1‑251 (MQ=255)                                                                          gACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGtactt                                                                                                                                                              >  1:57604/1‑250 (MQ=255)                                                                                 cACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTCCTAATa                                                                                                                                                       <  1:74531/250‑1 (MQ=255)                                                                                       cAAAGCCCCGGGCACCTCATCACATGACGCGTCAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCAt                                                                                                                                                >  2:13685/1‑251 (MQ=255)                                                                                         aaGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGt                                                                                                                                              >  1:99165/1‑251 (MQ=255)                                                                                            cccAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCcg                                                                                                                                           >  2:103205/1‑251 (MQ=255)                                                                                                            acaTGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATTTCCGCGATTCACCTGGTAAc                                                                                                                           <  2:233371/251‑1 (MQ=255)                                                                                                                    cgTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGCTCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCt                                                                                                                   >  2:53582/1‑251 (MQ=255)                                                                                                                        aaCAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCttttt                                                                                                               <  2:147692/251‑1 (MQ=255)                                                                                                                           aGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGCAACGACCGGCTTTTTGa                                                                                                             >  2:137843/1‑250 (MQ=255)                                                                                                                              ccaccaTCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAAcc                                                                                                         >  1:140498/1‑251 (MQ=255)                                                                                                                                                        aaTGACGACGCACAGGATGGTCGCGTTGTCCCACCGATCCATCTTTCAAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTCAT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTAATTACTAATACCGGCATGTCCGCGATTCAACTGGTACCGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGcgccgaccg                                                                               >  2:226650/1‑247 (MQ=255)                                                                                                                                                            acgacgAACAGTATGGTTGCGTTGTTCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCTGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTg                                                                           >  2:172570/1‑251 (MQ=255)                                                                                                                                                            acgacgAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTg                                                                           >  1:99491/1‑251 (MQ=255)                                                                                                                                                                   aCAGTATGGTTGCGTTGTCCCAGCGATCCAACTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGCAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTTcggc                                                                    <  2:8871/251‑1 (MQ=255)                                                                                                                                                                        aTGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAg                                                               <  2:42802/251‑1 (MQ=255)                                                                                                                                                                              gCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCg                                                         <  1:56689/251‑1 (MQ=255)                                                                                                                                                                                gTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGcc                                                       <  1:73674/251‑1 (MQ=255)                                                                                                                                                                                   gTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGGAGCCATCGCCTGt                                                    >  1:214802/1‑251 (MQ=255)                                                                                                                                                                                     cccACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTAcggcgggagactcaacctgtc                                                   >  2:136445/1‑235 (MQ=255)                                                                                                                                                                                             tCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTcc                                          >  1:88249/1‑250 (MQ=255)                                                                                                                                                                                                                                       cgcgcgcgCATGATTACT‑‑CGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGGGTTGTTTGTTGAttca  >  1:88812/1‑246 (MQ=255)                                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                                        CAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATAACTTTACCGGATTTAATGAACCGCGCGCGCATGATTACTCGCGTCGCGGCAACCCAACGCGCGATGTGGTTCAGCGTGCGCTGGCAGAACTGGAAGGTGGTGCTGGTGCAGTACTTACTAATACCGGCATGTCCGCGATTCACCTGGTAACGACCGTCTTTTTGAAACCTGGCGATCTGCTGGTTGCGCCGCACGACTGCTACGGCGGTAGCTATCGCCTGTTCGACAGTCTGGCGAAACGCGGTTGCTATCGCGTGTTGTTTGTTGATCAA  >  NC_000913/4128563‑4129043

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