Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	36,231	T→C	intergenic (‑69/+40)	caiD ← / ← caiC	carnitinyl‑CoA dehydratase/putative crotonobetaine/carnitine‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	36,231	0	T	C	100.0%	188.0 / NA	53	intergenic (‑69/+40)	caiD/caiC	carnitinyl‑CoA dehydratase/putative crotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/24);  total (29/24)

TGCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTC  >  NC_000913/36135‑36326                                                                                                 |                                                                                                tGCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGcgtg                                                                                              <  2:231338/100‑1 (MQ=255)  gCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGa                                                                                             <  2:283400/100‑1 (MQ=255)   cGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGat                                                                                            <  1:802813/100‑1 (MQ=255)    gggTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTTTTTTTTTTTTTACGACAACATCCCTGCAAAAAATTCATATTGTTTTAGAGCTTGatt                                                                                           >  2:791749/1‑100 (MQ=255)     ggTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGatta                                                                                          >  2:1304591/1‑100 (MQ=255)      gTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGattat                                                                                         <  1:78933/100‑1 (MQ=255)            tGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTg                                                                                   >  1:1021085/1‑100 (MQ=255)              tAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGc                                                                                 <  1:631115/100‑1 (MQ=255)               aaTGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCa                                                                                <  1:598109/100‑1 (MQ=255)                 tGATTCACTCATTCCTTTTCTCCATTGTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAgg                                                                              <  2:550644/100‑1 (MQ=255)                        cTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTc                                                                       >  1:795219/1‑100 (MQ=255)                        cTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTc                                                                       <  2:931373/100‑1 (MQ=255)                          cATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTccc                                                                     <  2:258446/100‑1 (MQ=255)                          cATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTccc                                                                     >  1:603229/1‑100 (MQ=255)                           aTTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCct                                                                    >  1:272752/1‑100 (MQ=255)                                  ttCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGttt                                                             <  2:1179376/100‑1 (MQ=255)                                    ctcCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCa                                                           >  2:282196/1‑100 (MQ=255)                                    ctcCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATGATTAGCTGTCAGGGTTGTTCCCTGCTGTTTCa                                                           >  1:844864/1‑100 (MQ=255)                                      ccATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCatt                                                         >  1:603195/1‑100 (MQ=255)                                      ccATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCatt                                                         >  2:429128/1‑100 (MQ=255)                                          ttttGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTCGCTGGCAGGGTAGTTATCTGCTGTTTCatttat                                                      >  2:1238008/1‑99 (MQ=255)                                             tGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCa                                                  <  2:206268/100‑1 (MQ=255)                                              gCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAg                                                 <  2:1171146/100‑1 (MQ=255)                                              gCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAg                                                 <  1:714167/100‑1 (MQ=255)                                               cTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGa                                                >  1:125158/1‑100 (MQ=255)                                                      aGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttctttc                                         <  1:496667/100‑1 (MQ=255)                                                         gacgacAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTaa                                      >  1:425391/1‑100 (MQ=255)                                                          acgacAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAat                                     <  2:811366/100‑1 (MQ=255)                                                             acaacaTCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAattat                                  >  2:51297/1‑100 (MQ=255)                                                             acaacaTCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAattat                                  >  1:60851/1‑100 (MQ=255)                                                                 caTCCCTGCAAAAAATGCATATTGTTTTAGCGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAattatt                                 <  2:272749/97‑1 (MQ=255)                                                                   tCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTccc                            <  1:5578/100‑1 (MQ=255)                                                                    cccTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTcccc                           >  2:1132066/1‑100 (MQ=255)                                                                     ccTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCg                          <  2:361796/100‑1 (MQ=255)                                                                     ccTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGCGTAGTTCCCTGCTGTTTTATTTATTTCAGTTTCTTTCTAATTATTTCcccc                           <  2:784042/99‑1 (MQ=255)                                                                      cTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGa                         >  1:478479/1‑100 (MQ=255)                                                                      cTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGa                         >  2:518255/1‑100 (MQ=255)                                                                       tGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCATTCTAATTATTTTCCCCGAg                        <  2:795214/100‑1 (MQ=255)                                                                         cAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCa                      <  1:1304602/100‑1 (MQ=255)                                                                         cAAAAAATGCATATTGGTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCa                      >  2:258858/1‑100 (MQ=255)                                                                           aaaaaTGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAAt                    >  1:1284804/1‑100 (MQ=255)                                                                            aaaaTGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAAtt                   >  1:323961/1‑100 (MQ=255)                                                                                  cATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGt               >  2:546034/1‑98 (MQ=255)                                                                                  cATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTgg             >  1:616655/1‑100 (MQ=255)                                                                                   atatTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGc            >  2:570550/1‑100 (MQ=255)                                                                                   atatTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGc            >  1:434993/1‑100 (MQ=255)                                                                                    tatTGTTTTAGAGCGTGATTATTAGCTGGCAGGTTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCa           <  2:402882/100‑1 (MQ=255)                                                                                    tatTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCa           >  1:369242/1‑100 (MQ=255)                                                                                     atTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAg          <  1:182339/100‑1 (MQ=255)                                                                                          tttAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCtt     >  2:1101419/1‑100 (MQ=255)                                                                                           ttAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCttt    >  1:1252115/1‑100 (MQ=255)                                                                                            tAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCtttt   <  1:756800/100‑1 (MQ=255)                                                                                             ngagCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATAtttt   >  2:783683/2‑99 (MQ=255)                                                                                             agagCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTc  <  1:280916/100‑1 (MQ=255)                                                                                                 |                                                                                                TGCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTC  >  NC_000913/36135‑36326

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