Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	37,903	C→T	D405N (GAC→AAC)	caiB ←	crotonobetainyl CoA:carnitine CoA transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	37,903	0	C	T	100.0%	117.5 / NA	38	D405N (GAC→AAC)	caiB	crotonobetainyl CoA:carnitine CoA transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/18);  total (20/18)

TCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTCCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATT  >  NC_000913/37781‑38018                                                                                                                           |                                                                                                                    tCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGa                                                                                                     <  1:108776/139‑1 (MQ=255) tCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGa                                                                                                     >  2:640073/1‑139 (MQ=255)  cGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGAc                                                                                                    <  2:626641/139‑1 (MQ=255)  cGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGAc                                                                                                    <  2:376624/139‑1 (MQ=255)            ttGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGAc                                                                                          <  2:18492/139‑1 (MQ=255)               aCGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGATGACCaa                                                                                       <  1:547649/139‑1 (MQ=255)                cGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAAc                                                                                      >  2:567263/1‑139 (MQ=255)                        ttgtCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAAt                                                                              <  2:465515/139‑1 (MQ=255)                             cGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGt                                                                         >  2:763641/1‑139 (MQ=255)                               ccAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCACCTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGttt                                                                       <  2:711986/139‑1 (MQ=255)                               ccAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGttt                                                                       >  1:303165/1‑139 (MQ=255)                                cAATGATATCCATTGCACCTCTATCCATTTTTGTTCGCTTGTTATTGGGCGGGCGCTACTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGtttt                                                                      >  1:656274/1‑139 (MQ=255)                                         ccATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGc                                                             >  1:698046/1‑139 (MQ=255)                                             tGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGat                                                         >  1:256688/1‑139 (MQ=255)                                             tGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGat                                                         >  1:54538/1‑139 (MQ=255)                                               cACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGatat                                                       >  1:549435/1‑139 (MQ=255)                                                    ctATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTc                                                  >  1:468949/1‑139 (MQ=255)                                                       tCCATTTTTGTTCGTTTGTTATTTGGCGGGCGCTAGTCAGGCAAGCCGACCGACTCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACAAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCaaa                                               >  1:423281/1‑139 (MQ=255)                                                       tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCaaa                                               >  2:292796/1‑139 (MQ=255)                                                         cATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAAt                                             >  1:34310/1‑139 (MQ=255)                                                           gttttGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATgg                                           <  1:451959/138‑1 (MQ=255)                                                                gTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGGCCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                      <  2:303165/139‑1 (MQ=255)                                                                gTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                      <  2:644455/139‑1 (MQ=255)                                                                gTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                      <  2:584684/139‑1 (MQ=255)                                                                gTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCCGTAGCCGATATTTTTCAAAATGGCAGcc                                      <  1:292796/139‑1 (MQ=255)                                                                 ttCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCg                                     >  1:683452/1‑139 (MQ=255)                                                                 ttCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCg                                     >  2:404901/1‑139 (MQ=255)                                                                      ttGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTcc                                <  1:640073/139‑1 (MQ=255)                                                                         ttATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTccatg                             <  2:25168/139‑1 (MQ=255)                                                                          tATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTccatgc                            <  1:567263/139‑1 (MQ=255)                                                                              gggcgggcGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTccatgccatg                        <  1:501256/139‑1 (MQ=255)                                                                                   ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCg                                     >  2:317091/1‑121 (MQ=255)                                                                                   ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCg                                     <  1:317091/121‑1 (MQ=255)                                                                                   ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTccgccctcacatggc                   <  2:468949/139‑14 (MQ=255)                                                                                   ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAgg                   >  2:727597/1‑139 (MQ=255)                                                                                                  aGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGa    >  2:785726/1‑139 (MQ=255)                                                                                                   gCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAGTGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTCAGGGCATTACGCGCCAGAt   >  1:519042/1‑139 (MQ=255)                                                                                                    ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGAtt  >  1:38519/1‑139 (MQ=255)                                                                                                                           |                                                                                                                    TCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTCCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATT  >  NC_000913/37781‑38018

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