Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	40,278	C→T	V37I (GTC→ATC)	caiA ←	crotonobetaine reductase subunit II, FAD‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	40,278	0	C	T	100.0%	221.3 / NA	68	V37I (GTC→ATC)	caiA	crotonobetaine reductase subunit II, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/33);  total (35/33)

ACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGG  >  NC_000913/40169‑40411                                                                                                              |                                                                                                                                      acGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAgg                                                                                                          >  2:983910/1‑139 (MQ=255) acGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAgg                                                                                                          <  1:983910/139‑1 (MQ=255)    gcggcgAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCt                                                                                                       >  1:523665/1‑139 (MQ=255)      ggcgAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTcc                                                                                                     <  2:1311344/139‑1 (MQ=255)       gcgAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTccc                                                                                                    >  1:1913869/1‑139 (MQ=255)          agagTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCa                                                                                                                        <  1:912377/116‑1 (MQ=255)          agagTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCa                                                                                                                        >  2:912377/1‑116 (MQ=255)                 cAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCt                                                                                          <  2:1933955/139‑1 (MQ=255)                          cgTCCAGACCACCGTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAg                                                                                 >  2:1252438/1‑139 (MQ=255)                                       gTGCTCTTCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCg                                                                    >  2:761699/1‑139 (MQ=255)                                             ttCAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACa                                                              <  1:761699/139‑1 (MQ=255)                                               cAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTcc                                                                                                     >  2:1553505/1‑98 (MQ=255)                                               cAGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTcc                                                                                                     <  1:1553505/98‑1 (MQ=255)                                                aGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACACAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAAc                                                           <  1:1252438/139‑1 (MQ=255)                                                aGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAAc                                                           <  2:188554/139‑1 (MQ=255)                                                aGGGATCAGCAGGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAAc                                                           <  2:523665/139‑1 (MQ=255)                                                       agcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCgg                                                                                                                   <  2:1134537/76‑1 (MQ=255)                                                       agcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCgg                                                                                                                   <  2:314247/76‑1 (MQ=255)                                                       agcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCgg                                                                                                                   >  1:314247/1‑76 (MQ=255)                                                       agcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCgg                                                                                                                   >  1:1134537/1‑76 (MQ=255)                                                       agcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGt                                                        <  1:1855083/135‑1 (MQ=255)                                                       agcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGt                                                        >  2:1855083/1‑135 (MQ=255)                                                       agcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCt                                                    >  2:1173300/1‑139 (MQ=255)                                                       agcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCt                                                    >  2:1078577/1‑139 (MQ=255)                                                       agcagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCt                                                    >  1:584178/1‑139 (MQ=255)                                                         cagGCTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGc                                                  >  1:661890/1‑139 (MQ=255)                                                                tCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCgg                                                                                                                   >  1:129148/1‑67 (MQ=255)                                                                tCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCgg                                                                                                                   <  2:129148/67‑1 (MQ=255)                                                                tCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTcatcat                                           >  1:109379/1‑139 (MQ=255)                                                                 cGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAAc                                                           >  2:1317504/1‑122 (MQ=255)                                                                 cGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAAc                                                           <  1:1317504/122‑1 (MQ=255)                                                                 cGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCAtt                                          >  1:1999380/1‑139 (MQ=255)                                                                   aTACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATtta                                        <  2:1913869/139‑1 (MQ=255)                                                                           atCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCgg                                                                                                                   >  1:745373/1‑56 (MQ=255)                                                                           atCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCgg                                                                                                                   <  2:745373/56‑1 (MQ=255)                                                                                 cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCt                                                                                          <  2:1461289/75‑1 (MQ=255)                                                                                 cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCt                                                                                          >  1:1461289/1‑75 (MQ=255)                                                                                 cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCtca                                               <  1:153322/118‑1 (MQ=255)                                                                                 cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCtca                                               >  2:153322/1‑118 (MQ=255)                                                                                 cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCt                                                 >  1:893420/1‑116 (MQ=255)                                                                                 cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCt                                                 <  2:893420/116‑1 (MQ=255)                                                                                 cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCAtt                                          >  1:1496492/1‑123 (MQ=255)                                                                                 cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCAtt                                          <  2:1496492/123‑1 (MQ=255)                                                                                 cAGTGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGACCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATc                          >  1:88183/1‑139 (MQ=255)                                                                                    tGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCAt                                                                                    >  1:14190/1‑78 (MQ=255)                                                                                    tGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCAt                                                                                    <  2:14190/78‑1 (MQ=255)                                                                                    tGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATcc                             >  1:980453/1‑133 (MQ=255)                                                                                    tGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATcc                             <  2:980453/133‑1 (MQ=255)                                                                                    tGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATc                              <  2:1166590/132‑1 (MQ=255)                                                                                    tGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATc                              <  1:516016/132‑1 (MQ=255)                                                                                    tGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATc                              >  1:1166590/1‑132 (MQ=255)                                                                                    tGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATc                              >  2:516016/1‑132 (MQ=255)                                                                                    tGCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCttt                       >  2:1125383/1‑139 (MQ=255)                                                                                     gCTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTc                      <  2:109379/139‑1 (MQ=255)                                                                                      cTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGa                                                                >  1:964582/1‑96 (MQ=255)                                                                                      cTTTGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGa                                                                <  2:964582/96‑1 (MQ=255)                                                                                       tttGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCaa                    <  1:1125383/139‑1 (MQ=255)                                                                                         tGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATAc                                                                      >  2:1942847/1‑87 (MQ=255)                                                                                         tGACAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATAc                                                                      <  1:1942847/87‑1 (MQ=255)                                                                                           aCAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAAc                                                           >  2:534651/1‑96 (MQ=255)                                                                                           aCAAAACGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAAc                                                           <  1:534651/96‑1 (MQ=255)                                                                                              aaaCGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAAcc                  <  2:1979630/134‑1 (MQ=255)                                                                                              aaaCGTTCCGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAAcc                  >  1:1979630/1‑134 (MQ=255)                                                                                                     ccGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCAt                                                                                    >  2:696024/1‑61 (MQ=255)                                                                                                     ccGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCAt                                                                                    <  1:696024/61‑1 (MQ=255)                                                                                                     ccGGGTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATAtttt      >  2:947716/1‑139 (MQ=255)                                                                                                         gTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTgggg  <  2:584178/139‑1 (MQ=255)                                                                                                         gTAGATGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGACCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTgggg  <  2:88183/139‑1 (MQ=255)                                                                                                              |                                                                                                                                      ACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCGCGGATACCGGCGACAAACAGTTCCTGCTCATCATTTAAATTAAAATCCATCTTTCAACCTCTTGATATTTTGGGG  >  NC_000913/40169‑40411

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