Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	39,967	A→G	Y140Y (TAT→TAC)	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	39,967	0	A	G	100.0%	190.8 / NA	54	Y140Y (TAT→TAC)	caiA	crotonobetaine reductase subunit II, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/23);  total (31/23)

TGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGAACCCGCC  >  NC_000913/39822‑40113                                                                                                                                                  |                                                                                                                                                   tGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGt                                                                                                                                               >  2:525065/1‑151 (MQ=255)  gCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGt                                                                                                                                               <  1:438430/150‑1 (MQ=255)     cATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGt                                                                                                                                            <  1:413329/150‑1 (MQ=255)          caacaaACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCa                                                                                                                                       <  1:619323/150‑1 (MQ=255)          caacaaACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGCCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTc                                                                                                                                        >  2:646909/1‑149 (MQ=255)            acaaACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCa                                                                                                                                       >  2:90226/1‑148 (MQ=255)              aaaCCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTa                                                                                                                                  <  2:559691/151‑1 (MQ=255)                          gtAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAg                                                                                                                      <  1:145329/151‑1 (MQ=255)                            aGACAGGATTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGCGTAGGCGCCGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTAGGAGcc                                                                                                                    >  2:518236/1‑151 (MQ=255)                              aCAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCcg                                                                                                                  <  2:504470/151‑1 (MQ=255)                                            agaAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTa                                                                                                    <  2:298826/151‑1 (MQ=255)                                              aaGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAAt                                                                                                  >  2:437120/1‑151 (MQ=255)                                                  cccGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGc                                                                                               >  2:619079/1‑150 (MQ=255)                                                    cGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTg                                                                                             >  1:83746/1‑150 (MQ=255)                                                     gTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGa                                                                                            <  2:532262/150‑1 (MQ=255)                                                       cgcgcgCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGtt                                                                                         >  1:127210/1‑151 (MQ=255)                                                       cgcgcgCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGtt                                                                                         >  1:420958/1‑151 (MQ=255)                                                        gcgcgcCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTc                                                                                        >  2:128107/1‑151 (MQ=255)                                                          gcgcCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCgcgc                                                                                                               <  1:187743/126‑1 (MQ=255)                                                          gcgcCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCgcgc                                                                                                               >  2:187738/1‑126 (MQ=255)                                                              cATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATc                                                                                  >  1:491083/1‑151 (MQ=255)                                                               aTCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCt                                                                                 >  2:642395/1‑151 (MQ=255)                                                                            aCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGcc                                                                    >  2:533349/1‑151 (MQ=255)                                                                             cGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCcg                                                                   <  1:590688/151‑1 (MQ=255)                                                                                         gctgctGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCAt                                                        >  1:423592/1‑150 (MQ=255)                                                                                            gctgGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAAtt                                                    >  2:350023/1‑151 (MQ=255)                                                                                               ggTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAAtttt                                                  <  1:457747/150‑1 (MQ=255)                                                                                                 taataaAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTc                                               >  1:101798/1‑151 (MQ=255)                                                                                                  aataaAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCg                                              <  1:480448/151‑1 (MQ=255)                                                                                                   ataaAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGa                                             <  2:405392/151‑1 (MQ=255)                                                                                                    taaAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGAt                                            >  2:223868/1‑151 (MQ=255)                                                                                                          acTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGctc                                      >  2:483168/1‑151 (MQ=255)                                                                                                           cTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGctct                                     >  1:602851/1‑151 (MQ=255)                                                                                                             tactacCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGctct                                     <  2:594712/149‑1 (MQ=255)                                                                                                              actacCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTtgt                                  <  1:421813/151‑1 (MQ=255)                                                                                                              actacCATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTtg                                   <  1:487978/150‑1 (MQ=255)                                                                                                                  ccATTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGc                              >  2:564902/1‑151 (MQ=255)                                                                                                                    aTTAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCt                            >  1:463107/1‑151 (MQ=255)                                                                                                                     ttAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCtt                           >  2:133201/1‑151 (MQ=255)                                                                                                                      tAAGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTc                          >  1:345088/1‑151 (MQ=255)                                                                                                                        aGATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTcgc                        >  2:473396/1‑151 (MQ=255)                                                                                                                         gATAAATCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTcgcg                       >  1:654074/1‑151 (MQ=255)                                                                                                                            aaaTCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCgcgc                      >  1:323840/1‑149 (MQ=255)                                                                                                                             aaTCTTACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAgg                   <  1:372710/151‑1 (MQ=255)                                                                                                                                  tACCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAgg               <  1:394998/150‑1 (MQ=255)                                                                                                                                   aCCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACTGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGgtg             >  1:384406/1‑151 (MQ=255)                                                                                                                                   aCCATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGgtg             >  2:409677/1‑151 (MQ=255)                                                                                                                                     cATTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGtt           <  2:422488/151‑1 (MQ=255)                                                                                                                                      aTTTCTACGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTg          >  1:5252/1‑151 (MQ=255)                                                                                                                                            aCGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGAAcccg    <  1:499459/151‑1 (MQ=255)                                                                                                                                            aCGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGAAccc     <  2:617851/150‑1 (MQ=255)                                                                                                                                            aCGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGAAcc      <  2:62058/149‑1 (MQ=255)                                                                                                                                             cGGGTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGAAcccgc   <  2:390769/151‑1 (MQ=255)                                                                                                                                              gggTGTAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGAAcccgcc  >  2:638438/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                   TGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCACATCTGCTTACCGGTGCCGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGAACCCGCC  >  NC_000913/39822‑40113

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