Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,873,782	G→A	G231G (GGC→GGT)	dgoK ←	2‑dehydro‑3‑deoxygalactonokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,873,782	0	G	A	100.0%	203.3 / NA	59	G231G (GGC→GGT)	dgoK	2‑dehydro‑3‑deoxygalactonokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/33);  total (26/33)

GTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAA  >  NC_000913/3873645‑3873916                                                                                                                                          |                                                                                                                                       gtgtCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCAcc                                                                                                                                      >  1:143925/1‑140 (MQ=255) gtgtCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCAcc                                                                                                                                      <  2:707959/140‑1 (MQ=255)         cgccACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATcaaca                                                                                                                              <  1:722451/140‑1 (MQ=255)                       cGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAAt                                                                                                                <  2:143925/140‑1 (MQ=255)                       cGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAAt                                                                                                                <  1:40970/140‑1 (MQ=255)                          cGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCg                                                                                                             <  1:247621/140‑1 (MQ=255)                               cccATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGAc                                                                                                        <  2:858926/140‑1 (MQ=255)                                  aTCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTg                                                                                                     >  2:715257/1‑140 (MQ=255)                                    cGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGtt                                                                                                   >  2:254839/1‑140 (MQ=255)                                     gCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTc                                                                                                  <  2:711531/140‑1 (MQ=255)                                          aaaGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCggg                                                                                             >  1:781517/1‑140 (MQ=255)                                             ggCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGaa                                                                                          <  2:564580/140‑1 (MQ=255)                                               cTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCggg                                                                                             <  1:346599/135‑1 (MQ=255)                                               cTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCggg                                                                                             >  2:346599/1‑135 (MQ=255)                                               cTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAgt                                                                                        <  1:715257/140‑1 (MQ=255)                                               cTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGT‑‑CGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAgtgt                                                                                      >  2:627700/1‑140 (MQ=255)                                                ttGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAgtg                                                                                       >  2:1070093/1‑140 (MQ=255)                                                  gCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGtt                                                                                     <  1:480309/140‑1 (MQ=255)                                                    tGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTcc                                                                                   <  2:136147/140‑1 (MQ=255)                                                     ggTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTccc                                                                                  <  1:266815/140‑1 (MQ=255)                                                         gcgcgcGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGca                                                                              <  1:31683/140‑1 (MQ=255)                                                            cgcgGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACAt                                                                           >  2:618708/1‑140 (MQ=255)                                                                        tGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGaa                                                               <  2:474308/140‑1 (MQ=255)                                                                           tCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACtt                                                            >  2:727360/1‑140 (MQ=255)                                                                               gCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCaaa                                                        <  1:254839/140‑1 (MQ=255)                                                                                    aaGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCt                                                   >  2:949116/1‑140 (MQ=255)                                                                                    aaGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCt                                                   >  1:965710/1‑140 (MQ=255)                                                                                      gggTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGc                                                 >  2:438365/1‑140 (MQ=255)                                                                                        gTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCgg                                               <  1:85234/140‑1 (MQ=255)                                                                                                 tgtTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggc                                      >  1:740309/1‑140 (MQ=255)                                                                                                 tgtTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACAAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggc                                      >  2:776863/1‑140 (MQ=255)                                                                                                  gtTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcg                                     >  2:603599/1‑140 (MQ=255)                                                                                                  gtTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcg                                     <  1:603599/140‑1 (MQ=255)                                                                                                  gtTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcg                                     <  2:504533/140‑1 (MQ=255)                                                                                                  gtTGATGGGCCACATAGT‑‑CGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcggg                                   <  1:627700/140‑1 (MQ=255)                                                                                                    tGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcggg                                   <  1:1070093/140‑1 (MQ=255)                                                                                                    tGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcggg                                   <  1:1015152/140‑1 (MQ=255)                                                                                                    tGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcggg                                   <  1:618708/140‑1 (MQ=255)                                                                                                    tGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcggg                                   <  1:727360/140‑1 (MQ=255)                                                                                                           ccACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATgg                                       <  2:487193/129‑1 (MQ=255)                                                                                                           ccACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATgg                                       >  1:487193/1‑129 (MQ=255)                                                                                                           ccACATAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGGATTa                            >  2:836928/1‑140 (MQ=255)                                                                                                             acaTAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcgggc                                  >  1:965619/1‑132 (MQ=255)                                                                                                             acaTAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcgggc                                  <  2:965619/132‑1 (MQ=255)                                                                                                              caTAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcggg                                   <  2:636668/130‑1 (MQ=255)                                                                                                              caTAGTCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcggg                                   >  1:636668/1‑130 (MQ=255)                                                                                                                   tCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcggg                                   >  1:812592/1‑125 (MQ=255)                                                                                                                   tCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcggg                                   <  1:975074/125‑1 (MQ=255)                                                                                                                   tCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcggg                                   >  2:975074/1‑125 (MQ=255)                                                                                                                   tCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcggg                                   <  2:812592/125‑1 (MQ=255)                                                                                                                   tCGCGCATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGGATTAAAACCAcc                    >  2:939002/1‑139 (MQ=255)                                                                                                                      cgcATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcggg                                   >  2:259722/1‑122 (MQ=255)                                                                                                                      cgcATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATggcggg                                   <  1:259722/122‑1 (MQ=255)                                                                                                                        cATACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGGATTTAAACCACCCCCaa               >  1:433500/1‑140 (MQ=255)                                                                                                                          tACTGGCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAgg             <  2:62216/140‑1 (MQ=255)                                                                                                                              ggCGACCTGTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAg         <  2:781517/140‑1 (MQ=255)                                                                                                                               gCGACCTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGGATTAAAACCACGCCCCAAGCCAGc        >  2:337228/1‑140 (MQ=255)                                                                                                                                   ccTCTGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAAACCACGCCCAAGGCCAGCtgtg    >  1:298469/1‑140 (MQ=255)                                                                                                                                     tctGCACCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGaa  <  2:965710/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                       GTGTCGCCCGCCACCGCCGTCACGTCGCAACCCATCGCCTGAAAGGCTTGCTGGTAGCGCGCGGTCAGCGATGTTCCGGCGACAAGGGTGATGGCGTGTTGATGGGCCACATAGTCGCGCATACTGGCGACCTCTGCGCCAATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAA  >  NC_000913/3873645‑3873916

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