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

CGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTCCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATT  >  NC_000913/37782‑38018                                                                                                                          |                                                                                                                    cGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGATTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGAc                                                                                                    >  1:362093/1‑139 (MQ=255)      cacaTTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTg                                                                                               <  1:1278738/139‑1 (MQ=255)            tGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGAcc                                                                                         <  2:994858/139‑1 (MQ=255)                    aTGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTg                                                                                 <  1:1124925/139‑1 (MQ=255)                              ccAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCt                                                                                                                   >  1:1292110/1‑95 (MQ=255)                              ccAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCt                                                                                                                   <  2:1292110/95‑1 (MQ=255)                              ccAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGc                                                                                                         <  2:116563/105‑1 (MQ=255)                              ccAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGc                                                                                                         >  1:116563/1‑105 (MQ=255)                              ccAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACAAACTCCTGAATGTAGttt                                                                       >  2:911693/1‑139 (MQ=255)                                    atatCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGtttt                                                                      >  1:1299952/1‑134 (MQ=255)                                    atatCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGtttt                                                                      <  2:1299952/134‑1 (MQ=255)                                            tGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGtttt                                                                      <  1:1068705/126‑1 (MQ=255)                                            tGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGtttt                                                                      >  2:1068705/1‑126 (MQ=255)                                                 ctctATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTg                                                                                            <  1:310582/99‑1 (MQ=255)                                                 ctctATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTg                                                                                            >  2:310582/1‑99 (MQ=255)                                                   ctATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                              >  2:585792/1‑79 (MQ=255)                                                   ctATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                              <  1:585792/79‑1 (MQ=255)                                                   ctATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTc                                                  >  2:1150064/1‑139 (MQ=255)                                                   ctATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTc                                                  >  2:1315621/1‑139 (MQ=255)                                                    tATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAACGTCGTTTTCGCTGTAGCCGATATTTTTCa                                                 >  1:980007/1‑139 (MQ=255)                                                      tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCt                                                                  <  1:327042/120‑1 (MQ=255)                                                      tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCt                                                                  >  2:327042/1‑120 (MQ=255)                                                           ttttGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGc                                          >  2:773447/1‑139 (MQ=255)                                                                ttCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGAc                                                                                                    <  1:882878/76‑1 (MQ=255)                                                                ttCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGAc                                                                                                    >  2:882878/1‑76 (MQ=255)                                                                  cGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGATTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCgtg                                   <  2:362093/139‑1 (MQ=255)                                                                    tttGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTc                                 <  1:1150064/139‑1 (MQ=255)                                                                               gcgggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCaa                                                                                       <  1:9853/74‑1 (MQ=255)                                                                               gcgggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCaa                                                                                       >  2:9853/1‑74 (MQ=255)                                                                                  ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAAc                                                                                                               >  1:320432/1‑47 (MQ=255)                                                                                  ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAAc                                                                                                               <  2:320432/47‑1 (MQ=255)                                                                                         gTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTcc            >  1:1180927/1‑139 (MQ=255)                                                                                             ggCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGc                                                                                                         <  1:1141088/42‑1 (MQ=255)                                                                                             ggCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGc                                                                                                         >  2:1141088/1‑42 (MQ=255)                                                                                                   ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCt                                                                  <  1:942813/75‑1 (MQ=255)                                                                                                   ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCt                                                                  >  2:942813/1‑75 (MQ=255)                                                                                                   ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGc                                       >  2:971881/1‑102 (MQ=255)                                                                                                   ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGc                                       <  1:971881/102‑1 (MQ=255)                                                                                                   ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGAtt  <  2:980007/139‑1 (MQ=255)                                                                                                                cgcgTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGAc                                                                                          >  2:417888/1‑38 (MQ=255)                                                                                                                cgcgTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGAc                                                                                          <  1:417888/38‑1 (MQ=255)                                                                                                                cgcgTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                      >  2:353324/1‑90 (MQ=255)                                                                                                                cgcgTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                      <  1:353324/90‑1 (MQ=255)                                                                                                                          |                                                                                                                    CGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTCCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATT  >  NC_000913/37782‑38018

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