Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	535,299	T→C	F462L (TTC→CTC)	gcl →	glyoxylate carboligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	535,299	0	T	C	97.4%	135.6 / ‑3.8	39	F462L (TTC→CTC)	gcl	glyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (19/19);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

TGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACG  >  NC_000913/535158‑535437                                                                                                                                              |                                                                                                                                           tgtgGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTc                                                                                                                                   >  2:362465/1‑151 (MQ=255)  gtgGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACAtt                                                                                                                                    <  1:551643/149‑1 (MQ=255)  gtgGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACAtt                                                                                                                                    <  1:89813/149‑1 (MQ=255)  gtgGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTcc                                                                                                                                  <  1:276416/151‑1 (MQ=255)     gTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTa                                                                                                                               <  1:33448/151‑1 (MQ=255)        aGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACAt                                                                                                                            >  2:332520/1‑151 (MQ=255)         ggCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATc                                                                                                                           >  2:427129/1‑151 (MQ=255)               tCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCAtgtg                                                                                                                     >  2:77645/1‑151 (MQ=255)                ccGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCAtgtg                                                                                                                     <  2:97654/150‑1 (MQ=255)                          ggACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTcaac                                                                                                            <  1:184118/149‑1 (MQ=255)                              gATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCt                                                                                                      <  2:479581/151‑1 (MQ=255)                               aTTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCtt                                                                                                     <  2:59982/151‑1 (MQ=255)                                             aGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGAtt                                                                                       <  1:216688/151‑1 (MQ=255)                                                ggTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgt                                                                                    >  2:135032/1‑151 (MQ=255)                                                 gTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgtc                                                                                   >  2:171984/1‑151 (MQ=255)                                                    tgtgCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCgtcagt                                                                                >  2:358775/1‑151 (MQ=255)                                                          gCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAAc                                                                          >  2:445027/1‑151 (MQ=255)                                                             gATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAAcgcg                                                                       <  2:210351/151‑1 (MQ=255)                                                                 cGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCttt                                                                   <  1:135035/151‑1 (MQ=255)                                                                      cGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGAc                                                               >  1:228999/1‑150 (MQ=255)                                                                      cGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACa                                                              <  2:350548/151‑1 (MQ=255)                                                                        cAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATg                                                            >  2:365137/1‑151 (MQ=255)                                                                              ggtggCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGacta                                                       <  2:352506/150‑1 (MQ=255)                                                                                    gATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTg                                                <  1:202345/151‑1 (MQ=255)                                                                                                cTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTc                                    >  2:363811/1‑151 (MQ=255)                                                                                                 tttGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGacta                                                       <  1:176706/131‑1 (MQ=255)                                                                                                 tttGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGacta                                                       >  2:176701/1‑131 (MQ=255)                                                                                                 tttGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCg                                   >  1:606493/1‑151 (MQ=255)                                                                                                     aCTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCgaga                                >  2:535430/1‑150 (MQ=255)                                                                                                     aCTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCgaga                                >  1:84815/1‑150 (MQ=255)                                                                                                            gTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAAc                        >  1:501860/1‑151 (MQ=255)                                                                                                              tCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACt                       <  1:405055/150‑1 (MQ=255)                                                                                                                  gATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCagt                  <  1:575567/151‑1 (MQ=255)                                                                                                                  gATTGAAGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCag                   <  2:431991/150‑1 (MQ=255)                                                                                                                       aaGAGTTAGCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCagtgaagt             >  2:473340/1‑151 (MQ=255)                                                                                                                               gCTGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCt     <  2:419590/151‑1 (MQ=255)                                                                                                                                 tGTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTAc   >  1:61014/1‑151 (MQ=255)                                                                                                                                  gTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTAc   <  1:332527/150‑1 (MQ=255)                                                                                                                                  gTTGGCGCGCAGCTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACg  >  1:72853/1‑151 (MQ=255)                                                                                                                                              |                                                                                                                                           TGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACG  >  NC_000913/535158‑535437

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