Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,937,802	A→G	D112G (GAC→GGC)	fucK →	L‑fuculokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,937,802	0	A	G	98.6%	245.9 / NA	69	D112G (GAC→GGC)	fucK	L‑fuculokinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (34/34);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCA  >  NC_000913/2937669‑2937931                                                                                                                                      |                                                                                                                                  gAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAata                                                                                                                              <  1:1775034/139‑1 (MQ=255)    ctgactgaATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTg                                                                                                                           <  1:1440549/139‑1 (MQ=255)     tgactgaATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGa                                                                                                                          >  2:178368/1‑139 (MQ=255)         tgaATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACgg                                                                                                                      >  2:1100/1‑139 (MQ=255)          gaATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCa                                                                                                                                 >  2:466038/1‑127 (MQ=255)          gaATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGt                                                                                                                     >  1:1360477/1‑139 (MQ=255)                   aTCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCg                                                                                                            <  2:1360477/139‑1 (MQ=255)                                 cGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGaa                                                                                                                         <  2:1384666/112‑1 (MQ=255)                                 cGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGaa                                                                                                                         >  1:1384666/1‑112 (MQ=255)                                   tCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGc                                                                                            >  2:735498/1‑139 (MQ=255)                                     accaccTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTa                                                                                          >  2:440749/1‑139 (MQ=255)                                               gtgtGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGc                                                                                                >  2:418463/1‑123 (MQ=255)                                               gtgtGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGc                                                                                                <  1:418463/123‑1 (MQ=255)                                                tgtgGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTc                                                                               <  1:914835/139‑1 (MQ=255)                                                  tgGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCgg                                                                             <  1:1100/139‑1 (MQ=255)                                                         cgcTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCttt                                                                      >  2:898648/1‑139 (MQ=255)                                                         cgcTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCttt                                                                      >  2:814804/1‑139 (MQ=255)                                                            tctGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGt                                                                   <  2:512294/139‑1 (MQ=255)                                                                gTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                                               <  1:178368/139‑1 (MQ=255)                                                                gTAGATAAGCAAGGCAATCTACTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGtttaa                                                               >  1:1055228/1‑137 (MQ=255)                                                                   gATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGAGGTTGCAGGCTATTTCTGGCGTCGGAGCCTGTAGTTTCAATa                                                            <  1:440749/139‑1 (MQ=255)                                                                     tAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACg                                                          <  2:1703888/139‑1 (MQ=255)                                                                     tAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACg                                                          <  1:814804/139‑1 (MQ=255)                                                                     tAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACg                                                          <  1:898648/139‑1 (MQ=255)                                                                       aGCAAGGCAATCTGCTCTATCTGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGtt                                                        <  2:1055228/139‑1 (MQ=255)                                                                        gCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAg                                                                                              <  1:682896/100‑1 (MQ=255)                                                                        gCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAg                                                                                              >  2:682896/1‑100 (MQ=255)                                                                        gCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACg                                                          >  1:617158/1‑136 (MQ=255)                                                                        gCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACg                                                          <  2:617158/136‑1 (MQ=255)                                                                          aaGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATa                                                            >  1:1689330/1‑132 (MQ=255)                                                                          aaGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATa                                                            <  2:1689330/132‑1 (MQ=255)                                                                          aaGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTtata                                                     >  2:1765427/1‑139 (MQ=255)                                                                          aaGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTtata                                                     >  2:698589/1‑139 (MQ=255)                                                                           aGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGtt                                                                                                                    >  1:97166/1‑75 (MQ=255)                                                                           aGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGtt                                                                                                                    <  2:97166/75‑1 (MQ=255)                                                                                 tCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAg                                                                           <  2:1072611/110‑1 (MQ=255)                                                                                 tCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAg                                                                           >  1:1072611/1‑110 (MQ=255)                                                                                 tCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTg                                              >  1:1117563/1‑139 (MQ=255)                                                                                               ttattaGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATa                                                            <  1:1798244/111‑1 (MQ=255)                                                                                               ttattaGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATa                                                            >  2:1798244/1‑111 (MQ=255)                                                                                               ttattaGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATAc                                                           >  1:1671338/1‑112 (MQ=255)                                                                                               ttattaGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATAc                                                           <  2:1671338/112‑1 (MQ=255)                                                                                               ttattaGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTgaaaga                                >  2:253052/1‑139 (MQ=255)                                                                                                    aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAg                                                                                                       <  2:1237871/63‑1 (MQ=255)                                                                                                    aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAg                                                                                                       >  1:1237871/1‑63 (MQ=255)                                                                                                    aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTg                                                                                                 <  2:1618201/69‑1 (MQ=255)                                                                                                    aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTg                                                                                                 >  1:1618201/1‑69 (MQ=255)                                                                                                    aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                                               >  2:1726331/1‑103 (MQ=255)                                                                                                    aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCa                                                               <  1:1726331/103‑1 (MQ=255)                                                                                                    aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATa                                                            >  1:550979/1‑106 (MQ=255)                                                                                                    aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATa                                                            <  2:550979/106‑1 (MQ=255)                                                                                                    aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACg                                                          <  2:1513050/108‑1 (MQ=255)                                                                                                    aGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACg                                                          >  1:1513050/1‑108 (MQ=255)                                                                                                         gAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCa                                                                                               >  1:733772/1‑66 (MQ=255)                                                                                                         gAAATGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCa                                                                                               <  2:733772/66‑1 (MQ=255)                                                                                                         gAAATGTCCGCGAACAGCAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCca                      >  2:1201145/1‑139 (MQ=255)                                                                                                          aaaTGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCcac                     <  1:1765427/139‑1 (MQ=255)                                                                                                           aaTGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGca                                                                                                          <  1:969073/53‑1 (MQ=255)                                                                                                           aaTGTCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGca                                                                                                          >  2:969073/1‑53 (MQ=255)                                                                                                               tCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATa                                                            <  1:954998/95‑1 (MQ=255)                                                                                                               tCCGCGAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATa                                                            >  2:954998/1‑95 (MQ=255)                                                                                                                 cgcgAACAGCAGCGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAActgct              >  2:186246/1‑139 (MQ=255)                                                                                                                        agcagcGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTgaaaga                                >  2:756855/1‑114 (MQ=255)                                                                                                                        agcagcGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTgaaaga                                <  1:756855/114‑1 (MQ=255)                                                                                                                             cGGTTATGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCa  <  1:698589/139‑1 (MQ=255)                                                                                                                                  aTGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCgg                                                                             >  1:165980/1‑59 (MQ=255)                                                                                                                                  aTGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCgg                                                                             <  2:165980/59‑1 (MQ=255)                                                                                                                                  aTGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCg                                                                              >  2:1170064/1‑58 (MQ=255)                                                                                                                                  aTGGGCAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCg                                                                              <  1:1170064/58‑1 (MQ=255)                                                                                                                                      |                                                                                                                                  GAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATGGCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCA  >  NC_000913/2937669‑2937931

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