Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,873,113	T→C	D167G (GAC→GGC)	dgoA ←	2‑dehydro‑3‑deoxy‑6‑phosphogalactonate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,873,113	0	T	C	100.0%	113.3 / NA	33	D167G (GAC→GGC)	dgoA	2‑dehydro‑3‑deoxy‑6‑phosphogalactonate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/12);  total (21/12)

GGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTG  >  NC_000913/3872978‑3873246                                                                                                                                        |                                                                                                                                      ggTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTAt                                                                                                                                   >  1:476860/1‑140 (MQ=255)    aaTTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATcga                                                                                                                                >  1:821375/1‑138 (MQ=255)                            gCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGGCCCTGCACAACCTGTGCCTATCCACGGCGCAAGGTTTTCTGGCGTCAc                                                                                                        >  2:350268/1‑140 (MQ=255)                             ccTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTTCCCCTATCCACTTCGCCGTGTTTTCTGGCGTCAcg                                                                                                       >  2:359531/1‑140 (MQ=255)                             ccTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCAcg                                                                                                       >  2:706457/1‑140 (MQ=255)                             ccTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCAcg                                                                                                       >  1:528270/1‑140 (MQ=255)                             ccTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCAcg                                                                                                       >  2:101719/1‑140 (MQ=255)                                  cGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACgccgcc                                                                                                  >  1:556131/1‑140 (MQ=255)                                       cGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACgg                                                                                             >  1:548521/1‑140 (MQ=255)                                        gCCTTAACAAATGCTCCTGCCTTCTGCGCTGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGc                                                                                            <  1:350268/140‑1 (MQ=255)                                          cTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCaa                                                                                          <  2:548521/140‑1 (MQ=255)                                                        ctgcctgctgcGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                                            <  1:101719/140‑1 (MQ=255)                                                        ctccctgctgcGCGGTGGGCTCTACGGATTGCCCGGCGCAATAGAGATGGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTcc                                                                            <  2:306362/137‑1 (MQ=255)                                                          gcctgctgcGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGa                                                                          <  1:706457/140‑1 (MQ=255)                                                              gctgcGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGGCAAAGGCAAAGACGGCGATGTCCGATGGc                                                                      >  1:374960/1‑140 (MQ=255)                                                                         cgcTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCAAACTGTATGGACTGCGATGTCTGATGGCAATACCGCttt                                                           >  2:399869/1‑140 (MQ=255)                                                                           ctctACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTa                                                         <  1:571655/140‑1 (MQ=255)                                                                                cGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCt                                                    <  2:528270/140‑1 (MQ=255)                                                                                 ggATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCtt                                                   >  2:52469/1‑140 (MQ=255)                                                                                   aTTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGCCCGTTGGCAATACAGCTTTTAACGCTTTg                                                 >  1:290489/1‑140 (MQ=255)                                                                                     tGCCCGGCCCGATAGAGATCGCCGCCTAAGCGCACACCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGAt                                               <  2:374960/140‑1 (MQ=255)                                                                                            cgcgATGGAGATCGCTGCCTAATCCCGCCTCTGCACAGCCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGc                                        <  1:359531/140‑1 (MQ=255)                                                                                                           tGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGAt                         >  2:144283/1‑140 (MQ=255)                                                                                                           tGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGAt                         >  2:532745/1‑140 (MQ=255)                                                                                                           tGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGAt                         >  2:634960/1‑140 (MQ=255)                                                                                                           tGCCTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGAt                         <  2:692066/140‑1 (MQ=255)                                                                                                             ccTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGa                       <  2:556131/140‑1 (MQ=255)                                                                                                             ccTAAGCCCGCCCCTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGa                       <  2:821375/140‑1 (MQ=255)                                                                                                                       ccccTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGTTGTATTGCGGACCAAAAACCGATGACGGAAATAtt             >  2:278277/1‑140 (MQ=255)                                                                                                                       ccccTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGa                       <  1:786685/130‑1 (MQ=255)                                                                                                                       ccccTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGa                       >  2:786685/1‑130 (MQ=255)                                                                                                                       ccccTGCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATAtt             >  2:516047/1‑140 (MQ=255)                                                                                                                            gCACAACCTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGAATTGCGGTCCATAAGCCGATGTCGGAAATATTTTCAg        >  1:340338/1‑140 (MQ=255)                                                                                                                                  ccTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTg  >  2:547496/1‑140 (MQ=255)                                                                                                                                  ccTGCGCCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTg  >  2:590104/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                      GGTAATTTTGGTGATTTTCATTGCACTGCCTCTCGATACGCCTTAACAAATGCTGCTGCCTGCTGCGCGGTGCGCTCTACGGATTGCCCGGCGCGATAGAGATCGCTGCCTAAGCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTG  >  NC_000913/3872978‑3873246

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