Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	535,337	T→C	A474A (GCT→GCC)	gcl →	glyoxylate carboligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	535,337	0	T	C	96.5%	96.1 / ‑6.1	29	A474A (GCT→GCC)	gcl	glyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (14/14);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTCTGTCACAAATCGCTGCGGCACAAATGCTGCATGTCTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCA  >  NC_000913/535094‑535546                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                  tGGTCTGTCACAAATCGCTGCGGCACAAATGCTGCATGTCTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTCGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTgg                                                                                                                                                                                                            >  1:151477/1‑251 (MQ=255)  ggTCTGTCACAAATCGCTGCGGCACAAATGCTGCATGTCTTTAAAGGCCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTggg                                                                                                                                                                                                           <  1:1721/251‑1 (MQ=255)         cacaAATCGCTGCGGCACAAATGCTGCATGTCTTTAAAGAACGCCACTGGATCAACTGTGGACAGGCTGGTCCGTCAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGCGCTGAtt                                                                                                                                                                                                    >  1:178614/1‑251 (MQ=255)                  cTGCGGCACAAATGCTGCATGTCTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCgtcagtc                                                                                                                                                                                            >  1:241538/1‑250 (MQ=255)                               gctgcATGTCTTTAAAGACAGCCACTGGATCACCTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCCGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCACTGTGGTGGCGATTTCTGGCGACTTTGACTTCCCGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACCTTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTg                                                                                                                                                                              >  2:297137/1‑251 (MQ=255)                                ctgcATGTCTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCCCAACGCGCTTTTg                                                                                                                                                                              >  2:89340/1‑250 (MQ=255)                                            aaaGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTg                                                                                                                                                                 >  1:145557/1‑251 (MQ=255)                                                cccGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACGTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACAGGGACTACTGCGTg                                                                                                                                                             <  2:177892/250‑1 (MQ=255)                                                   gCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCa                                                                                                                                                           >  1:239319/1‑250 (MQ=255)                                                     cACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACt                                                                                                                                                        >  1:244979/1‑251 (MQ=255)                                                     cACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACt                                                                                                                                                        <  1:163320/251‑1 (MQ=255)                                                          gATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCtt                                                                                                                                                   <  1:297137/251‑1 (MQ=255)                                                                                    tAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACAGTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTGCGAGAATATCAACTCCagtgaagt                                                                                                                          <  1:113900/250‑1 (MQ=255)                                                                                          ggACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGc                                                                                                                   >  1:227359/1‑251 (MQ=255)                                                                                                  ccGGCTGCGCTAGGGGTTTGTGCCGCTGGTCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTGTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATAGCTACGgtg                                                                                                            <  2:133520/250‑1 (MQ=255)                                                                                                           gtAGGGGTTTGTGCCGCTGGTCCGAAACGCAATGTGGGTGCGATTTCTGGCGACTGTGACTTCCAGGTCCTGATTGAAGAGTTAGCTGTGGGCGCGCAGGTCATCATTCCGTACATCCATGTGCTGGTCAACAACGCCTATTTGGGGCTGATTCGTCAGTCACAACGCGGTTTTGTCATTGACTACGTCTTGCAACTCGCTTTCGAGAATTTCAACTCCAGTGAAGGGAATTGCTACGGTGCTGACCAc                                                                                                    <  2:271241/248‑1 (MQ=255)                                                                                                             aGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTa                                                                                                 >  1:214085/1‑250 (MQ=255)                                                                                                             aGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTa                                                                                                 <  2:63548/250‑1 (MQ=255)                                                                                                                             gATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGttt                                                                                <  2:227359/251‑1 (MQ=255)                                                                                                                             gATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGGAGCGGAAGGttt                                                                                <  2:244979/251‑1 (MQ=255)                                                                                                                                  gAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGCTCCTGATTGAAGAGTTAGCTGTTGGCGAGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGtt                                                                           >  1:21326/1‑251 (MQ=255)                                                                                                                                   aaaCGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTg                                                                          >  2:2637/1‑251 (MQ=255)                                                                                                                                        cAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAg                                                                     <  1:86104/251‑1 (MQ=255)                                                                                                                                              ggtggCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGATCAACATTCCGTACCTCCCTGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGATATTc                                                               >  1:196967/1‑251 (MQ=255)                                                                                                                                                          tGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCCGTTCAACATTCCGTACATCCATGTGCTCGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAAc                                                   >  1:198791/1‑251 (MQ=255)                                                                                                                                                              gACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGa                                               <  2:170297/251‑1 (MQ=255)                                                                                                                                                              gACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACGCGCTTTCGAGAATGTCAACTCCAGTGAAGTGAATGGCTATGGTGGTGACCAGGTAAAAGTTGCTGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGa                                               <  1:182331/251‑1 (MQ=255)                                                                                                                                                                                        agagTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAgg                     >  1:174714/1‑251 (MQ=255)                                                                                                                                                                                                           cAGTTCAACATTCCGTACAGCCATGTTCTGGTCAACAACGCCTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCa  <  1:256834/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                  TGGTCTGTCACAAATCGCTGCGGCACAAATGCTGCATGTCTTTAAAGACCGCCACTGGATCAACTGTGGTCAGGCTGGTCCGTTAGGCTGGACGATTCCGGCTGCGCTAGGGGTTTGTGCCGCTGATCCGAAACGCAATGTGGTGGCGATTTCTGGCGACTTTGACTTCCAGTTCCTGATTGAAGAGTTAGCTGTTGGCGCGCAGTTCAACATTCCGTACATCCATGTGCTGGTCAACAACGCTTATCTGGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCA  >  NC_000913/535094‑535546

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