Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	39,838	C→T	W183* (TGG→TGA)	caiA ←	crotonobetainyl‑CoA reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	39,838	0	C	T	100.0%	128.3 / NA	39	W183* (TGG→TGA)	caiA	crotonobetainyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/26);  total (13/26)

CGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGT  >  NC_000913/39702‑39972                                                                                                                                         |                                                                                                                                       cGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCAt                                                                                                                                     <  1:364942/140‑1 (MQ=255)   ttACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTc                                                                                                                                   >  2:509457/1‑140 (MQ=255)              accgaACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAg                                                                                                                        >  1:263535/1‑140 (MQ=255)                    cATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGt                                                                                                                  <  1:21098/140‑1 (MQ=255)                      tGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTcc                                                                                                                <  2:448404/140‑1 (MQ=255)                          tttCTCGTTCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGgag                                                                                                            <  1:621988/140‑1 (MQ=255)                               cGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATACGGAGTTGACAGGATTGTCCGGAGAAGcc                                                                                                       >  1:222099/1‑140 (MQ=255)                                                tcAAAGGTGATTTCACAGCAGCTATGCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATcac                                                                                      >  2:453600/1‑140 (MQ=255)                                                  aaaGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCGTTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATcacca                                                                                    <  2:112100/140‑1 (MQ=255)                                                   aaGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATcaccac                                                                                   <  1:453600/140‑1 (MQ=255)                                                   aaGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATcaccac                                                                                   <  2:263535/140‑1 (MQ=255)                                                   aaGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATcaccac                                                                                   <  1:86022/140‑1 (MQ=255)                                                    aGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACg                                                                                  <  2:150825/140‑1 (MQ=255)                                                        gATTTCACAGCAGCTCTCCATACGCCGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGt                                                                              <  1:171208/140‑1 (MQ=255)                                                        gATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGt                                                                              >  1:603198/1‑140 (MQ=255)                                                         aTTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTa                                                                             >  1:77629/1‑140 (MQ=255)                                                         aTTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTa                                                                             <  2:133615/140‑1 (MQ=255)                                                         aTTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTa                                                                             <  2:603198/140‑1 (MQ=255)                                                         aTTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGACCGGAGATGCCGCGACGCGCGGCATCACCACGATGTa                                                                             >  1:176633/1‑140 (MQ=255)                                                           ttCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACg                                                                           <  1:141905/140‑1 (MQ=255)                                                           ttCACAGCAGCTATCCACACGCTGACCGAGCTTTTCAAGCTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACg                                                                           <  2:172950/140‑1 (MQ=255)                                                                           aTACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCgctgct                                                           <  2:606169/140‑1 (MQ=255)                                                                           aTACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCgctgct                                                           <  1:509457/140‑1 (MQ=255)                                                                             aCGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTgg                                                         <  1:508420/140‑1 (MQ=255)                                                                               gCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGta                                                       >  2:587376/1‑140 (MQ=255)                                                                                cAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGtaa                                                      >  2:401683/1‑140 (MQ=255)                                                                                   aCCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGtaataa                                                   >  2:105181/1‑140 (MQ=255)                                                                                   aCCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGtaataa                                                   >  1:542591/1‑140 (MQ=255)                                                                                     cGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATaaaa                                                 <  2:356417/140‑1 (MQ=255)                                                                                       aGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAaca                                               >  2:517470/1‑140 (MQ=255)                                                                                              cAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTtacta                                        <  1:22296/140‑1 (MQ=255)                                                                                                  tttGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCAt                                    <  2:149219/140‑1 (MQ=255)                                                                                                  tttGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCAt                                    <  1:26850/140‑1 (MQ=255)                                                                                                  tttGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCAt                                    <  1:105181/140‑1 (MQ=255)                                                                                                  tttGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCAt                                    <  2:77629/140‑1 (MQ=255)                                                                                                      gTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAg                                <  1:518575/140‑1 (MQ=255)                                                                                                      gTCACTTTGATGCCCGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAg                                >  2:518575/1‑140 (MQ=255)                                                                                                                  cccGGTTTGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTAcc                    <  1:587376/140‑1 (MQ=255)                                                                                                                         tGCTCATATCAACAAATCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTa             <  1:456370/140‑1 (MQ=255)                                                                                                                                    acaaATCATTCGGTGTATACAGCTTTGTCCGAAGACGCCTCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGt  <  2:208435/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       CGTTACCTTCCCGACCGAACATGTCTTTCTCGTCCAGTTCCACGTCGTCAAAGGTGATTTCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGT  >  NC_000913/39702‑39972

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