Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	354,592	G→A	*629* (TAG→TAA)	prpE →	propionate‑‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	354,592	0	G	A	98.5%	205.9 / NA	65	*629* (TAG→TAA)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (38/26);  minor base T (0/1);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACA  >  NC_000913/354485‑354721                                                                                                            |                                                                                                                                  ctgcGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATcc                                                                                                    >  1:47865/1‑139 (MQ=255)          aCGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACtgtg                                                                                           <  1:1661375/139‑1 (MQ=255)                        ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCgg                                                                                                  >  1:2002252/1‑118 (MQ=255)                        ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCgg                                                                                                  <  2:2002252/118‑1 (MQ=255)                        ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCa                                                                                                <  2:1676290/120‑1 (MQ=255)                        ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCa                                                                                                >  1:1676290/1‑120 (MQ=255)                        ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGAt                                                                             >  1:1833910/1‑139 (MQ=255)                        ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGAt                                                                             >  2:1419376/1‑139 (MQ=255)                        ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGAt                                                                             >  2:1443813/1‑139 (MQ=255)                        ttGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGAt                                                                             >  2:1723794/1‑139 (MQ=255)                                 aCGCGATTCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGATTAAGTAATTGTCGGATGCGTCGCACGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGCCGCTg                                                                    <  1:634091/139‑1 (MQ=255)                                  cgcgATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGa                                                                   <  1:1762923/139‑1 (MQ=255)                                  cgcgATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGa                                                                   <  2:1833910/139‑1 (MQ=255)                                   gcgATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAc                                                                  <  1:1443813/139‑1 (MQ=255)                                         cTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGttt                                                            <  2:47865/139‑1 (MQ=255)                                            gagaTCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTATGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCCATGCCTGATCCGACGCTGACGCGTTTTAt                                                         <  2:1373030/139‑1 (MQ=255)                                            gagaTCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTAt                                                         <  2:1762772/139‑1 (MQ=255)                                                 cTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCa                                                                                                >  2:1957228/1‑95 (MQ=255)                                                 cTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCa                                                                                                <  1:1957228/95‑1 (MQ=255)                                                 cTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGc                                                    >  1:2028278/1‑139 (MQ=255)                                                 cTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGc                                                    >  2:1173971/1‑139 (MQ=255)                                                 cTGACGACCATTGATAATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCGTCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGc                                                    >  1:1373030/1‑139 (MQ=255)                                                   gacgacCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCt                                                  >  2:682088/1‑139 (MQ=255)                                                       acCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGa                                                                   >  1:1712326/1‑118 (MQ=255)                                                       acCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGa                                                                   <  2:1712326/118‑1 (MQ=255)                                                           ttGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAcgcg                                                               >  2:571592/1‑118 (MQ=255)                                                           ttGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAcgcg                                                               <  1:571592/118‑1 (MQ=255)                                                           ttGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCt                                          >  1:1824686/1‑139 (MQ=255)                                                              atgatCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGaa                                       <  1:682088/139‑1 (MQ=255)                                                              atgatCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGaa                                       >  1:1933586/1‑139 (MQ=255)                                                                   ccTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCt                                          <  1:1126918/131‑1 (MQ=255)                                                                   ccTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCt                                          >  2:1126918/1‑131 (MQ=255)                                                                   ccTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTa                                  >  1:268511/1‑139 (MQ=255)                                                                    cTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTAtt                                                                                                                             <  1:1662324/47‑1 (MQ=255)                                                                    cTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTAtt                                                                                                                             >  2:1662324/1‑47 (MQ=255)                                                                      gCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGt                               >  1:114898/1‑139 (MQ=255)                                                                       cgtcgtTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTc                              >  1:1276083/1‑139 (MQ=255)                                                                        gtcgtTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCg                             >  1:1809404/1‑139 (MQ=255)                                                                        gtcgtTGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCg                             >  2:1386092/1‑139 (MQ=255)                                                                            ttGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACtg                                                                                             <  2:1125622/71‑1 (MQ=255)                                                                            ttGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACtg                                                                                             >  1:1125622/1‑71 (MQ=255)                                                                            ttGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATa                         >  1:112356/1‑139 (MQ=255)                                                                            ttGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATa                         >  1:1673163/1‑139 (MQ=255)                                                                             tGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGa                                                                   >  1:1755121/1‑96 (MQ=255)                                                                             tGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGa                                                                   <  2:1755121/96‑1 (MQ=255)                                                                             tGGATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATaa                        >  2:665510/1‑139 (MQ=255)                                                                              ggATCAGATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAg                       >  1:873549/1‑139 (MQ=255)                                                                                    gATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGTATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTTATGCGACGATg                                                                    <  1:1192973/88‑1 (MQ=255)                                                                                    gATCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTg                                                                    >  2:1192973/1‑88 (MQ=255)                                                                                     aTCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGTCCTGAACCGt                                   >  2:936111/1‑120 (MQ=255)                                                                                     aTCCGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGt                                   <  1:936111/120‑1 (MQ=255)                                                                                       ccGCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTcgcg              <  2:1673163/139‑1 (MQ=255)                                                                                         gCCAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTc            >  1:1529664/1‑139 (MQ=255)                                                                                          ccAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGa                                                                   >  1:92780/1‑83 (MQ=255)                                                                                          ccAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGa                                                                   <  2:92780/83‑1 (MQ=255)                                                                                          ccAGGCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAc                                                                  >  2:53845/1‑84 (MQ=255)                                                                                          ccAGGCGATGGAAGAGTAAATTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGAc                                                                  <  1:53845/84‑1 (MQ=255)                                                                                              gCGATGGAAGAGTCAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTa                                                 <  2:1095210/97‑1 (MQ=255)                                                                                              gCGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTa                                                 >  1:1095210/1‑97 (MQ=255)                                                                                               cGATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATcc      >  2:1861279/1‑139 (MQ=255)                                                                                                gATGGAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCg     <  1:1173971/139‑1 (MQ=255)                                                                                                   ggAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGTACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGaca  <  1:1419376/139‑1 (MQ=255)                                                                                                   ggAAGAGTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGaca  >  2:1150276/1‑139 (MQ=255)                                                                                                        agTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCa                                                                                                <  2:1809171/40‑1 (MQ=255)                                                                                                        agTAAGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCa                                                                                                >  1:1809171/1‑40 (MQ=255)                                                                                                            |                                                                                                                                  CTGCGCCGCACGATCCAGGCGATTTGCGAAGGACGCGATCCTGGAGATCTGACGACCATTGATGATCCTGCGTCGTTGGATCAGATCCGCCAGGCGATGGAAGAGTAGGTTATTGTCGGATGCGTCGCGCGGTGCATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACA  >  NC_000913/354485‑354721

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