Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,307,447	A→T	L118Q (CTG→CAG)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,307,447	0	A	T	100.0%	206.9 / NA	59	L118Q (CTG→CAG)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (30/29);  total (30/29)

CTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCG  >  NC_000913/4307316‑4307582                                                                                                                                    |                                                                                                                                        cTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACg                                                                                                                                 >  2:584381/1‑140 (MQ=255)  tGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGt                                                                                                                                >  2:528112/1‑140 (MQ=255)                 aCCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCaa                                                                                                                 <  1:283040/140‑1 (MQ=255)                 aCCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCaa                                                                                                                 <  2:519235/140‑1 (MQ=255)                 aCCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCaa                                                                                                                 >  2:727388/1‑140 (MQ=255)                      tCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGa                                                                                                            >  2:559158/1‑140 (MQ=255)                      tCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGa                                                                                                            <  2:152155/140‑1 (MQ=255)                      tCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGa                                                                                                            >  2:223983/1‑140 (MQ=255)                             cgcACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCa                                                                                                     >  1:608327/1‑140 (MQ=255)                             cgcACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCa                                                                                                     >  2:606488/1‑140 (MQ=255)                               cACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGc                                                                                                   <  2:445679/140‑1 (MQ=255)                               cACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGc                                                                                                   >  2:346939/1‑140 (MQ=255)                                            cacTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAAtt                                                                                                       <  2:536087/123‑1 (MQ=255)                                            cacTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAAtt                                                                                                       >  1:536087/1‑123 (MQ=255)                                            cacTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCt                                                                                      >  2:276000/1‑140 (MQ=255)                                            cacTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCt                                                                                      >  2:646806/1‑140 (MQ=255)                                              cTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTa                                                                                    >  2:268402/1‑140 (MQ=255)                                              cTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTa                                                                                    >  1:14971/1‑140 (MQ=255)                                                 cGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCg                                                                                 >  2:317574/1‑140 (MQ=255)                                                 cGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCg                                                                                 >  1:592381/1‑140 (MQ=255)                                                 cGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCg                                                                                 >  1:75643/1‑140 (MQ=255)                                                  gAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCgg                                                                                <  1:552328/140‑1 (MQ=255)                                                   aaCCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGc                                                                               <  2:14971/140‑1 (MQ=255)                                                    aCCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCg                                                                              <  1:346939/140‑1 (MQ=255)                                                     ccccATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCga                                                                             <  1:672819/140‑1 (MQ=255)                                                       ccATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCgaga                                                                           <  1:584381/140‑1 (MQ=255)                                                                tACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGGTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATaaa                                                                  >  2:13589/1‑140 (MQ=255)                                                                    gagaTAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATccg                                                              <  1:606488/140‑1 (MQ=255)                                                                     agaTAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATccgc                                                             >  1:764060/1‑140 (MQ=255)                                                                       aTAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATccgccg                                                           <  2:592381/140‑1 (MQ=255)                                                                       aTAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATccgccg                                                           <  2:75643/140‑1 (MQ=255)                                                                       aTAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATccgccg                                                           <  2:268422/140‑1 (MQ=255)                                                                               ccAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAAc                                                   >  2:752795/1‑140 (MQ=255)                                                                                 aGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACgg                                                 <  2:368575/140‑1 (MQ=255)                                                                                      cAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGt                                            <  2:699215/140‑1 (MQ=255)                                                                                      cAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGt                                            <  2:64530/140‑1 (MQ=255)                                                                                       aGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAgggt                                           >  1:50983/1‑138 (MQ=255)                                                                                       aGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGtt                                           <  1:98170/140‑1 (MQ=255)                                                                                       aGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGtt                                           <  1:523015/140‑1 (MQ=255)                                                                                       aGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGtt                                           <  1:444150/140‑1 (MQ=255)                                                                                             tgCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGt                                     >  1:622043/1‑140 (MQ=255)                                                                                             tgCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGt                                     <  2:585879/140‑1 (MQ=255)                                                                                             tgCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGt                                     <  2:14245/140‑1 (MQ=255)                                                                                              gCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGGTAGGCGTa                                    >  2:431868/1‑140 (MQ=255)                                                                                                  aaGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGaaa                                >  1:242225/1‑140 (MQ=255)                                                                                                        gTTTTCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATgg                          <  1:729014/140‑1 (MQ=255)                                                                                                          tttCTACTACGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTg                        <  2:622043/140‑1 (MQ=255)                                                                                                              tactacGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGtt                    <  1:559158/140‑1 (MQ=255)                                                                                                                 tacGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTAt                 <  2:296954/140‑1 (MQ=255)                                                                                                                 tacGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTAt                 <  1:727388/140‑1 (MQ=255)                                                                                                                 tacGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTAt                 <  1:528112/140‑1 (MQ=255)                                                                                                                 tacGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTAt                 <  2:764060/140‑1 (MQ=255)                                                                                                                   cGTCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGTGGAGAAATAAATGGGGGTTTATCt               >  2:365788/1‑140 (MQ=255)                                                                                                                     tCCCAGGAGAGTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGGAGAAATAATGGGGCGGTTATtttt             >  2:108563/1‑140 (MQ=255)                                                                                                                              agTTGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAg    >  1:697485/1‑140 (MQ=255)                                                                                                                                ttGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAgcg  >  2:638916/1‑140 (MQ=255)                                                                                                                                ttGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAgcg  >  1:368853/1‑140 (MQ=255)                                                                                                                                ttGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAAGCGTAGAAATAATGGGGCGTTTATCTTTACTGACCAgcc  >  1:530257/1‑139 (MQ=255)                                                                                                                                ttGCTGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGGTAGGGGGAGAAAAAATGGGGGGTTTATCTTTACTGACCaccc  >  2:181038/1‑137 (MQ=255)                                                                                                                                    |                                                                                                                                        CTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTCTACTACGTCCCAGGAGAGTTGCAGGTTAACGTCGCGGGAAAACTCAACCGGACAATTCAGCGTATTTTCGAGCTTATCGGCGAGATCATATAAATCCGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCG  >  NC_000913/4307316‑4307582

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: ‑