Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,070,858	T→G	K129N (AAA→AAC)	fbaA ←	fructose‑bisphosphate aldolase class II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,070,858	0	T	G	96.8%	112.0 / NA	31	K129N (AAA→AAC)	fbaA	fructose‑bisphosphate aldolase class II
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (17/13);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCG  >  NC_000913/3070622‑3071098                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 gtgtATGCGTAATCAACGTCTTCCGGCTGGGAGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAAAAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTc                                                                                                                                                                                                                                    >  1:241539/1‑251 (MQ=255)    tATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGAAGAAGCGTCCATGTGACTGTTGTCAACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTACAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGAGGTTTTCAACCGCGTCCa                                                                                                                                                                                                                                  >  1:160601/1‑250 (MQ=255)                      tCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCGTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCAc                                                                                                                                                                                                                 <  2:43297/249‑1 (MQ=255)                                                                        gTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCa                                                                                                                                                                               <  2:240177/233‑1 (MQ=255)                                                                        gTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCa                                                                                                                                                                               >  1:240177/1‑233 (MQ=255)                                                                          cttcttCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATa                                                                                                                                                            >  2:2618/1‑250 (MQ=255)                                                                               tcACCACCGGTGCAACACAGTTCGATTTACAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTc                                                                                                                                                      >  1:73003/1‑251 (MQ=255)                                                                                 accaccGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTc                                                                                                                                                      >  1:35313/1‑249 (MQ=255)                                                                                          gCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGt                                                                                                                                           <  2:239081/251‑1 (MQ=255)                                                                                               ccAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACg                                                                                                                                      <  2:31303/251‑1 (MQ=255)                                                                                                      gATTTCCATAGTCATTCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATgcg                                                                                                                               >  1:155694/1‑251 (MQ=255)                                                                                                         ttCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCAc                                                                                                                            >  1:39728/1‑251 (MQ=255)                                                                                                          tCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCTTAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCAc                                                                                                                            <  1:38452/250‑1 (MQ=255)                                                                                                              gagTCATGCCGATTTTGTAAATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCGTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCCCCagag                                                                                                                       <  2:158352/251‑1 (MQ=255)                                                                                                                  cATGCCGATGTTGGACAGGCGCTCCAGGTATTTAGAGCATATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACATCGGTTTACCTGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCATGAGAACCGGAACACCATAATGTTCAGCCATCTGGTGAACTTGATGCGCTCCAGAGATc                                                                                                                    <  2:22189/250‑1 (MQ=255)                                                                                                                   aTGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATcgc                                                                                                                  <  2:324/251‑1 (MQ=255)                                                                                                                   aTGCCGATGTTGGACAGGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATcgc                                                                                                                  <  2:22563/251‑1 (MQ=255)                                                                                                                     gCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCgcgc                                                                                                                >  2:17945/1‑251 (MQ=255)                                                                                                                      ccGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCTCACCAGAGATCGCGcc                                                                                                               <  1:142340/251‑1 (MQ=255)                                                                                                                                           aGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCgg                                                                                          >  1:140803/1‑251 (MQ=255)                                                                                                                                                                                      cAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGcaca                                                >  2:167965/1‑249 (MQ=255)                                                                                                                                                                                         aCAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGAATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCg                                            <  2:25639/251‑1 (MQ=255)                                                                                                                                                                                                gATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCAAGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTgg                                        >  1:213384/1‑248 (MQ=255)                                                                                                                                                                                                   cATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACt                                  <  2:199747/251‑1 (MQ=255)                                                                                                                                                                                                     tgtgAGAAGAAAACAGAGGACTACCGGTCGCTGCGAAGAAGATTTCACCCGCGTCCAACAGAACGTCGATCCACGGCAGCAGTTTCTTAGCGCAGTCGTCAGCGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTAAGATTTACCGCCTTTACCAGCGCTAAAGGAAGCACCACCGTTGGAGAACTg                                 >  2:147083/1‑250 (MQ=255)                                                                                                                                                                                                         agaagagaaCAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACg                             >  2:180924/1‑250 (MQ=255)                                                                                                                                                                                                           aagagaaCAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGAt                           >  1:159893/1‑250 (MQ=255)                                                                                                                                                                                                                aaCAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGCCCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCgg                     >  1:98729/1‑251 (MQ=255)                                                                                                                                                                                                                  cAGCGGCTTACCGGTAGCTGCGAAGTGGTTTTCCCCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTGCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATGTCACGCCTTTTCCAGCGATAAAGGAAGCACCTCCGTTGGAGAACTGAACGATAACCGgcg                   <  2:235498/251‑1 (MQ=255)                                                                                                                                                                                                                              ggTAGCTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTag       <  2:200143/251‑1 (MQ=255)                                                                                                                                                                                                                                   cTGCGAAGTGGTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCg  >  1:227966/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 GTGTATGCGTAATCAACGTCTTCCGGCTGGGTGTACAGTGCAGAAGCGTCCATGTGGCTGTTGTCCACGCCGTCTTCTTCACCACCGGTGCAACCCAGTTCGATTTCCAGAGTCATGCCGATTTTGGACATGCGCTCCAGGTATTTAGAGCAGATTTCGATGTTCTCTTGCAGAGATTCTTCAGACAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGGAACACCATAATGTTCAGCCATCTGGTGAACGTGATGCGCACCAGAGATCGCGCCCAGGATAGCAGCACCCTGCGGAACGTCAGATTTCACGCCTTTACCAGCGATAAAGGAAGCACCACCGTTGGAGAACTGAACGATAACCGGCGCTTTAACTTTAGCAGCG  >  NC_000913/3070622‑3071098

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