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

TCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTCCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTATTT  >  NC_000913/37784‑38033                                                                                                                        |                                                                                                                                   tCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGc                                                                                                                        >  2:1340280/1‑132 (MQ=255) tCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGc                                                                                                                        <  1:1340280/132‑1 (MQ=255)                            ccAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGttt                                                                                      >  1:1430983/1‑139 (MQ=255)                                         ttGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGa                                                                         >  2:1391856/1‑139 (MQ=255)                                          tGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCt                                                                                                                                  <  2:20886/81‑1 (MQ=255)                                          tGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCt                                                                                                                                  >  1:20886/1‑81 (MQ=255)                                          tGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCa                                                                                                       >  2:306516/1‑108 (MQ=255)                                          tGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCa                                                                                                       <  1:306516/108‑1 (MQ=255)                                            cACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                                             <  2:1427220/84‑1 (MQ=255)                                            cACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                                             >  1:1427220/1‑84 (MQ=255)                                                 ctATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGa                                                                                                          <  2:188176/98‑1 (MQ=255)                                                 ctATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGa                                                                                                          >  1:188176/1‑98 (MQ=255)                                                 ctATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTc                                                                 >  1:1040579/1‑139 (MQ=255)                                                    tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGAcc                                                                                                        >  1:915028/1‑97 (MQ=255)                                                    tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGAcc                                                                                                        <  2:915028/97‑1 (MQ=255)                                                    tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATAtt                                                                     <  1:1484715/132‑1 (MQ=255)                                                    tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATAtt                                                                     >  2:1484715/1‑132 (MQ=255)                                                    tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCaaa                                                              >  1:216794/1‑139 (MQ=255)                                                            tGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGc                                                      <  2:216794/139‑1 (MQ=255)                                                             gTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                                     <  1:406304/139‑1 (MQ=255)                                                             gTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                                     <  2:1040579/139‑1 (MQ=255)                                                              ttCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCg                                                    >  1:610139/1‑139 (MQ=255)                                                              ttCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCg                                                    <  1:1425575/139‑1 (MQ=255)                                                                    tGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTcca                                              <  2:610139/139‑1 (MQ=255)                                                                               gggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAg                                   >  1:957653/1‑139 (MQ=255)                                                                                ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGAcc                                                                                                        <  1:1080795/69‑1 (MQ=255)                                                                                ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGAcc                                                                                                        >  2:1080795/1‑69 (MQ=255)                                                                                ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCt                                                                                 <  2:1392807/92‑1 (MQ=255)                                                                                ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCt                                                                                 >  1:1392807/1‑92 (MQ=255)                                                                                ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGc                                                      >  2:922437/1‑119 (MQ=255)                                                                                ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGc                                                      <  1:922437/119‑1 (MQ=255)                                                                                      aGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTc                            <  1:1391856/139‑1 (MQ=255)                                                                                            gCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                                             >  2:1083082/1‑36 (MQ=255)                                                                                            gCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                                             <  1:1083082/36‑1 (MQ=255)                                                                                            gCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGAcc                                                                                                        >  2:749006/1‑57 (MQ=255)                                                                                            gCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGAcc                                                                                                        <  1:749006/57‑1 (MQ=255)                                                                                                 ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGAc                                                                                                         >  2:963151/1‑51 (MQ=255)                                                                                                 ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGAc                                                                                                         <  1:963151/51‑1 (MQ=255)                                                                                                 ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGAtt                 >  2:1395775/1‑139 (MQ=255)                                                                                                     cTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTc             <  2:957653/139‑1 (MQ=255)                                                                                                      tGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTg                                                                                                           <  1:1490166/44‑1 (MQ=255)                                                                                                      tGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTg                                                                                                           >  2:1490166/1‑44 (MQ=255)                                                                                                      tGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAAGGGCAGCCGTGTCCATGCCATGTGAGGGCATGCCGCGCCAGGATTGTcc            >  2:342587/1‑139 (MQ=255)                                                                                                              cgcgTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGCTTTGTCCGGGGTTAt    >  2:1391661/1‑139 (MQ=255)                                                                                                                cgTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTAttt  >  2:644568/1‑139 (MQ=255)                                                                                                                        |                                                                                                                                   TCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTCCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTATTT  >  NC_000913/37784‑38033

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