Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	535,582	A→T	E556V (GAG→GTG)	gcl →	glyoxylate carboligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	535,582	0	A	T	95.2%	47.1 / NA	22	E556V (GAG→GTG)	gcl	glyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base T (4/16);  minor base G (0/1);  total (4/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGAGCGTGTGACCAATATTTCGATGGGCAGCGAACTGGATAACGTCATGGAATTTGAAGATATCGCCGATAACGCAGCGGACGCACCGACTGAAACCTGCTTCATGCACTATGAATAAGGGAGATAAATAATGTTACGTTTCTCTGCTAATTTATCGATGTTATTTGGAGAATATGATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGA  >  NC_000913/535343‑535825                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    ggggCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCTGAAGGGTTAGGTTGTAAATCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGTGCGTGTGACCa                                                                                                                                                                                                                                          <  2:267642/251‑1 (MQ=255)  ggtcTGATTCGTCATTCAACACGCGCTTTTGACATGGACTTCTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGTCTACGGTGTTTACCACTTAAAAGTAGCGGAAGGTTTAGTTTGTAAAGCTATTCGGGTCTTCAAACCTGAAGATATTGCGCCAGTCTTTGACCAGGCGAAAGCCTTAATGGCTCAATATCGGGTACCGGTATTCGTGGAAGTTATTCTCGTGCGTGTGACCaa                                                                                                                                                                                                                                         <  1:289949/248‑1 (MQ=255)              cagtcaCAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGTGCGTGTGACCAATATTTCGATgg                                                                                                                                                                                                                              <  1:119558/250‑1 (MQ=255)                       cgcgCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGTGCGTGTCACCAATATTTCGATGGGCAGCGAACt                                                                                                                                                                                                                    <  1:191578/251‑1 (MQ=255)                         cgcTTTTGACATGGACTACTGCGTGCAACTCGCGTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTAGTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGTGCGTGTGACCAATATTTCGATGGGCAGAGAACTgg                                                                                                                                                                                                                  <  1:249565/251‑1 (MQ=255)                                     ggACTACTGCGTGCAACGCGCTTTCGAGAATATCAACTCCAGTGAAGTTGATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATGCTCGTGCGTGTGACCAATATTACGATGGGCAGCGAACTGGATAACGTCATgg                                                                                                                                                                                                      <  1:173459/251‑1 (MQ=255)                                     ggACGACTGCTTGCAACTCGCTTTCGAGAATATCAAATCCAGTTAAGTGAATGGCGACGTTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGTGCGTGTGACCAATATTGCGATGGGCAGCGAACTGGATAACGTCATg                                                                                                                                                                                                       <  2:110264/250‑1 (MQ=255)                                                     cTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGTGCGTGTGACCAATATTTCGATGGGCAGCGAACTGGATAACGTCATGGAATTTGAAGATATCGc                                                                                                                                                                                      >  1:37034/1‑251 (MQ=255)                                                      tCGCTTTCGAGATTCTCAACTCCAGTGAAGTGAATGGCTACGGTGTTCACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAAATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATTTCGGGTACCGGTTGTCGTGGAAGATATTCTCGTGCTTTTGACCAGTGTTTCGATGGGCAGCGAACTGGATAACGTCATGGAGTTTGatggtactgcc                                                                                                                                                                                     >  2:135338/1‑241 (MQ=255)                                                       cGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGTGCGTGTGACCAATATTTCGATGGGCAGCGAACTGGATAACGTCATGGAATTTGAAGATATCGcc                                                                                                                                                                                     <  1:296908/250‑1 (MQ=255)                                                                           ccAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGTGCGTGTGACCAATATTTCGATGGGCCGCGAACTGGATAACGTCATGGAATTTGGAGATATCGCCGATAAAGCAGCGGAAGCACCg                                                                                                                                                                >  1:23261/1‑251 (MQ=255)                                                                                     gaaTGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGTGCGTGTGACCAATATTTCGATGGGCAGCGAACTGGATAACGTCATGGAATTTGAAGATATCGCCGATAACGCAGCGGACGCACCGACTGAAACCt                                                                                                                                                      <  1:73986/251‑1 (MQ=255)                                                                                         ggCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTGTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGTGCGTGTGACCAATATTTCGATGGGCAGCGAACTGGATAACGTCATGGAATTTGAAGATATCGCCGATAACGCAGCGGACGCACCGACTGAAACCTGCt                                                                                                                                                   <  1:106973/250‑1 (MQ=255)                                                                                                  ggtggccACGAAAAAGTAGCTGAAGTTTTTGTTGGTAAAGCTAATCGGGTCTTCAAACCGGGAGATATTGCGCCTGTCTTTGAACAGGCTAAAGCCTTAATGGCGCAATATCGGGTACCGTGAGTCGTGGAAGTTATTCTCGGGCGTGTGACCAATATTTCGATGGGCAGCGAACTGGATAACGTCATGGAATTTGAAGATATCGCCGATAACGCAGCGGACGCACCGACTGAATCCTGCTTCATGCACt                                                                                                                                          <  2:222703/245‑1 (MQ=255)                                                                                                       ccACGTAAAAGTAGCGGAAGGTGTAGGTTGTAAAGCGAGTCGGGGCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAGCAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGTTCGTGTGACCAATAGTGCGATGGGCAGCGAACTGGATAACGTCATGGAATTTGAAGATATCGCCGATAACGCAGCGGACGCACCGACTGAAACCTTCTTCATGCACTATGAAt                                                                                                                                    <  2:161047/251‑1 (MQ=255)                                                                                                       ccACGAAAAAGTAGGGGAAGGTTAGGGTTGTAAAGTTATCGGGGTCTCCAAACGGTAAAATATGGGGCCATCCTTGGACCAGG‑GTAAGTCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGTGCGTGTGACCAATATTTCGATGGGCAGCGAACTGGATAACGTCATGGAATTTGAAGATATCGCCGATAACGCAGCGGACGCACCGACTGAAACCTGCTTCATGCACTATGAATa                                                                                                                                   <  2:200049/251‑1 (MQ=255)                                                                                                                          ggTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCTGTAGTCGTGGAAGTTATTCTCGTGCGTGTGACCAAGATTTCGATGGGCAGCGAACTGGATAACGTCATGGAATTTTAAGATATCGCCGATAACTCAGCGGACGCACCGACTGAAACCTGCTTCATGCACTATGAATAAGGGAGATAAATAATGtt                                                                                                                 <  1:196991/251‑1 (MQ=255)                                                                                                                                               cGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGTGCGTGTGACCAATATTTCGATGGGCAGCGAACTGGATAACGTCATGGAATTTGAAGATATCGCCGATAACGCAGCGGACGCACCGACTGAAACCTGCGTCATGCTCTATGAATAAGGGAGATAAATAATGTTACGTTTCTCTGCTAATTTAt                                                                                             <  1:215316/250‑1 (MQ=255)                                                                                                                                                    catCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGTGCGTGTGACCAATATTTCGATGGGCAGCGAACTGGATAACGTCATGGAATTTGAAGATATCGCCGATAACGCAGCGGACGCACCGACTGAAACCTGCTTCATGCACTATGAATAAGGGAGATAAATAATGTTACGTTTCTCTGCTAATTTATCGATg                                                                                        <  1:2997/248‑1 (MQ=255)                                                                                                                                                               aaGATATTGCGCCAGCCTTGGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGTGCGTGTGACCAAGAGTTCGATGGGCAGCGAACTGGATAACGTCATGGAATTTGAAGATATCGCCGATAACGCAGCGGACGCAACGACTGAAACCTGCTTCATGCACTATGAATAAGGGAGATAAATAATGTTACGTTTCTCTGCTAATTTATCGATGTTATTTGGAGaa                                                                            <  2:179130/251‑1 (MQ=255)                                                                                                                                                                                                               atagggTACCGGTAGTCGTGGAAGTTATTCTCGTGCGTGTGACCAATATTTCGATGGGCAGCGAACTGGATAACGTCATGGAATTTGAAGATATAGCCGATAACGCAGCGGAAGCACCGACTGAAAACTGCTTCATGCACTATGAATAAGGAAGATAAAAAATGTAAGTTTTCGCTGCAAATTTATCGATGGTCTTTAGAGAATTTGATTTTCTCGCCGGATTTTAGAAAGCTGCGCGGTGTGGtttttgg                            >  2:79418/4‑248 (MQ=255)                                                                                                                                                                                                                                         attcgagAGCGTGTGATCAATATTGCGATGGGCAGAGAACTGGATAACGTCAGGGAATTTGAAGATATCGCCGATAAAGCAGCGGACGCACCGACTGAAACCTGCTTCATGTACTATGAATAAGGGAGATAAATAATGTTACGTTTCTCTGCTAATTTATCGAGGTTATTTGGAGAATAGGATTTTCTCTCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTGCGCGGCGTTGAATTTATGTTTCCTTATGa  <  1:112852/245‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    GGGGCTGATTCGTCAGTCACAACGCGCTTTTGACATGGACTACTGCGTGCAACTCGCTTTCGAGAATATCAACTCCAGTGAAGTGAATGGCTACGGTGTTGACCACGTAAAAGTAGCGGAAGGTTTAGGTTGTAAAGCTATTCGGGTCTTCAAACCGGAAGATATTGCGCCAGCCTTTGAACAGGCGAAAGCCTTAATGGCGCAATATCGGGTACCGGTAGTCGTGGAAGTTATTCTCGAGCGTGTGACCAATATTTCGATGGGCAGCGAACTGGATAACGTCATGGAATTTGAAGATATCGCCGATAACGCAGCGGACGCACCGACTGAAACCTGCTTCATGCACTATGAATAAGGGAGATAAATAATGTTACGTTTCTCTGCTAATTTATCGATGTTATTTGGAGAATATGATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGA  >  NC_000913/535343‑535825

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