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	93.2%	231.4 / 3.0	74	intergenic (‑69/+40)	caiD/caiC	carnitinyl‑CoA dehydratase/putative crotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base C (35/34);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGAT  >  NC_000913/36134‑36330                                                                                                  |                                                                                                    ttGCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAAAAAAAAAAAAAAAAATGCATATTGTTTTAGAgcgt                                                                                                   <  1:455360/100‑1 (MQ=255)  tGCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGcgtg                                                                                                  >  2:376657/1‑100 (MQ=255)   gCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTTCTTTTTTAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGa                                                                                                 >  1:1700078/1‑100 (MQ=255)    cGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGat                                                                                                >  2:1582594/1‑100 (MQ=255)     gggTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGatt                                                                                               >  2:1764788/1‑100 (MQ=255)      ggTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGatta                                                                                              <  2:1313422/100‑1 (MQ=255)          aGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGGGAttatta                                                                                           >  2:1526745/1‑99 (MQ=255)          aGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAg                                                                                          <  1:1840053/100‑1 (MQ=255)               tAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGTGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGc                                                                                     >  1:1898856/1‑100 (MQ=255)                aaTGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCa                                                                                    <  2:466779/100‑1 (MQ=255)                 aTGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAg                                                                                   >  1:1753856/1‑100 (MQ=255)                 aTGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAg                                                                                   <  1:1285728/100‑1 (MQ=255)                 aTGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAg                                                                                   <  2:1140773/100‑1 (MQ=255)                   gATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAggg                                                                                 >  2:521634/1‑100 (MQ=255)                    aTTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGt                                                                                <  2:1482655/100‑1 (MQ=255)                      tCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAg                                                                              <  1:658969/100‑1 (MQ=255)                      tCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATAACTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAg                                                                              <  1:1432177/100‑1 (MQ=255)                       cACTCATTCCTTTTCTCCATTTTTTTTTTTTTTTTTACGACAACATCCCTGCAAAAAATGCATTTTGTTTTAGAGCGTGATTATTAGCTGTCTGGGTAGt                                                                             >  1:118769/1‑100 (MQ=255)                       cACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGt                                                                             >  2:712042/1‑100 (MQ=255)                       cACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGt                                                                             >  1:1331554/1‑100 (MQ=255)                           cATTCCTTTTCTCCATTTTTTTTTTTTTTTTTTCGACAACATCCCTGCAAATAATTTATATTGTTTTTGTGCGTGATTATTAGCTGTCTTGGTAGtt                                                                            >  1:222748/1‑97 (MQ=255)                           cATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTccc                                                                         >  1:1193987/1‑100 (MQ=255)                               ccTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCctgct                                                                     >  1:996582/1‑100 (MQ=255)                                cTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCctgctg                                                                    <  1:1427807/100‑1 (MQ=255)                                 ttttCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGt                                                                   <  1:1752081/100‑1 (MQ=255)                                     ctcCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCa                                                               >  1:1108406/1‑100 (MQ=255)                                     ctcCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCa                                                               <  1:1843931/100‑1 (MQ=255)                                        cATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCattt                                                            <  2:387113/100‑1 (MQ=255)                                        cATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCattt                                                            >  1:946467/1‑100 (MQ=255)                                             ttGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTc                                                       <  1:993037/100‑1 (MQ=255)                                              tGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCa                                                      <  2:1671596/100‑1 (MQ=255)                                              tGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCa                                                      <  1:319795/100‑1 (MQ=255)                                               gCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAg                                                     <  1:946090/100‑1 (MQ=255)                                                cTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGa                                                    >  2:640403/1‑100 (MQ=255)                                                cTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGa                                                    >  1:1936479/1‑100 (MQ=255)                                                cTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGa                                                    >  1:1643469/1‑100 (MQ=255)                                                cTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGa                                                    <  2:112489/100‑1 (MQ=255)                                                cTTTTTCAGGGACGACAAAATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGa                                                    >  2:1343682/1‑100 (MQ=255)                                                  ttttCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAtt                                                  >  2:1730396/1‑100 (MQ=255)                                                     tCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttctt                                               >  2:1697169/1‑100 (MQ=255)                                                        gggACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttctttct                                            <  2:44486/100‑1 (MQ=255)                                                        gggACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttctttct                                            <  1:487032/100‑1 (MQ=255)                                                        gggACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGAttctttct                                            >  1:506111/1‑100 (MQ=255)                                                         ggACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCGGATTCTTTCTa                                           <  2:497834/100‑1 (MQ=255)                                                          gacgacAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTaa                                          >  1:60143/1‑100 (MQ=255)                                                          gacgacAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTaa                                          >  2:1926506/1‑100 (MQ=255)                                                           acgacAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTGTTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAat                                         <  2:1233399/100‑1 (MQ=255)                                                           acgacAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAat                                         >  2:1422888/1‑100 (MQ=255)                                                           acgacAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAat                                         <  1:801909/100‑1 (MQ=255)                                                            cgacAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAatt                                        >  2:1045588/1‑100 (MQ=255)                                                            cgacAACATCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAatt                                        >  2:705516/1‑100 (MQ=255)                                                               caacaTCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAattatt                                     <  1:1113757/100‑1 (MQ=255)                                                                    tCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTccc                                <  2:1666804/100‑1 (MQ=255)                                                                    tCCCTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTccc                                >  2:987290/1‑100 (MQ=255)                                                                     cccTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTcccc                               >  1:1158798/1‑100 (MQ=255)                                                                     cccTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTcccc                               <  1:381342/100‑1 (MQ=255)                                                                     cccTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTcccc                               <  2:247147/100‑1 (MQ=255)                                                                     cccTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTcccc                               <  2:246851/100‑1 (MQ=255)                                                                       cTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGa                             >  2:1380142/1‑100 (MQ=255)                                                                       cTGCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGa                             <  1:1723795/100‑1 (MQ=255)                                                                         gCAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGc                           <  2:410971/100‑1 (MQ=255)                                                                          cAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCa                          >  1:1774507/1‑100 (MQ=255)                                                                          cAAAAAATGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCa                          >  1:434245/1‑100 (MQ=255)                                                                            aaaaaTGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAAt                        <  2:1700089/100‑1 (MQ=255)                                                                             aaaaTGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAAtt                       <  2:1108407/100‑1 (MQ=255)                                                                              aaaTGCATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTa                      <  2:105813/100‑1 (MQ=255)                                                                                   cATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTTCCCGAGCAATTACGTgg                 >  1:88654/1‑100 (MQ=255)                                                                                   cATATTGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTgg                 <  1:1343674/100‑1 (MQ=255)                                                                                   cATATTGTTTTAGAGCATGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTgg                 <  2:705706/100‑1 (MQ=255)                                                                                        tGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGAt            >  2:252611/1‑100 (MQ=255)                                                                                        tGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGAt            >  2:716921/1‑100 (MQ=255)                                                                                        tGTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGAt            >  2:883374/1‑100 (MQ=255)                                                                                         gTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATc           >  2:387593/1‑100 (MQ=255)                                                                                         gTTTTAGAGCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATc           >  1:248325/1‑100 (MQ=255)                                                                                              agagCGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTc      <  1:712041/100‑1 (MQ=255)                                                                                               gagTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCt     <  2:1753867/100‑1 (MQ=255)                                                                                                  tgtgATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGAt  >  1:118295/1‑100 (MQ=255)                                                                                                  |                                                                                                    TTGCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGAT  >  NC_000913/36134‑36330

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