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%	94.5 / NA	31	D405N (GAC→AAC)	caiB	crotonobetainyl CoA:carnitine CoA transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/12);  total (19/12)

AACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTCCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTC  >  NC_000913/37768‑38022                                                                                                                                        |                                                                                                                        aaCGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTCGTTCGTTTGTTATTGGGAGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTtctc                                                                                                                      >  1:516456/1‑139 (MQ=255)                 tCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCt                                                                                                      >  1:561484/1‑139 (MQ=255)                     acaTTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGc                                                                                                  <  2:561484/139‑1 (MQ=255)                      caTTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCt                                                                                                 >  2:552283/1‑139 (MQ=255)                          ggACGTAGATGTTGTCCTCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCCACTTTGGCCAGACCTTTGCTGAcc                                                                                             <  1:292000/138‑1 (MQ=255)                           gACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCa                                                                                            <  1:16353/139‑1 (MQ=255)                                        tCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGc                                                                                                             >  1:165971/1‑109 (MQ=255)                                        tCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGc                                                                                                             <  2:165971/109‑1 (MQ=255)                                            ccAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGttt                                                                           >  1:694129/1‑139 (MQ=255)                                                          tGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGat                                                             >  2:741376/1‑139 (MQ=255)                                                           gCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGata                                                            >  1:556496/1‑139 (MQ=255)                                                           gCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGata                                                            >  1:732019/1‑139 (MQ=255)                                                                 ctATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTc                                                      >  2:71141/1‑139 (MQ=255)                                                                    tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCaaa                                                   >  2:358345/1‑139 (MQ=255)                                                                    tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCaaa                                                   >  2:219210/1‑139 (MQ=255)                                                                    tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCaaa                                                   >  1:720378/1‑139 (MQ=255)                                                                         ttttGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGc                                              >  1:401334/1‑139 (MQ=255)                                                                         ttttGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGc                                              >  2:656184/1‑139 (MQ=255)                                                                           ttGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAg                                            >  1:495701/1‑139 (MQ=255)                                                                            tGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGc                                           <  2:495701/139‑1 (MQ=255)                                                                            tGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGc                                           <  2:516456/139‑1 (MQ=255)                                                                             gTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                          <  2:720378/139‑1 (MQ=255)                                                                             gTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                          <  1:552283/139‑1 (MQ=255)                                                                             gTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                          <  2:61846/139‑1 (MQ=255)                                                                              ttCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCg                                         >  1:519815/1‑139 (MQ=255)                                                                                      ttATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTccatg                                 <  2:401334/139‑1 (MQ=255)                                                                                        aTTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCg                                                                        >  1:149456/1‑98 (MQ=255)                                                                                        aTTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCg                                                                        <  2:149456/98‑1 (MQ=255)                                                                                                               aGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGa        >  2:316782/1‑139 (MQ=255)                                                                                                                 ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGAtt      >  1:35323/1‑139 (MQ=255)                                                                                                                 ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGAtt      >  2:25960/1‑139 (MQ=255)                                                                                                                     cTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTc  <  1:358345/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                        AACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTCCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTC  >  NC_000913/37768‑38022

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