Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,096,517	T→C	Q311R (CAA→CGA)	rhaB ←	rhamnulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,096,517	0	T	C	100.0%	158.1 / NA	45	Q311R (CAA→CGA)	rhaB	rhamnulokinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/24);  total (21/24)

GTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTTGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGCCAGTGCCGTGTCATTGGTAAATGGCGTC  >  NC_000913/4096382‑4096645                                                                                                                                        |                                                                                                                                 gTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTcgct                                                                                                                               <  2:1010620/139‑1 (MQ=255) gTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTcgct                                                                                                                               <  1:59726/139‑1 (MQ=255)             gCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCa                                                                                                                   <  2:180830/139‑1 (MQ=255)                        ttCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAagc                                                                                                        <  1:422268/139‑1 (MQ=255)                             tGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCa                                                                                                   >  1:132502/1‑139 (MQ=255)                                     gTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCa                                                                                           <  2:874635/139‑1 (MQ=255)                                     gTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCa                                                                                           <  2:896588/139‑1 (MQ=255)                                                aataaAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATaa                                                                                                >  2:841487/1‑123 (MQ=255)                                                aataaAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATaa                                                                                                <  1:841487/123‑1 (MQ=255)                                                aataaAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTc                                                                                >  1:21409/1‑139 (MQ=255)                                                 ataaAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCa                                                                               >  2:602965/1‑139 (MQ=255)                                                      gCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCAt                                                                                                  <  1:62519/115‑1 (MQ=255)                                                      gCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCAt                                                                                                  >  2:62519/1‑115 (MQ=255)                                                      gCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCAcc                                                                          >  1:564888/1‑139 (MQ=255)                                                       cGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCAccc                                                                         >  1:478786/1‑139 (MQ=255)                                                         aTCGTAATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGTAGCCATAAGCCCATAATATTTTTCAGCACCCGa                                                                       >  1:744250/1‑139 (MQ=255)                                                                    gATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCg                                                                  <  2:909577/133‑1 (MQ=255)                                                                    gATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCg                                                                  >  1:909577/1‑133 (MQ=255)                                                                       tGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTccgc                                                         <  1:442206/139‑1 (MQ=255)                                                                        gATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTccgcc                                                        >  2:663322/1‑139 (MQ=255)                                                                        gATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTccgcc                                                        >  1:859848/1‑139 (MQ=255)                                                                        gATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTccgcc                                                        >  1:838033/1‑139 (MQ=255)                                                                         aTAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTccgccc                                                       >  2:163621/1‑139 (MQ=255)                                                                             tGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTccgccccgcc                                                   <  2:21409/139‑1 (MQ=255)                                                                                   ggCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTccgcccc                                                      <  2:769236/130‑1 (MQ=255)                                                                                   ggCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTccgcccc                                                      >  1:769236/1‑130 (MQ=255)                                                                                   ggCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCAtt                                             >  2:164407/1‑139 (MQ=255)                                                                                    gCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTc                                                <  2:516523/135‑1 (MQ=255)                                                                                    gCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTc                                                >  1:516523/1‑135 (MQ=255)                                                                                            gAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGtt                                    <  1:397894/139‑1 (MQ=255)                                                                                              aGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCATTCGCTGAAGCAGCCACAAGCCCATAATATTTTTCAGCACCCGATAGCGTCCTTCCGCCACGCCTTCATTGGTGATGTTgg                                  <  2:132502/139‑1 (MQ=255)                                                                                                  ccTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGc                              >  1:171104/1‑139 (MQ=255)                                                                                                   cTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCAcc                                                                          <  2:962416/94‑1 (MQ=255)                                                                                                   cTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCAcc                                                                          >  1:962416/1‑94 (MQ=255)                                                                                                   cTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGcc                             <  2:859848/139‑1 (MQ=255)                                                                                                   cTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGcc                             <  1:163621/139‑1 (MQ=255)                                                                                                    tgtgTCGCGGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGCAGTAGCCATAAGCCCTTAATATTTTTCAGCACCCGACAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGCCa                            <  2:744250/139‑1 (MQ=255)                                                                                                         cgcgGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGCCAGTGcc                       <  1:663322/139‑1 (MQ=255)                                                                                                          gcgGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGCCAGTGCCg                      <  2:838033/139‑1 (MQ=255)                                                                                                           cgGAGATAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGCCAGTGCCgt                     <  2:898889/139‑1 (MQ=255)                                                                                                               gaTAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGCCAGTGCCGTGTCa                 <  2:1988/139‑1 (MQ=255)                                                                                                                aTAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGCCAGTGcc                       >  2:334052/1‑132 (MQ=255)                                                                                                                aTAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGCCAGTGcc                       <  1:334052/132‑1 (MQ=255)                                                                                                                aTAAGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGCCAGTGCCGTGTCAt                <  1:90209/139‑1 (MQ=255)                                                                                                                  aaGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGCCAGTGCCgtgt                   >  2:527934/1‑134 (MQ=255)                                                                                                                  aaGCGCCGGAAGATCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGCCAGTGCCgtgt                   <  1:527934/134‑1 (MQ=255)                                                                                                                              aTCGTTGATTCGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGCCAGTGCCGTGTCATTGGTAAATGGCGTc  >  2:873079/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                 GTTTCCCGACACGCAGCCTGAATTTCGCTGCACATCGTCTCAGGATTAATAAAGCGATCGTCATTGGGATTGATAATGAAGCGGCAAGCCGGAAGTGCCTGTGTCGCGGAGATAAGCGCCGGAAGATCGTTGATTTGCTGCTCCTGAAGCACTCGCTGAAGCAGCCATAAGCCCATAATATTTTTCAGCACCCGATAGCGACCTTCCGCCCCGCCTTCATTGGTGATGTTGGCTGCCAGTGCCGTGTCATTGGTAAATGGCGTC  >  NC_000913/4096382‑4096645

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